On the Role of Short-Term Animal Movements on the Persistence of Brucellosis

Short-term animal movements play an integral role in the transmission and control of zoonotic infections such as brucellosis, in communal farming zones where animal movements are highly uncontrolled. Such movements need to be incorporated in models that aim at informing animal managers effective ways to control the spread of zoonotic diseases. We developed, analyzed and simulated a two-patch mathematical model for brucellosis transmission that incorporates short-term animal mobility. We computed the basic reproduction number and demonstrated that it is a sharp threshold for disease dynamics. In particular, we demonstrated that, when the basic reproduction number is less than unity, then the disease dies out. However, if the basic reproduction number is greater than unity, the disease persists. Meanwhile, we applied optimal control theory to the proposed model with the aim of exploring the cost-effectiveness of different culling strategies. The results demonstrate that animal mobility plays an important role in shaping optimal control strategy.


Introduction
Brucellosis, a highly contagious zoonotic disease, remains a significant public health threat worldwide.It is estimated that more than 500,000 new cases of the disease are reported annually [1], with incidence as high as 200 cases per 100,000 in most endemic countries [2].Majority of brucellosis infections occur in: Sub-Sahara Africa in countries such as Ethiopia, Chad, Tanzania, Nigeria, Uganda, Kenya, Zimbabwe and Somalia due to high level of pastoralism; and the Middle East, Spain, Latin America and Asia-in particular Southeast Asia-where factors such as pastoral farming practices, beliefs and lack of bio-security have been attributed to persistence of the disease [3].Since human transmission of brucellosis is considered to be negligible [4], measures to effectively control brucellosis in humans ultimately require a thorough control of the disease among domestic cattle, camels, goats and sheep.
Transmission and control of brucellosis in both human and animal population remains a complex phenomenon that possibly involves the type of farming practiced in the area; economic, geographic and environmental structures; and the intrinsic disease biology and ecology.In particular, animal movement plays crucial role on transmission and control of the disease.For example, in communal farming zones, animal movements are highly uncontrolled compared to private farming.Prior studies have demonstrated that, on a daily basis, a single cattle herd in a communal farming zone has the potential to mix with at least five heterogeneous herds at both the communal grazing and watering points.Since livestock management varies from one farmer to another, it is evident that understanding the volume of these movements and the risks associated with them is fundamental in elucidating the epidemiology and control of animal diseases.
Mathematical models have proven to be important tools that can aid our understanding as well as provide solutions to phenomena which are complex to measure in the field.Recently, mathematical models have been proposed to explore brucellosis transmission and control (see, e.g., [5][6][7][8][9][10][11]).For example, Dobson and Meagherin [6] used nonlinear ordinary differential equations to describe brucellosis transmission among the bison population in the Yellowstone National Park (YNP).Abatih et al. [7] mathematically analyzed the brucellosis model proposed in [6].Lolika et al. [9] applied a non-autonomous model to discuss the effects of optimal vaccination and environmental decontamination on long-term brucellosis dynamics among cattle in periodic environments.Yang et al. [5] developed a two-patch model with risk heterogeneity in which animals immigrated between two different risk environments.Their work utilized a Eulerian approach for mobility.However, the Eulerian approach has some limitations, for instance it neither incorporates the concept of residence times nor the effective population size.Here, the term residence times refers to the average proportion of daily time an animal spends in a given patch.Therefore, to gain a better and more comprehensive understanding of effects of animals movements on brucellosis dynamics, a model should incorporate a Lagrangian approach that can account for the effects of residence time and the effective population size per patch.
In this paper, we consider a dynamical model to describe the role of short-term animal movements on the persistence of brucellosis.The proposed two-patch model incorporates all the relevant biological and ecological factors as well as short-term animal movements which are modeled using the Lagrangian approach.For the purpose of distinction between the hosts, we assumed that Patch 1 is a high risk environment, that is, brucellosis control measures in this patch are poorly managed.The reverse is assumed for Patch 2. Thus, disease transmission in Patch 1 is assumed to be higher relative to Patch 2. Further, disease transmission is assumed to occur through direct contact and vertical transmission.In addition, since vaccines are often unavailable or expensive to farmers in communal farming zones, we assumed that a more sensible approach to control the spread of the disease is culling of infected animals.

Modeling Framework
We developed a mathematical model to study the transmission and control of brucellosis within an environment defined by two patches of heterogeneous risk.Our model is a modification of the one developed in [7].Precisely, the model in [7] is a single-patch framework.
Let N i (t) represent the total population of animals in Patch i at time t, i = 1, 2. We assume that animals of Patch i spend Consequently, the effective proportion of infectious individual in Patch 2 is The following system of ordinary differential equations (ODES) account for the brucellosis dynamics in two patches: ( where the variables S i (t), I i (t) and R i (t) represent the susceptible, infectious and recovered population, respectively; µ i is recruitment rate of animals and it is assumed to be equal to natural death rate of animals, thus µ −1 i represents the animal's commercial lifespan; e i (0 ≤ e i ≤ 1) denotes a proportion of new recruits that are infected with brucellosis and the complementary proportion (1 − e i ) represents those that are susceptible to infection; β i denotes the disease transmission; α i is the recovery rate; and δ i denotes immunity waning rate.Disease related mortality is considered negligible.Thus, the total population is constant and is given by N i (t) = S i (t) + I i (t) + R i (t).The Parameters and values are shown in Table 1.

Positivity and Boundedness of Solutions
It can easily be verified that the domain of biological interest is positively invariant and attracting with respect to the model in Equation (1).

Disease Dynamics for a Single Patch
If only a single patch, that is, i = 1, is considered, then the system in Equation (1) reduces to (3) The system in Equation ( 3) is isomorphic to the model proposed by Dobson and Meagherin [6] and analyzed by Abatih et al. [7].As highlighted in [7], the model in Equation ( 1) is well defined, supporting a sharp threshold property, namely, the disease dies out if the basic reproduction number R 01 is less than unity, persisting whenever R 01 > 1 where R 01 = (β 1 + e 1 µ 1 ) (α 1 + µ 1 ) .

The Reproduction Number
The disease-free equilibrium E 0 of the system in Equation ( 1) is E 0 : (S 0 1 , S 0 2 , I 0 1 , I 0 2 , R 0 1 , R 0 2 ) = (N 1 , N 2 , 0, 0, 0, 0).The basic reproduction number, denoted by R 0 , is integral quantity in epidemiological model.It accounts for the average number of secondary infections generated by a single infectious animal introduced in a fully susceptible population during its average infectious period [12].We utilized the next generation matrix approach [12] to determine R 0 .See Appendix A for the derivation.The basic reproduction number for the system in Equation ( 1) is with We can write Equation (4) as follows where R 0i (i = 1, 2) represents the disease risks for Patches 1 and 2 in the absence of animal mobility.From Equation (4), we can observe that the basic reproduction number is influenced by short-term animal dispersal.
To investigate the effects of short-term animal dispersal on the generation of new infections, we compute the values of the basic reproduction number using a residence-time matrix in Table 2.
Table 2. Association between the basic reproduction number and the residence-time matrix.

Description
R 0 More precisely, the residence-time matrix configuration incorporates the coupling intensity and mobility patterns.For instance, weak coupling implies that most animals stay in their own patch while strong coupling implies that certain proportions of animals move to the other patch.Mobility patterns represent the symmetry of animal movement between the two patches.For example, symmetric mobility represents a scenario when an equal ratio of animals move from Patch 1 to Patch 2 and vice versa.However, if the ratio of animals that move between the two patches is not equal, then the mobility pattern is asymmetric.Note that the total population of animals in the two patches is assumed to be the same.
Results in Table 2 demonstrate that the basic reproduction number will always be high when coupling intensity is weak, that is, when most animals stay in their patch.Further, the highest value of the basic reproduction number occurs when the mobility pattern is symmetric.Using parameters and values in defined in Tables 1 and 2, we calculated the reproduction numbers for Patches 1 and 2 in the absence of animal dispersal and we obtained R 01 = 1.4 and R 02 = 0.05.We can observe that, based on our assumption that Patch 1 is high risk, the highest reproductive number came from this patch.In addition, we can observe that, whenever there is animal mobility, the disease transmission risk increases globally compared to locally; for instance, in the absence of animal mobility, we expect brucellosis to die off in Patch 2. It is worth noting that the results in Table 2 shows that, when animal mobility increases, the basic reproduction number decreases; however, for all cases demonstrated in Table 2, it will never drop below 1.Hence, under our assumption, we can conclude that effective brucellosis control will always be difficult to attain whenever there is animal mobility.

Disease Invasion and Persistence
From the work in [12], we know that the DFE is locally asymptotically stable when R 0 < 1, and unstable when R 0 > 1.Indeed, we can establish a stronger result regarding the global dynamics of the DFE.
The system in Equation ( 1) is said to be uniformly persistent in the interior Ω if there exists a constant η 0 > 0 such that lim inf Biologically, a uniform persistent system indicates that the infection persists for a long period of time.Thus, we have the following result.
Theorem 2. If R 0 > 1, then the DFE is unstable and the system in Equation ( 1) is uniformly persistent in Ω.
Theorem 3. If R 0 > 1, the system in Equation ( 1) has a unique equilibrium E * , which is globally asymptotically stable.
The proof of Theorems 2 and 3 in Appendixes C and D, respectively.In addition, due to the complexity of the non-trivial solutions of the system in Equation (1), we utilized numerical simulations to demonstrate that the the model in Equation ( 1) has a globally stable endemic equilibrium point (see Figure A1).

Optimal Culling
Vaccination and culling of infected animals are the only feasible ways to control brucellosis transmission.Vaccinating animals prevents susceptibility to the disease and culling of infectious animals reduces the density of infected animals thereby reducing the contact between susceptible and infected animals.However, in many brucellosis endemic countries, farmers cannot afford the cost of vaccines, and this leaves culling as the only disease intervention strategy.In this section, we explore the impact of culling on controlling the spread of the disease.Thus, we modify the model in Equation ( 1) to include culling control u i (t), i = 1, 2. The controls, u i (t), are represented as functions of time and assigned reasonable upper and lower bounds.The modified model is given by ( The control set is defined as where U i denotes the upper bound for the culling effort in Patch i. In the following, we introduce an objective functional J to formulate the optimization problem of interest, namely, that of identifying the most effective strategies over the admissible set of (u 1 (t), u 2 (t)).The overall objective is to minimize the numbers of infectious animals over a finite time interval [0, T] at minimal costs.The objective functional J is thus defined as where J 1 and J 2 represent objective functions for Patches 1 and 2, respectively.A i , B i and C i are positive balancing coefficients transferring the integrals into monetary quantity over a finite period of T years.Precisely, A i represents the cost (due to the loss of animals) associated with the number of infected animals in Patch i and B i represents the cost associated with the number of infected animals culled in Patch i.The objective function in Equation ( 6) also includes quadratic terms with coefficients C i to indicate potential non-linearities in the costs.
The existence and uniqueness of optimal control can be proven by applying a standard results in optimal control theory [13,14].The necessary conditions that optimal controls must satisfy are derived using Pontryagin's Maximum Principle [15].Thus, the system in Equation ( 5) is converted into an equivalent problem, namely the problem of minimizing the Hamiltonian H given by: where λ g i (t), g = S, I, R, i = 1, 2, are the adjoint functions to be determined.Thus, given an optimal control pair (u * 1 , u * 2 ) and corresponding states (S i , I i , R i ), there exist adjoint functions [13] satisfying From Equation ( 7), we have with transversality conditions λ g i (T) = 0. Furthermore, the optimal controls are characterized by the optimality conditions: In the following, we utilize the forward-backward sweep method [13] together with parameter values in Table 1 and the residence-matrix defined in Table 2 to determine numerical solutions of our optimality system.Our main goal is to explore the effects of optimal culling on the transmission and control of brucellosis under the following cases: In the above scenarios, we assumed that culling intensity in low risk population is always above average and we fixed it at u 2 = 0.8.Scenario 1 is assumed to apply to farmers who rear livestock near game reserves.Prior studies highlighted that livestocks reared in proximity to game reserves mix with wildlife on almost daily basis [16], even though, in many countries where brucellosis is endemic, intervention measures to control the spread of zoonotic infections among wildlife are not available.Scenario 2 represents heterogeneity on culling intensity.This scenario may exist in communal farming zones where one farmer, X, may have resources (knowledge and financial capacity) to perform culling at the high intensity while another farmer, Y, does not have enough resources to perform culling at an intensity that does not exceed the average.
In all simulation results presented in this section, we used parameter and initial values in Table 1 as well as the residence matrix in Table 2.For simplicity, in our numerical simulation, we set A 1 = A 2 = 1 so that the minimization of the infectious animal population has the same importance/weight in all patches.Further, we set The values of the weight constants B i and C i were determined through numerical simulations, precisely for these values the cost are low and the control efforts can be applied at maximum intensity in all scenarios suggested above.
For each strategy and coupling intensity described in Table 2, we found the total number of new infections given by the following formula where Γ i represent the total number of new infections for path i and the total cost associated with infected animals and the controls J, which is given by Equation ( 6).In the following, we determine the effects of optimal culling under different coupling intensity and mobility patterns (see Table 2).
In Table 3, we present the values of the total number of new infections and J for Scenario 1.We can clearly observe that the highest total number of new infections recorded in Patch 1 over a ten-year period under all possible coupling cases is Γ 1 = 7.12 × 10 3 and this occurs when the coupling intensity is weak and the mobility pattern is symmetric.Moreover, when the coupling intensity is weak and the mobility pattern is symmetric, Patch 2 records the lowest total number of new infections as Γ 2 = 950.321under all possible coupling cases over the same period.However, this coupling case (weak and symmetric) is associated with the lowest total number of new infections Γ = 8.06 × 10 3 as well as the total cost J = 162.15.We surmise that, due to weak animal mobility, the spread of the disease will be highly confined in independent patches, with more infections being observed in the high risk patch (Patch 1).
Table 3.The total number of newly infected animals over a ten-year period and the total cost J with respect to the control strategy under Scenario 1.In Table 3, we can also observe that strong symmetric coupling gives the lowest total number of new infections for Patch 1 only, Γ = 5.81 × 10 3 , while Patch 2 will record the highest total of new infections, Γ 2 = 3.95 × 10 3 , and overall this will yield the highest total of new infections, Γ = 9.76 × 10 3 , in the community.This clearly demonstrates that increased short-term dispersal of animals strongly influences the transmission and control of brucellosis.
Next, we compare the impact of presence and absence of time dependent culling on brucellosis transmission dynamics under Scenario 1 (Figures 1-4) over a ten-year period.Figures 1-4 show the number of infected animals per patch, with and without optimal culling under weak symmetric coupling, strong symmetric coupling, weak asymmetric coupling and strong asymmetric coupling, respectively.As we can observe, whenever the coupling is weak, despite its skewness, the optimal control policy will not have a significant impact in Patch 1 compared to Patch 2 where the number of infections decrease with time.However, whenever the coupling is strong, the number of infected animals in both patches decrease with time but with more effect being noticed in Patch 2 where there is disease control.Figure 5 shows the optimal control profile for u 2 (t): (a) when the costs of culling are low; and (b) when the costs of culling are high (we set B 2 = C 2 = 2).Recall that, due to the absence of control in Patch 1, u 1 (t) = 0.As shown, when the costs of culling are either low or high, the control profile starts from the maximum initially and stays there for more than half of the entire period before it switches to its minimum.Precisely, when the costs of culling are low, the control profile stays at its maximum for a longer period compared to when the costs are high.This clearly demonstrates that the control is highly sensitive cost parameters, thus under low costs optimal culling can be implemented at maximum intensity for a long period.We further investigate the impact of low intensity optimal culling in the risk patch (Patch 1); we set u 1 = 0.45 while u 2 remains fixed at 0.8.Results for this scenario are depicted in Table 4 and Figures 6-10.As observed earlier (Table 3), the highest total number of new infections occurs when the coupling intensity is weak and symmetric.We also observe that the presence of control in Patch 1 leads to a reduction in the total number of new infections by 30.1%, 21.4% and 28.9% in Patch 1 only, Patch 2 only and overall (Patch 1 and Patch 2 combined), respectively.In Table 4, it is also evident that the lowest total number of new infections occurs when we have strong asymmetric coupling, Γ = 5.12 × 10 3 .As observed in Table 3, the highest total number of new infections in the community will occur under strong symmetric coupling, Γ = 5.79 × 10 3 .Table 4.The total number of newly infected animals over a ten-year period and the total cost J with respect to the control strategy under Scenario 2.    Figures 6-9 demonstrate the impact of optimal culling under all possible coupling cases.As shown in Figures 6-9, the total number of infected animals per patch decreases as a result of the optimal policy.Figure 10 shows the optimal control profiles for controls u 1 (t) and u 2 (t) with low cost parameters.As we can observe, both u 1 and u 2 start from the maximum initially, and stay there for a long time before they switch to the minimum just before the final time horizon.

Discussion
We have provided a mathematical framework to investigate the role of short-term animal dispersal on transmission and control of brucellosis in a heterogeneous population.The proposed model comprises two-patches and animal dispersal has been modeled using a Lagrangian approach.Our study is applicable in communal lands where animal mobility is highly uncontrolled.Hence, it is well known that a single herd of livestock in these communities can be exposed to a highly variable number of contacts with others herds of livestock for a short time frame.This heterogeneity in animal contacts may contribute significantly to the transmission and control of brucellosis.
The basic reproduction number R 0 of the proposed model was computed and analyzed.We observed that it is a function of several factors such as the transmission rates, natural mortality rate, proportions of vertical transmission and the proportion of time that animals of each patch spend in their patch and the other patch.Precisely, we found that R 0 depends on the characteristics of both patches.However, in the absence of animal mobility we observed that each patch has its own reproduction number R 0i i = 1, 2, which depends entirely on the characteristics of that patch.With the aid of model parameter values and initial population levels in [7], we demonstrated numerically that, whenever there is no animal mobility R 01 = 1.4 and R 02 = 0.04, which implies that the disease dies out in low risk patch (Patch 2) and persists in high risk patch (Patch 1).However, with animal mobility incorporated, we noted that R 0 will always be greater than 2 demonstrating that animal mobility will increase the spread of the disease in the community.In particular, we observed that R 0 will be highest when the coupling intensity is weak and the mobility pattern is symmetric, R 0 = 3.03.Analytical methods were also used to demonstrate that, when R 0 ≤ 1, the brucellosis dies out in the community; and when R 0 > 1, a unique endemic equilibrium exists and the disease is uniformly persistent.
Meanwhile, we applied optimal control theory to the proposed model to identify optimal culling strategies that can lead to effective control of brucellosis in the community.Two controls representing culling of infectious animals in each patch were incorporated into the original model.Two possible scenarios that characterize disease control in developing nations were evaluated.Scenario 1 entails no control (we set u 1 = 0) in high risk patch while control is above average (we set u 2 = 0.8) among the low risk population.We hypothesized that this scenario mirrors livestock farming in areas that are in proximity to wildlife.Due to the unavailability of resources in most developing nations, it follows that control of brucellosis among wildlife is less prioritized.In Scenario 2, we set u 1 = 0.45 and u 2 = 0.8.We also suggested that this scenario may represent two herds of livestock that belong to two different farmers who share grazing lands.One farmer may have some financial resources to maintain culling at an intensity above average while the other does not have enough financial capacity to do so.
Under Scenario 1, we observed that the lowest and highest total number of new infections will be recorded in the community under weak symmetric coupling and strong symmetric coupling, respectively.Meanwhile we observed that by introducing a control in high risk patch, the total number of new infections decreases by 30.1%, 21.4% and 28.9% in Patch 1 only, Patch 2 only and overall (Patch 1 and Patch 2 combined), respectively.The numerical results provided evidence that, as expected, controlling the two patches gives the best reduction in brucellosis prevalence.Our result show that animal mobility plays an important role in shaping the long term dynamics of brucellosis, which subsequently impacts the design of its optimal control strategies.
Several avenues for future research arise from this work.First, future research should assess the role of seasonal variations and short-term animal mobility on the persistence of brucellosis.Seasonal availability of water and pastures have a significant influence on pastoral farming, hence there is need to investigate its impact on the persistence of brucellosis.Second, although we were able to establish the uniqueness and uniform persistence result for the endemic equilibrium, we did not resolve the stability of this equilibrium point analytically and that remains an interesting topic for our future research.
Using the next-generation matrix notations in [12], the non-negative matrix F that represents the generation of new infection and the non-singular matrix V that denotes the disease transfer among compartments, are, respectively, given by Then, R 0 , which corresponds to the dominant eigenvalue of the matrix F V −1 , is given by where F and V are defined in Equation (A2).Motivated by [17], we define a Lyapunov function as follows U = w T V −1 Y.
Differentiating U along solutions of Equation (1), we have It can be easily verified that the largest invariant subset of Ω where U = 0 is the singleton {E 0 }.Therefore, by LaSalle's invariance principle [18], E 0 is globally asymptotically stable in Ω when R 0 ≤ 1.
If R 0 > 1, then by continuity, U > 0 in a neighbourhood of E 0 in Ω. Solutions in Ω sufficiently close to E 0 move away from the DFE, implying that the DFE is unstable.In the following, we demonstrate that, if R 0 > 1, then the disease persists and a unique endemic equilibrium point exists.

Appendix C. Uniform Persistence
Proof of Theorem 2.
Therefore, a global attractor for ψ t exists.The disease-free equilibrium is the unique equilibrium on the manifold ∂X 0 and is globally asymptotically stable on ∂X 0 .Moreover, ∪ x∈M ∂ ω(x) = {E 0 } and no subset of M forms a cycle in ∂X 0 .Finally, since the disease-free equilibrium is unstable on X 0 if R 0 > 1, we deduce that the system in Equation ( 1) is uniformly persistent by using a result from [19] (Theorem 1.3.1 and Remark 1.3.1).This completes the proof of Theorem 2.

Appendix D. Existence of a Unique Endemic Equilibrium Point
Proof of Theorem 3. We can reduce the system in Equation ( 1) into a four-dimensional system by setting R We use a result by Hethcote and Thieme in [20] to prove the uniqueness of the endemic equilibrium.An endemic equilibrium (S * i , I * i ) satisfies: The first part of Equation (A5) gives Hence, from the last part of Equation (A5), we deduce that Since J 2 (x) ≥ 0 and J 3 (x) ≥ 0, all off-diagonal entries of the Jacobian matrix are nonnegative, thus the function H(x) is monotone.Therefore, monoticity of a matrix H(x) implies that the model in Equation (1) has a unique positive fixed point if and only if R 0 > 1.This completes the first part of the proof for Theorem 3 and due to less traceability of our model we will utilizing numerical simulations to demonstrate the global stability of the endemic equilibrium (see Figure A1).= 1000) over time for each i.Thus, animals of Patch 1 spends, on the average, the proportion p 11 of their time in residency in Patch 1 and the proportion p 12 of their time in Patch 2 such that p 11 + p 12 = 1.Similarly, animals of Patch 2 spend the proportion p 22 of their time in Patch 2 and p 21 = 1 − p 22 in Patch 1.Therefore, at time t, the effective population in Patch 1 is p 11 N 1 + p 21 N 2 while the effective population of Patch 2, at time t is p 12 N 1 + p 22 N 2 .Susceptible animals of Patch 1 (S 1 ) could be infected contagiously, in Patch 1 (if currently in Patch 1, that is, p 11 S 1 ) or in Patch 2 (if currently in Patch 2, that is, p 12 S 1 ).It follows from the above discussion that the effective proportion of infectious individual in Patch 1 is

Figure 1 .
Figure 1.Simulation results of the proposed two patch brucellosis model for scenario 1 under weak symmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 2 .Figure 3 .
Figure 2. Simulation results of the proposed two patch brucellosis model for scenario 1 under strong symmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 4 .
Figure 4. Numerical illustrations depicting the effects of optimal intervention strategies on controlling the long-term brucellosis dynamics for scenario 1 under strong asymmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 5 .
Figure 5.The control profile for Scenario 1: (a) low cost of culling; and (b) high cost of culling.

Figure 6 .Figure 7 .
Figure 6.Simulation results of the proposed brucellosis model demonstrating the disease dynamics for scenario 2 under weak symmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 8 .
Figure 8. Numerical results highlighting the impact of optimal intervention strategies on controlling the long-term brucellosis dynamics for scenario 2 under weak asymmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 9 .
Figure 9. Numerical illustrations demonstrating the impact of optimal intervention strategies on controlling the long-term brucellosis dynamics for scenario 2 under strong asymmetric coupling (a) the numbers of infected animals in patch 1 (b) the numbers of infected animals in patch 2. In all the figures the dotted blue and solid black curves represent the infected population, without and with control, respectively.

Figure 10 .
Figure 10.The control profile for Scenario 2.

Figure A1 .
Figure A1.Phase portrait illustrating the global stability of E * for the system in Equation (1) in the S 1 -I 1 plane with R 0 = 2.84 (we set β 1 = β 2 = 1.5).Each curve in the plot corresponds to a different initial condition, and all these curves converge to the equilibrium E * (where S 1 = S 2 .=1500, I 1 = I 2 .=1000) over time

Table 1 .
Parameters and values.