Asymptotic Periodicity in the Fecally-Orally Epidemic Model in a Heterogeneous Environment

Abdelrazig K. Tarboush^1,2, Zhengdi Zhang¹

¹Faculty of Science, Jiangsu University, Zhenjiang, China.
²Department of Mathematics, Faculty of Education, University of Khartoum, Khartoum, Sudan.

DOI: 10.4236/jamp.2019.75069                

Abstract

To understand the influence of seasonal periodicity and environmental heterogeneity on the transmission dynamics of an infectious disease, we consider asymptotic periodicity in the fecally-orally epidemic model in a heterogeneous environment. By using the next generation operator and the related eigenvalue problems, the basic reproduction number is introduced and shows that it plays an important role in the existence and non-existence of a positive T-periodic solution. The sufficient conditions for the existence and non-existence of a positive T-periodic solution are provided by applying upper and lower solutions method. Our results showed that the fecally-orally epidemic model in a heterogeneous environment admits at least one positive T-periodic solution if the basic reproduction number is greater than one, while no T-periodic solution exists if the basic reproduction number is less than or equal to one. By means of monotone iterative schemes, we construct the true positive solutions. The asymptotic behavior of periodic solutions is presented. To illustrate our theoretical results, some numerical simulations are given. The paper ends with some conclusions and future considerations.

Keywords

Fecally-Orally Epidemic Model, Basic Reproduction Number, Time Periodicity, Asymptotic Behavior

1. Introduction

The geographic transmission of infectious diseases is an important issue in mathematical epidemiology and various epidemic models have been proposed and analyzed by many researchers [1] [2] [3] . For instance, Capasso and Paveri-Fontana [2] formulated and studied a non-spatial model to investigate the cholera epidemic which spread in the European Mediterranean regions. To discover the impact of spatial spread of a class of bacterial and viral diseases, Capasso and Maddalena [3] proposed the following reaction-diffusion system

$\{\begin{array}{ll}{u}_t\left(x,t\right)-d_1\Delta u\left(x,t\right)=-\mu u\left(x\mathrm{,}t\right)+\alpha v\left(x\mathrm{,}t\right),\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \\ v_t\left(x,t\right)-d_2\Delta v\left(x\mathrm{,}t\right)=-\gamma v\left(x\mathrm{,}t\right)+g\left(u\left(x\mathrm{,}t\right)\right),\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \end{array}$ (1)

where $u\left(x\mathrm{,}t\right)$ and $v\left(x\mathrm{,}t\right)$ represent the spatial density of an infectious agent and an infected human population at point x in the habitat $\Omega$ and time $t\ge 0$ ,

respectively. Here, $d_1$ , $d_2$ , $\mu$ , $\gamma$ and $\alpha$ are positive constants, $\frac{1}{\mu }$ means the lifetime of the agent in the environment, $\frac{1}{\gamma }$ denotes the mean infectious period of the human infections, $\frac{1}{\alpha }$ is the multiplicative factor of the mean

infectious agent to the human population, and $g\left(u\right)$ is the force of infection on the human population due to a concentration u of the infectious agent.

The dynamics of the spatially dependent model (1) and its corresponding Cauchy problem have been considered by many scholars. For example, the traveling waves’ solutions were studied in [4] . Wang et al. [5] studied entire solutions in a time-delayed and diffusive epidemic model. In addition, Ahn et al. [1] discussed the corresponding free boundary problem. Recently, the corresponding strongly coupled elliptic system has been proposed in [6] to investigate the existence and non-existence of coexistence states.

In recent years, a great deal of mathematical models has been developed to study the impact of seasonal periodicity and environmental heterogeneity on the dynamics of infectious diseases [7] [8] [9] . On the other hand, the periodicity of environmental factors is realistic and highly important for the dynamics of infectious diseases.

Taking this fact into account, we extend model (1) to the following reaction-diffusion system

$\{\begin{array}{ll}\frac{\partial u}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta u=-\mu \left(x\mathrm{,}t\right)u+\alpha \left(x\mathrm{,}t\right)v,\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial v}{\partial t}-d_2\left(x\mathrm{,}t\right)\Delta v=-\gamma \left(x\mathrm{,}t\right)v+g\left(u\right),\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=0,\hfill & \left(x\mathrm{,}t\right)\in \partial \Omega \times \left(0,\infty \right),\hfill \end{array}$ (2)

with the periodic condition

$u\left(x,0\right)=u\left(x,T\right),v\left(x,0\right)=v\left(x,T\right),x\in \Omega ,$ (3)

where $\Omega$ is a bounded domain in ${\mathcal{R}}^N$ with smooth boundary $\partial \Omega$ , $\nu$ denotes the outward unit vector on $\partial \Omega$ . Here, we assumed that the boundary $\partial \Omega$ is of a class $C^{2+\alpha }$ for some $0<\alpha <1$ ; and $d_1\left(x\mathrm{,}t\right)$ , $d_2\left(x\mathrm{,}t\right)$ , $\mu \left(x\mathrm{,}t\right)$ , $\alpha \left(x\mathrm{,}t\right)$ and $\gamma \left(x\mathrm{,}t\right)$ are all sufficiently smooth and strictly positive functions defined in $\bar{\Omega }\times \left[\mathrm{0,}\infty \right)$ and are T-periodic.

Moreover, in order to study the attractivity of T-periodic solutions of problem (2), (3), we consider model (2) under the initial condition

$u\left(x,0\right)={\eta }_1\left(x\right),v\left(x,0\right)={\eta }_2\left(x\right),x\in \Omega ,$ (4)

and the initial functions ${\eta }_i\left(x\right)\in C^2\left(\bar{\Omega }\right)$ $\left(i=1,2\right)$ are nonnegative.

Furthermore, we make the following hypothesis on the function g

(H₁) $g\in C^1\left(\left[0,\infty \right)\right),g\left(0\right),g^{\prime }\left(z\right)>0,\forall z\ge 0$ ;

(H₂) $\frac{g\left(z\right)}{z}$ is decreasing and

$\lim {\sup }_{z\to \infty }\frac{g\left(z\right)}{z}<{\min }_{\bar{\Omega }\times \left(0,T\right]}\left\{\frac{\mu \left(x,t\right)}{\alpha \left(x,t\right)}\right\}\cdot {{\displaystyle \min }}_{\bar{\Omega }\times \left(0,T\right]}\left\{\gamma \left(x,t\right)\right\}$

Interest in problems (2), (3) and (2), (4) is motivated by valuable results about the periodic solution of a weakly-parabolic systems [10] [11] [12] [13] , all those researchers have used the upper and lower solutions method developed by Pao [14] , which is powerful and effective to derive the periodic solutions. For instance, various computation algorithms for numerical solutions of the periodic boundary problem were studied in [11] . In [12] [13] , the stability and attractivity analysis, which are for quasimonotone nondecreasing and mixed quasimonotone reaction functions by the monotone iteration scheme, were given. It is important to mentioned that there are other standard approaches to derive the periodic solution of a weakly-parabolic systems, Hopf bifurcation theorem and Lyapunov functional [7] [15] , numerical methods [11] [16] , Poincaré index [17] , Schauder fixed point theorem [18] [19] , etc.

The remainder of this paper is organized as follows. In the next section, we deal with the basic reproduction number of problem (2), (3) and its properties. The existence and non-existence of T-periodic solutions of problem (2), (3) and the long time behavior of problem (2), (4) are given in Section 3. Section 4 is devoted to numerical simulations and a brief discussion. Finally, we give some conclusions and future considerations in Section 5.

2. Basic Reproduction Numbers

The focus of this section is to present the basic reproduction number and its properties for the corresponding system in $\Omega$ . It is worth emphasis that the basic reproduction number is defined as the expected number of secondary cases produced, in a completely susceptible population, by a typical infected individual during its entire interval of infectiousness [20] . For spatially-independent epidemic models, which are described by ordinary differential systems, the numbers are usually calculated by the next generation matrix method [21] , while for the models constructed by reaction-diffusion systems, the numbers are formulated as the spectral radius of next infection operator induced by a new infection rate matrix and an evolution operator of an infective distribution [22] , and the numbers could be expressed in the term of the principal eigenvalues of relevant eigenvalue problems [17] [23] .

Considering the linearized problem of (2), (3), we have

$\{\begin{array}{ll}\frac{\partial u}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta u=-\mu \left(x\mathrm{,}t\right)u+\alpha \left(x\mathrm{,}t\right)v,\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial v}{\partial t}-d_2\left(x\mathrm{,}t\right)\Delta v=-\gamma \left(x\mathrm{,}t\right)v+g^{\prime }\left(0\right)u,\hfill & \left(x\mathrm{,}t\right)\in \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=0,\hfill & \left(x\mathrm{,}t\right)\in \partial \Omega \times \left(0,\infty \right),\hfill \\ u\left(x,0\right)=u\left(x,T\right),v\left(x,0\right)=v\left(x,T\right),\hfill & \Omega .\hfill \end{array}$ (5)

Letting

$w=\left(\begin{array}{c}u\\ v\end{array}\right),D=\left(\begin{array}{cc}{d}_1\left(x,t\right)& 0\\ 0& d_2\left(x,t\right)\end{array}\right),$

$F\left(x,t\right)=\left(\begin{array}{cc}0& \alpha \left(x,t\right)\\ g^{\prime }\left(0\right)& 0\end{array}\right),V\left(x,t\right)=\left(\begin{array}{cc}\mu \left(x,t\right)& 0\\ 0& \gamma \left(x,t\right)\end{array}\right),$

then we have

$\{\begin{array}{ll}\frac{\partial w}{\partial t}-D\Delta w=F\left(x,t\right)w-V\left(x,t\right)w,\hfill & \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial w}{\partial \nu }=0,\hfill & \partial \Omega \times \left(0,+\infty \right).\hfill \end{array}$ (6)

Consider the following problem

$\{\begin{array}{ll}\frac{\partial w}{\partial t}-D\Delta w=-V\left(x,t\right)w,\hfill & \Omega \times \left(0,+\infty \right),\hfill \\ \frac{\partial w}{\partial \nu }=0,\hfill & \partial \Omega \times \left(0,+\infty \right).\hfill \end{array}$ (7)

and let $E\left(t\mathrm{,}s\right)$ be the evolution operator of (7), then according to the standard semigroup theory there exist positive constants $K,c_0>0$ such that

$\Vert E\left(t\mathrm{,}s\right)\Vert \le K{\text{e}}^{-c_0\left(t-s\right)}$

for any $t\ge s\mathrm{,}t\mathrm{,}s\in \mathcal{R}$ .

Suppose that $\hslash =\left(\varphi ,\psi \right)$ is the density distribution of u at the spatial location $x\in \Omega$ and time s. We use the same idea in [22] [23] to introduce the linear operator

$L\left(\hslash \right)\left(t\right):={\displaystyle {\int }_0^{\infty }}E\left(t,t-s\right)F\left(\cdot ,t-s\right)\hslash \left(\cdot ,t-s\right)\text{d}s,$

where $\hslash \in {\mathcal{C}}_T:=\left\{\hslash :t\to w\left(t\right)\in {\left(\mathcal{C}\left(\Omega \right)\right)}^2,w\left(0\right)\left(x\right)=w\left(T\right)\left(x\right)\right\}$ . According to the definition, we can easily know that L is continuous, strong positive and compact on ${\mathcal{C}}_T$ . We now define the basic reproduction number of system (5) by the spectral radius of L as follows

$R_0^N:=\rho \left(L\right).$

To ensure the existence of the basic reproduction numbers, we consider now the following linear periodic-parabolic eigenvalue problem

$\{\begin{array}{ll}\frac{\partial \varphi }{\partial t}-d_1\left(x,t\right)\Delta \varphi =\frac{\alpha \left(x,t\right)}{R}\psi -\mu \left(x,t\right)\varphi +\lambda \varphi ,\hfill & \Omega \times \left(0,T\right],\hfill \\ \frac{\partial \psi }{\partial t}-d_2\left(x,t\right)\Delta \psi =g^{\prime }\left(0\right)\varphi -\gamma \left(x,t\right)\psi +\lambda \psi ,\hfill & \Omega \times \left(0,T\right],\hfill \\ \frac{\partial \varphi }{\partial \nu }=\frac{\partial \psi }{\partial \nu }=0,\hfill & \partial \Omega \times \left(0,T\right],\hfill \\ \varphi \left(x,0\right)=\varphi \left(x,T\right),\psi \left(x,0\right)=\psi \left(x,T\right),\hfill & \Omega ,\hfill \end{array}$ (8)

where $R>0$ . Letting

$L_R=\left(\begin{array}{cc}\frac{\partial }{\partial t}-d_1\left(x,t\right)\Delta +\mu \left(x,t\right)& -\frac{\alpha \left(x,t\right)}{R}\\ -\frac{g^{\prime }\left(0\right)}{R}& \frac{\partial }{\partial t}-d_2\left(x,t\right)\Delta +\gamma \left(x,t\right)\end{array}\right)$

then the corresponding abstract eigenvalue problem of (8) becomes

$L_R\left(\begin{array}{c}\varphi \\ \psi \end{array}\right)=\lambda \left(\begin{array}{c}\varphi \\ \psi \end{array}\right)$ (9)

in the space

$X=\left\{\left(\varphi \mathrm{,}\psi \right)\in {\left({\mathcal{C}}^{\alpha \mathrm{,}\frac{\alpha }{2}}\left(\bar{\Omega }\times \left(\mathrm{0,}T\right]\right)\right)}^2\mathrm{<mo>\; :\; </mo>}\varphi \mathrm{,}\psi \text{are}T\text{-periodicwithrespectto}t\right\}\mathrm{,}$

and the domain of the operator $\text{dom}\left(L_R\right)=X_1$ is denoted by

$\begin{array}{l}{X}_1=\{\left(\varphi ,\psi \right)\in {\left({\mathcal{C}}^{2+\alpha ,1+\frac{\alpha }{2}}\left(\bar{\Omega }\times \left(0,T\right]\right)\right)}^2:\varphi =\psi =0\text{on}\partial \Omega \times \left(0,T\right],\\ \underset{\begin{array}{l}\\ \end{array}}{\overset{}{}}\varphi ,\psi \text{are}T\text{-periodicin}t\}.\end{array}$

Since the system (9) is strongly cooperative, it follows from [24] [25] [26] that for any $R>0$ , there are a unique value $\lambda :={\lambda }_1\left(R\right)$ , and called the principal eigenvalue, such that problem (8) admits a unique solution pair $\left({\varphi }_R\mathrm{,}{\psi }_R\right)$ (subject to constant multiples) with ${\varphi }_R>0$ and ${\psi }_R>0$ in $\Omega \times \left(\mathrm{0,}T\right]$ . That is, the solution pair $\left({\varphi }_R\mathrm{,}{\psi }_R\right)\in X_1$ is called the principal eigenfunction. Besides, ${\lambda }_1\left(R\right)$ is algebraically simple and dominant, and the following result hold.

Lemma 2.1 ${\lambda }_1\left(R\right)$ is continuous and strictly increasing with respect to R.

In light of the above discussion, we can further have the following relation between the two eigenvalues.

Theorem 2.1 $sign\left(1-R_0^N\right)=sign\left({\lambda }_0\right)$ , and ${\lambda }_0$ is the principal eigenvalue of the eigenvalue problem

$\{\begin{array}{ll}\frac{\partial \varphi }{\partial t}-d_1\left(x,t\right)\Delta \varphi =\alpha \left(x,t\right)\psi -\mu \left(x,t\right)\varphi +{\lambda }_0\varphi ,\hfill & \Omega \times (\left(0,T\right],\hfill \\ \frac{\partial \psi }{\partial t}-d_2\left(x,t\right)\Delta \psi =g^{\prime }\left(0\right)\varphi -\gamma \left(x,t\right)\psi +{\lambda }_0\psi ,\hfill & \Omega \times \left(0,T\right],\hfill \\ \frac{\partial \varphi }{\partial \nu }=\frac{\partial \psi }{\partial \nu }=0,\hfill & \partial \Omega \times \left(0,T\right],\hfill \\ \varphi \left(x,0\right)=\varphi \left(x,T\right),\psi \left(x,0\right)=\psi \left(x,T\right),\hfill & \Omega ,\hfill \end{array}$ (10)

and $R_0^N$ is the unique principal eigenvalue of problem

$\{\begin{array}{ll}\frac{\partial \varphi }{\partial t}-d_1\left(x,t\right)\Delta \varphi =\frac{\alpha \left(x,t\right)}{R_0^N}\psi -\mu \left(x,t\right)\varphi ,\hfill & \Omega \times \left(0,T\right],\hfill \\ \frac{\partial \psi }{\partial t}-d_2\left(x,t\right)\Delta \psi =\frac{g^{\prime }\left(0\right)}{R_0^N}\varphi -\gamma \left(x,t\right)\psi ,\hfill & \Omega \times \left(0,T\right],\hfill \\ \frac{\partial \varphi }{\partial \nu }=\frac{\partial \psi }{\partial \nu }=0,\hfill & \partial \Omega \times \left(0,T\right],\hfill \\ \varphi \left(x,0\right)=\varphi \left(x,T\right),\psi \left(x,0\right)=\psi \left(x,T\right),\hfill & \Omega .\hfill \end{array}$ (11)

Proof: According to (8) and (10), it is obvious that ${\lambda }_0={\lambda }_1\left(1\right)$ . Meanwhile, it follows from the monotonicity of coefficients in (8), one can easily deduce that ${\lim }_{R\to 0^{+}}{\lambda }_1\left(R\right)<0$ and ${\lim }_{R\to +\infty }{\lambda }_1\left(R\right)>0$ , therefore $R_0^N$ is the unique positive root of the equation ${\lambda }_1\left(R\right)=0$ . One finally gets the desired result by the monotonicity of ${\lambda }_1\left(R\right)$ with respect to R. □

Remark 2.1 Recalling the monotone non-increasing of ${\lambda }_1$ with respect to $\alpha \left(x\mathrm{,}t\right)$ in the sense that ${\lambda }_1\left({\alpha }_1\left(x,t\right)\right)<{\lambda }_1\left({\alpha }_2\left(x,t\right)\right)$ , if ${\alpha }_1\left(x,t\right)\ge {\alpha }_2\left(x,t\right)$ and ${\alpha }_1\left(x,t\right)\cancel{\equiv }{\alpha }_2\left(x,t\right)$ in $\Omega$ , one can deduce from Theorem 2.1 that $R_0^N$ is monotonically nondecreasing with respect to $\alpha \left(x\mathrm{,}t\right)$ and $R_0^N>1$ if $\alpha \left(x\mathrm{,}t\right)$ is big enough.

In what follows, we provide an explicit formula for $R_0^N$ if all coefficients are constants.

Theorem 2.2 Suppose that $d_1\left(x,t\right)=d_1^{*},d_2\left(x,t\right)=d_2^{*}$ , $\alpha \left(x,t\right)={\alpha }^{*}$ and $\mu \left(x,t\right)={\mu }^{*}$ . Then the principal eigenvalue $R_0^N$ for (11), or a threshold parameter for model (2), (4), is expressed by

$R_0^N\left(\Omega \right)=\sqrt{\frac{{\alpha }^{*}{g}^{\prime }\left(0\right)}{{\mu }^{*}{\gamma }^{*}}}.$ (12)

Proof: Let ${\psi }^{\mathrm{<mo>\; *\; </mo>}}\left(x\mathrm{,}t\right)\equiv 1$ in $\bar{\Omega }\times \left[\mathrm{0,}+\infty \right)$ and

$C^{*}=\frac{{\alpha }^{*}{g}^{\prime }\left(0\right)}{{\mu }^{*}{\gamma }^{*}},$

${\varphi }^{*}=\frac{{\alpha }^{*}}{\sqrt{C^{*}}{\mu }^{*}}{\psi }^{*}\text{in}\bar{\Omega }\times \left[0,+\infty \right).$

Then we can check that $\left({\varphi }^{\mathrm{<mo>\; *\; </mo>}}\mathrm{,}{\psi }^{\mathrm{<mo>\; *\; </mo>}}\right)$ is a positive solution of problem (11) with $R_0^N=\sqrt{R^{*}}$ , and directly lead to (12) from the uniqueness of the principal eigenvalue of problem (11).

Define $f^i={\min }_{\bar{\Omega }\times \left(0,T\right]}f\left(x,t\right)$ and $f^M={\max }_{\bar{\Omega }\times \left(0,T\right]}f\left(x,t\right)$ for any given continuous T-periodic function f, we state the following estimate after combining Theorem 2.2 and the monotonicity of $R_0^N$ with respect to all coefficients in (11).

Corollary 2.1 The principal eigenvalue $R_0^N$ for (11) satisfies

$\sqrt{\frac{g^{\prime }\left(0\right){\alpha }^i}{{\gamma }^M{\mu }^M}}\le R_0^N\left(\Omega \right)\le \sqrt{\frac{g^{\prime }\left(0\right){\alpha }^M}{{\gamma }^i{\mu }^i}}\mathrm{.}$ (13)

3. T-Periodic Solution

In this section, we discuss the existence and non-existence of a positive T-periodic solution of problem (2), (3) and the attractivity of the initial-boundary value problem (2), (4) in relation to the maximal and minimal T-periodic solutions of problem (2), (4). To begin with, the following theorem gives the non-existence of a T-periodic solution for problem (2), (3).

Theorem 3.1 If $R_0^N\le 1$ , problem (2), (3) has no positive T-periodic solution.

Proof: If the assertion is not hold, then one can suppose that $\left(u^{\mathrm{<mo>\; *\; </mo>}}\left(x,t\right)\mathrm{,}{v}^{\mathrm{<mo>\; *\; </mo>}}\left(x,t\right)\right)$ is a positive T-periodic solution of problem (2), (3), that is, $u^{*}\left(x,t\right),v^{*}\left(x,t\right)>\left(0,0\right)$ in $\bar{\Omega }\times \left[\mathrm{0,}+\infty \right)$ and satisfy

$\{\begin{array}{ll}\frac{\partial U^{*}}{\partial t}-d_1\left(x,t\right)\Delta u^{*}=\alpha \left(x,t\right)v^{*}-\mu \left(x,t\right)u^{*},\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial v^{*}}{\partial t}-d_2\left(x,t\right)\Delta v^{*}=g\left(u^{*}\right)-\gamma \left(x,t\right)v^{*},\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=0,\hfill & \left(x,t\right)\in \partial \Omega \times \left[0,\infty \right),\hfill \\ u^{*}\left(x,0\right)=u^{*}\left(x,T\right),v^{*}\left(x,0\right)=v^{*}\left(x,T\right),\hfill & x\in \bar{\Omega }.\hfill \end{array}$ (14)

From (14), we can obtain that

$\{\begin{array}{ll}\frac{\partial U^{*}}{\partial t}-d_1\left(x,t\right)\Delta u^{*}<\alpha \left(x,t\right)v^{*}-\mu \left(x,t\right)u^{*},\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial v^{*}}{\partial t}-d_2\left(x,t\right)\Delta v^{*}<g^{\prime }\left(0\right)u^{*}-\gamma \left(x,t\right)v^{*},\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=0,\hfill & \left(x,t\right)\in \partial \Omega \times \left[0,\infty \right),\hfill \\ u^{*}\left(x,0\right)=u^{*}\left(x,T\right),v^{*}\left(x,0\right)=v^{*}\left(x,T\right),\hfill & x\in \bar{\Omega }.\hfill \end{array}$ (15)

On the other hand, the principal eigenvalue ${\lambda }_0$ in problem (10) meets

$\{\begin{array}{ll}\frac{\partial \varphi }{\partial t}-d_1\left(x,t\right)\Delta \varphi =\alpha \left(x,t\right)\psi -\mu \left(x,t\right)\varphi +{\lambda }_0\varphi ,\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial \psi }{\partial t}-d_2\left(x,t\right)\Delta \psi =g^{\prime }\left(0\right)\varphi -\gamma \left(x,t\right)\psi +{\lambda }_0\psi ,\hfill & \left(x,t\right)\in \Omega \times \left[0,\infty \right),\hfill \\ \frac{\partial \varphi }{\partial \nu }=\frac{\partial \psi }{\partial \nu }=0,\hfill & \left(x,t\right)\in \partial \Omega \times \left[0,\infty \right),\hfill \\ \varphi \left(x,0\right)=\varphi \left(x,T\right),\psi \left(x,0\right)=\psi \left(x,T\right),\hfill & x\in \bar{\Omega },\hfill \end{array}$ (16)

It is follows from Lemma 2.1 that ${\lambda }_0$ is monotone non-increasing with respect to $\alpha$ , then by comparing (15) and (16) we can deduce that ${\lambda }_0<0$ . Consequently, $R_0^N>1$ by Theorem 2.1. This leads to a contradiction. □

Before approaching the existence result of a positive T-periodic solution to problem (2), (3), we need some preliminaries. Set $D=\Omega \times \left[0,\infty \right)$ , $\bar{D}=\bar{\Omega }\times \left[0,\infty \right)$ , $\Gamma =\partial \Omega \times \left[0,\infty \right)$ , $\bar{\Omega }=\Omega \cup \partial \Omega$ and

$S=\left\{\left(u\mathrm{,}v\right)\in \mathcal{C}\left(\bar{D}\right)\mathrm{<mo>\; ;\; </mo>}\left(\hat{u}\mathrm{,}\hat{v}\right)\le \left(u\mathrm{,}v\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)\mathrm{,}\left(x\mathrm{,}t\right)\in \bar{D}\right\}\mathrm{,}$

where $\left(\hat{u}\mathrm{,}\hat{v}\right)$ and $\left(\tilde{u}\mathrm{,}\tilde{v}\right)$ are given in the following definition. Then we have an equivalent form of problem (2), (3) as

$\{\begin{array}{ll}\frac{\partial u}{\partial t}-d_1\left(x,t\right)\Delta u=f_1\left(x,t,u,v\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial v}{\partial t}-d_1\left(x,t\right)\Delta v=f_2\left(x,t,u,v\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial h}{\partial \nu }=\frac{\partial m}{\partial \nu }=0,\hfill & \left(x,t\right)\in \Gamma ,\hfill \\ u\left(x,0\right)=u\left(x,T\right),v\left(x,0\right)=v\left(x,T\right),\hfill & x\in \bar{\Omega },\hfill \end{array}$ (17)

where,

$\begin{array}{l}{f}_1\left(x,t,u,v\right)=-\mu \left(x,t\right)u+\alpha \left(x,t\right)v,\\ f_2\left(x,t,u,v\right)=-\gamma \left(x,t\right)v+g\left(u\right).\end{array}$

With regard to problem (17) we state the following definition about the upper and lower solutions.

Definition 3.1 $\left(\tilde{u}\mathrm{,}\tilde{v}\right)$ , $\left(\hat{u}\mathrm{,}\hat{v}\right)$ are ordered upper and lower solutions of problem (17), if $\left(\hat{u}\mathrm{,}\hat{v}\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)$ and

$\{\begin{array}{ll}\frac{\partial u}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta u\ge f_1\left(x\mathrm{,}t\mathrm{,}u\mathrm{,}v\right)\mathrm{,}\hfill & \left(x\mathrm{,}t\right)\in D\mathrm{,}\hfill \\ \frac{\partial v}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta v\ge f_2\left(x\mathrm{,}t\mathrm{,}u\mathrm{,}v\right)\mathrm{,}\hfill & \left(x\mathrm{,}t\right)\in D\mathrm{,}\hfill \\ \frac{\partial u}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta u\le f_1\left(x\mathrm{,}t\mathrm{,}u\mathrm{,}v\right)\mathrm{,}\hfill & \left(x\mathrm{,}t\right)\in D\mathrm{,}\hfill \\ \frac{\partial v}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta v\le f_2\left(x\mathrm{,}t\mathrm{,}u\mathrm{,}v\right)\mathrm{,}\hfill & \left(x\mathrm{,}t\right)\in D\mathrm{,}\hfill \\ \frac{\partial \tilde{u}}{\partial \nu }\ge 0\ge \frac{\partial \hat{u}}{\partial \nu }\mathrm{,}\frac{\partial \tilde{v}}{\partial \nu }\ge 0\ge \frac{\partial \hat{v}}{\partial \nu }\mathrm{,}\hfill & \left(x\mathrm{,}t\right)\in \Gamma \mathrm{,}\hfill \\ \tilde{u}\left(x\mathrm{,0}\right)\ge \tilde{u}\left(x\mathrm{,}T\right)\mathrm{,}\hat{u}\left(x\mathrm{,0}\right)\le \hat{u}\left(x\mathrm{,}T\right)\mathrm{,}\hfill & x\in \bar{\Omega }\mathrm{,}\hfill \\ \tilde{v}\left(x\mathrm{,0}\right)\ge \tilde{v}\left(x\mathrm{,}T\right)\mathrm{,}\hat{v}\left(x\mathrm{,0}\right)\le \hat{v}\left(x\mathrm{,}T\right)\mathrm{,}\hfill & x\in \bar{\Omega }\mathrm{.}\hfill \end{array}$ (18)

It is worth emphasizing that the upper and lower solutions in the above definition are not required to be T-periodic with respect to t.

Now to verify the existence of a positive T-periodic solution to problem (2), (3), it suffices to find a pair of ordered upper and lower solutions to problem (2), (3). One can seek such as in the form $\left(\tilde{u},\tilde{v}\right)=\left(M_1,M_2\right)$ , $\left(\hat{u}\mathrm{,}\hat{v}\right)=\left(\delta \varphi \mathrm{,}\delta \psi \right)$ where $M_i$ $\left(i=1,2\right)$ and $\delta$ are positive constant with $\delta$ sufficiently small, $\left(\varphi \mathrm{,}\psi \right)\equiv \left(\varphi \left(x\mathrm{,}t\right)\mathrm{,}\psi \left(x\mathrm{,}t\right)\right)$ is (normalized) positive eigenfunction corresponding to ${\lambda }_0$ , and ${\lambda }_0$ is the principal eigenvalue of periodic-parabolic eigenvalue problem (10). Then it is easy to verify that $\left(\tilde{u}\mathrm{,}\tilde{v}\right)$ and $\left(\hat{u}\mathrm{,}\hat{v}\right)$ satisfy (18) if

$\{\begin{array}{l}\frac{\partial M_1}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta M_1\ge \alpha \left(x\mathrm{,}t\right)M_2-\mu \left(x\mathrm{,}t\right)M_1\mathrm{,}\hfill \\ \frac{\partial M_2}{\partial t}-d_2\left(x\mathrm{,}t\right)\Delta M_2\ge g\left(M_1\right)-\gamma \left(x\mathrm{,}t\right)M_2\mathrm{,}\hfill \\ \frac{\partial \left(\delta \varphi \right)}{\partial t}-d_1\left(x\mathrm{,}t\right)\Delta \left(\delta \varphi \right)\le \alpha \left(x\mathrm{,}t\right)\left(\delta \psi \right)-{\gamma }_b\left(x\mathrm{,}t\right)\left(\delta \varphi \right)\mathrm{,}\hfill \\ \frac{\partial \left(\delta \psi \right)}{\partial t}-D_2\left(x\mathrm{,}t\right)\Delta \left(\delta \psi \right)\le g^{\prime }\left(\wp \right)\left(\delta \varphi \right)-\gamma \left(x\mathrm{,}t\right)\left(\delta \psi \right)\hfill \end{array}$ (19)

hold, where $\wp \in \left(\mathrm{0,}\hat{u}\right)$ .

Recalling that $\lim {\sup }_{z\to \infty }\frac{g\left(z\right)}{z}<{\min }_{\bar{\Omega }\times \left(0,T\right]}\left\{\frac{\mu \left(x,t\right)}{\alpha \left(x,t\right)}\right\}\cdot {\min }_{\bar{\Omega }\times \left(0,T\right]}\left\{\gamma \left(x,t\right)\right\}$ , there exits constant $M_0$ such that $\frac{g\left(z\right)}{z}<{\min }_{\bar{\Omega }\times \left(0,T\right]}\left\{\frac{\mu \left(x,t\right)}{\alpha \left(x,t\right)}\right\}\cdot {\min }_{\bar{\Omega }\times \left(0,T\right]}\left\{\gamma \left(x,t\right)\right\}$ , for $z\ge M_0$ . This implies that the first two inequalities of (19) will hold, if we take $\left(\tilde{u},\tilde{v}\right)=\left(M_1,M_2\right)$ , where $M_1=\max \left\{M_0,{\max }_{\bar{\Omega }}u\left(x,0\right),{\max }_{\bar{\Omega }}v\left(x,0\right)\right\}$ , $M_2=M_1\cdot {\max }_{\bar{\Omega }\times \left(0,T\right]}\left\{\frac{\mu \left(x,t\right)}{\alpha \left(x,t\right)}\right\}$ .

Since $R_0^N>1$ , the principal eigenvalue of problem (10) is ${\lambda }_0<0$ , so one can select $\delta$ small enough such that the last two inequalities of (19) hold.

Therefore, the function pair $\left(\tilde{u},\tilde{v}\right)=\left(M_1,M_2\right)$ , $\left(\hat{u},\hat{v}\right)=\left(\delta \varphi ,\delta \psi \right)$ are ordered upper and lower solutions of problem (17), respectively.

To summarize the above conclusions, we have the following result.

Theorem 3.2 If $R_0^N>1$ , then problem (2), (3) admits at least one positive T-periodic solution $\left(u\left(x\mathrm{,}t\right)\mathrm{,}v\left(x\mathrm{,}t\right)\right)$ .

Remark 3.1 Suppose that all coefficients of (2) are constants. The corresponding basic reproduction ratio $R_0^N$ is represented by (12). If ${\alpha }_b^{\mathrm{<mo>\; *\; </mo>}}$ is big enough, then problem (2), (3) admits at least one positive T-periodic solution. On the other hand, if ${\alpha }_b^{\mathrm{<mo>\; *\; </mo>}}$ small enough, then problem (2), (3) has no positive T-periodic solution.

We can now construct the true solutions of problem (17) by applying the monotone iterative scheme associated with the method of upper and lower solutions. Due to

$F_1=k_{1}u+f_1,F_2=k_{2}v+f_2,$

we denote

$k_1=\underset{\left(x,t\right)\in \bar{\Omega }\times \left(0,T\right]}{\max }\mu \left(x,t\right),k_2=\underset{\left(x,t\right)\in \bar{\Omega }\times \left(0,T\right]}{\max }\gamma \left(x,t\right)$

then problem (17) is equivalent to

$\{\begin{array}{ll}\frac{\partial u}{\partial t}-d_1\left(x,t\right)\Delta u+k_{1}u=F_1\left(x,t,u,v\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial v}{\partial t}-d_2\left(x,t\right)\Delta v+k_{2}v=F_2\left(x,t,u,v\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=0,\hfill & \left(x,t\right)\in \Gamma ,\hfill \\ u\left(x,0\right)=u\left(x,T\right),v\left(x,0\right)=v\left(x,T\right),\hfill & x\in \bar{\Omega }.\hfill \end{array}$ (20)

With respect to (20), $F_1$ and $F_2$ are quasimonotone nondecreasing with respect to u and v, respectively. Choose $\left({\bar{u}}^{\left(0\right)}\mathrm{,}{\bar{v}}^{\left(0\right)}\right)=\left(M_1\mathrm{,}{M}_2\right)$ and $\left({\underset{\_}{u}}^{\left(0\right)},{\underset{\_}{v}}^{\left(0\right)}\right)=\left(\delta \varphi ,\delta \psi \right)$ as an initial iteration, we can construct a sequence $\left\{\left(u^{\left(n\right)}\mathrm{,}{v}^{\left(n\right)}\right)\right\}$ from the iteration process

$\{\begin{array}{ll}\frac{\partial {\bar{u}}^{\left(n\right)}}{\partial t}-d_1\left(x,t\right)\Delta {\bar{u}}^{\left(n\right)}+k_1{\bar{u}}^{\left(n\right)}=F_1\left(x,t,{\bar{u}}^{\left(n-1\right)},{\bar{v}}^{\left(n-1\right)}\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial {\bar{v}}^{\left(n\right)}}{\partial t}-d_2\left(x,t\right)\Delta {\bar{v}}^{\left(n\right)}+k_2{\bar{v}}^{\left(n\right)}=F_2\left(x,t,{\bar{u}}^{\left(n-1\right)},{\bar{v}}^{\left(n-1\right)}\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial {\underset{\_}{u}}^{\left(n\right)}}{\partial t}-d_1\left(x,t\right)\Delta {\underset{\_}{u}}^{\left(n\right)}+k_1{\underset{\_}{u}}^{\left(n\right)}=F_1\left(x,t,{\underset{\_}{u}}^{\left(n-1\right)},{\underset{\_}{v}}^{\left(n-1\right)}\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial {\underset{\_}{v}}^{\left(n\right)}}{\partial t}-d_2\left(x,t\right)\Delta {\underset{\_}{v}}^{\left(n\right)}+k_2{\underset{\_}{v}}^{\left(n\right)}=F_2\left(x,t,{\underset{\_}{u}}^{\left(n-1\right)},{\underset{\_}{v}}^{\left(n-1\right)}\right),\hfill & \left(x,t\right)\in D,\hfill \\ \frac{\partial {\bar{u}}^{\left(n\right)}}{\partial \nu }=\frac{\partial {\bar{v}}^{\left(n\right)}}{\partial \nu }=0,\frac{\partial {\underset{\_}{u}}^{\left(n\right)}}{\partial \nu }=\frac{\partial {\underset{\_}{v}}^{\left(n\right)}}{\partial \nu }=0,\hfill & \left(x,t\right)\in \Gamma ,\hfill \\ {\bar{u}}^{\left(n\right)}\left(x,0\right)={\bar{u}}^{\left(n-1\right)}\left(x,T\right),{\bar{v}}^{\left(n\right)}\left(x,0\right)={\bar{v}}^{\left(n-1\right)}\left(x,T\right),\hfill & x\in \bar{\Omega },\hfill \\ {\underset{\_}{u}}^{\left(n\right)}\left(x,0\right)={\underset{\_}{u}}^{\left(n-1\right)}\left(x,T\right),{\underset{\_}{v}}^{\left(n\right)}\left(x,0\right)={\underset{\_}{v}}^{\left(n-1\right)}\left(x,T\right),\hfill & x\in \bar{\Omega }.\hfill \end{array}$ (21)

where $n=1,2,\cdots$ .

Then we can easily see that the sequences $\left\{\left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\right\}$ and $\left\{\left({\underset{\_}{u}}^{\left(n\right)},{\underset{\_}{v}}^{\left(n\right)}\right)\right\}$ governed by (21) are well-defined. Consequently, to show the monotone property of these sequences we have the following result.

Lemma 3.1. The maximal and minimal sequences $\left\{\left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\right\}$ and $\left\{\left({\underset{\_}{u}}^{\left(n\right)},{\underset{\_}{v}}^{\left(n\right)}\right)\right\}$ are well-defined and possess the monotone property

$\left(\hat{u}\mathrm{,}\hat{v}\right)\le \left({\underset{\_}{u}}^{\left(n\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n\right)}\right)\le \left({\underset{\_}{u}}^{\left(n+1\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n+1\right)}\right)\le \left({\bar{v}}^{\left(n+1\right)}\mathrm{,}{\bar{v}}^{\left(n+1\right)}\right)\le \left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)\mathrm{,}$

for every $n=1,2,\cdots$ .

According to above lemma, the pointwise limits

$\underset{n\to \infty }{\lim }\left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)=\left(\bar{u}\mathrm{,}\bar{v}\right)\mathrm{,}\underset{n\to \infty }{\lim }\left({\underset{\_}{u}}^{\left(n\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n\right)}\right)=\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$

exist and satisfy the relation

$\left(\hat{u}\mathrm{,}\hat{v}\right)\le \left({\underset{\_}{u}}^{\left(n\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n\right)}\right)\le \left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)\le \left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)$ (22)

for every $n=1,2,\cdots$ .

Therefore, the maximal and minimal property of $\left(\bar{u}\mathrm{,}\bar{v}\right)$ and $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$ is in the sense that $\left(u\mathrm{,}v\right)$ is any other T-periodic solution of (17) in $\langle \left(\hat{u}\mathrm{,}\hat{v}\right)\mathrm{,}\left(\tilde{u}\mathrm{,}\tilde{v}\right)\rangle$ , then $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\le \left(u\mathrm{,}v\right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)$ for $\left(x\mathrm{,}t\right)\in \bar{D}$ .

Next, to shows that $\left(\bar{u}\mathrm{,}\bar{v}\right)$ and $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$ are the maximal and minimal solutions of (2), (3), respectively, we have the following theorem.

Theorem 3.3 Let $\left(\tilde{u}\mathrm{,}\tilde{v}\right)$ and $\left(\hat{u}\mathrm{,}\hat{v}\right)$ be a pair of ordered upper and lower solutions of (2), (3), respectively, then the sequences $\left\{\left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\right\}$ and $\left\{\left({\underset{\_}{u}}^{\left(n\right)},{\underset{\_}{v}}^{\left(n\right)}\right)\right\}$ provided from (21) converge monotonically from above to a maximal T-periodic solution $\left(\bar{u}\mathrm{,}\bar{v}\right)$ and from below to a minimal T-periodic solution $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$ in S, respectively, and satisfy relation

$\begin{array}{c}\left(\bar{u}\mathrm{,}\bar{v}\right)\le \left({\underset{\_}{u}}^{\left(n\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n\right)}\right)\le \left({\underset{\_}{u}}^{\left(n+1\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n+1\right)}\right)\le \left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)\le \left({\bar{u}}^{\left(n+1\right)}\mathrm{,}{\bar{v}}^{\left(n+1\right)}\right)\\ \le \left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)\mathrm{,\; <mi>\; </mi>}\text{for}\text{every}n=\mathrm{1,2,}\cdots \end{array}$

Moreover, if $\bar{u}=\underset{\_}{u}$ or $\bar{v}=\underset{\_}{v}$ then $\left(\bar{u}\mathrm{,}\bar{v}\right)=\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\left(\equiv \left(u^{\mathrm{<mo>\; *\; </mo>}}\mathrm{,}{v}^{\mathrm{<mo>\; *\; </mo>}}\right)\right)$ and $\left(u^{\mathrm{<mo>\; *\; </mo>}}\mathrm{,}{v}^{\mathrm{<mo>\; *\; </mo>}}\right)$ is the unique T-periodic solution of (2), (3) in S.

In what follows, we will study the attractivity of problem (2), (4) in relation to the maximal and minimal T-periodic solutions of problem (2), (3). The following lemma plays an important role in the establishment of attractivity result and its proof is similar to that in [12] , so we omit it here.

Lemma 3.2 let $\left(u\mathrm{,}v\right)\left(x\mathrm{,}t\mathrm{<mo>\; ;\; </mo>}\eta \right)$ be the solution of (2), (4). Then

$\left({\underset{\_}{u}}^{\left(n\right)}\mathrm{,}{\underset{\_}{v}}^{\left(n\right)}\right)\left(x\mathrm{,}t\right)\le \left(u\mathrm{,}v\right)\left(x\mathrm{,}t+nT\mathrm{<mo>\; ;\; </mo>}\eta \right)\le \left({\bar{u}}^{\left(n\right)}\mathrm{,}{\bar{v}}^{\left(n\right)}\right)\left(x\mathrm{,}t\right)$

on $\bar{D}$ for every $n=1,2,\cdots$ .

With the help of above lemma and Theorem 3.1 of [12] , the solution $\left(u\mathrm{,}v\right)$ of the problem (2) under the initial condition (4) possesses the following convergence

$\underset{n\to \infty }{\lim }\left(u\mathrm{,}v\right)\left(x\mathrm{,}t+nT\mathrm{<mo>\; ;\; </mo>}\eta \right)=\{\begin{array}{l}\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\left(x\mathrm{,}t\right)\mathrm{<mi>\; </mi>\; <mi>\; </mi>}\text{if}\mathrm{<mi>\; </mi>}\left(\hat{u}\mathrm{,}\hat{v}\right)\le \eta \left(x\right)\le \left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\text{in}\bar{\Omega }\mathrm{,}\\ \left(\bar{u}\mathrm{,}\bar{v}\right)\left(x\mathrm{,}t\right)\mathrm{<mi>\; </mi>}\text{if}\mathrm{<mi>\; </mi>}\left(\bar{u}\mathrm{,}\bar{v}\right)\le \eta \left(x\right)\le \left(\tilde{u}\mathrm{,}\tilde{v}\right)\text{in}\bar{\Omega }\mathrm{.}\end{array}$ (23)

and

$\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\left(x\mathrm{,}t\right)\le \left(u\mathrm{,}v\right)\left(x\mathrm{,}t+nT\mathrm{<mo>\; ;\; </mo>}\eta \right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)\left(x\mathrm{,}t\right)\text{on}\bar{D}\text{as}n\to \infty \mathrm{.}$ (24)

Consequently, the function pair $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$ and $\left(\bar{u}\mathrm{,}\bar{v}\right)$ are positive T-periodic solutions of problem (2), (3) and satisfy the following relation $\left(\delta \varphi \mathrm{,}\delta \psi \right)\le \left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)\le \left(M_1\mathrm{,}{M}_2\right)$ on $\bar{D}$ .

On one hand, for any nonnegative and nontrivial $\left({\eta }_1\left(x\right)\mathrm{,}{\eta }_2\left(x\right)\right)$ the solution $\left(u\mathrm{,}v\right)$ with $\left({\eta }_1\left(x\right)\mathrm{,}{\eta }_2\left(x\right)\right)$ is positive in $\bar{\Omega }$ for any $t>0$ . So for any $t_0>0$ , we can find $\delta >0$ , such that

$\left(\delta \varphi \mathrm{,}\delta \psi \right)\le \left(u\left(x\mathrm{,}{t}_0\right)\mathrm{,}v\left(x\mathrm{,}{t}_0\right)\right)\le \left(M_1\mathrm{,}{M}_2\right)\text{on}\bar{D}\mathrm{.}$

The above derivations lead to the following theorem.

Theorem 3.4 Let $\left(u\left(x\mathrm{,}t\mathrm{<mo>\; ;\; </mo>}{\eta }_1\right)\mathrm{,}v\left(x\mathrm{,}t\mathrm{<mo>\; ;\; </mo>}{\eta }_2\right)\right)$ be the solution of (2), (4) for $\left({\eta }_1\mathrm{,}{\eta }_2\right)$ with $0<{\eta }_i\le M_i$ , i = 1, 2 and let condition $R_0^N>1$ be satisfied. Then we have

1) problem (2), (3) has maximal and minimal positive T-periodic solutions $\left(\bar{u}\mathrm{,}\bar{v}\right)$ and $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)$ such that $\left(\underset{\_}{u}\mathrm{,}\underset{\_}{v}\right)\le \left(\bar{u}\mathrm{,}\bar{v}\right)$ on $\bar{D}$ ;

2) the solution $\left(u\left(x\mathrm{,}t\mathrm{<mo>\; ;\; </mo>}{\eta }_1\right)\mathrm{,}v\left(x\mathrm{,}t\mathrm{<mo>\; ;\; </mo>}{\eta }_2\right)\right)$ of (2), (4) possesses the convergence properties (23) and (24);

3) if $\left(\bar{u},\bar{v}\right)=\left(\underset{\_}{u},\underset{\_}{v}\right)=\left(u^{*},v^{*}\right)$ , then

$\underset{n\to \infty }{\lim }\left(u\left(x\mathrm{,}t+nT\mathrm{<mo>\; ;\; </mo>}{\eta }_1\right)\mathrm{,}v\left(x\mathrm{,}t+nT\mathrm{<mo>\; ;\; </mo>}{\eta }_2\right)\right)=\left(u^{\mathrm{<mo>\; *\; </mo>}}\left(x\mathrm{,}t\right)\mathrm{,}{v}^{\mathrm{<mo>\; *\; </mo>}}\left(x\mathrm{,}t\right)\right)\mathrm{,}t>\mathrm{0,}x\in \bar{\Omega }\mathrm{.}$ .

4. Numerical Simulation and Discussion

In this section, we provide some simulations for problem (2), (4) to illustrate our analytical results. We first fix the following diffusion coefficients and parameters:

$d_1\left(x,t\right)=1.5+0.5\left(\cos \left(x\right)+\cos \left(10t\right)\right),d_2\left(x,t\right)=1.4+0.5\left(\sin \left(x\right)+\sin \left(10t\right)\right),$

$\mu \left(x,t\right)=2+0.5\left(\sin \left(x\right)+\sin \left(10t\right)\right),\gamma \left(x,t\right)=1.5+0.5\left(\cos \left(x\right)+\cos \left(10t\right)\right)$

and $g\left(u\right)=\sqrt{u+1}-1$ .

Therefore, now we will change the value of parameter $\alpha \left(x\mathrm{,}t\right)$ and then observe the long time behavior of problem (2), (4).

Example 1: Chose big $\alpha =20+0.2\left(\sin \left(x\right)+\sin \left(10t\right)\right)$ so that

$R_0^N\ge \sqrt{\frac{{\alpha }^i{g}^{\prime }\left(0\right)}{{\gamma }^M{\mu }^M}}=\sqrt{\frac{19.6\times \frac{1}{2}}{3\times 2.5}}>1,$

and from Figure 1 it is easy to see that the solutions $u\left(x\mathrm{,}t\right)$ and $v\left(x\mathrm{,}t\right)$ of (2), (4) tends to a positive T-periodic solution.

Example 2: Take small $\alpha =0.05+0.02\left(\sin \left(x\right)+\sin \left(10t\right)\right)$ , then

$R_0^N\le \sqrt{\frac{{\alpha }^M{g}^{\prime }\left(0\right)}{{\gamma }^i{\mu }^i}}=\sqrt{\frac{0.09\times \frac{1}{2}}{1\times 0.5}}<1,$

hence the solutions $u\left(x\mathrm{,}t\right)$ and $v\left(x\mathrm{,}t\right)$ of (2), (4) decays to zero and no positive periodic solution exists to problem (2), (3) (see Figure 2).

To understand the impact of seasonal periodicity and environmental heterogeneity on the fecally-orally epidemic model, we consider a T-periodic

Figure 1. The solution of (2), (4) with big $\alpha$ presents a state of periodic oscillation.

Figure 2. The solution of (2), (4) with small $\alpha$ decays quickly to zero.

solution of problem (2), (3) and the attractivity of problem (2), (4) in relation to the maximal and minimal T-periodic solutions of problem (2), (3). Firstly, we introduce the basic reproduction number $R_0^N$ by using the next generation operator and associated eigenvalue problems, which is known as a threshold parameter of problem (2), (3) (Theorem 2.1). In the case that all coefficients are constants, we provide an explicit formula for $R_0^N$ (Theorem 2.2). Secondly, the existence of a T-periodic solution of problem (2), (3) is investigated by combining the method of upper and lower solutions, the eigenvalue problems and associated monotone iterative schemes. Our results indicate that if $R_0^N>1$ and $\alpha \left(x\mathrm{,}t\right)$ big enough problem (2), (3) admits at least one positive T-periodic solution (Theorem 3.4, Theorem 3.5 and Figure 1), while if $R_0^N\le 1$ and $\alpha \left(x\mathrm{,}t\right)$ sufficiently small, problem (2), (3) has no positive T-periodic solution (Theorem 3.3 and Figure 2).

Moreover, we discuss the attractivity of problem (2), (4) in relation to the maximal and minimal T-periodic solutions of problem (2), (3) by applying the monotone convergence property (Lemma 3.3 and Theorem 3.6). As result, it shown that the long time state of the solution to problem (2), (4) is between the maximal and minimal T-periodic solutions of problem (2), (3).

5. Conclusion

To better understand the impact of environmental heterogeneity and seasonal periodicity on the spatial spread of a class of bacterial and viral diseases, the fecally-orally epidemic model in heterogeneous environment has been considered. By means of next generation infection operator and associated eigenvalue problems, the spatial-temporal basic reproduction number is defined. We use the number to study the existence and non-existence of T-periodic solutions. The attractivity of our problem in relation to the maximal and minimal T-periodic solutions is discussed. We conclude that if the environment is T-periodic, then the solution of any initial-boundary system will present T-periodic phenomena gradually. Moreover, we believe that the investigation of the transmission dynamics of infectious diseases in heterogeneous environment is more close to reality than in homogeneous environment. However, here we considered the periodic solutions on specified period T. In the future, we will consider the time-periodic solutions for various possible periods. Moreover, the uniqueness of T-periodic solutions is still unclear.

Acknowledgements

We thank the Editor and the referee for their comments. The work is supported by People’s Republic of China grant National Natural Science Foundation (11872189, 11472116).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Ahn, I.K., Baek, S.Y. and Lin, Z.G. (2016) The Spreading Fronts of an Infective Environment in a Man-Environment-Man Epidemic Model. Applied Mathematical Modelling, 40, 7082-7101. https://doi.org/10.1016/j.apm.2016.02.038
[2]	Capasso, V. and Paveri-Fontana, S.L. (1980) A Mathematical Model for 1973 Cholera Epidemic in European Mediterranean Region. Revue d’épidémiologie et de Santé Publique, 27, 121-132.
[3]	Capasso, V. and Maddalena, L. (1981) Convergence to Equilibrium States for a Reaction-Diffusion System Modeling the Spatial Spread of a Class of Bacterial and Viral Diseases. Journal of Mathematical Biology, 13, 173-184. https://doi.org/10.1007/BF00275212
[4]	Xu, D.S. and Zhao, X.Q. (2005) Asymptotic Speeds of Spread and Traveling Waves for a Nonlocal Epidemic Model. Discrete and Continuous Dynamical Systems—Series B, 5, 1043-1056. https://doi.org/10.3934/dcdsb.2005.5.1043
[5]	Wang, Y.X. and Wang, Z.C. (2013) Entire Solutions in a Time-Delayed and Diffusive Epidemic Model. Applied Mathematics and Computation, 219, 5033-5041. https://doi.org/10.1016/j.amc.2012.11.029
[6]	Ge, J., Lin, Z. and Zhang, Q. (2017) Influence of Cross-Diffusion on the Fecally-Orally Epidemic Model with Spatial Heterogeneity. Advances in Difference Equations, 2017, 371. https://doi.org/10.1186/s13662-017-1422-9
[7]	Li, Y. and Zhang, T. (2016) Existence and Multiplicity of Positive Almost Periodic Solutions for a Non-Autonomous SIR Epidemic Model. Bulletin of the Malaysian Mathematical Sciences Society, 39, 359-379. https://doi.org/10.1007/s40840-015-0176-3
[8]	Tarboush, A.K., Ge, J. and Lin, Z.G. (2017) Asymptotic Periodicity in a Diffusive West Nile Virus Model in a Heterogeneous Environment. International Journal of Biomathematics, 10, Article ID: 1750110. https://doi.org/10.1142/S1793524517501108
[9]	Wang, B.G., Li, W.T. and Wang, Z.C. (2015) A Reaction-Diffusion SIS Epidemic Model in an Almost Periodic Environment. Zeitschrift für Angewandte Mathematik und Physik, 66, 3085-3108. https://doi.org/10.1007/s00033-015-0585-z
[10]	Pao, C.V. (1999) Periodic Solutions of Parabolic Systems with Nonlinear Boundary Conditions. Journal of Mathematical Analysis and Applications, 234, 695-716. https://doi.org/10.1006/jmaa.1999.6412
[11]	Pao, C.V. (2001) Numerical Methods of Time-Periodic Solutions for Nonlinear Parabolic Boundary Value Problems. SIAM Journal on Numerical Analysis, 39, 647-667. https://doi.org/10.1137/S0036142999361396
[12]	Pao, C.V. (2005) Stability and Attractivity of Periodic Solutions of Parabolic Systems with Time Delays. Journal of Mathematical Analysis and Applications, 304, 423-450. https://doi.org/10.1016/j.jmaa.2004.09.014
[13]	Zhou, L. and Fu, Y.P. (2000) Periodic Quasimonotone Global Attractor of Nonlinear Parabolic Systems with Discrete Delays. Journal of Mathematical Analysis and Applications, 250, 139-161. https://doi.org/10.1006/jmaa.2000.6986
[14]	Pao, C.V. (1992) Nonlinear Parabolic and Elliptic Equations. Springer Science and Business Media, Berlin. https://doi.org/10.1007/978-1-4615-3034-3
[15]	Xu, C., Zhang, Q., Liao, M., et al. (2013) Existence and Global Attractivity of Positive Periodic Solutions for a Delayed Competitive System with the Effect of Toxic Substances and Impulses. Applied Mathematics, 58, 309-328. https://doi.org/10.1007/s10492-013-0015-5
[16]	Yang, J. (2015) A Numerical Method for Computing Time-Periodic Solutions in Dissipative Wave Systems. Studies in Applied Mathematics, 134, 420-455. https://doi.org/10.1111/sapm.12071
[17]	Alexander, M.E. and Moghadas, S.M. (2004) Periodicity in an Epidemic Model with a Generalized Non-Linear Incidence. Mathematical Biosciences, 189, 75-96. https://doi.org/10.1016/j.mbs.2004.01.003
[18]	Kuniya, T. (2014) Existence of a Nontrival Periodic Solution in an Age-Structured SIR Epidemic Model with Time Periodic Coefficients. Applied Mathematics Letters, 27, 15-20. https://doi.org/10.1016/j.aml.2013.08.008
[19]	Wang, C. (2008) Existence and Stability of Periodic Solutions for Parabolic Systems with Time Delays. Journal of Mathematical Analysis and Applications, 339, 1354-1361. https://doi.org/10.1016/j.jmaa.2007.07.082
[20]	Diekmann, O., Heesterbeek, J.A.P. and Metz, J.A.J. (1990) On the Definition and the Computation of the Basic Reproduction Ratio R0 in Models for Infectious Diseases in Heterogeneous Populations. Journal of Mathematical Biology, 28, 365-382. https://doi.org/10.1007/BF00178324
[21]	van den Driessche, P. and Watmough, J. (2002) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48. https://doi.org/10.1016/S0025-5564(02)00108-6
[22]	Wang, W.D. and Zhao, X.Q. (2012) Basic Reproduction Numbers for Reaction-Diffusion Epidemic Models. SIAM Journal on Applied Dynamical Systems, 11, 1652-1673. https://doi.org/10.1137/120872942
[23]	Zhao, X.Q. (2017) Dynamical Systems in Population Biology. Second Edition, CMS Books in Mathematics/Ouvrages de Mathématiques de la SMC, Springer, Cham.
[24]	álvarez-Caudevilla, P. and López-Gómez, J. (2008) Asymptotic Behaviour of Principal Eigenvalues for a Class of Cooperative Systems. Journal of Differential Equations, 244, 1093-1113. https://doi.org/10.1016/j.jde.2007.10.004
[25]	de Figueiredo, D.G. and Mitidieri, E. (1986) A Maximum Principle for an Elliptic System and Applications to Semilinear Problems. SIAM Journal on Mathematical Analysis, 17, 836-849. https://doi.org/10.1137/0517060
[26]	Sweers, G. (1992) Strong Positivity in for Elliptic Systems. Mathematische Zeitschrift, 209, 251-271. https://doi.org/10.1007/BF02570833