Abstract
This paper presents a differential equation model which describes a possible transmission route for Q fever dynamics in cattle herds. The model seeks to ascertain epidemiological and theoretical inferences in understanding how to avert an outbreak of Q fever in dairy cattle herds (livestock). To prove the stability of the model’s equilibria, we use a matrix-theoretic method and a Lyapunov function which establishes the local and global asymptotic behaviour of the model. We introduce time-dependent vaccination, environmental hygiene, and culling and then solve for optimal strategies. The optimal control strategies are necessary management practices that may increase animal health in a Coxiella burnetii-induced environment and may also reduce the transmission of the disease from livestock into the human population. The sensitivity analysis presents the relative importance of the various generic parameters in the model. We hope that the description of the results and the optimality trajectories provides some guidelines for animal owners and veterinary officers on how to effectively minimize the bacteria and control cost before/during an outbreak.
1. Introduction
Q fever is a bacterial disease which is caused by Coxiella burnetii (C. burnetii) [1, 2]. Q fever is a potential biological warfare agent being very infectious and very durable in the environment as well as capable of windborne spread [3]. C. burnetii affects a wide range of animals and also causes illness in humans [1, 2]. The bacteria are mostly found in farm animals such as goats, sheep, and cattle, but can also be found in cats, dogs, rodents, birds, and other wildlife [4, 5]. Coxiella burnetii in ruminants result in abortion, stillbirth, mastitis, infertility, premature delivery, and weak offspring [6–8]. Other clinical signs in animals may include, fever, mild coughing, anorexia, and rhinitis [9]. The bacteria are shed from an infected animal into the environment through urine, faeces, milk, and vaginal fluids, but most commonly, the bacterium is in the amniotic fluids and the placenta discharge during abortion or parturition of an infected animal [5]. Q fever is noted as the second most commonly reported laboratory infection with several recorded outbreaks involving 15 or more persons [3]. The disease dynamics in both animals and humans start primarily through inhalation of contaminated dust, contact with placenta discharge during an abortion, parturition of an infected animal, drinking unpasteurised infected milk, ingestion meat containing Coxiella burnetii, or contact with contaminated wool [3, 5, 10, 11]. Outbreaks typically occur following birth or abortion where the environment becomes contaminated with birthing fluids of an infected animal [3]. The bacteria contaminate dust and then are spread by the wind for long distances [3].
Q fever is a reemerging zoonosis in most parts of Europe, which has seen a sharp rise in recent years, especially in the Netherlands with a large number of human cases, attributed to livestock [12–14]. An outbreak of Q fever in the Netherlands between 2007 and 2010 resulted in about 4000 reported human cases, and efforts to end the epidemic resulted in the culling of more than 50,000 small ruminants and a temporary ban on breeding of animals [10]. A survey indicates that 82% of cows in some California dairies were seropositive, as well as 78% of coyotes, 55% of foxes, 53% of brush rabbits, and 22% of deer in Northern California [9]. Australia meat industry makes about $1 million losses annually as a result of Q fever [15]. Geographically, Q fever infection in animals has been detected in most countries, except, New Zealand, Iceland, French Polynesia, and Norway [3, 10].
In recent time, Q fever has received an increase in publicity on the potential source of transmission and possible preventative measures in both livestock and humans (see, for example, [11, 16] and the references therein). The research by De Rooij et al. [17] suggested that during any future epidemics, attention should be placed on a rapid source of identification, quantification of emission, accurate data collection, and smooth data exchange amongst relevant actors to enable effective risk assessment and risk management. Mori and Roest [18] pointed out that in achieving an effective control and management, there should be collective expertise from other scientific disciplines and stakeholder so that the expectations of farmers and the wider community can be met on the spread of the disease. Courcoul et al. [19] used the Bayesian approach to obtain key epidemiological parameters from field data and suggested that the values of the parameters obtained can be used to provide information for calibrating simulation models to assess control strategies for C. burnetii. Courcoul et al. [20] further proposed a model to assess the effectiveness of vaccination in dairy cattle and suggested that their model should be modified to simulate various control strategies such as environmental and hygienic measures.
In this sense, our goal is to use the parameter values presented in [19, 20] and some assumptions to study the impact of vaccination, environmental hygiene (cleaning/burying of placenta discharge after birth or abortion), and culling/isolation using ordinary differential equation model, sensitivity analysis, and optimal control theory. Thus, we systematically analysed the environmental transmission of the disease and then obtained mathematical and epidemiological properties, such as basic reproduction number, equilibrium points, local and global stability analysis, sensitivity analysis, and optimality conditions. The optimal control introduced is to obtain optimal trajectories that depict the most effective control measures and also account for the costs involved. Optimal control is a useful mathematical tool which has been recently used to determine optimal strategies for other infectious diseases, although not for Q fever.
The rest of the paper is arranged as follows: Section 2 presents a vaccination-induced model and mathematical analysis, which establishes the stability of the proposed model; in Section 3, we study the local and global sensitivity analysis of the model’s parameters; in Section 4, we set up an optimal control problem; finally, in Section 5, we draw conclusions from this paper and assess what new insights this work gives to the body of knowledge on Q fever transmission dynamics and the way forward.
2. Basic Model for Indirect Transmission
The model is divided into susceptible , asymptomatic (expose) , vaccinated , and symptomatic cattle, respectively, with an assumed bacteria- (C. burnetii-) induced environment , where t represents time. The rate at which new cattle enter the susceptible population is denoted as , μ is the constant rate of death, and thus is the average lifetime. The asymptomatic cattle become symptomatic at a constant rate α, and so is the average asymptomatic period. The susceptible cattle get vaccinated at the rate ν, θ is the rate at which vaccinated cattle lose immunity and regain susceptibility. The symptomatic cattle are affected by an additional exit which may serve as a control measure at a constant rate c (culling/isolation rate). The rate of natural decay of the bacteria from the environment is denoted by δ and ε represents the rate of environmental hygiene, thus cleaning/burying of placenta discharge after birth or abortion. Asymptomatic cattle who develop temporal resistance to the bacteria become susceptible again at the rate ρ, if means no temporal resistance to the bacteria, and means resistant to the bacteria (thus ), and hence throughout this paper, we assume . We assume a mass action for the transmission of the disease, and thus susceptible cattle become asymptomatic at a proportional rate , where β is the effective rate of contracting the bacteria through indirect means (thus environmental transmission rate through inhaling contaminated dust/through grazing). The rate of inflow of the bacteria into the environment by both asymptomatic and symptomatic cattle is η (thus an assumed shedding rate). Here, the total cattle population for the vaccination-induced model is denoted as . The following assumptions govern the vaccination-induced model:(i)The parameters are all positive, and thus and (ii)Newly purchased or new birth goes into the susceptible class only(iii)No permanent recovery(iv)A successfully vaccinated animal stays immune for a period of time
Figure 1 depicts the compartmental diagrams of the disease in a community with control measures; model (1) gives the dynamical nature of the individual compartments with respect to time.
2.1. Basic Model and Analysis
with initial conditionsand total differential population
2.2. Positivity and Feasible Region
Lemma 1. For any nonnegative initial data,model (1) has a unique global solution for all .
Lemma 2. From the initial conditions and for epidemiological meaningfulness, we assume that all parameters are positive. Hence, the proposed model (1) will then be analysed in a feasible region given as
Now, we need to show that the region is positively invariant, to support the assertion that the model is epidemiologically useful and that all its state variables are nonnegative as time increases and stays in .
Proof. Let where denotes transposition, also if we replace with , then model (1) can be written as , whereNow, from the second equation in model (1), we haveapplying the integration factor method and change of variables, we haveFrom the third equation of model (1), we also obtainIt follows from standard property that , , and , provided , with similar properties being valid for , , which ensures that the state variable remains positive during the entire scope of the study (see [21] and Appendix A of [22]). Based on this, we see that matrix D has all its off-diagonal entries to be nonnegative and that matrix , which proves the properties of Metzler matrix [23]. Therefore, it implies that model (1) is positively invariant in . Finally, using equation (3), it can be shown that all solution sets of model (1) remain in the feasible region , with the notation that and . Hence, the proposed model is considered to be epidemiologically meaningful (mathematically well posed) and it is sufficient to say that all solutions of the Q fever model with initial conditions remain in the feasible region .
2.3. Equilibrium Points and Basic Reproduction Number
The equilibrium points are the constant solutions to the differential equations in model (1), and these equilibrium points and the basic reproduction number play an important role in the long-term behaviour of the solution to the epidemiological model (1). The equilibrium points of model (1) are obtained by setting the right-hand side of model (1) to zero, which gives
The two main equilibria considered here are the disease-free equilibrium and endemic equilibrium. We obtained the disease-free equilibrium by assuming that the disease-induced states and a free C. burnetii environment and thus . Therefore, by solving for the nonzero state variables in equations (10) and (11), a unique disease-free equilibrium is obtained as follows:
2.3.1. Basic Reproduction Numbers and
The basic reproduction number is obtained by using the disease-free equilibrium points in terms of the original parameters of the model; epidemiologically, the reproduction number informs us on the number of secondary infections produced when a single infectious individual is introduced into a completely susceptible population. Here, we define the basic reproduction number as the number of secondary infections produced in a completely susceptible cattle population when the organism is released into the environment. Also, the analytic nature of the basic reproduction number gives clues on parameters in the model that contributes to the spread of the disease. Hence, with the aid of the next-generation matrix as introduced by [24, 25] and the disease-free equilibrium (15), the epidemiological reproduction numbers for model (1) are as follows:where equation (16) is called the control reproduction number and equation (17) is the basic reproduction number without controls (thus ). From the basic reproduction numbers, it is evidential that there exists a linear relationship between the rate of environmental transmission β and shedding rate η, and hence, this shows that an increase in bacterial load in the environment may have a corresponding increase in the transmission process. This assertion has been numerically varied using Figure 2(b) in Section 3.2.1.
(a)
(b)
(c)
2.3.2. Endemic Equilibrium Point
Now, taking in the feasible set , and using equations (10)–(14) together with (16), we obtain the following endemic equilibrium point . Thus
The endemic equilibrium has a unique equilibrium point if the vaccination-induced reproduction number is greater than unity . In the next subsection, we investigate the stability of the equilibrium points.
2.4. Stability Analysis
In this section, we focus on establishing the local and global stability of the model equilibria. To obtain the local stability of the disease-free equilibrium, we constructed a linearized Jacobian matrix , evaluated at the disease-free equilibrium, which giveswhere , is irreducible and nonnegative. The simple eigenvalue of is denoted as which is obtained from the relation . The following equivalence holds:
Using Theorem 2 of van den Driessche and Watmough [24], we see that assumptions are satisfied, and also is satisfied if all the eigenvalues of matrix have negative real parts, and thuswhere , and .
The simple eigenvalue of is . Hence, model (1) fully satisfies the assumptions of van den Driessche and Watmough [24], and therefore, the disease-free equilibrium of the Q fever model is locally asymptotically stable if and unstable if . The below theorem addresses the global stability of the disease-free equilibrium.
Theorem 1. If , then the disease-free equilibrium of model (1) is globally asymptotically stable.
Proof. Applying Theorem 2.1 of the matrix-theoretic method in [26], the global stability of the disease-free equilibrium can be established by constructing a Lyapunov function , since is reducible (the third column is the only nonzero column), is the left Perron eigenvector of the matrix , and is the disease-induced classes. Algebraic operation of giveswhere . Therefore, by Theorem 5.1 of [26], the left Perron eigenvector of for the largest eigenvalue is defined asThis giveswhich implies that ; therefore, . Now simplifying yieldsUsing [26], the derivative of along the solution of model (1) is illustrated as follows:where . Therefore, after some algebra, we haveNow using equation (26) and some algebraic simplification leads toTherefore, provided and . Furthermore, implies that and or . It can be shown that the only invariant set where or is a singleton of the set . Since the four state variables used in constructing the Lyapunov functional lead to global stability, it then means that as . Hence, by [27] (LaSalle’s invariance principle), the disease-free equilibrium for the vaccination model is globally asymptotically stable. The ramification of this global stability in the presence of vaccination implies that no matter how many initial asymptomatic and symptomatic cattle in the population, the disease cannot be established (persist) and eventually dies out. It also indicates that no matter how many initial bacterial loads in the environment, the disease will die out from the population in the long round.
2.4.1. Global Stability of Endemic Equilibrium
The nature of the endemic equilibrium indicates that for , then there exists a unique endemic equilibrium resulting in an unstable disease-free equilibrium. To establish the global stability of , here we explore a special case (simpler case) of the model (1), where . Let be the stable manifold of . Therefore, we claim the following result.
Theorem 2. The unique endemic equilibrium of the model (1) is globally asymptotically stable in whenever and .
Proof. At the endemic equilibrium state of model (1) with , the following equalitieshold, which further leads toSubstituting equation (34) into model (1) result in model (35), in which its stability is the same as model (1):Now, we define a Lyapunov function aswhere and are positive constants to be determined, and thus . The time derivative of the Lyapunov function along with the solution of model (35) yields the following equation:Simplifying the above expression and rearranging lead toWe now letMaking and ν the subject in equations (29) and (30), respectively, and substituting into (39) with series of simplification giveFurther elaboration on results inTo obtain , we make the following assertion by equating the constant coefficients of the terms and in to zero, thusSolving (44) and (45), we obtain and . Replacing into (42) and simplifying yieldHence, adding equations (41) and (46), we havewhereThus, if and from the fact that the arithmetic mean is greater than or equal to the geometric mean, then it follows that , provided the coefficients are nonnegative. With the claim thatit leads toTo verify the above expression, the definition of and properties of logarithms were used to simplify the expression on the right hand of . Then applying the proposition of [28] givesIt is evidential that only holds when and that is the only equilibrium state of model (1) on its line. Therefore, by Lyapunov–LaSalle asymptotic stability theorem [27], the positive equilibrium state is globally asymptotically stable in the positive region whenever and .
3. Sensitivity Analysis
Here, we apply the concept of sensitivity analysis to ascertain the relative importance of each generic parameter in the basic reproduction number and the disease-induced state variables, . Secondly, with parameter values from published papers and some realistic assumptions, we obtain both analytic and numerical values of the various parameters in . The analytic expressions obtained can be used to shed some light on how to control the onset of the disease in variant localities, if and only if the dynamics follow model (1). Parameter values were obtained from [19, 20], which were then converted into days; this is to enable us to see the daily trajectories of the model. The estimated parameters were calculated based on the life expectancy of dairy cattle and the natural decay rate of the bacteria from dust/soil (which ranges from 120 to 150 days [9]), and finally, we used the gestation period of cattle to denote the rate of new inflow of susceptible cattle/recruited cattle. All parameters used are given in Tables 1 and 2, respectively.
3.1. The Impact of Transition Rate ρ on the Reproduction Number
Taking the partial derivative of equation (17) with respect to ρ gives the relative importance of the discontinuation rate ρ (the feedback term) and thus the rate at which asymptomatic cattle build antibiotics and reenter the susceptible class:
Hence, the basic reproduction number without culling, vaccination, and environmental hygiene is a decreasing factor of the parameter ρ as shown in Subsection 2.4.1. This implies that the basic reproduction number decreases with increase in the value of ρ. Consequently, taking the limit of as ρ becomes large gives
This means that any attempt to increase ρ will result in reducing the number of secondary infections, thereby eradicating the disease in the absence of other control measures, but a low rate of ρ will produce a high rate of asymptomatic cattle which may require effective control measures as shown in Section 4. From this, we infer that the development of any daily antibody boosters for dairy cattle (livestock) will help curtail the reproduction losses caused by C. burnetii and other bacterial diseases.
3.2. Calculation of Sensitivity/Elasticity Indices
Fundamentally, the onset of disease transmission is directly associated with the basic/control reproduction number , and the prevalence of the disease is associated with the endemic equilibrium point [29]. The sensitivity indices for the state variables are calculated as follows:where represent the equilibrium points at the endemic state and are the corresponding parameters in [30]. For simplicity, the numerical values of the endemic state variables are shown in Table 2 without analytic expression due to the complex nature of the analytic result. Similarly, the sensitivity/elasticity indices of to the parameters in the model are also defined as follows:where the meaning of and remains unchanged. Therefore
Remark 1. whenever , this implies that an increase in culling/isolation of symptomatic cattle from farms during an outbreak will reduce the impact of α, thereby reducing the number of secondary infections and bacterial load in the environment as shown in Table 2. This further suggests that in choosing between culling and antibiotic booster during an outbreak, it is prudent to focus on culling/isolation.
In Table 2, we obtain numerical values which indicate the relative importance of parameters in and the state variables at the endemic stage when other parameters are kept constant (depicting the local sensitivity of parameters in model (1)). The result reveals that and θ influence the stability of the disease and are also closely related. Thus, a increase in , and θ will have a corresponding percentage increase in the basic reproduction number. Among the control measures implemented, the elasticity index in Table 2 indicates that the most efficient method in the infected classes will be vaccination ν and the most efficient method of reducing the bacteria from the environment may be a good environmental hygiene ε. The total efficacy of these control measures will be shown in Section 4. From Table 2, it also shows that a 0.1 effort in keeping the environment clean/burying of placenta discharge after birth or abortion may result in a 2.70 reduction in bacterial load in the affected herds’ environment and a corresponding reduction in the number of asymptomatic and symptomatic by 1.78, respectively. Interestingly, we notice that a change in natural death may reduce the reproduction number by 0.95, but natural death cannot be used as a control measure; it, therefore, presupposes that the disease could be established as long as cattle live longer in an environment that is induced by the bacteria when there are no control measures; hence, vaccinating cattle and good environmental hygiene in farms will help mitigate the intensity of the organism on farms. It is well noticing that the analysis done so far represents the local sensitivity analysis of parameters in the model, which shows a close relationship between some of the parameters; hence, to see the global effect of these parameters on the model, we adopt the Latin hypercube sampling (LHS) and partial rank correlation coefficient (PRCC) which is demonstrated in the next Subsection 3.2.1.
3.2.1. Global Sensitivity Analysis
In carrying out the global sensitivity analysis on the control reproduction number , we use the analytic expression in equation (16) and the range of values in Table 1. The objective is to ascertain the most influential parameters in and also give some insightful ideas from the PRCC plot, by using a sample size of 1000.
The scatter plot in Figure 2(a) shows a positive relation between and the parameters and a negative relation between and the parameters ν and The positive relation indicates that a high rate of either of these parameters will produce a high transmission rate during an outbreak, while the negative relation for ν and ρ indicates that an increase in parameter value of ν and ρ will help reduce the severity of the disease-induced rate. Figure 2(b) gives sharply defined simulation results of the scatter plots in Figure 2(a) and that of the numerical signs in Table 2, but we notice a reverse sign for α, which is as a result of a decrease in the culling/isolation of symptomatic dairy cattle, as indicated in Remark 1; therefore, as shown in Figure 2(b) that, to have a positive impact of c in combination of other control measures so as to reduce α, c should be greater than 0.12, and thus it is represented by the red line as . Figure 2(b) also indicates that the most influential parameters that promote the disease spread are and α (thus with low culling/isolation rate) and the most influential control parameter is ε (environmental hygiene). Therefore, any measures to reduce the relative relevance of and α (with an increase in culling/isolation rate) will reduce the spreading rate of the bacteria in dairy cattle herds. Figure 2(c) shows the saturation effect of vaccination and environmental hygiene, and thus having and may lead to a Coxiella-free environment.
4. Optimal Control Problem
Whenever the number of unknown parameters in a model is larger than the number of state equations, then the model in question usually allows for more than one solution. The optimization problem for such situations in continuous models is known as an optimal control problem. Here, we investigate the impact of time-variant controls: culling, vaccination, and environmental hygiene (cleaning/burying of placenta discharge after birth or abortion). This is to obtain the best control measure in preventing an outbreak of the disease in a situation where the cause of new infections is solely environmental. Thus, we seek to find optimal trajectories that show the relative importance of applying vaccination, culling, and environmental hygiene in controlling Q fever in an endemic setting. To see the impact of these controls, we modify model (1), thus setting (no natural resistance). The objective function is to reduce the numbers of asymptomatic cattle, symptomatic cattle, and bacterial load through vaccination ν, culling c, and environmental hygiene ε, respectively, and the cost associated with the implementation of the controls. Thus is the percentage of susceptible and asymptomatic cattle who are vaccinated per unit time, is the percentage of symptomatic cattle who are culled per unit time, and is admissible measure for keeping the environment clean. Here , where and is the final time for running the control measures. Now the objective function that minimizes our desirable control problem is given assubject towith initial conditionswhere is the reduction in transmission due to environmental hygiene, are small positive constant to keep balance in the size of , , , and , respectively [31]. The square term of the controls reflects the nonlinearity in the cost of controls, and the half-term minimizes the effect of applying the controls [32]. The positive constants which are associated with are weight values such that , and . The relative weights given to the positive constants associated with the control terms indicate greater or lower importance placed on minimizing the cost of a control measure [33]. Now to obtain the optimal solutions, we defined a Lagrangian L for the control problem, given as follows:
We now seek the minimal value of the Lagrangian function. Hence, we obtain a Hamiltonian H for the control problem using Pontryagin’s maximum principle:where are the adjoint variables associated with the states and B. Finally, taking the partial derivatives of H with respect to the state variables and the controls, we obtain the following relation for the adjoint and the optimal controls:with transversality conditions (or boundary conditions) [31], . Therefore, the adjoint system is given as
By the optimality conditions, we have
Using the optimal solution in the control space [31], we have
Therefore, the compact notation (characterization) of the admissible control is then written as
The control terms in the bound of 0.9 and 0.8 satisfy the inability to vaccinate all susceptible cattle at a time and a 0.2 failure in eliminating all bacteria in a particular farm at a given time, note that the control bounds are assumed for a realistic illustrative purposes and may not be practically practised, but the variation of these bounds will not narrow the essence of the studies. Substituting (67) into (59), (62), and (64), respectively, we obtain an optimality system with the corresponding Hamiltonian system which is then solved numerically to obtain optimal trajectories for our dynamical system.
4.1. Optimal Control Trajectories
In order to obtain the optimal trajectories for model (59), we employ the forward-backward sweep method, as described in [34]. We begin with an initial guess for the controls. The state equation (59) is solved forward in time while the adjoint function (64) is solved backward in time. The controls are then updated for a period of time until convergence is achieved. For better numerical output, we firstly made use of the following initial state values: and the stated parameter values in Table 1. The initial bacterial load is assumed to be small, the reason being that per the numerical simulations, the bacterial population grows more and more rapidly as depicted in Figures 3(f) and 4(f), respectively, by the red line, hence, the need to start with small initial values for the bacterial load so as to get a suitable numerical output for that of the initial conditions chosen for the animal population. The corresponding weight in the objective function (58) is taken to be . The figures below depict the relative importance of finding optimal strategy/strategies in controlling the disease spread in cattle (animal population). Table 3 shows the total cost associated with the objective function (58) after a simulation time of 100 days. Secondly, we look at the control trajectories and the accompanying cost by changing the associated weights of control, cost of control implantation, and initial state values and thus and as depicted in Figure 4 and Table 3, item 3.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
Figure 3(a) shows the control profile for the various bounds employed in model (59), which indicates that a bound of 0.8 for environmental hygiene is maintained for the entire simulation time, but a drop in the control bound for both vaccination and culling. Thus, the bound for culling and vaccination switches to 0.32 after 10 days. This indicates during the initial onset of an outbreak, effort should be made to cull all symptomatic cattle or isolate them from other susceptible (healthy) cattle for the first 10 days, after which less culling approach can be adopted, and also for first 10 days, at least of cattle should receive vaccination during the initial stage of the disease outbreak. Figures 3(b) and 3(c) show that applying vaccination, culling and environmental hygiene simultaneously increases the rate of healthy cattle, thereby reducing the risk of infection in the susceptible cattle population and also reducing the number of asymptomatic cattle, symptomatic cattle, and bacterial load in the environment as shown in Figures 3(d)–3(f), respectively. We also notice that the gradual use of these controls in the long run is more cost-effective as shown in Table 3, item 2.
In Figure 4(a), we changed the initial population size, which indicates that an increase in population size requires that the initial control bound for vaccination be kept for 80 days, and thus an increase in population size of dairy cattle herds requires high control effort and high resources as depicted in item 3 of Table 3. The control trajectory for the vaccination bound and Figure 4(c) further indicate that for optimal eradication or reduction of the disease in a population of 1000 cattle (livestock), at least of the animals should be vaccinated within the first 20–80 days. We also noticed that the use of vaccination, culling, and environmental hygiene has similar control efficacy in the disease compartments irrespective of the population size as shown in Figures 4(b)–4(f), respectively.
In Figure 5, we considered variations that may arise in the control measures, thus sticking to one or more controls. In Figure 5(a), we assumed that , that is a situation where culling is not applied, and it reveals that symptomatic cattle will always be in the population but stays below the trajectory when all controls are zero (see Figure 5(b)). Therefore, we notice that in a situation where a farm owner chooses to ignore culling and decide to stick to only vaccination and environmental hygiene during an outbreak of the disease may not lead to the eradication of the disease; hence, the combination of vaccination and environmental hygiene is not a good control measure in reducing the number of symptomatic cattle (see Figure 6(d)) for the implementation of culling only as a control measure on the symptomatic compartment. Figure 5(c) shows that in the absence of vaccination, the control bound for culling should remain at 1, which means that during the cause of culling, the aim should be focused on destroying all infected farms or animals which is sometimes practically hard to achieve and also comes with high control cost as shown in Table 3, item 7. Therefore, it clearly shows in Figures 5(d)–5(f) that effective vaccination plays a major role in reducing the disease-induced rate in dairy cattle herds depending only on culling. It can also be seen that without vaccination, the trajectories of Figures 5(d)–5(f) of the optimal control stay very close to the trajectories without optimal control.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
Figure 6(a) shows the control profile in the absence of environmental hygiene; the figure depicts that with only vaccination and culling as control measures, the bound of vaccination is kept at 0.9 for 78 days before switching to 0.8, and the bound of culling switches to 0.45 as compared to Figure 3(a). Hence, Figure 6(b) depicts the noneffectiveness of vaccination and culling on the bacterial load in the environment until day 47, but due to vaccination, the number of asymptomatic cattle reduces as shown in Figure 6(c). Figure 6(d) shows the effect of culling on the number of symptomatic cattle, which indicates that early detection of symptomatic cattle on the farm will help in mitigating the disease within the first 5 days.
5. Concluding Remarks
We presented a theoretical analysis of Q fever within dairy cattle herds, which may lead to new research on the disease dynamics, such as the effect of direct transmission, seasonal transmission, the dynamical role of ticks in maintaining the circle of Q fever in livestock, the role of media in averting the disease in humans, and the dynamical spread of the disease between livestock and humans. Here we analytically obtained the control states and numerically simulated it for different outcomes. We notice that culling only as a control measure produces a higher control cost, which also does not completely eradicate the number of symptomatic cattle due to the asymptomatic nature of the disease. We obtained a cost profile for the various control measures as shown in Table 3 which suggests that the most effective method in getting a lower control cost, in the long run, is the continuous application of vaccination and good environmental hygiene together with culling or isolation of symptomatic cattle. The numerical simulation also reveals that irrespective of the control employed, the disease will always be present in the animal population, but vaccinating and cleaning/burying of placenta discharge after birth or abortion will help to mitigate the disease in the animal population. In totality, our work has revealed some global dynamics of Q fever in dairy cattle herds and the impact of sensitivity analysis and optimal control strategies in managing Q fever disease in dairy cattle herds (livestock). We also analytically obtained the local and global stability analysis of the disease by using matrix-theoretic method and a Lyapunov functional. The model studied here did not assume a complex model, but the results can provide some ideas for further studies and the prevention of the disease.
Data Availability
The parameter values (data) used to support the findings of this study have been described in Section 3, and other data (parameter values) used were previously reported data and are cited at relevant places within the text as references [19, 20].
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (11331009), National Natural Science Foundation of China (61873154), Shanxi Key Laboratory (201705D111006), and Shanxi Scientific and Technology Innovation Team (201705D15111172). The first author is grateful to the Chinese Government and the Complex Systems Research Center, Shanxi University, for their support. The first author also thanks Ming-Tao Li for his advice on the global stability of the endemic equilibrium.