# Modeling the Impact of Seasonal Weather Variations on the Infectiology of Brucellosis

Authors: Nkuba Nyerere; Livingstone S. Luboobi; Saul C. Mpeshe; Gabriel M. Shirima
Journal: Computational and Mathematical Methods in Medicine
 (2020)

---

## Abstract

A deterministic mathematical model for brucellosis that incorporates seasonality on direct and indirect transmission parameters for domestic ruminants, wild animals, humans, and the environment was formulated and analyzed in this paper. Both analytical and numerical simulations are presented. From this study, the findings show that variations in seasonal weather have the great impact on the transmission dynamics of brucellosis in humans, livestock, and wild animals. Thus, in order for the disease to be controlled or eliminated, measures should be timely implemented upon the fluctuation in the transmission of the disease.


---

## Body

## 1. Introduction

Brucellosis is a bacterial zoonosis that causes potential loss of production in livestock and undulant fever in humans in many countries all over the world [1]. The infection is caused by the genus Brucella with B. melitensis, B. suis, and B. abortus being predominant in domestic animals and also infecting humans [2–4]. International organizations like the World Organisation for Animal Health (Office International des Epizooties (OIE)), the World Health Organization (WHO), and the Food and Agriculture Organization (FAO) identify brucellosis as one of the most prevalent zoonoses in the world alongside bovine tuberculosis and rabies [5].In most parts of the developing world, brucellosis is endemic and leads to devastating losses in the livestock industry especially to smallholder keepers and to international market [6]. The disease results in huge financial losses by causing abortions, sterility, decreased milk production, veterinary fees, and cost of replacing animals. In many countries of sub-Saharan Africa, the control of the disease had proven to be a challenge because of different farming systems, low community awareness about the disease, poor health network systems, weak surveillance programmes, and limited vaccinations [7]. In animals, brucellosis is transmitted when a susceptible animal ingests contaminated materials such as pastures or discharges from infected animals while in humans, the bacteria are transmitted through ingestion of contaminated raw blood, meat, dairy products, and unpasteurized milk. Brucellosis is an occupational disease to abattoir workers, farmers, veterinarians, and laboratory personnel through direct contact with aborted materials and discharges, handling of suspected samples, and handling of livestock during deliveries [8]. Although traditionally Brucella species are host specific, recent studies revealed that cattle are also susceptible to B. melitensis [9–11].Infected animals exhibit clinical signs like reduced fertility, late-term abortion, considerable drop in milk production, retained placenta, metritis, and hygromas in chronic cases in cattle [6, 12]. Symptoms in humans include headache, weakness, continuous or intermittent fever, chills, joint pains, profuse sweats, weight loss, aches, and devastating complications that may lead to miscarriage during the first trimester in pregnant women. Endocarditis, bone abscesses, or testicular and neurological complications can also occur [1, 13]. Human brucellosis is debilitating and needs prolonged treatment using a combination of antibiotics [14]. Furthermore, the clinical signs of the disease in humans are not pathognomonic; hence, patients were clinically misdiagnosed with malaria, rheumatic fever, typhoid fever, elapsing fever, and joint diseases [15].Globally, the burden of human brucellosis remains huge; more than 500,000 new cases per year are reported [8]. Brucellosis exists throughout the sub-Saharan African region, it is poorly understood with fluctuating records from one country to another, and its prevalence is still unclear [16]. In many parts of Tanzania, brucellosis is a highly prevalent disease. However, very limited data is available regarding its distribution, affected host species, and impact. In addition, it has been demonstrated that the cattle seroprevalence level in various production systems, zones, and regions varies from 1 to 30% while in humans, the average prevalence is from 1 to 5% [17]. A study by Carugati et al. [18] shows that brucellosis incidences are moderate in the northern part of Tanzania and that it is a common human health problem since it is endemic in the region. Human brucellosis cases have also been reported in parts of eastern, lake, and western regions of Tanzania with seroprevalences varying from 0.7 to 20.5% [19, 20].The incidence and prevalence of most infectious diseases are directly linked to seasonal weather variations. The understanding of seasonal patterns in infectious disease occurrences dates back to the Hippocratic era [21]. The seasonal weather variations influence the dynamics of infectious diseases by affecting the host-pathogen interactions which alters the components of the reproduction number [22]. In particular, cold or wet seasons are associated with high disease incidences due to the abundance, survival, and virulence of pathogens and the fact that most people spend their time in poorly ventilated houses. On the other hand, warm or dry seasons are associated with decreased disease incidences due to increased outdoor activities and exposure of the pathogens to UV light. In addition, the survival of pathogens outside their hosts depends on other environmental factors such as humidity, salinity, temperature, and soil pH, abundance of vectors and nonhuman hosts, host immune function, and host behavior [23].Mathematical models can give insight into how the mechanisms and strength of seasonality affect the persistence and spreading of communicable diseases. In this view, understanding the impact of seasonality and timing offers important intuitions on parasite-host system operation, how and when the parasite control measures should be applied, and the response of disease risks to anthropogenic climate change and patterns of seasonality.Seasonal variations are exhibited in brucellosis incidences where a large number of new cases are expected in months with wet or dry seasons of the year in both developing and developed countries [19]. The disease incidence is higher during the wet season; breeding is synchronized for animals to give birth during the wet season when pastures are available. Pastoral and agropastoral settings depend on natural pastures. During this time, infected animals shed pathogens into the environment through birth fluids and tissues that contaminate pastures and the surroundings. In addition, during the wet season, it is anticipated that the cold weather favours survival of Brucella pathogens in the environment compared to the hot dry season hence influencing the transmission rate [24]. For instance, high transmission rates between domestic and wild animals are expected during the dry season due to sharing of pastures and water points, while the within-herd transmission is expected during the wet season due to a high birth rate and abortion storms [25]. According to the WHO [8], in countries with cold or temperate climates, there are notable seasonal variations in brucellosis incidences with most occurring cases in the summer and spring. This concurs with the peak period for parturitions and abortions in animals and consequently the highest level of exposure to other animals and people consuming their products or attending the animals. Seasonality in transmission dynamics of the disease is also attributed to seasonal livestock movements due to the availability of water and grasslands. This is the common practice in sub-Saharan Africa countries; for instance, during the dry seasons, 83.1% of the cattle owners in Northern Tanzania move their cattle away from homes for pasture and water needs [25]. This changes the disease dynamics since the concentration of animals is expected near water bodies and wildlife parks and increases the contact rates between susceptible and infected animals.Brucella is a robust pathogen, and it can persevere outside and inside the mammalian hosts for a long time despite the unfriendly conditions; it remains in food for up to 15 months given adverse conditions such as acidity and temperature between 14°C and 11°C or for two to three days under 37°C. When Brucella is exposed directly to sunlight, it may survive for few hours while its survival in contaminated manure and aborted foeti is more than 2 months during the winter season [26]. Furthermore, in an ideal environment, the survival of Brucella spp. is reported to last up to 135 days [27]. Therefore, to estimate the impact of seasonality on brucellosis transmission in animals and humans using mathematical modeling becomes imperative to device timely interventions. Despite the fact that the WHO, FAO, and OIE efforts and interventions are available, brucellosis continues to pose great economic threats and it affects livelihoods and food security mostly in developing countries. Thus, there is need to assess the impact of the current control strategies if we are to control or eradicate the disease. So far, a few studies [28–34] analyzed the dynamics and spread of brucellosis in homogeneous/heterogeneous populations. However, none of these studies have considered the mathematical approach to analyze the impact of seasonal weather variations on the transmission of brucellosis in human, livestock, and wildlife populations. In this paper, the impacts of seasonal weather parameters on the transmission of brucellosis are studied using a mathematical model.


## 2. Model Formulation

A deterministic mathematical model that illustrates the transmission of brucellosis in humans and domestic and wild animals is formulated and analyzed under this section. More importantly, in incorporating the variations on seasonal weather in both direct and indirect transmission routes of the disease, we follow the approach presented in [33, 35, 36]. The stimuli of seasonal variations on the direct transmission of brucellosis in domestic ruminants, humans, and wild animals are, respectively, modeled by the periodic continuous functions βat=b11+a1sinωt, βht=b21+a2sinωt, and βwt=b31+a3sinωt while the indirect transmission in the three populations is captured by αat=c11+r1sinωt, αht=c21+r2sinωt, and αwt=c31+r3sinωt, respectively.Furthermore, we consider the pathogen shedding rate by the infective livestock and wild animals to be represented by the periodic functions of the formρt=ρ01+ρ1sinωt and ρwt=ρ21+ρ3sinωt, respectively. The decaying rate of the pathogens in the environment is also represented by the periodic continuous function εt=ε01+ε1sinωt. The constants b1, b2, b3, c1, c2, c3, ρ0, ρ2, and ε0 are the baseline values of the parameters βa, βh, βw, αa, αh, αw, ρa, ρw, and ε, respectively, whereas 0<a1,a2,a3,r1,r2,r3,ρ1,ρ3,ε1<1 are the strength of seasonal forcing in transmission (amplitudes of seasonal variations) for each of the seasonal parameters, and ω=π/6 corresponds to a one-year period of time.


### 2.1. Model Assumptions

The following assumptions are considered in the formulation of the brucellosis model:(i)
Mixing of individuals in each population is homogeneous(ii)
Infected animals shedBrucella in the environment(iii)
Domestic and wild animals’ seropositivity is lifelong(iv)
Immunized livestock cannot be infected unless their resistance to infection wanes(v)
The natural mortality rate in each of the species is constant(vi)
The birth rate for each population is greater than the natural mortality rateThe variables and parameter values per year incorporated in this model are summarized in Tables1 and 2, respectively.Table 1
Model variables.
VariableDescriptionShtNumber of susceptible humans at timetIhtNumber of infected humans at timetRhtNumber of recovered humans at timetSatNumber of susceptible animals at timetIatNumber of infected animals at timetVatNumber of vaccinated animals at timetBtNumber ofBrucella load per unit volume in the environment at time tTable 2
Parameters of the model and their description.
ParameterDescriptionValueSourceπaPer-capita livestock birth rate0.1[37]ϕaLivestock vaccination rate0.7[37]πhPer-capita human birth rate0.02[38]σHuman recovery rate0.25[37]μhPer-capita human natural death rate0.02[38]ψLivestock vaccine efficacy waning rate0.4[31]βaWithin-livestock transmission rate0.0011[31]dGradual culling of seropositive livestock0.35[31]μaPer-capita livestock natural mortality rate0.25[31]πwPer-capita wild animal birth rate0.08[39]βwWithin-wild animal transmission rate0.05[39]αwBrucella from B to wild animal transmission rate0.00035[3]μwPer-capita natural death rate of wild animals0.07[39]αBrucella from B to livestock transmission rate0.00035[3]αhBrucella from B to human transmission rate0.002[37]ρBrucella shedding rate of infected livestock0.5[37]ρwBrucella shedding rate of infected wild animals15[30]βhLivestock to human transmission rate0.0002[37]εDecaying rate ofBrucella in the environment8[31]τEnvironmental hygiene and sanitation rate12[3]The interactions between humans, animals, and pathogens in the environment are shown in Figure1, and the resulting model system is shown by equation (1).
(1)dVadt=ϕSa−ψ+μaVa,dSadt=πaNa+ψVa−βatIa+αatB+ϕ+μaSa,dIadt=βatIa+αatBSa−μa+dIa,dShdt=πhNh+γRh−βhtIa+βhIh+αhtB+μhSh,dIhdt=βhtIa+αhtBSh−σ+μhIh,dRhdt=σIh−γ+μhRh,dSwdt=πwNw−βwtIw+αwtB+μwSw,dIwdt=βwIw+αwtB−μwIw,dBdt=ρtIa+ρwtIw−τ+εtB.Figure 1
Flow diagram for brucellosis dynamics in animals, environment, and humans.


### 2.2. Model Properties

In this section, we use the box invariance method proposed by [40] to assess the well-posedness of the model (1) (existence and feasibility of its solution). In other words, we investigate whether the solutions of system (1) that have nonnegative initial values remain nonnegative for all times t≥0. The compact form of system (1) can be expressed as
(2)dXdt=AX+F,where X=Va,Sa,Ia,Sh,Ih,Rh,Sw,Iw,B and F is a column vector given by
(3)F=0,πcNc,0,πhNh,0,0,πwNw,0,0T,A=−ψ+μaϕ0000000ψ−λ100000000λ1−μa+dt000000000λ2+μh0γ000000λ1−σ+μh00000000σ−γ+μh000000000−λ3+μw00000000λ3−μw000ρ0000ρw−λ,where
(4)λa=βatIa+αtB+ϕ+μa,λ2=βhtIa+βhtIh+αhB,λ3=βwtIw+αwtB,λ=τ+εt.It can be noticed thatA is the Metzler matrix for all X∈ℝ+9. Therefore, based on the fact that F≥0, model (1) is positively invariant in ℝ+9. This implies that an arbitrary trajectory of the system starting in ℝ+9 forever remains in ℝ+9. In addition, F is Lipschitz continuous. Thus, a unique maximal solution exists, and so
(5)D=Va,Sa,Ia,Sh,Ih,Rh,Sw,Iw,B≥0∈ℝ+9is the feasible region for the model (1). Thus, model (1) is epidemiologically and mathematically well-posed in the region D.


### 2.3. Brucellosis-Free Equilibrium

The brucellosis-free equilibrium solution for system (1) is computed and found to be
(6)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0=ϕπaNaμaϕ+ψ+μa,ψ+μaπaNaμaϕ+ψ+μa,0,πhNhμh,0,0,πwNwμw,0,0,where Na, Nh, and Nw are, respectively, the initial total populations of the livestock, humans, and wild animals.


### 2.4. The Reproduction Number

A heterogeneous population with individuals which can be grouped inton homogeneous compartments is considered in this section. Let x=x1,⋯,xnT, with xi≥0, be the state of individuals in each compartment. It is assumed that the compartments can be divided into the following: infected designated as i=1,⋯,m and uninfected designated as i=m+1,⋯,n. We also define Xs to be the set of all disease-free states:
(7)Xs=x≥0:xi=0,∀i=1,⋯,m.LetFit,x be the input rate of newly infected individuals in the ith compartment, Vi+t,x be the input rate of individuals by other means (for example, births and immigrations), and Vi−t,x be the rate of transfer of individuals out of compartment i (for example, deaths, recovery, and emigrations). Henceforth, the disease transmission model is governed by a nonautonomous ordinary differential system:
(8)dxidt=Fit,x−Vit,x≜fit,x,i=1,⋯,n,where Vit,x=Vi−t,x−Vi+t,x.Succeeding the approach by [41] and that of [42] for epidemic models, we look at conditions (A1)–(A7) for the brucellosis model. The model (1) is equivalent to periodic ordinary differential system (8), we can easily see that conditions (A1)–(A5) stated below are satisfied.(A1) For each1≤i≤n, the functions Fit,x, Vi+t,x, and Vi−t,x are nonnegative and continuous on ℝ×ℝ+n and continuously differential with respect to x. This is based on the fact that each function denotes a directed nonnegative transfer of individuals(A2) There is a real numberω>0 such that for each 1≤i≤n, the functions Fit,x, Vi+t,x, and Vi−t,x are ω-periodic in t. This biologically describes a periodic environment due to seasonality(A3) Ifxi=0, then Vi−t,x=0. In particular, if x∈Xs, then Vi−t,x=0 for i=1,⋯,m. That is, if a compartment is empty, then there is no transfer of individuals out of it(A4)Fi=0 for i>m. This means that the infection incidence for uninfected compartments is zero(A5) Ifx∈Xs, then Fi=Vi+=0 for i=1,⋯,m. This implies that if the population is disease-free in the beginning, it will remain soWe know that model (8) has a disease-free periodic solution, so we define a 5×5 matrix for the nontransmitting compartments as
(9)Mt=−ψ+μaϕ000ψ−ϕ+μa00000−μhγ0000−γ+μh00000−μw.LetΦMt be the monodromy matrix of the linear ω-periodic system dz/dt=Mtz. Then, ρΦMω<1 implying that E0t is linearly asymptotically stable in the disease-free subspace Xs; that is,(A6)ρΦMω<1, where ρΦMω is the spectral radius of ΦMω, is satisfiedFor convenience purposes and easy presentation of the results, we letC denote all continuous functions on the real line. If f is a periodic function in C, then we use f¯ for the average value of the time interval 0T defined by
(10)f¯=1T∫0Tftdt,for continuous T periodic function ft. Inspired by the approach of [41, 43], we obtain
(11)F=ψ+μaβ¯atπaNaμaϕ+ψ+μa00ψ+μaα¯atπaNaμaϕ+ψ+μaβ¯htπhNhμhβ¯htπhNhμh0α¯htπhNhμh00β¯wtπwNwμwα¯wtπwNwμwρ¯at0ρ¯wt0,(12)V=μa+d0000σ+μh0000μω0000τ+ε¯tand observe that F is nonnegative and −V is cooperative because its off-diagonal elements are nonnegative.It follows that the effective reproductive number of the time-averaged autonomous system is(13)Re=R11+R33+R11−R332+4R13R312,where
(14)R11=β¯atτ+ε¯t+α¯atρ¯tψ+μaπaNaμaϕ+ψ+μaμa+dτ+ε¯t,R33=β¯wtτ+ε¯t+α¯wtρ¯wtψ+μaπwNwμw2τ+ε¯t,R13=α¯atρ¯wtψ+μaπaNaμaμwτ+ε¯tϕ+ψ+μa,R31=α¯wtρtπwNwμwμa+dτ+ε¯t.Generally, the time-averaged effective reproduction number is computed as the dominant eigenvalue ofFV−1 using the Maple package and is found to be
(15)ρFV−1=Re=1T∫0TR11+R33+R11−R332+4R13R312ds.If no interventions are administered, the time-averaged basic reproductive number for model system (1) is found to be
(16)R0=1T∫0TR110+R330+R110−R3302+4R130R3102ds,where
(17)R11=β¯atε¯t+α¯atρ¯tπaNaμa2ε¯t,R33=β¯wtε¯t+α¯wtρ¯wtμaπwNwμw2ε¯t,R13=α¯atρ¯wtπaNaμaμwε¯t,R31=α¯wtρtπwNwμwμaε¯t.R0 may be interpreted as the average number of secondary cases arising from the introduction of a single infected person into a completely susceptible population at a random time of the year. The condition R0<1 is sufficient and necessary for long-term disease extinction. Furthermore, let Yt,s,t≥s, be the evolution operator of the linear ω-periodic system:
(18)dydt=−Vty.That is, for eachs∈ℝ, the 4×4 matrix Yt,s satisfies
(19)ddtYt,s=−VtYt,s,∀t≥s,Ys,s=I,where I is a 4×4 identity matrix. Therefore, the monodromy matrix ΦVt of (18) equals Yt,0,t≥0. Thus, condition (A7) below is satisfied.(A7) The internal evolution of individuals in the infectious compartments due to deaths and movements is dissipative and decays exponentially in many cases. This is because of loss of infective members from natural and disease-induced mortality. Thus,ρΦVω<1Based on the assumptions (A1)–(A7), we are now able to analyze the reproduction ratios for the epidemic model system (1). For this purpose, we always assume that the population is near the disease-free periodic state E0t. By the standard theory of linear periodic systems [44], there exist K>0 and α>0 such that
(20)Yt,s≤Ke−αt−s,∀t≥s,s∈ℝ.Consequently,(21)Yt,t−aFt−a≤KFt−ae−αa,∀t∈ℝ,a∈0,∞.In the computation of the basic reproduction number for the nonautonomous model system (1), we follow the method by [42]. Suppose Γs is the initial distribution of infectious individuals in this periodic environment; then, FsΓs is the rate of new infectious individuals produced by the infected individuals who were introduced at time s. Yt,sFsΓs represents the distribution of the newly infected at time s and remains in the infected compartment at time t≥s. It follows that the cumulative distribution of new infections at t produced by all infected Γt individuals introduced prior to t=s is given by
(22)Ψt=∫−∞tYt,sFsΓsds=∫0∞Yt,t−aFt,t−aΓt−ada,∀t∈ℝ,Γ∈Cω.LetCω be the ordered Banach space of all ω-periodic functions from ℝ to ℝn, which is equipped with the maximum norm, .∞, and the positive cone Cω+=Γ∈CωΓt≥0,t∈ℝ. We define the linear operator L:CωCω by
(23)LΓt=∫0∞Yt,t−aFt,t−aΓt−ada,∀t∈ℝ,Γ∈Cω,where L is the next infection operator. Then, the basic reproduction number is given by
(24)Rω=ρL,where ρL is the spectral radius of L. By direct calculation, the evolution operator Yt,s for the system (1) is found to be
(25)Yt,s=e−μa+dt−s0000e−σ+μht−s0000e−μωt−s0000Y¯t,s,with
(26)Y¯t,s=e−τ+ε0t−s+6ε0ε1/πcosπt/6−cosπs/6.Motivated by [45], the next infection operator can be numerically evaluated as
(27)Lφt=∫0∞Yt,t−aFt,t−aΓt−ada=∫0ωGt,aΓt−ada,where
(28)Gt,s≈∑k=0MYt,t−s−kωFt−s≈∑k=0Mm1100m14m21m220m2400m33m33m410m430,for positive integers M which are large enough, and
(29)m11=βat−sψ+μaπaNaμaϕ+ψ+μae−μa+dt−s,m14=αat−sψ+μaπaNaμaϕ+ψ+μae−μa+dt−s,m21=βht−sπhNhμhe−σ+μat−s,m22=βht−sπhNhμhe−σ+μat−s,m24=αht−sπhNhμhe−σ+μat−s,m33=βwt−sπhNhμhe−ωt−s,m34=αwt−sπwNwμwe−ωt−s,m41=ρt−se−τ+ε0t−s+6ε0ε1/πcosπt6−cosπs6,m43=ρwt−se−τ+ε0t−s+6ε0ε1/πcosπt6−cosπs6.


### 2.5. Global Stability of the Brucellosis-Free Solution

In this section, we establish the conditions for global stability of a disease-free periodic solution.Theorem 1.
The disease-free solution of system (1) is globally asymptotically stable if the basic reproduction number in D is less than one.Proof.
Consider the matrix function:(30)Ft−Vt=βatSa0−μa+d00αatSa0βhtπhNh0μhβhtSh0−σ+μh0αhtπhNh0μh00βwtSw0−μwαwtπwNwμwρt0ρwt−τ+εt.
We verify that matrix function (30) is continuous, cooperative, irreducible, and ω-periodic. Let ΦF−V.t be the fundamental solution matrix of the linear ordinary differential system:
(31)ẋ=Ft−Vtx,and ρΦF−V.ω be the dominant eigenvalue of ΦF−V.ω. From Theorem 2.2 in [42], we have R0<1 if and only if ρΦF−V.ω<1.Lemma 2.
Letv=1/ωlnρΦF−V.ω. Then, there exists a positive ω-periodic function vt such that evtvt is a solution to equation (31).
From the nondisease transmitting equations of system (1), we obtain the following:
(32)Vat≤ϕπaNaμaϕ+ψ+μa≜Va0,Sat≤ψ+μaπaNaμaϕ+ψ+μa≜Sa0,Sht≤πhNhμh≜Sh0,Sωt≤πwNwμw≜Sw0.
Again, from the infectious and recovered classes of system (1), we have the following:
(33)ddtIatIhtIwtBt≤F−VIatIhtIwtBt.Based on Lemma2, there exists vt such that xt=I¯at,I¯ht,I¯wt,B¯t=vtevt is a solution to equation (31) with v=1/ωlnρΦF−V..Based on the fact thatR0<1, we have ρΦF−V.<1 and v<0. Thus,
(34)Iat,Iht,Iwt,Bt≤I¯at,I¯ht,I¯wt,B¯t,when t is very large which would imply that
(35)limt→∞Iat=limt→∞Iht=limt→∞Iwt=limt→∞Bt=0.Moreover, ast⟶∞, we have
(36)ddtVa+Sa⟶πaNa−μaVa+Sa,which implies
(37)dVadt⟶ϕπaNaμa−Va−ψ+μaVa=ψπaNaμa−ϕ+ψ+μaVa,(38)orϕπaNaμaϕ+ψ+μa=Va0,which leads to
(39)Sat⟶πaNaμa−Va0=ψ+μaπaNaμaπ+ψ+μa=Sa0.Again,(40)dShdt⟶πhNh−μhSh,dSwdt⟶πwNw−μwSw,which gives
(41)Sh0=πhNhμh,Sw0=πwNwμw.Therefore,(42)limt→∞xt=Va0,Sa0,0,Sh0,0,0,Sw0,0,0,for each solution xt in system (1).


### 2.6. Endemic Equilibrium Solution

This section is aimed at investigating the behavior of model system (1) when R0>1. We show that if R0>1, brucellosis infection persists in the animal and human populations and there exists a positive periodic solution. Following the approach in [46, 47], we define
(43)X=ℝ+9;X0=ℝ+4×Intℝ+5;∂X0=X\X0.LetL:X⟶X be the Poncaré map associated with model system (1) such that Px0=uω,x0∀x0∈X, where ut,x0 denotes a unique solution of the system with u0,x0=x0.Definition 3.
The solutions of the model system (1) are said to be uniformly persistent if there exists some ξ>0 such that
(44)limt→∞InfVat>ξ,limt→∞InfSat>ξ,limt→∞InfIat>ξ,limt→∞InfSht>ξ,limt→∞InfIht>ξ,limt→∞InfRht>ξ,limt→∞InfSwt>ξ,limt→∞InfIwt>ξ,limt→∞InfBt>ξ,whenever
(45)Va0>0,Sa0>0,Ia0>0,Sh0>0,Ih0>0,Rh0>0,Sw0>0,Iw0>0,B0>0.Theorem 4.
The solutions of the model system (1) are uniformly persistent, and the system admits at least one positive ω-periodic solution if R0>1.Proof.
We define(46)H∂=Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0:PmVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0,∀m≥0,H~=Va0,Sa0,0,Sh0,0,0,Sw0,0,0:Va0≥0,Sa0≥0,Sh0≥0,Sw0≥0.
It is evident thatH~∂⊆H∂.
We first show thatH∂=H~∂. Consider the initial values:
(47)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0H~.
IfIa0=0, Ih0, Iw0=0, and B0>0, then based on the fact that there is a recruitment rate for susceptible individuals, we have Ia′>0. Similarly, if Iw0=0, Ih0, B0=0, and Ia0>0, then B′0>0, Ia0=0, Ih0, Iw0=0, and B0>0, and if Ia0=0, Ih0, B0=0, and Iw0>0, then B′0>0. It follows that Vat,Sat,Iat,Sht,Iht,Rht,Swt,Iwt,Bt∉∂X0 for 0<t≪1. The positive invariance of X0 implies that H∂=H~∂.
Again, if we consider the fixed point:(48)H0=ϕπaNaμaϕ+ψ+μa,ψ+μaπaNaμaϕ+ψ+μa,0,πhNhμh,0,πwNwμw,0,0,we define
(49)WSH0=x0:Lmx0⟶H0,x⟶∞.
It can be deduced from system (1) that if Ia=Ih=Iw=B=0 and t⟶∞,
(50)Vat⟶Va0=ϕπaNaμaϕ+ψ+μa,Sat⟶Sa0ψ+μaπaNaμaϕ+ψ+μa,Sht⟶Sh0=πhNhμh,Swt⟶Sw0=πwNwμw.
We prove thatWSH0∩X0=∅.
Let. denote a norm on ℝ+9. Based on the continuity of solutions with respect to the initial conditions, for every ε>0, there exists δ>0 but small such that for all
(51)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0,with
(52)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0≤δ,we have
(53)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0≤ε,∀t∈0,ω.
So we claim that(54)limt→∞supVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0≥δ,∀Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈X0and prove by contradiction as follows:
Suppose(55)limt→∞supVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0<δ,for some
(56)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈X0.
In addition, we assume without loss of generality that(57)PmVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0<δ,∀m≥0.
Therefore,∀t∈0,ω,m≥0, we have
(58)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0≤ε.
Furthermore, for any nonnegativet, we can write t=t0+nω with t0∈0,ω and n being the greatest integer less than or equal to t/ω. Then, we get
(59)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0=ut0,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut0,H0≤ε,for any t>0.
Let(60)Vat,Sat,Iat,Sht,Iht,Rht,Swt,Iwt,Bt=Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0.
It follows that(61)ϕπaNaμaϕ+ψ+μa−ε<Vat<ϕπaNaμaϕ+ψ+μa+ε,ψ+μaπaNaμaϕ+ψ+μa−ε<Sat<ψ+μaπaNaμaϕ+ψ+μa+ε,πhNhμh−ε<Sht<πhNhμh,πwNwμw−ε<Swt<πwNwμw,0<Iat<ε,0<Iht<ε,0<Iwt<ε,0<Bt<ε.
Then, we have(62)dIadt=β1tIa+α1tBSa−μa+dIa≥β1tIa+α1tBψ+μaπaNaμaϕ+ψ+μa−ε−μa+dIa=β1tIa+α1tBψ+μaπaNaμaϕ+ψ+μa−μa+dIa−εβ1tIa+α1tB.
Similarly,(63)dIhdt≥β2tIh+α2tBπhNhμh−σ+μhIh−εβ2tIh+α2tB,dIwdt≥βwtIw+αwtBπwNwμw−μwIw−εβwtIw+αwtB.
Thus, we obtain(64)ddtIatIhtIwtBt≥F−V−εKIatIhtIwtBt.
ButR0>1 if and only if ρΦF−V.>1. Thus, for ε>0 whenever small, we have ρΦF−V.>1. Using Lemma 2 and the comparison principle, we get
(65)limt→∞Iat=limt→∞Iht=limt→∞Iwt=limt→∞Bt=∞,which contradicts our original assumption.
Thus,H0 is acyclic in H∂, and P is uniformly persistent with respect to X0,∂X0, which implies the uniform persistence of the solutions to the original system [47]. Consequently, the Poincaré map p has a fixed point:
(66)V¯a0,S¯a0,I¯a0,S¯h0,I¯h0,S¯w0,I¯w0,B¯0∈X0,with Va0,Sa0,Sh0,Sw0≠0. Thus,
(67)V¯a0,S¯a0,I¯a0,S¯h0,I¯h0,S¯w0,I¯w0,B¯0∈Intℝ+9,and
(68)V~a0,S~a0,I~a0,S~h0,I~h0,S~w0,I¯w0,B~0=ut,V~a0,S~a0,I~a0,S~h0,I~h0,S~w0,I¯w0,B~0is a positive ω-periodic solution of the system.


## 2.1. Model Assumptions

The following assumptions are considered in the formulation of the brucellosis model:(i)
Mixing of individuals in each population is homogeneous(ii)
Infected animals shedBrucella in the environment(iii)
Domestic and wild animals’ seropositivity is lifelong(iv)
Immunized livestock cannot be infected unless their resistance to infection wanes(v)
The natural mortality rate in each of the species is constant(vi)
The birth rate for each population is greater than the natural mortality rateThe variables and parameter values per year incorporated in this model are summarized in Tables1 and 2, respectively.Table 1
Model variables.
VariableDescriptionShtNumber of susceptible humans at timetIhtNumber of infected humans at timetRhtNumber of recovered humans at timetSatNumber of susceptible animals at timetIatNumber of infected animals at timetVatNumber of vaccinated animals at timetBtNumber ofBrucella load per unit volume in the environment at time tTable 2
Parameters of the model and their description.
ParameterDescriptionValueSourceπaPer-capita livestock birth rate0.1[37]ϕaLivestock vaccination rate0.7[37]πhPer-capita human birth rate0.02[38]σHuman recovery rate0.25[37]μhPer-capita human natural death rate0.02[38]ψLivestock vaccine efficacy waning rate0.4[31]βaWithin-livestock transmission rate0.0011[31]dGradual culling of seropositive livestock0.35[31]μaPer-capita livestock natural mortality rate0.25[31]πwPer-capita wild animal birth rate0.08[39]βwWithin-wild animal transmission rate0.05[39]αwBrucella from B to wild animal transmission rate0.00035[3]μwPer-capita natural death rate of wild animals0.07[39]αBrucella from B to livestock transmission rate0.00035[3]αhBrucella from B to human transmission rate0.002[37]ρBrucella shedding rate of infected livestock0.5[37]ρwBrucella shedding rate of infected wild animals15[30]βhLivestock to human transmission rate0.0002[37]εDecaying rate ofBrucella in the environment8[31]τEnvironmental hygiene and sanitation rate12[3]The interactions between humans, animals, and pathogens in the environment are shown in Figure1, and the resulting model system is shown by equation (1).
(1)dVadt=ϕSa−ψ+μaVa,dSadt=πaNa+ψVa−βatIa+αatB+ϕ+μaSa,dIadt=βatIa+αatBSa−μa+dIa,dShdt=πhNh+γRh−βhtIa+βhIh+αhtB+μhSh,dIhdt=βhtIa+αhtBSh−σ+μhIh,dRhdt=σIh−γ+μhRh,dSwdt=πwNw−βwtIw+αwtB+μwSw,dIwdt=βwIw+αwtB−μwIw,dBdt=ρtIa+ρwtIw−τ+εtB.Figure 1
Flow diagram for brucellosis dynamics in animals, environment, and humans.


## 2.2. Model Properties

In this section, we use the box invariance method proposed by [40] to assess the well-posedness of the model (1) (existence and feasibility of its solution). In other words, we investigate whether the solutions of system (1) that have nonnegative initial values remain nonnegative for all times t≥0. The compact form of system (1) can be expressed as
(2)dXdt=AX+F,where X=Va,Sa,Ia,Sh,Ih,Rh,Sw,Iw,B and F is a column vector given by
(3)F=0,πcNc,0,πhNh,0,0,πwNw,0,0T,A=−ψ+μaϕ0000000ψ−λ100000000λ1−μa+dt000000000λ2+μh0γ000000λ1−σ+μh00000000σ−γ+μh000000000−λ3+μw00000000λ3−μw000ρ0000ρw−λ,where
(4)λa=βatIa+αtB+ϕ+μa,λ2=βhtIa+βhtIh+αhB,λ3=βwtIw+αwtB,λ=τ+εt.It can be noticed thatA is the Metzler matrix for all X∈ℝ+9. Therefore, based on the fact that F≥0, model (1) is positively invariant in ℝ+9. This implies that an arbitrary trajectory of the system starting in ℝ+9 forever remains in ℝ+9. In addition, F is Lipschitz continuous. Thus, a unique maximal solution exists, and so
(5)D=Va,Sa,Ia,Sh,Ih,Rh,Sw,Iw,B≥0∈ℝ+9is the feasible region for the model (1). Thus, model (1) is epidemiologically and mathematically well-posed in the region D.


## 2.3. Brucellosis-Free Equilibrium

The brucellosis-free equilibrium solution for system (1) is computed and found to be
(6)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0=ϕπaNaμaϕ+ψ+μa,ψ+μaπaNaμaϕ+ψ+μa,0,πhNhμh,0,0,πwNwμw,0,0,where Na, Nh, and Nw are, respectively, the initial total populations of the livestock, humans, and wild animals.


## 2.4. The Reproduction Number

A heterogeneous population with individuals which can be grouped inton homogeneous compartments is considered in this section. Let x=x1,⋯,xnT, with xi≥0, be the state of individuals in each compartment. It is assumed that the compartments can be divided into the following: infected designated as i=1,⋯,m and uninfected designated as i=m+1,⋯,n. We also define Xs to be the set of all disease-free states:
(7)Xs=x≥0:xi=0,∀i=1,⋯,m.LetFit,x be the input rate of newly infected individuals in the ith compartment, Vi+t,x be the input rate of individuals by other means (for example, births and immigrations), and Vi−t,x be the rate of transfer of individuals out of compartment i (for example, deaths, recovery, and emigrations). Henceforth, the disease transmission model is governed by a nonautonomous ordinary differential system:
(8)dxidt=Fit,x−Vit,x≜fit,x,i=1,⋯,n,where Vit,x=Vi−t,x−Vi+t,x.Succeeding the approach by [41] and that of [42] for epidemic models, we look at conditions (A1)–(A7) for the brucellosis model. The model (1) is equivalent to periodic ordinary differential system (8), we can easily see that conditions (A1)–(A5) stated below are satisfied.(A1) For each1≤i≤n, the functions Fit,x, Vi+t,x, and Vi−t,x are nonnegative and continuous on ℝ×ℝ+n and continuously differential with respect to x. This is based on the fact that each function denotes a directed nonnegative transfer of individuals(A2) There is a real numberω>0 such that for each 1≤i≤n, the functions Fit,x, Vi+t,x, and Vi−t,x are ω-periodic in t. This biologically describes a periodic environment due to seasonality(A3) Ifxi=0, then Vi−t,x=0. In particular, if x∈Xs, then Vi−t,x=0 for i=1,⋯,m. That is, if a compartment is empty, then there is no transfer of individuals out of it(A4)Fi=0 for i>m. This means that the infection incidence for uninfected compartments is zero(A5) Ifx∈Xs, then Fi=Vi+=0 for i=1,⋯,m. This implies that if the population is disease-free in the beginning, it will remain soWe know that model (8) has a disease-free periodic solution, so we define a 5×5 matrix for the nontransmitting compartments as
(9)Mt=−ψ+μaϕ000ψ−ϕ+μa00000−μhγ0000−γ+μh00000−μw.LetΦMt be the monodromy matrix of the linear ω-periodic system dz/dt=Mtz. Then, ρΦMω<1 implying that E0t is linearly asymptotically stable in the disease-free subspace Xs; that is,(A6)ρΦMω<1, where ρΦMω is the spectral radius of ΦMω, is satisfiedFor convenience purposes and easy presentation of the results, we letC denote all continuous functions on the real line. If f is a periodic function in C, then we use f¯ for the average value of the time interval 0T defined by
(10)f¯=1T∫0Tftdt,for continuous T periodic function ft. Inspired by the approach of [41, 43], we obtain
(11)F=ψ+μaβ¯atπaNaμaϕ+ψ+μa00ψ+μaα¯atπaNaμaϕ+ψ+μaβ¯htπhNhμhβ¯htπhNhμh0α¯htπhNhμh00β¯wtπwNwμwα¯wtπwNwμwρ¯at0ρ¯wt0,(12)V=μa+d0000σ+μh0000μω0000τ+ε¯tand observe that F is nonnegative and −V is cooperative because its off-diagonal elements are nonnegative.It follows that the effective reproductive number of the time-averaged autonomous system is(13)Re=R11+R33+R11−R332+4R13R312,where
(14)R11=β¯atτ+ε¯t+α¯atρ¯tψ+μaπaNaμaϕ+ψ+μaμa+dτ+ε¯t,R33=β¯wtτ+ε¯t+α¯wtρ¯wtψ+μaπwNwμw2τ+ε¯t,R13=α¯atρ¯wtψ+μaπaNaμaμwτ+ε¯tϕ+ψ+μa,R31=α¯wtρtπwNwμwμa+dτ+ε¯t.Generally, the time-averaged effective reproduction number is computed as the dominant eigenvalue ofFV−1 using the Maple package and is found to be
(15)ρFV−1=Re=1T∫0TR11+R33+R11−R332+4R13R312ds.If no interventions are administered, the time-averaged basic reproductive number for model system (1) is found to be
(16)R0=1T∫0TR110+R330+R110−R3302+4R130R3102ds,where
(17)R11=β¯atε¯t+α¯atρ¯tπaNaμa2ε¯t,R33=β¯wtε¯t+α¯wtρ¯wtμaπwNwμw2ε¯t,R13=α¯atρ¯wtπaNaμaμwε¯t,R31=α¯wtρtπwNwμwμaε¯t.R0 may be interpreted as the average number of secondary cases arising from the introduction of a single infected person into a completely susceptible population at a random time of the year. The condition R0<1 is sufficient and necessary for long-term disease extinction. Furthermore, let Yt,s,t≥s, be the evolution operator of the linear ω-periodic system:
(18)dydt=−Vty.That is, for eachs∈ℝ, the 4×4 matrix Yt,s satisfies
(19)ddtYt,s=−VtYt,s,∀t≥s,Ys,s=I,where I is a 4×4 identity matrix. Therefore, the monodromy matrix ΦVt of (18) equals Yt,0,t≥0. Thus, condition (A7) below is satisfied.(A7) The internal evolution of individuals in the infectious compartments due to deaths and movements is dissipative and decays exponentially in many cases. This is because of loss of infective members from natural and disease-induced mortality. Thus,ρΦVω<1Based on the assumptions (A1)–(A7), we are now able to analyze the reproduction ratios for the epidemic model system (1). For this purpose, we always assume that the population is near the disease-free periodic state E0t. By the standard theory of linear periodic systems [44], there exist K>0 and α>0 such that
(20)Yt,s≤Ke−αt−s,∀t≥s,s∈ℝ.Consequently,(21)Yt,t−aFt−a≤KFt−ae−αa,∀t∈ℝ,a∈0,∞.In the computation of the basic reproduction number for the nonautonomous model system (1), we follow the method by [42]. Suppose Γs is the initial distribution of infectious individuals in this periodic environment; then, FsΓs is the rate of new infectious individuals produced by the infected individuals who were introduced at time s. Yt,sFsΓs represents the distribution of the newly infected at time s and remains in the infected compartment at time t≥s. It follows that the cumulative distribution of new infections at t produced by all infected Γt individuals introduced prior to t=s is given by
(22)Ψt=∫−∞tYt,sFsΓsds=∫0∞Yt,t−aFt,t−aΓt−ada,∀t∈ℝ,Γ∈Cω.LetCω be the ordered Banach space of all ω-periodic functions from ℝ to ℝn, which is equipped with the maximum norm, .∞, and the positive cone Cω+=Γ∈CωΓt≥0,t∈ℝ. We define the linear operator L:CωCω by
(23)LΓt=∫0∞Yt,t−aFt,t−aΓt−ada,∀t∈ℝ,Γ∈Cω,where L is the next infection operator. Then, the basic reproduction number is given by
(24)Rω=ρL,where ρL is the spectral radius of L. By direct calculation, the evolution operator Yt,s for the system (1) is found to be
(25)Yt,s=e−μa+dt−s0000e−σ+μht−s0000e−μωt−s0000Y¯t,s,with
(26)Y¯t,s=e−τ+ε0t−s+6ε0ε1/πcosπt/6−cosπs/6.Motivated by [45], the next infection operator can be numerically evaluated as
(27)Lφt=∫0∞Yt,t−aFt,t−aΓt−ada=∫0ωGt,aΓt−ada,where
(28)Gt,s≈∑k=0MYt,t−s−kωFt−s≈∑k=0Mm1100m14m21m220m2400m33m33m410m430,for positive integers M which are large enough, and
(29)m11=βat−sψ+μaπaNaμaϕ+ψ+μae−μa+dt−s,m14=αat−sψ+μaπaNaμaϕ+ψ+μae−μa+dt−s,m21=βht−sπhNhμhe−σ+μat−s,m22=βht−sπhNhμhe−σ+μat−s,m24=αht−sπhNhμhe−σ+μat−s,m33=βwt−sπhNhμhe−ωt−s,m34=αwt−sπwNwμwe−ωt−s,m41=ρt−se−τ+ε0t−s+6ε0ε1/πcosπt6−cosπs6,m43=ρwt−se−τ+ε0t−s+6ε0ε1/πcosπt6−cosπs6.


## 2.5. Global Stability of the Brucellosis-Free Solution

In this section, we establish the conditions for global stability of a disease-free periodic solution.Theorem 1.
The disease-free solution of system (1) is globally asymptotically stable if the basic reproduction number in D is less than one.Proof.
Consider the matrix function:(30)Ft−Vt=βatSa0−μa+d00αatSa0βhtπhNh0μhβhtSh0−σ+μh0αhtπhNh0μh00βwtSw0−μwαwtπwNwμwρt0ρwt−τ+εt.
We verify that matrix function (30) is continuous, cooperative, irreducible, and ω-periodic. Let ΦF−V.t be the fundamental solution matrix of the linear ordinary differential system:
(31)ẋ=Ft−Vtx,and ρΦF−V.ω be the dominant eigenvalue of ΦF−V.ω. From Theorem 2.2 in [42], we have R0<1 if and only if ρΦF−V.ω<1.Lemma 2.
Letv=1/ωlnρΦF−V.ω. Then, there exists a positive ω-periodic function vt such that evtvt is a solution to equation (31).
From the nondisease transmitting equations of system (1), we obtain the following:
(32)Vat≤ϕπaNaμaϕ+ψ+μa≜Va0,Sat≤ψ+μaπaNaμaϕ+ψ+μa≜Sa0,Sht≤πhNhμh≜Sh0,Sωt≤πwNwμw≜Sw0.
Again, from the infectious and recovered classes of system (1), we have the following:
(33)ddtIatIhtIwtBt≤F−VIatIhtIwtBt.Based on Lemma2, there exists vt such that xt=I¯at,I¯ht,I¯wt,B¯t=vtevt is a solution to equation (31) with v=1/ωlnρΦF−V..Based on the fact thatR0<1, we have ρΦF−V.<1 and v<0. Thus,
(34)Iat,Iht,Iwt,Bt≤I¯at,I¯ht,I¯wt,B¯t,when t is very large which would imply that
(35)limt→∞Iat=limt→∞Iht=limt→∞Iwt=limt→∞Bt=0.Moreover, ast⟶∞, we have
(36)ddtVa+Sa⟶πaNa−μaVa+Sa,which implies
(37)dVadt⟶ϕπaNaμa−Va−ψ+μaVa=ψπaNaμa−ϕ+ψ+μaVa,(38)orϕπaNaμaϕ+ψ+μa=Va0,which leads to
(39)Sat⟶πaNaμa−Va0=ψ+μaπaNaμaπ+ψ+μa=Sa0.Again,(40)dShdt⟶πhNh−μhSh,dSwdt⟶πwNw−μwSw,which gives
(41)Sh0=πhNhμh,Sw0=πwNwμw.Therefore,(42)limt→∞xt=Va0,Sa0,0,Sh0,0,0,Sw0,0,0,for each solution xt in system (1).


## 2.6. Endemic Equilibrium Solution

This section is aimed at investigating the behavior of model system (1) when R0>1. We show that if R0>1, brucellosis infection persists in the animal and human populations and there exists a positive periodic solution. Following the approach in [46, 47], we define
(43)X=ℝ+9;X0=ℝ+4×Intℝ+5;∂X0=X\X0.LetL:X⟶X be the Poncaré map associated with model system (1) such that Px0=uω,x0∀x0∈X, where ut,x0 denotes a unique solution of the system with u0,x0=x0.Definition 3.
The solutions of the model system (1) are said to be uniformly persistent if there exists some ξ>0 such that
(44)limt→∞InfVat>ξ,limt→∞InfSat>ξ,limt→∞InfIat>ξ,limt→∞InfSht>ξ,limt→∞InfIht>ξ,limt→∞InfRht>ξ,limt→∞InfSwt>ξ,limt→∞InfIwt>ξ,limt→∞InfBt>ξ,whenever
(45)Va0>0,Sa0>0,Ia0>0,Sh0>0,Ih0>0,Rh0>0,Sw0>0,Iw0>0,B0>0.Theorem 4.
The solutions of the model system (1) are uniformly persistent, and the system admits at least one positive ω-periodic solution if R0>1.Proof.
We define(46)H∂=Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0:PmVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0,∀m≥0,H~=Va0,Sa0,0,Sh0,0,0,Sw0,0,0:Va0≥0,Sa0≥0,Sh0≥0,Sw0≥0.
It is evident thatH~∂⊆H∂.
We first show thatH∂=H~∂. Consider the initial values:
(47)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0H~.
IfIa0=0, Ih0, Iw0=0, and B0>0, then based on the fact that there is a recruitment rate for susceptible individuals, we have Ia′>0. Similarly, if Iw0=0, Ih0, B0=0, and Ia0>0, then B′0>0, Ia0=0, Ih0, Iw0=0, and B0>0, and if Ia0=0, Ih0, B0=0, and Iw0>0, then B′0>0. It follows that Vat,Sat,Iat,Sht,Iht,Rht,Swt,Iwt,Bt∉∂X0 for 0<t≪1. The positive invariance of X0 implies that H∂=H~∂.
Again, if we consider the fixed point:(48)H0=ϕπaNaμaϕ+ψ+μa,ψ+μaπaNaμaϕ+ψ+μa,0,πhNhμh,0,πwNwμw,0,0,we define
(49)WSH0=x0:Lmx0⟶H0,x⟶∞.
It can be deduced from system (1) that if Ia=Ih=Iw=B=0 and t⟶∞,
(50)Vat⟶Va0=ϕπaNaμaϕ+ψ+μa,Sat⟶Sa0ψ+μaπaNaμaϕ+ψ+μa,Sht⟶Sh0=πhNhμh,Swt⟶Sw0=πwNwμw.
We prove thatWSH0∩X0=∅.
Let. denote a norm on ℝ+9. Based on the continuity of solutions with respect to the initial conditions, for every ε>0, there exists δ>0 but small such that for all
(51)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈∂X0,with
(52)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0≤δ,we have
(53)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0≤ε,∀t∈0,ω.
So we claim that(54)limt→∞supVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0≥δ,∀Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈X0and prove by contradiction as follows:
Suppose(55)limt→∞supVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0<δ,for some
(56)Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0∈X0.
In addition, we assume without loss of generality that(57)PmVa0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−H0<δ,∀m≥0.
Therefore,∀t∈0,ω,m≥0, we have
(58)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0≤ε.
Furthermore, for any nonnegativet, we can write t=t0+nω with t0∈0,ω and n being the greatest integer less than or equal to t/ω. Then, we get
(59)ut,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut,H0=ut0,Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0−ut0,H0≤ε,for any t>0.
Let(60)Vat,Sat,Iat,Sht,Iht,Rht,Swt,Iwt,Bt=Va0,Sa0,Ia0,Sh0,Ih0,Rh0,Sw0,Iw0,B0.
It follows that(61)ϕπaNaμaϕ+ψ+μa−ε<Vat<ϕπaNaμaϕ+ψ+μa+ε,ψ+μaπaNaμaϕ+ψ+μa−ε<Sat<ψ+μaπaNaμaϕ+ψ+μa+ε,πhNhμh−ε<Sht<πhNhμh,πwNwμw−ε<Swt<πwNwμw,0<Iat<ε,0<Iht<ε,0<Iwt<ε,0<Bt<ε.
Then, we have(62)dIadt=β1tIa+α1tBSa−μa+dIa≥β1tIa+α1tBψ+μaπaNaμaϕ+ψ+μa−ε−μa+dIa=β1tIa+α1tBψ+μaπaNaμaϕ+ψ+μa−μa+dIa−εβ1tIa+α1tB.
Similarly,(63)dIhdt≥β2tIh+α2tBπhNhμh−σ+μhIh−εβ2tIh+α2tB,dIwdt≥βwtIw+αwtBπwNwμw−μwIw−εβwtIw+αwtB.
Thus, we obtain(64)ddtIatIhtIwtBt≥F−V−εKIatIhtIwtBt.
ButR0>1 if and only if ρΦF−V.>1. Thus, for ε>0 whenever small, we have ρΦF−V.>1. Using Lemma 2 and the comparison principle, we get
(65)limt→∞Iat=limt→∞Iht=limt→∞Iwt=limt→∞Bt=∞,which contradicts our original assumption.
Thus,H0 is acyclic in H∂, and P is uniformly persistent with respect to X0,∂X0, which implies the uniform persistence of the solutions to the original system [47]. Consequently, the Poincaré map p has a fixed point:
(66)V¯a0,S¯a0,I¯a0,S¯h0,I¯h0,S¯w0,I¯w0,B¯0∈X0,with Va0,Sa0,Sh0,Sw0≠0. Thus,
(67)V¯a0,S¯a0,I¯a0,S¯h0,I¯h0,S¯w0,I¯w0,B¯0∈Intℝ+9,and
(68)V~a0,S~a0,I~a0,S~h0,I~h0,S~w0,I¯w0,B~0=ut,V~a0,S~a0,I~a0,S~h0,I~h0,S~w0,I¯w0,B~0is a positive ω-periodic solution of the system.


## 3. Numerical Simulations

In this part, we perform numerical simulations for model system (1) for the purpose of verifying some of the analytical findings. The baseline parameter values used in our computations are mainly from literature similar to this work, and unavailable parameter values are assumed for illustration. The parameter descriptions and values per year are shown in Table 2. Figures 2–5 illustrate the variations in human, wild animal, and livestock subpopulations while Figure 6 shows the existence of a globally stable disease-free periodic solution. Additionally, Figures 7–9 highlight the impact of temperature variations on the transmission dynamics of brucellosis. Figure 2 shows that the number of infective livestock decreases seasonally with an increase in time while Figure 2(a) illustrates a decrease in the susceptible animal subpopulation as time increases. The decrease in the number of infective livestock is due to proper implementation of vaccination and gradual culling of seropositive animals as control strategies. On the other hand, the sharp decrease in the susceptible animal subpopulation can be associated with the large number of infective animals and consequently high transmission rate in less than a one-year period of time while the gradual decrease in the next two years is due to vaccination programmes and decreased infection rate. Figure 3 shows a strong relationship between the number of infective and susceptible humans. For instance, at t=0, Sa=5000 and Ia=0 while at t=3, Sa=2555 and Ia=1850. The seasonal increase in the individuals in Figure 3(a) is associated with the low human treatment rate and poor control of the disease from infective livestock as well as contaminated environment. Besides, the decrease in the number of susceptible humans in Figure 3(b) is due to the high transmission rate from both infective animals and their products while the increase may be associated with proper implementation of the control strategies such as environmental hygiene, animal vaccination, and gradual culling of seropositive animals [24]. Figure 4 shows that the number of susceptible wild animals decreases with the increase in infective wild animals. In particular, the introduction of 200 susceptible wild animals in the contaminated environment produces more than 200 infective wild animals. This is based on the fact that both infective and susceptible animals have free movements and interactions within their parks. Besides, lack of wild animal brucellosis control measures and the fact that the disease does not kill keep the number of infected wild animals seasonally increasing. This implies that, in order to control the transmission dynamics of brucellosis in livestock and humans, interactions between domestic and wild animals should be restricted. Figure 5(a) shows that the number of Brucella bacteria in the environment decreases seasonally as the time increases while Figure 5(b) illustrates the variations in the number of recovered humans with respect to increase in time. These variations are associated with the regular implementation of the control strategies like environment hygiene and sanitation, human treatment, and gradual culling of infective animals. Furthermore, the recovered human population in the first six years increases due to effective treatment of the infective animals, and its decrease is associated with the decrease in the number of infected humans as well as proper control of the disease from livestock and their products. Figure 6 shows the existence of a stable periodic solution between the animal subpopulations and the number of Brucella bacteria in the environment. Figure 7(a) shows the seasonal variations in the effective reproductive number with respect to maximum daily temperature while Figure 7(b) illustrates the changes in the effective reproduction number with respect to seasonal variations in minimum daily temperature. Figure 8(a) illustrates the variations in the effective reproduction number versus maximum daily temperature while Figure 8(b) depicts the changes in the effective reproduction number with respect to seasonal variations in minimum daily temperature. Figure 9 presents the comparison between direct and indirect routes of brucellosis transmission. In particular, high strength of seasonal forcing shown in Figure 9(a) is due to seasonality in both direct and indirect routes of disease transmission while the curve with low amplitude shows the impact of lack of seasonality on the direct disease transmission. Moreover, Figure 9(b) indicates that seasonality in direct transmission has a significant contribution to the brucellosis transmission than that in indirect transmission; the graph in red is for seasonality in both direct and indirect transmission while the one in blue is for seasonality in both.Figure 2
Seasonal variations in the number of infective and susceptible animals.
(a)(b)Figure 3
Seasonal variations in the number of infective and susceptible humans.
(a)(b)Figure 4
Seasonal variations in the number of infective and susceptible wild animals.
(a)(b)Figure 5
Variations in the effective reproduction number with respect to changes in environmental hygiene and human treatment.
(a)(b)Figure 6
Relationship betweenBrucella spp. and susceptible and infected subpopulations.
(a)(b)Figure 7
Variations in the effective reproduction number with seasonal changes in temperature for the year 1979 in Mpwapwa District, Dodoma.
(a)(b)Figure 8
Variations in the effective reproduction number with seasonal changes in temperature for the year 2014 in Ngorongoro District, Arusha.
(a)(b)Figure 9
Variations in the effective reproduction number with seasonal changes in temperature for the year 2014 in Ngorongoro District, Arusha.
(a)(b)Generally, findings from this study advocate that, when the weather condition favours the increase in the transmission rates of brucellosis in livestock, humans, wild animals, and the environment, the incidence of the disease increases significantly and vice versa. This implies that in order to effectively prevent, control, eliminate, or eradicate brucellosis from the community, measures should be timely taken in accordance with the fluctuation in the disease transmission rates as a result of daily temperature variations. Thus, to avoid underestimation or overestimation of the resources when dealing with brucellosis, the aspect of seasonal weather variation should be taken into account when planning for prevention, control, elimination, or eradication of brucellosis infections.



---

*Source: 8972063-2020-10-17.xml*