Loading [MathJax]/jax/output/SVG/jax.js
Research article

Discrete Dirac reduction of implicit Lagrangian systems with abelian symmetry groups

  • This paper develops the theory of discrete Dirac reduction of discrete Lagrange–Dirac systems with an abelian symmetry group acting on the configuration space. We begin with the linear theory and, then, we extend it to the nonlinear setting using retraction compatible charts. We consider the reduction of both the discrete Dirac structure and the discrete Lagrange–Pontryagin principle, and show that they both lead to the same discrete Lagrange–Poincaré–Dirac equations. The coordinatization of the discrete reduced spaces relies on the notion of discrete connections on principal bundles. At last, we demonstrate the method obtained by applying it to a charged particle in a magnetic field, and to the double spherical pendulum.

    Citation: Álvaro Rodríguez Abella, Melvin Leok. Discrete Dirac reduction of implicit Lagrangian systems with abelian symmetry groups[J]. Journal of Geometric Mechanics, 2023, 15(1): 319-356. doi: 10.3934/jgm.2023013

    Related Papers:

    [1] Bibi Fatima, Mehmet Yavuz, Mati ur Rahman, Fuad S. Al-Duais . Modeling the epidemic trend of middle eastern respiratory syndrome coronavirus with optimal control. Mathematical Biosciences and Engineering, 2023, 20(7): 11847-11874. doi: 10.3934/mbe.2023527
    [2] Chentong Li, Jinhu Xu, Jiawei Liu, Yicang Zhou . The within-host viral kinetics of SARS-CoV-2. Mathematical Biosciences and Engineering, 2020, 17(4): 2853-2861. doi: 10.3934/mbe.2020159
    [3] Ugo Avila-Ponce de León, Angel G. C. Pérez, Eric Avila-Vales . Modeling the SARS-CoV-2 Omicron variant dynamics in the United States with booster dose vaccination and waning immunity. Mathematical Biosciences and Engineering, 2023, 20(6): 10909-10953. doi: 10.3934/mbe.2023484
    [4] Tahir Khan, Roman Ullah, Gul Zaman, Jehad Alzabut . A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator. Mathematical Biosciences and Engineering, 2021, 18(5): 6095-6116. doi: 10.3934/mbe.2021305
    [5] Salman Safdar, Calistus N. Ngonghala, Abba B. Gumel . Mathematical assessment of the role of waning and boosting immunity against the BA.1 Omicron variant in the United States. Mathematical Biosciences and Engineering, 2023, 20(1): 179-212. doi: 10.3934/mbe.2023009
    [6] Mohanad A. Deif, Ahmed A. A. Solyman, Mehrdad Ahmadi Kamarposhti, Shahab S. Band, Rania E. Hammam . A deep bidirectional recurrent neural network for identification of SARS-CoV-2 from viral genome sequences. Mathematical Biosciences and Engineering, 2021, 18(6): 8933-8950. doi: 10.3934/mbe.2021440
    [7] Hamdy M. Youssef, Najat A. Alghamdi, Magdy A. Ezzat, Alaa A. El-Bary, Ahmed M. Shawky . A new dynamical modeling SEIR with global analysis applied to the real data of spreading COVID-19 in Saudi Arabia. Mathematical Biosciences and Engineering, 2020, 17(6): 7018-7044. doi: 10.3934/mbe.2020362
    [8] Junyuan Yang, Guoqiang Wang, Shuo Zhang . Impact of household quarantine on SARS-Cov-2 infection in mainland China: A mean-field modelling approach. Mathematical Biosciences and Engineering, 2020, 17(5): 4500-4512. doi: 10.3934/mbe.2020248
    [9] Yuta Okada, Taishi Kayano, Asami Anzai, Tong Zhang, Hiroshi Nishiura . Protection against SARS-CoV-2 BA.4 and BA.5 subvariants via vaccination and natural infection: A modeling study. Mathematical Biosciences and Engineering, 2023, 20(2): 2530-2543. doi: 10.3934/mbe.2023118
    [10] Fang Wang, Lianying Cao, Xiaoji Song . Mathematical modeling of mutated COVID-19 transmission with quarantine, isolation and vaccination. Mathematical Biosciences and Engineering, 2022, 19(8): 8035-8056. doi: 10.3934/mbe.2022376
  • This paper develops the theory of discrete Dirac reduction of discrete Lagrange–Dirac systems with an abelian symmetry group acting on the configuration space. We begin with the linear theory and, then, we extend it to the nonlinear setting using retraction compatible charts. We consider the reduction of both the discrete Dirac structure and the discrete Lagrange–Pontryagin principle, and show that they both lead to the same discrete Lagrange–Poincaré–Dirac equations. The coordinatization of the discrete reduced spaces relies on the notion of discrete connections on principal bundles. At last, we demonstrate the method obtained by applying it to a charged particle in a magnetic field, and to the double spherical pendulum.



    Middle Eastern respiratory syndrome (MERS), a respiratory virus caused by MERS-CoV coronavirus, first emerged in Jeddah, Saudi Arabia, on June 13, 2012 [1,2,3,4]. The disease, which is of the same genealogy as SARS-CoV-1 (which emerged in 2002–2003) and SARS-CoV-2 (which emerged in 2019) (see Table 1 for a comparison of the three human coronaviruses), rapidly spread to numerous countries (including Algeria, Austria, Bahrain, China, Egypt, France, Germany, Greece, Iran, Italy, Jordan, Kuwait, Lebanon, Malaysia, Netherlands, Oman, Philippines, Qatar, Republic of Korea, Thailand, Tunisia, Turkey, United Kingdom, United Arab Emirates, United States, and Yemen) [2,5]. Data from the World Health Organization report for May 2023 showed that MERS caused 2604 confirmed cases and 936 deaths (representing a 36% case-fatality ratio) globally (as of June 5, 2023). Furthermore, the Kingdom of Saudi Arabia (KSA) and the Republic of Korea recorded the majority of the global MERS-CoV cases (in particular, the KSA reported 2196 cases and 855 deaths) [6,7] (Figure 1 depicts the number of weekly new cases of MERS-CoV for the KSA). The KSA is vulnerable to MERS-CoV-2 outbreaks due to the numerous large gathering events it hosts, notably the annual Hajj pilgrimage (once a year) and Umrah year [8], and the camel festival in Riyadh (which runs for a two-month period) [9,10].

    Table 1.  Comparison of epidemiological and clinical characteristics of novel coronaviruses [34,38,39].
    Characteristics SARS-CoV-1 MERS-CoV SARS-CoV-2
    Origin location Guandong, China Jeddah, Saudi Arabia Wuhan, China
    Duration 2002–2003 2012–to date Dec. 2019–to date
    Reservoir Bats & civet cats Bats & camels Bats (Pangolins)
    Countries 29 27 220
    Incubation period 2–7 days 5 days 2–14 days
    Confirmed cases 8098 2604 676,609,955
    Global deaths 774 936 6,881,955
    Cases-fatality ratio 9.5% 36% 1%

     | Show Table
    DownLoad: CSV
    Figure 1.  Weekly new cases of MERS-CoV in the Kingdom of Saudi Arabia for the period from the second week of June 2012 to the fourth week of June 2019 [6].

    Several empirical studies have shown that bats and dromedary camels are the main reservoir hosts of MERS-CoV [11,12,13,14,15] (see Figure 2). Transmission from camel to camel or camel to human could occur by contact with the animal's nasal secretions, saliva, or respiratory droplets [16,17]. MERS-CoV spreads among dromedary camels via respiratory routes [18], and dromedary camels infected with MERS-CoV exhibit symptoms (typically manifesting within two days and continue for up to two weeks [19]) which include mucopurulent nasal and lachrymal discharge, appetite loss, fever, rhinorrhea (runny nose), coughing, and sneezing [19]. The disease is primarily transmitted to humans by infected dromedary camels [2,11,17,20,21,22] through inhalation of respiratory droplets and airborne particles. The virus could also be transmitted to humans after drinking the raw milk or eating the meat of infected camels [17,19,23]. Other potential sources of MERS-CoV transmission from camels to humans include crowding (i.e., camel-to-camel or camel-to-human), mixing camels from different regions, transportation of camels, and markets for trade of live animals [17,21,24]. Furthermore, camel racing and camel show festivals (where camels are kept in a large space for 2 months [18]) increases the likelihood of MERS-CoV transmission between camels and humans [18,19,25] (MERS-CoV can be detected in viral RNA in camels for up to 35 days post-inoculation [26,27]). Furthermore, although the mechanisms for camel-to-human transmission of MERS-CoV are well-understood, there is a general lack of knowledge of the mechanism resulting in human-to-camel transmission of MERS-CoV [17,18,28,29]. MERS-CoV transmission from human to human is well-documented in family clusters, community settings, and more frequently in healthcare settings [19,30,31,32,33]. Finally, MERS-CoV can also be transmitted through direct contact with contaminated objects and surfaces [34,35,36,37].

    Figure 2.  Zoonotic origins and transmission pathways of novel human coronaviruses [43].

    The main non-pharmaceutical interventions (NPIs) implemented against MERS-CoV in the KSA were focused on preventing human-to-human transmission (notably reducing the likelihood of inhaling infected respiratory droplets and minimizing contacts with infected humans [40]) and minimizing or avoiding contacts with infected camels [29,41]. Consequently, the main NPIs implemented against MERS-CoV in the KSA include the use of face masks by humans who are within close proximity (typically a meter) of MERS-infected individuals and avoiding contact with MERS-infected patients during visitations in a healthcare setting [40]. Furthermore, individuals who have close contact with camels are encouraged to wear proper personal protective equipment (PPE), especially N95 masks, surgical gowns, and goggles [18,42].

    To reduce the likelihood of MERS-CoV transmission through the consumption of infected milk, the Saudi Ministry of Health strongly recommended that camel milk should be thoroughly boiled or pasteurized before use. Similarly, it is strongly advised that camel meat must be thoroughly cooked before consumption [18,42]. Healthcare workers are at a high risk of acquiring a MERS-CoV infection (they account for 31% of transmission [44,45]). Thus, healthcare workers are strongly encouraged to wear N95 respirators (data showed a marked reduction, by at least 56%, in sero-positivity, in comparison to healthcare workers who do not, or rarely, wear N95 masks [42]). The importation and transportation of dromedary camels between nations is another major risk factor for MERS-CoV transmission [19,42]. Camels destined for exportation are tested in the originating country for the presence of MERS-CoV antibodies and antigens [19,42]. Upon arrival at the destination country, the camels are typically placed in quarantine while being sampled and tested for MERS-CoV. Before being allowed to be imported into the country, camels that have tested positive for MERS-CoV or are sero-positive are isolated at border crossings. Camels brought for other objectives, such as breeding and racing competitions, also represent a high risk of MERS-CoV transmission to the indigenous camel breeds (because they will eventually be mixed with the indigenous herds). Camels suspected of being infected with MERS-CoV are transported internally in properly equipped, closed, and ventilated vehicles. Furthermore, in general, camels are not crammed too closely together in a single vehicle [42].

    Although there is currently no authorized antiviral medication or vaccine to protect humans from MERS-CoV [2,27,46,47,48,49,50,51], there, however, are two candidate anti-MERS human vaccines (produced by Novavax and Greffex Inc.) currently undergoing advanced stages of clinical trials [2,48,52]. Additionally, some potentially promising vaccines are being developed for use in camels [53]. For example, the biotechnology company INOVIO, in collaboration with the Coalition for Epidemic Preparedness Innovations (CEPI), developed a DNA-based human MERS-CoV vaccine that is currently at the second stage of clinical trials [54]. Furthermore, INOVIO, in collaboration with GeneOne Life Science Inc., developed a DNA-based vaccine [53] that is also undergoing the second stage of clinical trials in humans [51,55]. Finally, Johnson & Johnson, in collaboration with CEPI and the Jenner Institute at the University of Oxford, is currently designing an adenovirus-based vaccine against MERS-CoV in humans [56]. Similarly, although there are currently no vaccines against MERS-CoV in camels [2,27,46,47,48,49,50,51], some promising candidate vaccines are undergoing clinical trials. These include: recombinant DNA and protein-based vaccines, subunit vaccines, viral vector-based vaccines (including vaccines designed based on using the measles viral vector, Poxviral vector, human adenoviral vectors, and Chimpanzee adenoviral vectors) [53]. INOVIO, in collaboration with GeneOne Life Science Inc., developed a DNA-based vaccine [53] that was tested in dromedary camels, showing encouraging protective outcomes against the disease [51,53]. An adenovirus-based vaccine developed by researchers from the University of Oxford is being tested in dromedary camels and humans[47,51,53,57].

    Mathematical models of varying types (such as compartmental [2,28,58,59], statistical [60,61], and stochastic [15,62,63,64]) have been developed and used to study the transmission dynamics and control of MERS-CoV in camel, human, or human and camel populations. For example, compartmental models were formulated and used to study the transmission dynamics of MERS-CoV among humans, camels, and between humans and camels [2,28,31,59,65,66]. One of the earliest MERS-CoV models was a two-group compartmental model developed by Malik [2] to study the dynamics of the disease in a mass gathering of humans (specifically, during the Hajj and Umrah gatherings in the Kingdom of Saudi Arabia [2]). The study proposed optimal strategies, based on the quarantine of suspected cases, the isolation of symptomatic individuals, and the use of a potential future vaccine, for effectively combating or mitigating the disease during those mass gathering events. Ghosh and Nadim [28] used a compartmental model of two co-circulating zoonotic MERS-CoV strains to study the spread of the disease in both humans and camels. This study emphasized preventive measures for susceptible humans and the isolation of infected camels. Chowell et al.[15] evaluated the multiple MERS-CoV transmission routes in humans using a stochastic model. The study, which used epidemiological data for the MERS-CoV outbreak that occurred in the KSA during April to October of 2013 [15], discussed the uncertainty in epidemic risk while taking observation bias in their consideration. They match the evaluation of MERS-CoV cases during that period with a dynamic transmission model that included community and hospital compartments.

    Since the KSA attracts such mass gatherings from around the world every year (and year-around), the main motivation of the current study is that effectively controlling outbreaks of MERS-CoV in the KSA could prevent the likelihood of a major global pandemic starting from the KSA. The objective of this study is to use mathematical modeling approaches, together with data analytics and computation, to study the transmission dynamics of the MERS-CoV endemic in the KSA. Specifically, the objective will be achieved by, first, using a basic deterministic model that captures the essential epidemiological and biological features of the disease. The basic model, which will be parameterized using observed cumulative MERS-CoV case data for the KSA, will be rigorously analyzed to gain insight into its qualitative features (with respect to the existence and asymptotic stability of its associated equilibria). Numerical simulations will be carried out to assess the population-level impacts of various control and mitigation strategies. The basic model will then be extended to allow for the assessment of the population-level impact of public health intervention (pharmaceutical and non-pharmaceutical), notably face mask usage and the vaccination of both humans and camels. The extended model will also be rigorously analyzed and simulated.

    The paper is organized as follows. The basic model (2.1) is formulated, and its basic qualitative properties analyzed, in Section 2.1. The basic model is also fitted with the cumulative cases data for the Kingdom of Saudi Arabia in Section 3.2. The basic model is rigorously analyzed, with respect to the existence and asymptotic stability of its associated equilibria, in Section 3. Furthermore, the basic model is extended in Section 4.1 to incorporate the aforementioned intervention and control strategies. The extended model is rigorously analyzed and simulated in Sections 4.2 and 5, respectively.

    To develop the basic mathematical model for the transmission dynamics of MERS-CoV within and across human and dromedary camel populations in the KSA, the total human population in the KSA at time t, denoted by NH(t), is split into the mutually exclusive compartments of susceptible (SH(t)), exposed (EH(t)), asymptomatically infectious (IHA(t)), symptomatically infectious (IHS(t)), hospitalized (IHH(t)), and recovered (RH(t)) humans, so that:

    NH(t)=SH(t)+EH(t)+IHA(t)+IHS(t)+IHH(t)+RH(t).

    Similarly, the total population of dromedary camels in the KSA at time t, denoted by NC(t), is subdivided into the sub-populations of susceptible (SC(t)), exposed (EC(t)), infectious (IC(t)), and recovered (RC(t)) camels. Hence,

    NC(t)=SC(t)+EC(t)+IC(t)+RC(t).

    The equations for the dynamics of MERS-CoV in a humans-camels community is given by the following deterministic system of nonlinear differential equations [4]:

    ˙SH=ΠHλHSHλCHSHμHSH,˙EH=λHSH+λCHSH(σH+μH)EH,˙IHA=rσHEH(γHA+μH+δHA)IHA,˙IHS=(1r)σHEH(ϕHS+γHS+μH+δHS)IHS,˙IHH=ϕHSIHS(γHH+μH+δHH)IHH,˙RH=γHAIHA+γHSIHS+γHHIHHμHRH,˙SC=ΠC+ψCRCλCSC(λHCA+λHCS)SCμCSC,˙EC=λCSC+(λHCA+λHCS)SC(σC+μC)EC,˙IC=σCEC(γC+μC+δC)IC,˙RC=γCIC(ψC+μC)RC, (2.1)

    where the infection rates for MERS CoV-2 transmission from human to human (λH), camel to camel (λC), camel to human (λCH), asymptomatically infectious human to camel (λHCA), and symptomatically infectious humans to camel (λHCS) are given, respectively, by:

    λH=βHH[ηHAIHA+IHS+ηHHIHH]NH, λC=βCCICNC, λCH=βCHICNC, λHCA=βHC ηHAIHANH, λHCS=βHCIHSNH. (2.2)

    In the model (2.1), ΠH (ΠC) is the recruitment of humans (camels) into the population (all newly recruited humans and camels are assumed to be susceptible to MERS-CoV), and βHH (βCC) represents the rate at which MERS-CoV is transmitted to susceptible humans (camels). βHC (βCH) represents the rate at which MERS-CoV is transmitted from humans to camels (camels to humans). The parameter μH (μC) represents the natural death rate of humans (camels). We differentiate the transmission probability of symptomatic and hospitalized individual with the respective weights, 0<ηHA,ηHH<1. Exposed humans (camels) become infectious at a rate of σH (σC) (i.e., 1/σH (1/σC) is the average duration from the initial infection to infectiousness for humans (camels)). The proportion r of these humans who become infectious after surviving the exposed class show no clinical symptoms of MERS-CoV, and move to the asymptomatically infectious class IHA (at a rate r σH), and the remaining proportion, 1r, show clinical symptoms of MERS-CoV and move to the symptomatic class IHS (at a rate (1r)σH). The parameter γHA(γHS)(γHH)(γC) represents the rate at which asymptomatic (symptomatic)(hospitalized) humans (camels) recover from MERS-CoV infection. The parameter δHA(δHS)(δHH)(δC) is the disease induced mortality rate for asymptomatic (symptomatic) (hospitalized) humans (camels). Let DH(t) and DC(t) represent the total number of humans and camels, respectively, who died of MERS-CoV at time t. Hence, the rates of change of the deceased human and camel populations are given, respectively, by:

    ˙DH=δHAIHA+δHSIHS+δHHIHH, (2.3)

    and

    ˙DC=δCIC. (2.4)

    A flow diagram of the model (2.1) is depicted in Figure 3, and the state variables and parameters of the model are described in Tables 5 and 6, respectively [4].

    Figure 3.  Flow diagram of the model (2.1), where the infection rates λH, λC, λCH, λHCA, and λHCS are as given in (2.2) [4].

    The main assumptions made in the formulation of the model (2.1) are:

    1) Exponentially distributed waiting times in each epidemiological compartment for humans and camels [15,67].

    2) Well-mixed homogeneous human and camel populations: it is assumed that the populations for humans and camels are well-mixed, and that every member of the community (human or camel) has equal likelihood of mixing with every other member of the community. Furthermore, every member has equal probability of acquiring and/or transmitting infection to hosts of the same or different types (i.e., human to camel and camel to human) [2,15].

    3) It is assumed that between-species transmission occurs at a reduced rate, in comparison to within-species transmission (i.e., humans (camels) are assumed to transmit to camels (humans) at a lower rate, in comparison to human-to-human (camel-to-camel) transmission [28,66].

    4) Hospitalized humans do not transmit infection to camels.

    5) Humans who have recovered from MERS-CoV acquire lasting (and permanent) natural immunity against re-infection [2,28]. However, recovered camels can lose their infection-acquired immunity and become susceptible again [63].

    6) Since MERS-CoV has been in circulation for more than ten years in the Kingdom of Saudi Arabia, it is assumed to be endemic (so that the demographic (birth and natural death parameters) processes are accounted for in the model).

    The basic qualitative properties of the model (2.1) will now be assessed.

    Since the model (2.1) monitors the temporal dynamics of the populations of humans and camels, all its associated parameters and state variables are non-negative for all t0. We claim the following non-negativity result for the model:

    Theorem 2.1. Consider the closed sets

    DH={(SH,EH,IHA,IHS,IHH,RH)R6+}andDC={(SC,EC,IC,RC)R4+}.

    The region

    D=DHDC

    is positively invariant and attracting with respect to the model (2.1).

    Proof. Adding the first six and the last four equations of the model (2.1) gives, respectively,

    ˙NH=ΠHμHNH(δHAIHA+δHSIHS+δHHIHH),˙NC=ΠCμCNCδCIC. (2.5)

    Since all of the parameters of the model (2.1) are non-negative, it follows from (2.5) that

    ˙NH(t)ΠHμHNH(t), (2.6)

    and

    ˙NC(t)ΠCμCNC(t). (2.7)

    NC(t)ΠC/μC, respectively. Consequently, it follows from (2.6), using the comparison theorem [68], that:

    NH(t)NH(0)eμH(t)+ΠH/μH[1eμH(t)].

    Similarly, it follows from (2.7) that:

    NC(t)NC(0)eμC(t)+ΠC/μC[1eμC(t)].

    In particular, NH(t)ΠH/μH if NH(0)ΠH/μH and NC(t)ΠC/μC if NC(0)ΠC/μC, respectively. Thus, if NH(0)>ΠH/μH, then either the solution of the human component of the model enters DH in finite time, or NH(t) approaches ΠH/μH and the infected variables EH(t),IHA(t),IHS(t), and IHH(t) approach zero as t [69]. Hence, the feasible region DH is invariant and attracting (i.e., all solutions of the human component of the model lie in, or eventually enter, R6+). Similarly, it can be seen that if NC(0)>ΠC/μC, then either the solution of the camel component of the model enters DC in finite time, or NC(t) approaches ΠC/μC (and the infected variables EC(t) and IC(t) approach zero as t [69]). Hence, DC is also positively invariant and attracting. Since the union and intersection of invariant sets is also invariant [70,71], it follows that D is positively invariant and attracting with respect to the model.

    The epidemiological consequence of Theorem 2.1 is that the model (2.1) is well-posed and mathematically and epidemiologically in the region D (hence, it is sufficient to consider the dynamics of the flow generated by the model (2.1) in D [72]).

    The model (2.1) has a unique disease-free equilibrium (DFE) given by:

    E0=(SH,EH,IHA,IHS,IHH,RH,SC,EC,IC,RC)=(ΠH/μH,0,0,0,0,0,ΠC/μC,0,0,0). (3.1)

    It follows, using the next-generation operator method [73,74], that the DFE of the model is locally asymptotically stable if the spectral radius of the matrix FV1 is less than one, where matrices F (of the new infection terms) and V (of the linear transition terms in the infected compartments) are given, respectively, by:

    F=[ 0βHHηHA(SHNH)βHH(SHNH)βHHηHH(SHNH)0βCH(SHNC) 000000 000000 000000 0βHCηHA(SCNH)βHC(SCNH)00βCC(SCNC) 000000]

    and

    V=[K100000rσHK20000(1r)σH0K300000ϕHSK4000000K500000σCK6],

    where K1=σH+μH, K2=γHA+μH+δHA, K3=ϕHS+γHS+μH+δHS, K4=γHH+μH+δHH, K5=σC+μC, K6=γC+μC+δC, and NH=ΠH/μH. It is convenient to define R0M, the reproduction number of the model (2.1).

    R0M=ρ(FV1) =12{RC+RH+(RCRH)2+4 (RHAC +RHSC)}, (3.2)

    where

    RH=βHH(K3K4rσHηHA+K2(1r)σH(K4+ηHHϕHS)K1K2K3K4),RC=βCC(σCK5K6), RHAC=βHCβCHrσHηHAσCK1K2K5K6, RHSC=βHCβCH(1r)σHσCK1K3K5K6. (3.3)

    The basic reproduction number R0M, of model (2.1), represents the average number of new cases (in humans or camels) generated by a typical infectious human or camel if introduced in a completely susceptible population of both humans and camels. The result below follows from Theorem 2 of [73].

    Theorem 3.1. The disease-free equilibrium of the model (2.1) is locally asymptotically stable if R0M <1, and unstable if R0M >1.

    The epidemiological implication of Theorem 3.1 is that a small influx of infected humans or camels in the well-mixed human-camel population (when R0M <1) will not generate a significant MERS-CoV outbreak in the human-camel community if the initial number of infected humans or camels is small enough (i.e., in the basin of attraction of the aforementioned disease-free equilibrium). The terms in the expression for the basic reproduction number (R0M) of the model are epidemiologically interpreted below.

    The terms in the expression for the basic reproduction number (R0M) for the model (2.1), given by (3.2), will be epidemiologically interpreted below. The expression for R0 contains four terms, namely RH, RC, RHAC, and RHSC. These terms will now be epidemiologically interpreted (it should be noted that SH=NH=ΠH/μHand SC=NC=ΠC/μC).

    Interpretation of RH

    The quantity RH is associated with the transmission dynamics of MERS-CoV between humans (only) in the community. It is convenient to re-write RH in the following form:

    RH=RHA+RHS+RHH,

    where

    RHA=βHHηHArσHK1K2, RHS=βHH(1r)σHK1K3, RHH=βHHηHH(1r)σHϕHSK1K3K4.  (3.4)

    That is, the quantity RH is now expressed as the sum of the three constituent reproduction numbers for the average number of new cases generated by infectious humans in the asymptomatic (RHA), symptomatic (RHS), and hospitalized (RHH) compartments.

    The terms in the constituent reproduction number RHA can be expressed as follows. They are given by the product of the infection rate of susceptible humans with asymptomatically infectious humans near the DFE, given by (βHHηHASHNH), the proportion that survived the exposed class and moved to the asymptomatically infectious class (rσH/K1), and the average duration in the asymptomatically infectious class (1/K2). Similarly, RHS is the product of the infection rate of symptomatic humans near the DFE (βHHSHNH), the proportion that survived the exposed class and moved to the symptomatically infectious class ((1r) σHK1), and the average duration in the symptomatically infectious class (1/K3). Finally, RHH is the product of the infection rate of hospitalized infectious humans near the DFE (βHHηHHSHNH), the proportion that survived the exposed class and moved to the symptomatically infectious class ((1r)σHK1), the proportion that survived the symptomatically infectious class and moved to the hospitalized class (ϕHS/K3), and the average duration in the hospitalized class (1/K4). The sum of the constituent reproduction numbers, RHA, RHS, and RHH, gives RH.

    Interpretation of RC

    The quantity RC is associated with the transmission dynamics of MERS-CoV among camels (only) in the community. The terms in RC represent the product of the infection rate of susceptible camels with infected camels near the DFE given by (βCCSCNC), the proportion that survived the exposed class and moved to the infected camel class (σC/K5), and the average duration in the infectious class (1/K6).

    Interpretation of RHAC

    This quantity is associated with the asymptomatically infectious humans camels asymptomatically infectious humans MERS-CoV transmission cycle. It is convenient to re-write it as:

    RHAC=RhacRCH,

    where

    Rhac=(βHCηHArσHSCK1K2NH)andRCH=(βCHSHK5K6NC)

    with Rhac representing a measure of MERS-CoV transmission from asymptomatically infectious humans to susceptible camels and RCH accounting for MERS-CoV transmission from infectious camels to susceptible humans. The quantity Rhac represents the product of the infection rate of susceptible camels by asymptomatically infected humans (βHCηHASCNH), the proportion of infected humans that survived the exposed class and moved to the asymptomatically infectious IHA class, (rσHK1), and the average duration in the IHA class (1K2). Similarly, the quantity RCH is the product of the infection rate of susceptible humans by infectious camels, (βCHSHNC), the proportion of infected camels that survived the exposed class and moved to the symptomatic class IC, (σCK5), and the average duration in the IC class (1K6).

    Interpretation of RHSC

    This threshold quantity is associated with the symptomatic humans-camels-symptomatic humans MERS-CoV transmission cycle. It can be re-written as:

    RHSC=RhscRCH,

    where Rhsc is a measure of MERS-CoV transmission from symptomatically infectious humans to susceptible camels and RCH accounts for the MERS-CoV transmission from infectious camels to susceptible humans. Specifically, the quantity Rhsc is the product of the infection rate of susceptible camels by MERS-CoV symptomatic-infected humans (βHCSCNH), the proportion of infected humans who survived the exposed class and moved to the symptomatically infectious IHS class, ((1r)σHK1), and the average duration in the symptomatically infectious IHS class, (1K3). Similarly, the quantity RCH is the product of the infection rate of susceptible humans by infectious camels, (βCHSHNC), the proportion of infected camels that survived the exposed class and moved to the IC class, (σCK5), and the average duration in the IC class (1K6). The combinations of RH, RC, RHAC, and RHSC, in the way described in (3.2), gives R0M.

    It should be mentioned that the square root in the expression for R0M, given by (3.2), is a consequence of the two generations (humans-camels-humans) needed to complete the cycles for generating RHAC and RHSC. Furthermore, the numbers 1/2 and 4 in the expression (3.3) appeared because of the weighted averaging of the two modes (human-to-camel and camel-to-human) of MERS-CoV transmission [75]. It is also worth mentioning that, using the baseline values of the fixed and fitted parameters of the model (2.1) tabulated in Table 7, the value of the reproduction number (R0M) for MERS-CoV transmission in the Kingdom of Saudi Arabia (using observed cumulative data from Week 13, 2013 to Week 9, 2014) is R0M0.860. Thus, since R0M is less than one, it follows from the result of Theorem 3.1 that MERS-CoV can be effectively controlled in the Kingdom, and that such effective control is dependent on whether or not the initial sizes of the populations of infected humans or camels are small enough (i.e., if the initial numbers of infected hosts are in the basin of attraction of the disease-free equilibrium). Since R0M is less than one for the Kingdom, our study shows that the prospect for MERS-CoV elimination in the Kingdom is promising provided the initial numbers of infected humans and camels are small enough.

    The values of the fitted parameters of the model (2.1) can be used to determine the main drivers of the MERS-CoV pandemic in the Kingdom of Saudi Arabia, by computing the relative proportion of each of the fitted parameters in comparison to the total sum of these parameters). For instance, since the fitted value of the transmission rate for human to human is βHH=0.114 per week, and the sum of the values of all four fitted transmission rates is 0.842, it follows that, during the fitting period (i.e., from Week 13, 2013 to Week 42, 2013), human-to-human transmission accounts for about 14% (i.e., 0.114/0.842 expressed as a percentage) of all new MERS-CoV cases in the human-camel population of the Kingdom of Saudi Arabia. Table 2 depicts the proportion of new MERS-CoV cases generated by the various transmission pathways. It follows from this table that, during the fitting period, human-to-camel transmission accounts for 67% of all new MERS-CoV infections in the human-camel community within the Kingdom. Camel-to-camel transmission accounts for 16% of the new cases, while the camel-to-human transmission pathway accounts for only 3% of new cases.

    Table 2.  Contributions of various transmission pathways on the number of new MERS-CoV cases in the human-camel community of the Kingdom of Saudi Arabia, as measured by the relative proportion of the fitted values of the transmission rate parameters (βHH, βCC, βCH, and βHC) tabulated in Table 7.
    Transmission rate parameter Proportion of new infected cases
    βHH 14%
    βCC 16%
    βCH 3%
    βHC 67%

     | Show Table
    DownLoad: CSV

    The model (2.1) contains 23 parameters. While realistic values of most of them are available in the modeling literature (see Table 7), realistic values of the four contact rate parameters (namely βHH, βCC, βHC, and βCH) are unknown, and need to be estimated from observed data. To do so, we specifically fit the model (2.1) with the observed cumulative case data for MERS-CoV in the Kingdom of Saudi Arabia for the period from Week 13 of 2013 to Week 42 of 2013, and the segment of the data for the period from Week 43 of 2013 to Week 9 of 2014 is used to cross-validate the fitted model [6]. The least squares regression method is used to fit the model with the data [76,77,78]. This method entails finding the optimal parameter region that minimizes the root mean square difference between the predicted cumulative confirmed cases (generated from the model (2.1)) and the observed confirmed data [76,77,78]. The results obtained from the data fitting, depicted in Figure 4, show a reasonably good fit (and cross-validation). The estimated values of the four parameters (generated from the model fitting), and the associated 95% confidence intervals, are tabulated in Table 7.

    Figure 4.  Data fitting of the model (2.1) using cumulative MERS-CoV new cases data for the Kingdom of Saudi Arabia from Week 13 of 2013 to Week 42 of 2013 (blue curve) and cross-validated using the segment of the data from Week 43 of 2013 to Week 9 of 2014 (green curve) [6]. The baseline values and the values of the four fitted parameters in Table 7 and Figure 1 were used in generating this figure.

    The model (2.1) contains 23 parameters, and sensitivity analysis will be conducted to determine the parameters that have the most influence on the transmission dynamics of the disease. The global parameter sensitivity analysis entails using the Latin hypercube sampling (LHS) approach and partial rank correlation coefficients (PRCCs) [79,80,81,82,83] to identify the parameters of the model that have the most influence on the chosen response function, R0 [79,80,84,85]. In order to carry out the sensitivity analysis, a range (lower and upper bound) and distribution must first be established for each parameter of the model. Since the values of the demographic parameters for humans and camels, namely ΠH, ΠC, μH, and μC, are realistically known from census data [28,47,86,87,88], these four parameters are not included in the sensitivity analysis. Parameters with high PRCC values close to -1 or +1 are said to be highly correlated with the response function [12,89]. Those with negative (positive) PRCC values are said to be negatively (positively) correlated with the response function (R0M) [12,89]. The baseline values and ranges of the parameters of the model (2.1) needed for the sensitivity analysis are tabulated in Table 4. The control basic reproduction number, R0M, is chosen as the response function. The results obtained for the global sensitivity analysis of the model (2.1), with respect to the response function R0M, are depicted in Figure 5 (and also tabulated in Table 4). It follows from Figure 5 and Table 4 that the parameters that have the most effect on the response function (R0M) are:

    Table 3.  Baseline values of the four fitted parameters (and their confidence intervals) of the model (2.1).
    Parameter Fitted value (week1) 95% Confidence Interval (week1)
    βHH 0.114 [0.000-0.204]
    βCC 0.135 [0.126-0.143]
    βCH 0.023 [0.023-0.025]
    βHC 0.568 [0.000-0.999]

     | Show Table
    DownLoad: CSV
    Table 4.  PRCC values of the parameters (19 of 23) of the model (2.1), with R0M chosen as the response function. Parameters with PRCCs 0.5 in magnitude are considered to be highly correlated with the response function, and are highlighted in bold font.
    Parameter Baseline value Range PRCC
    βHH 0.114 0.091–0.137 0.216
    βCC 0.135 0.108–0.163 0.921
    βCH 0.023 0.018–0.027 0.119
    βHC 0.568 0.454–0.682 0.569
    ηHA 0.500 0.400–0.600 0.238
    ηHH 0.100 0.080–0.120 0.144
    σH 1.400 1.120–1.680 -0.067
    σC 2 1.600–2.400 -0.160
    ϕHS 1.750 1.400–2.100 0.098
    r 0.251 0.200–0.301 0.000
    ψC 0.200 0.160–0.240 -0.299
    γHA 1 0.800–1.200 0.065
    γHS 0.583 0.466–0.699 -0.053
    γHH 0.388 0.311–0.466 -0.162
    γC 0.166 0.133–0.200 -0.980
    δHA 3.50×107 2.80×1074.20×107 0.125
    δHS 3.37×106 2.69×1064.04×106 0.024
    δHH 3.37×106 2.69×1064.04×106 -0.104
    δC 3.60×107 2.88×1074.32×107 0.065

     | Show Table
    DownLoad: CSV
    Figure 5.  Partial rank correlation coefficients (PRCCs) of the parameters in the expression of the basic reproduction number R0M of the model (2.1). The baseline values of the parameters used in generating the PRCC values are as given in Table 7.

    (a) the rate at which camels recover from infection (γC, with a PRCC value 0.980),

    (b) the transmission rate from infected camels to susceptible camels (βCC, with a PRCC value +0.921),

    (c) the transmission rate from infected humans to susceptible camels (βHC, with a PRCC value +0.569).

    Thus, based on the parameter values used in the sensitivity analysis and the data used to fit the model (2.1), we identified the top-three parameters that have the highest influence on the dynamics and burden of MERS-CoV (as measured by the changes in the basic reproduction number, R0M) in the Kingdom of Saudi Arabia. Since the recovery rate parameter for camels (γC) has the highest negative correlation with the response function, it follows that control and mitigation strategies (such as isolation and treatment of infectious camels with highly effective antiviral drugs) that increase the recovery rates strategy will reduce the response function (and, consequently, reduce the burden of the disease in the community). Furthermore, the transmission rate parameter from infected camels to susceptible camels (βCC) and from infected humans to susceptible camels (βHC) are highly positively correlated with the response function. Hence, the implementation of control and mitigation strategies that reduce the values of these transmission parameters, such as the quarantine of susceptible camels suspected of being exposed to MERS-CoV and isolation of camels with confirmed symptoms of the disease, will reduce the disease burden in the human and camel populations. Furthermore, these parameters can be reduced by the use of effective face masks by humans who have frequent contacts with camels or the quarantine and isolation of humans suspected of being exposed or those with confirmed symptoms, respectively. In summary, the global sensitivity analysis carried out in this section identified three parameters that have the highest influence on the MERS-CoV transmission dynamics in the KSA. Consequently, control strategies that target these parameters (as described above) will be most effective in reducing the disease burden (and/or eliminating the disease) in the Kingdom.

    As discussed in Section 1, numerous public health control and mitigation strategies are being implemented or proposed in an attempt to effectively curtail the spread of the MERS-CoV pandemic in the human and camel populations. In particular, non-pharmaceutical interventions, such as the use of face masks (to prevent human-to-human or human-to-camel or camel-to-human transmission), quarantine and isolation of suspected humans (camels) and those with clinical symptoms of the disease, are deployed or proposed in an attempt to effectively control the spread and burden of MERS-CoV in the human and camel populations [18,40,42]. Similarly, pharmaceutical interventions, such as the vaccination of humans and/or camels are being proposed as possible future additions to the armory of interventions against the MERS-CoV pandemic [49,50,51,53,57]. In this section, the basic model (2.1) developed in Section 2.1 will be extended to allow for the assessment of the population-level impact of public health intervention (pharmaceutical and non-pharmaceutical), notably face mask usage and the vaccination of both humans and camels. The extended model is formulated below.

    The extended model to be formulated in this section is based on extending the basic model (2.1) to include the effect of face mask usage (by humans) and vaccination (of both humans and camels). In order to incorporate vaccination and the use of face masks, two new state variables will be added to the basic model (2.1). Specifically, we let VH(t) represent the total number of susceptible humans vaccinated against MERS-CoV at time t. Similarly, let VC(t) be the total number of susceptible camels vaccinated against MERS-CoV at time t. Furthermore, in order to incorporate the use of the face mask among humans, we re-scaled the transmission rate parameter from infectious humans to susceptible humans, βHH, by βHH(1εmhcmh), where 0εmh1 represents the efficacy (effectiveness) of the face mask to prevent transmission or acquisition of MERS-CoV infection (from infectious humans to susceptible humans) and 0cmh1 accounts for the compliance (or coverage) of face mask usage by the members of the public. Similarly, in order to incorporate the use of the face mask by humans who have close contact with camels, we re-scaled the transmission rate parameter from infectious camels to susceptible humans, βCH, by βCH(1εmccmc); where 0εmc1 represents the efficacy of the face mask to prevent transmission or acquisition of MERS-CoV infection by infectious camels to susceptible humans, and 0cmc1 measures the masking compliance by individuals who have close contact with camels. Finally, in order to incorporate the use of the face mask by infected humans who have close contact with susceptible camels, we re-scaled the transmission rate parameter from infectious humans to susceptible camels, βHC, by βHC(1εmccmc), where 0εmc1 represents the efficacy of the face mask to prevent transmission of infection and 0cmc1 accounts for the compliance in its usage.

    It follows, based on the above definitions and assumptions, that the equation for the rate of change of the susceptible human population is now given by:

    ˙SH=ΠH+ωVHVH(λH+λCH+ξVH+μH)SH, (4.1)

    where the parameter ΠH is the recruitment rate of humans, ωVH is the waning rate of the vaccine, ξVH is the vaccination rate of humans, μH is the natural death rate of humans, and λH and λCH are the human-to-human transmission rate and the camel-to-human transmission rate now defined, respectively (note that the parameters βHH, ηHA, ηHH, and βCH are as defined above and in Section 2.1):

    λH=(1εmhcmh)[βHH(ηHAIHA+IHS+ηHHIHH)NH], λCH=[βCH(1εmccmc)](ICNC).  (4.2)

    Similarly, the rate of change of the vaccinated human population is given by:

    ˙VH=ξVHSHωVHVH(1εVH)(λH+λCH)VHμHVH, (4.3)

    where ξVH is the vaccination rate, ωVH is the vaccine waning rate, and 0<εVH<1 is the average protective efficacy of the vaccine to prevent vaccinated susceptible humans from acquiring MERS-CoV infection.

    Furthermore, the rate of change of the population of vaccinated camels is given by:

    ˙VC=ξVCSCωVCVC(1εVC)(λC+λHCA+λHCS)VCμCVC, (4.4)

    where ξVC is the vaccination rate of camels, ωVC is the vaccine waning rate for vaccinated camels, 0<εVC<1 is the average protective efficacy of the vaccine to prevent susceptible camels from acquiring MERS-CoV infection, μC is the natural death rate of camels, and the infection rates λC, λHCA, and λHCS (for the camel-to-camel transmission, asymptomatic-human-to-camel transmission, and symptomatic-human-to-camel transmission, respectively) are now given, respectively, by:

    λC=(βCC)(ICNC), λHCA=[βHC(1εmccmc)](ηHAIHANH), λHCS=[βHC(1εmccmc)](IHSNH). (4.5)

    It follows, based on the above derivations and assumptions, that the extended MERS-CoV model with public health interventions is given by the original (basic) model (2.1), but with the equation for the rate of change of the susceptible population now replaced by Eq (4.1), together with Eqs (4.3) and (4.4). In other words, the extended model for the transmission dynamic of MERS-CoV in a human-camel population, that incorporates vaccination of humans and camels and the use of face masks by humans, is given by the following deterministic system of nonlinear differential equations. The flow diagram of the extended model is depicted in Figure 6, and the state variables and parameters of the extended model (4.6) are described in Tables 5 and 6, respectively [4].

    ˙SH=ΠH+ωVHVH(λH+λCH+ξVH+μH)SH,˙VH=ξVHSHωVHVH(1εVH)(λH+λCH)VHμHVH,˙EH=(λH+λCH)SH+(1εVH)(λH+λCH)VH(σH+μH)EH,˙IHA=rσHEH(γHA+μH+δHA)IHA,˙IHS=(1r)σHEH(ϕHS+γHS+μH+δHS)IHS,˙IHH=ϕHSIHS(γHH+μH+δHH)IHH,˙RH=γHAIHA+γHSIHS+γHHIHHμHRH,˙SC=ΠC+ψCRC+ωVCVC(λC+λHCA+λHCS+ξVC+μC)SC,˙VC=ξVCSCωVCVC(1εVC)(λC+λHCA+λHCS)VCμCVC,˙EC=(λC+λHCA+λHCS)SC+(1εVC)(λC+λHCA+λHCS)VC(σC+μC)EC,˙IC=σCEC(γC+μC+δC)IC,˙RC=γCIC(ψC+μC)RC, (4.6)

    where the rates for the human-to-human (λH), camel-to-camel (λC), camel-to-human (λCH), asymptomatically-infectious-human-to-camel (λHCA), and symptomatically-infectious-human-to-camel (λHCS) transmission are given in (4.2) and (4.5).

    Figure 6.  Flow diagram of the extended model (4.6), with the infection rates λH,λC,λCH,λHCA, and λHCS as defined in Eqs (4.2) and (4.5) [4].
    Table 5.  A Description of the state variables of the original model (2.1) and the extended model (4.6).
    State variable Description
    SH Population of susceptible humans
    VH Population of vaccinated humans
    EH Population of exposed (newly infected) humans
    IHA Population of asymptomatically infectious humans
    IHS Population of symptomatically infectious humans
    IHH Population of hospitalized infectious humans
    RH Population of recovered humans
    SC Population of susceptible camels
    VC Population of vaccinated camels
    EC Population of exposed (newly infected) camels
    IC Population of symptomatically infectious camels
    RC Population of recovered camels

     | Show Table
    DownLoad: CSV
    Table 6.  Description of the parameters of the original model (2.1) and the extended model (4.6).
    Parameter Description
    ΠH Recruitment (birth or immigration) rate of humans into the community
    ΠC Recruitment (birth or importation) rate of camels into the community
    μH Natural death rate for humans
    μC Natural death rate for camels
    βHH Transmission rate from infected humans to susceptible humans
    βCC Transmission rate from infected camels to susceptible camels
    βHC Transmission rate from infected humans to susceptible camels
    βCH Transmission rate from infected camels to susceptible humans
    ηHA Modification parameter for the infectiousness of asymptomatically infectious humans, in comparison to the infectiousness of symptomatic humans
    ηHH Modification parameter for the infectiousness of hospitalized humans, in comparison to the infectiousness of symptomatic humans
    σH Rate at which exposed (newly infected) humans become infectious
    σC Rate at which exposed (newly infected) camels become infectious
    ϕHS Hospitalization rate of symptomatic infectious humans
    r Proportion of newly infectious humans who do not show clinical symptoms of the disease
    1r Proportion of newly infectious humans who show clinical symptoms of the disease
    ψC Rate of loss of natural immunity for camels
    γHA Recovery rate for asymptomatically infectious humans
    γHS Recovery rate for symptomatically infectious humans
    γHH Recovery rate for hospitalized infectious humans
    γC Recovery rate for infectious camels
    δHA Disease-induced mortality rate for asymptomatically-infectious humans
    δHS Disease-induced mortality rate for symptomatically infectious humans
    δHH Disease-induced mortality rate for hospitalized humans
    δC Disease-induced mortality rate for camels
    ξVH Vaccination rate for humans
    ωVH Vaccine waning rate for vaccinated humans
    εVH Efficacy of the vaccine to prevent susceptible humans from acquiring infection
    εmh Efficacy of face mask to prevent disease transmission from infectious humans to susceptible humans
    εmc Efficacy of face mask to prevent disease transmission from infectious camels to susceptible humans
    cmh Face mask coverage (compliance) by humans to prevent human-to-human transmission
    cmc Face mask coverage (compliance) for humans who have frequent close contact with camels
    ξVC Vaccination rate for camels
    ωVC Vaccine waning rate for vaccinated camels
    εVC Efficacy of the vaccine to prevent susceptible camels from acquiring infection

     | Show Table
    DownLoad: CSV

    Some of the additional assumptions made in the formulation of the extended model (4.6) are:

    (ⅰ) Vaccinated susceptible humans (in the VH class) are assumed to have received the full required vaccine doses and that enough time has elapsed for the body to develop vaccine-derived immunity. Similarly, vaccinated susceptible camels (in the VC class) are assumed to have received the full required vaccine doses.

    (ⅱ) The vaccines (for both humans and camels) are assumed to be imperfect. Further, the protective efficacy of each of the two vaccines is assumed to wane with time.

    The following result can also be proved using the approach used to prove Theorem (2.1).

    Theorem 4.1. Consider the closed sets

    DHE={(SH,VH,EH,IHA,IHS,IHH,RH)R7+}andDCE={(SC,VC,EC,IC,RC)R5+}.

    The region DE=DHEDCE is positively invariant and attracting with respect to the extended model (4.6).

    The unique DFE of the extended model (4.6) is given by:

    E0v=(SH,VH,EH,IHA,IHS,IHH,RH,SC,VC,EC,IC,RC) =(ΠH(ωVH+μH)μH(ωVH+ξVH+μH),ΠHξVHμH(ωVH+ξVH+μH),0,0,0,0,0,ΠC(ωVC+μC)μC(ωVC+ξVC+μC),ΠCξVCμC(ωVC+ξVC+μC),0,0,0). (4.7)

    Here, too, the linear asymptotic stability of the DFE will be explored using the next generation operator method (used in Section 3). Specifically, the associated non-negative matrix (Fv) of the new infection terms is given by:

    Fv=[ 0f1f2f30f4 000000 000000 000000 0g1g200g3 000000]

    where

    f1=(1εmhcmh)βHHηHA[SH+(1εVH)VHNH],f2=(1εmhcmh)βHH[SH+(1εVH)VHNH],f3=(1εmhcmh)βHHηHH[SH+(1εVH)VHNH],f4=(1εmccmc)βCH[SH+(1εVH)VHNC],g1=(1εmccmc)βHCηHA[SC+(1εVC)VCNH],g2=(1εmccmc)βHC[SC+(1εVC)VCNH]

    and g3=βCC[SC+(1εVC)VCNC]. It should be noted that the proportion of humans vaccinated at the disease-free steady-state is given by VHNH=ξVHωVH+ξVH+μH=fVH. Similarly, the proportion of camels vaccinated at the disease-free steady-state is VCNC=ξVCωVC+ξVC+μC=fVC.

    Similarly, the matrix (V) of the linear transition terms (as defined in Section 3) is given by:

    V=[K100000rσHK20000(1r)σH0K300000ϕHSK4000000K500000σCK6]

    where Ki(i=1,,6) are as defined in Section 3. It follows that the control reproduction number of the extended model (4.6), denoted by RCM, is given by:

    RCM=ρ(FV1) =12[RCC+RHH+(RCCRHH)2+4 RCHRHC ] (4.8)

    where (following the approach in [95]), RCC,RHH,RCH,andRHC are given by:

    RCC=K1K2K3K4σCg36i=1Ki,RHH=K5K6σH[rf1K3K4+(1r)K2(f2K4+f3ϕHS)]6i=1Ki,RCH=σCK1K2K3K4f46i=1Ki,RHC=σHK4K5K6[rg1K3+(1r)g2K2]6i=1Ki. (4.9)

    The quantity RCM is the control reproduction number of the extended model (4.6), and it measures the average number of new MERS-CoV cases (in humans or camels) generated by a typical infectious human or camel if introduced in a human-camel population where a vaccination program (of both humans and camels, using imperfect human and camel vaccines) is implemented, and a certain proportion of humans wear face masks in public (to prevent human-to-human transmission) and when having direct contact with camels (to prevent human-to-camel or camel-to-human transmission). The result below follows from Theorem 2 of [73].

    Theorem 4.2. The disease-free equilibrium (E0v) of the vaccination model (4.6) is locally asymptotically stable if RCM <1, and unstable if RCM >1.

    The epidemiological implication of Theorem 4.2 is that a small influx of infected humans or camels in the well-mixed human-camel population (when RCM <1) will not generate a significant MERS-CoV outbreak in the human-camel community if the initial number of infected humans or camels is small enough (i.e., in the basin of attraction of the aforementioned disease-free equilibrium).

    The extended model (4.6) will now be simulated to assess the population-level impact of pharmaceutical and non-pharmaceutical intervention strategies, such as the use of face masks (by humans in public or when in close proximity with camels) and vaccination (of both humans and camels), on the transmission dynamics and control of MERS-CoV in the Kingdom of Saudi Arabia. Unless otherwise noted, the simulations will be carried out using the baseline values of the fixed and fitted parameters given in Table 7. The following three intervention scenarios will be simulated:

    Table 7.  Baseline values of the fixed and estimated parameters of the original model (2.1) and the extended model (4.6).
    Parameter Baseline value (week1) Reference
    ΠH 191 [86,87]
    ΠC 14 [47,88]
    μH 5.21×106 [86,87]
    μC 9.78×106 [28,88]
    βHH 0.114 Fitted
    βCC 0.135 Fitted
    βCH 0.023 Fitted
    βHC 0.568 Fitted
    ηHA 0.5 (dimensionless) Assumed
    ηHH 0.1 (dimensionless) Assumed
    σH 1.400 [59,90,91]
    σC 2 [63]
    ϕHS 1.750 [15,90]
    r 0.251 (dimensionless) [45]
    1r 0.749 (dimensionless) [45]
    ψC 0.2 [26]
    γHA 1 [15]
    γHS 0.583 [90]
    γHH 0.388 [28]
    γC 0.166 [92]
    δHA 3.5×107 Assumed
    δHS 3.37×106 Adapted from [7,87]
    δHH 3.37×106 Adapted from [7,87]
    δC 3.6×107 Assumed
    ξVH 2.8×106 [93]
    ωVH 0.027 [94]
    ξVC 1×103 Assumed
    ωVC 0.027 [94]

     | Show Table
    DownLoad: CSV

    (a) Scenario 1 (face mask as a sole intervention): Thr use of face masks worn by individuals having close contact with camels is implemented as the sole intervention in the community.

    (b) Scenario 2 (vaccinating humans or camels as a sole intervention):

    (2a) The vaccination of humans as the sole intervention is implemented in the community.

    (2b) The vaccination of camels as the sole intervention is implemented in the community.

    (c) Scenario 3 (combined mask-vaccination (hybrid) strategy): The vaccination (of both humans and camels) and the use of face masks by humans in public or when in close proximity with camels is administered as the main intervention in the community.

    Most of the simulations to be carried out are for the period between May 1, 2023 (considered as Day 1 of our simulations) and April 30, 2024. Further, for simulation purposes, the initial conditions of the model will be chosen based on the initial cumulative number of MERS-CoV cases (2196) and mortality (855) in humans, as well as the cumulative number of cases in camels (977), in the Kingdom of Saudi Arabia recorded on May 1, 2023 (considered as Day 1 of the simulations) [7,96]). The aforementioned three simulation scenarios are briefly discussed below.

    To assess the impact of face mask use as a sole intervention (where the masks are worn by humans when they come in close proximity with camels only) on the transmission dynamics and control of MERS-CoV in the KSA, the extended model (4.6) is simulated for the special case without vaccination (i.e., ξVH=εVH=ωVH=ξVC=εVC=ωVC=0) and humans do not wear face masks in public (i.e., εmh=cmh=0). The simulation results obtained are illustrated in the form of contour plots of the associated control reproduction number, denoted by ˜RCMmc (˜RCMmc = RCM|ξVH=εVH=ωVH=ξVC=εVC=ωVC=εmh=cmh=0), as a function of face mask efficacy (εmc) and coverage (cmc), depicted in Figure 7. This figure shows that the associated control reproduction number ˜RCMmc decreases with increasing values of the mask efficacy and compliance, and that ˜RCMmc lies in the range ˜RCMmc[0.85,1.1] (suggesting that the use of face masks as a sole intervention can lead to the elimination of the disease in the population, depending on the efficacy of, and the compliance in using, face masks in the community. For example, this figure shows that if surgical masks (with an estimated efficacy of about 70% [93,97]) is prioritized, the disease can be effectively controlled or eliminated (in both the human and camel populations) if at least 45% of the human population consistently wear such masks whenever they are within close proximity to camels. This result is illustrated, by simulating the extended model with mask efficacy and compliance set at 70% and 45%, respectively (and other parameter values chosen such that the associated reproduction number of the model, denoted by ˜RCM<, is less than one), showing convergence of multiple initial conditions to the disease-free equilibrium (see Figure (8); which suggests disease elimination in about 300 days).

    Figure 7.  Simulations of the extended model (4.6) to assess the effect of the use of face masks by humans within close proximity to camels as the sole public health intervention (Scenario 1), depicting contour plots of the associated control reproduction number (˜RCM) of the extended model (4.6), as a function of mask efficacy (εmc) and compliance (cmc). The other parameter values used to generate the contour plots are as given in Table 7, with ξVH=εVH=ωVH=ξVC=εVC=ωVC=εmh=cmh=0and, now, βCH=0.223.
    Figure 8.  Simulations of the extended model (4.6), showing a time series of the total number of new cases as a function of time using the parameter values given in Table 7, with εmh=εmc=0.7andcmh=cmc=0.45 (so that, ˜RCM = 0.836 < 1), showing convergence of all of the initial conditions to the disease-free equilibrium.

    In this section, the extended model (4.6) will be simulated using the baseline values of the parameters in Table 7 (unless otherwise stated) to assess the population-level impact of the vaccination program (in both the human and camel population) implemented as the sole public health intervention in the human-camel community within the Kingdom of Saudi Arabia. In other words, in this section, the extended model (4.6) will be simulated for the scenario where masking is not implemented in the Kingdom (i.e., we set εmh=εmc=cmh=cmc=0), and vaccination is the only intervention implemented in the community. We consider the following two cases of vaccinating humans only and vaccinating camels only.

    Here, the extended model (4.6) is simulated to assess the impact of vaccinating humans only (i.e., camels are not vaccinated and humans do not wear face masks). For these simulations, camels are not vaccinated (so that, ξVC=εVC=ωVC=0) and humans do not wear masks in public or when in close contact with camels (i.e., εmh=εmc=cmh=cmc=0). The results obtained are depicted by the contour plots, of the human vaccination reproduction number ˜RCMH = RCM|ξVC=εVC=ωVC=εmh=εmc=cmh=cmc=0 of the model (4.6), as a function of vaccine efficacy (εVH) and vaccine coverage in human population (fVH) at the steady-state, in Figure 9. This figure shows, first of all, that, during the simulation period, the value of the human vaccination reproduction number of the extended model, ˜RCMH, decreases with increasing values of the vaccine efficacy and coverage in the human population. Furthermore, ˜RCMH lies in the range ˜RCMH[0.816,0.827], suggesting that (since ˜RCMH<1) the disease can be effectively controlled (in line with Theorem 4.2) in the Kingdom. Thus, the results depicted in this figure show that even if vaccination (of humans) is not implemented in the Kingdom, the human vaccination reproduction number is below one (with the maximum value of 0.827) during the simulation period. Furthermore, this figure shows that the effect of vaccinating humans only seems to generally have a marginal impact in reducing the human vaccination reproduction number (hence, vaccinating humans seems to have only a marginal impact in reducing the disease burden in the population). For instance, using a perfect vaccine with 100% coverage could only reduce the control reproduction from the maximum value of 0.827 to 0.816 (i.e., reduction by about 10%).

    Figure 9.  Simulations of the extended model (4.6) to assess the effect of vaccinating humans as a sole intervention (Scenario 2a), depicting contour plots of the human vaccination reproduction number (˜RCMH) of the extended model (4.6), as a function of the vaccination coverage for humans at steady-state (fVH=VH/NH) and efficacy (εVH). The other parameter values used to generate these contour plots are as given in Table 7 with ξVC=εVC=ωVC=εmh=εmc=cmh=cmc=0.

    The results for this scenario (with vaccination of humans as the only intervention implemented) is also presented in terms of the cumulative number of new cases for the following three levels of the human vaccination-only strategy:

    (a) The low effectiveness level of the human vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 25% (i.e., εVH = 0.25) and the vaccination rate to be ξVH = 5×103.

    (b) The moderate effectiveness level of the human vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 50% (i.e., εVH = 0.50) and the vaccination rate to be ξVH = 15×103.

    (c) The high effectiveness level of human the vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 75% (i.e., εVH = 0.75) and the weekly vaccination rate to be ξVH = 5×102.

    The results obtained, depicted in Figure 10, show a marginal decrease in the cumulative number of new cases in the human population, in comparison to the baseline scenario (where all parameters are maintained at their baseline values given in Table 7). For instance, the cumulative number of new cases at the end of the one-year simulation period under the baseline scenario is 2410. This number decreases to 2403 and 2378 under the low and moderate levels of the human vaccination-only strategy, respectively (this represents less than 1% and 1.3% reduction in the cumulative new cases, in comparison to the baseline, respectively). Similarly, the cumulative number of new cases obtained using the high effectiveness level of this intervention is 2316 (which represents about 4% reduction, in comparison to the cumulative new cases generated under the baseline scenario). Thus, these simulations also confirm that, based on the parameter values used in the simulations, vaccinating humans as a sole intervention strategy would have only a marginal impact in reducing the disease burden in the Kingdom of Saudi Arabia.

    Figure 10.  Simulation of the extended model (4.6) assessing the impact of vaccinating humans as a sole intervention (i.e., no vaccination of camels and no masking). Time series of the cumulative number of new MERS-CoV cases in humans, as a function of time, for various levels of effectiveness and efficacy of vaccination. Simulations are carried out for a one-year duration. (ⅰ) The red curve represents the results obtained for the case where the parameters of the model are maintained at their baseline values given in Table 7. (ⅱ) The black curve represents the results obtained for the low effectiveness level of this intervention (εVH = 0.25 and ξVH = 5×103). (ⅲ) The magenta curve represents the results obtained for the moderate effectiveness level of this intervention (εVH = 0.50 and ξVH = 15×103), and (iv) The blue curve represents the results obtained for the high effectiveness level of this intervention (εVH = 0.75 and ξVH = 5×102). Parameter values used in these simulations are as given in Table 7, with ξVC=εVC=ωVC=εmh=εmc=cmh=cmc=0 and various values of ξVH and εVH.

    Here, the extended model (4.6) is simulated to assess the impact of vaccinating camels only (i.e., humans are not vaccinated and humans do not wear face masks). The simulations are also carried out for the one-year period from May 1, 2023 to April 30, 2024. Since humans are not vaccinated, we set (for these simulations) ξVH=εVH=ωVH=0. Furthermore, since humans do not wear face masks under this scenario where vaccination of camels is the only intervention strategy implemented, we set the mask-related parameters to zero (i.e., we set εmh=εmc=cmh=cmc=0). Figure 11 depicts the results obtained for this camel vaccination-only intervention scenario, showing contour plots of the camel vaccination reproduction number (˜RCMC = RCM|ξVH=εVH=ωVH=εmh=εmc=cmh=cmc=0), as a function of the vaccine efficacy (εVC) and vaccine coverage in the camel population (fVC) at the steady-state. This figure also shows a decrease in the camel vaccination reproduction number with increasing efficacy and coverage of the camel vaccine. Furthermore, the camel vaccination reproduction number lies in the range ˜RCMC[0.1,0.8], suggesting that the disease can be effectively controlled even if the vaccination of camels (as a sole intervention) is not implemented. Unlike for the case where vaccinating humans was the sole intervention strategy implemented (which only had a marginal impact in reducing the associated reproduction number), the simulation results depicted in Figure 11 show that vaccinating camels as a sole intervention can significantly reduce the camel vaccination reproduction number (hence, reduce the disease burden in the camel population). For instance, using a camel vaccine with perfect efficacy (εVC=1) and perfect coverage (fVC=1) reduces the maximum value of the control reproduction number (0.8) to 0.1 (representing about 88% reduction). These results are also presented in terms of the cumulative number of new cases of MERS-CoV infection in camels, as a function of time, using the following effectiveness levels of the camel vaccination-only intervention:

    Figure 11.  Simulations of the extended model (4.6) to assess the effect of vaccinating camels as a sole intervention (Scenario 2b), depicting contour plots of the camel vaccination reproduction number (˜RCMC) of the extended model (4.6), as a function of vaccination coverage for camels at steady-state (fVC=VC/NC) and efficacy (εVC). The other parameter values used to generate the contour plots are as given in Table 7 with ξVH=εVH=ωVH=εmh=εmc=cmh=cmc=0.

    (a) The low effectiveness level of the camel vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 25% (i.e., εVC = 0.25) and the vaccination rate to be ξVC = 2×103.

    (b) The moderate effectiveness level of the camel vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 50% (i.e., εVC = 0.50) and the vaccination rate to be ξVC = 6×103.

    (c) The high effectiveness level of the camel vaccination-only strategy: for this strategy, we set the vaccination efficacy to be 75% (i.e., εVC = 0.75) and the weekly vaccination rate to be ξVC = 8×103.

    The results obtained are depicted in Figure 12, showing a significant reduction in the baseline cumulative number of new cases with an increasing effectiveness level of the camel vaccination-only intervention. For instance, the cumulative number of new cases under the baseline scenario at the end of the simulation period is 7154. It reduces to 6719 and 5413, respectively, under the low and moderate effectiveness levels of this intervention (representing, respectively, a 6% and 24% reduction, in comparison to the baseline value). Furthermore, the cumulative number of new cases recorded under the high effectiveness level of this intervention is 4730 (which represents a 34% reduction in comparison to the baseline).

    Figure 12.  Simulation of the extended model (4.6) assessing the impact of vaccinating camels as a sole intervention (i.e., no vaccination of humans and no masking). Time series of the cumulative number of new MERS-CoV cases in camels, as a function of time, for various levels of effectiveness and efficacy of vaccination. Simulations are carried out for a one-year duration. (ⅰ) The red curve represents the results obtained for the case where the parameters of the model are maintained at their baseline values given in Table 7. (ⅱ) The black curve represents the results obtained for the low effectiveness level of this intervention (εVC = 0.25 and ξVC = 2×103). (ⅲ) The magenta curve represents the results obtained for the moderate effectiveness level of this intervention (εVC = 0.50 and ξVC = 6×103), and (ⅳ) The blue curve represents the results obtained for the high effectiveness level of this intervention (εVC = 0.75 and ξVC = 8×102). Parameter values used in these simulations are as given in Table 7, with ξVH=εVH=ωVH=εmh=εmc=cmh=cmc=0 and various values of ξVC and εVC.

    The extended model (4.6) will now be simulated, using the baseline values of the parameters in Table 7 (unless otherwise stated), to assess the combined population-level impacts of human and camel vaccination and mask usage interventions in the Kingdom of Saudi Arabia. This intervention is referred to as a hybrid strategy. The simulations for the hybrid strategy will be carried out for a one-year duration, from May 1, 2023 to April 30, 2024. Furthermore, we consider the following effectiveness levels of the hybrid vaccination and mask usage strategy:

    (a) The low effectiveness level of the hybrid strategy: here, we set the protective efficacies of the vaccines (for humans and camels) and masks to be 30% (i.e., εVH = εmh = εmc = εVC = 0.3), mask compliance to be 25% (i.e., cmh=cmc = 0.25), and the vaccination rates for humans and camels set at ξVH=5×103 and ξVC=2×103, respectively.

    (b) The moderate effectiveness level of the hybrid strategy: for this effectiveness level, we set the efficacies of the vaccines and mask to be 60% (i.e., εVH = εmh = εmc = εVC = 0.6), mask compliance to be 50% (i.e., cmh=cmc = 0.50), and the vaccination rates for humans and camels set at ξVH=15×103 and ξVC=6×103, respectively.

    (c) The high effectiveness level of the hybrid strategy: for the high effectiveness level of the hybrid strategy, we set the protective efficacies of the vaccines and mask to be 80% (i.e., εVH = εmh = εmc = εVC = 0.8), mask compliance to be 75% (i.e., cmh=cmc = 0.75), and the vaccination rates for humans and camels set, respectively, at ξVH=5×102 and ξVC=8×103.

    The results obtained for a cumulative number of new cases in humans is depicted in Figure 13. This figure shows a marginal decrease in the cumulative number of new cases with increasing effectiveness levels of the hybrid strategy, in comparison to the baseline scenario. The results obtained for the cumulative number of new MERS-CoV cases in camels are depicted in Figure 14. Unlike in the case of the cumulative new cases in humans (depicted in Figure 13), this figure shows a significant decrease in the cumulative number of new MERS-CoV cases in camels with increasing effectiveness levels of the hybrid strategy, in comparison to the baseline scenario. For instance, the cumulative number of new cases at the end of the simulation period under the baseline scenario is 7154. This figure reduces to 6447 and 4902, respectively, under the low effectiveness and the moderate effectiveness level of the hybrid strategy (representing about 9% and 32% reduction, in comparison to the baseline). Similarly, the cumulative number of new cases generated using the high effectiveness level of the hybrid strategy is 4158 (representing about 42% reduction, in comparison to the baseline value). Thus, this study shows that while the hybrid strategy only resulted in a marginal decrease in the MERS-CoV burden in humans, it resulted in a very significant reduction in the cumulative number of new cases in camels (with increasing effectiveness levels of the hybrid strategy), in comparison to the baseline scenario. In other words, the hybrid strategy will be more effective in curtailing the disease burden in camels than in the human population (albeit there is also an indirect benefit in reducing the burden in humans, since reducing cases in camels also implies some reduction of cases in the human population). This study suggests that strategies that focus on reducing the disease burden in the camel population (notably vaccinating camels with an effective vaccine) would significantly contribute to the effort to effectively control (or eliminate) the MERS-CoV pandemic in the human population. The extended model was simulated to illustrate the high effectiveness level of the hybrid strategy (using parameter values such that the associated reproduction number, ˜RCM, is less than one) and the results obtained, depicted in Figure 15, show convergence of all of the initial conditions to the disease-free equilibrium and disease elimination is reached within about 150 days; which is faster than the 300 days to achieve elimination predicted by simulating the model for the human vaccination-only strategy, shown in Figure 8.

    Figure 13.  Simulation of the extended model (4.6) to assess the impact of the hybrid strategy (which entails the vaccination of humans and camels, and the use of face masks by humans in public and when having close contact with camels). Time series of the cumulative number of new MERS-CoV cases in humans, as a function of time, for the period from May 1, 2023 to April 30, 2024, for various effectiveness levels of the hybrid strategy. (ⅰ) The baseline strategy (red curve) where parameters are maintained at their baseline values; (ⅱ) the low effectiveness level of the hybrid strategy (black curve); (ⅲ) the moderate effectiveness level of the hybrid strategy (magenta curve); and (ⅳ) the high effectiveness level of the hybrid strategy (blue curve). Parameter values used are as given in Table 7 with various values of ξVH,εVH,εmh,εmc,cmh,cmc,ξVC,andεVC.
    Figure 14.  Simulation of the extended model (4.6) to assess the impact of the hybrid strategy (which entails the vaccination of humans and camels, and the use of face masks by humans in public and when having close contact with camels). Time series of the cumulative number of new MERS-CoV cases in camels, as a function of time, for the period from May 1, 2023 to April 30, 2024, for various effectiveness levels of the hybrid strategy. (ⅰ) The baseline strategy (red curve) where parameters are maintained at their baseline values; (ⅱ) the low effectiveness level of the hybrid strategy (black curve); (ⅲ) the moderate effectiveness level of the hybrid strategy (magenta curve); and (ⅳ) the high effectiveness level of the hybrid strategy (blue curve). Parameter values used are as given in Table 7 with various values of ξVH,εVH,εmh,εmc,cmh,cmc,ξVC,andεVC.
    Figure 15.  Simulations of the extended model (4.6), showing a time series of the total number of new cases as a function of time using the parameter values given in Table 7, with εVH = εmh = εmc = εVC = 0.8, cmh=cmc = 0.75, ξVH=5×102, and ξVC=8×103.

    Furthermore, it is worth stating that the cumulative number of new cases (at the end of the simulation period) for humans generated using the hybrid strategy is much lower than the corresponding cumulative number of new cases generated using any of the vaccination-only strategies discussed above. For instance, while the cumulative new cases generated using the high effectiveness level of the humans-only vaccination was about 2315 (see the blue curve in Figure 10), the cumulative number of new human cases generated using the corresponding high effectiveness level of the hybrid strategy (Figure 13, blue curve) is 2227 (which represents about 4% reduction, in comparison to the cumulative number of new cases generated using the corresponding high effectiveness level of the humans-only vaccination strategy). Similarly, the cumulative number of new cases in camels generated using the high effectiveness level of the camels-only vaccination strategy was 4700 (Figure 12, blue curve), whereas that for the corresponding high effectiveness level of the hybrid strategy is 4000 (Figure 14, blue curve), which represents about 15% reduction, in comparison to the cumulative number of new cases recorded under the camels-only vaccination strategy). The numerical simulations depicted in this section (Figures 10, 1214) consistently highlight the fact that the intervention strategies implemented caused greater reduction in new cases in camels than in humans.

    The Middle Eastern respiratory syndrome (MERS-CoV), a coronavirus transmitted to humans by dromedary camels, emerged out of the Kingdom of Saudi Arabia (KSA) in 2012 causing unprecedented burden in about 30 countries (with the KSA suffering the brunt of the burden). The KSA is most vulnerable to MERS-CoV-2 outbreaks due to its hosting of several large gathering events of humans and camels, notably the annual Hajj pilgrimage (once a year) and Umrah year [8], and the camel festival in the capital city of Riyadh (which typically runs for a two-month period) [9,10]. This study presented a new mathematical model for the transmission dynamics of MERS-CoV in the KSA. The basic model, which takes the form of a deterministic system of nonlinear differential equations, was, first of all, fitted using observed MERS-CoV data for the KSA (for the period from Week 13 to Week 42 of 2013), and the fitted model was used to estimate the unknown parameters of the model. The data for Week 43 of 2013 to Week 9 of 2014 was used to cross-validate the model (4.6). This cross-validation, together with simulations involving the cumulative cases of MERS-CoV, presented a good match to the observed data in the Kingdom of Saudi Arabia.

    The basic model was rigorously analyzed to gain insight into its qualitative properties. Specifically, results with respect to the boundedness and invariance of its solutions, as well as the asymptotic stability properties of its associated disease-free equilibrium, were provided. Specifically, it was shown that the disease-free equilibrium of the model was locally asymptotically stable whenever its reproduction number (denoted by R0M) was less than unity. The epidemiological implication of this result is that the effective control of the disease in the human-camel population within the KSA (when R0M) depends on the initial sizes of the sub-populations of the model (in other words, disease control or elimination depends on whether or not the initial number of infected humans and camels lie within the basin of attraction of the disease-free equilibrium). Using the fixed and estimated parameters of the basic model, this study showed that the value of R0M for the Kingdom of Saudi Arabia was approximately 0.84, suggesting that the prospects for MERS-CoV elimination in the KSA are highly promising.

    Global sensitivity analysis was carried out, using Latin hypercube sampling and partial rank correlation coefficients (PRCCs), to determine the parameters of the model (2.1) that have the greatest influence on the value of the chosen response function (R0M, which is a measure of the burden of the disease in the Kingdom). The top five PRCC-ranked parameters identified in this study were: (a) the parameters for the rates of camels recovered from the MERS-CoV infection (γC), (b) the transmission rate from infected to susceptible camels (βCC), and (c) the transmission rate from infected humans to susceptible camels (βHC). The numerical size and (negative) sign of the computed PRCC values of the parameters of the model (4.6) (tabulated in Table 4) indicate that reduction of the infection rate of infectious camels to susceptible camels (βCC) and the infection rate from infectious humans to susceptible camels (βHC), which can be achieved, for instance, by the quarantine and isolation of suspected and symptomatic humans and camels, will lead to a corresponding reduction of R0M (hence, a reduction in disease burden). The use of effective face masks by humans in the public or when having contact with camels, social-distancing between humans, and community lockdowns in both populations will also decrease the transmission rate parameters (hence, decrease the response function, R0M; and, consequently, the disease burden in both host populations). Similarly, the PRCC values obtained (tabulated in Table 4) show that the implementation of control and mitigation measures that increase the value of the recovery rate parameter for camels (γC), such as by the isolation and treatment of infectious camels with highly effective antiviral drugs, will correspondingly reduce the response function (and, consequently, reduce the burden of the disease in the Kingdom).

    The basic model was extended to incorporate non-pharmaceutical and pharmaceutical public health interventions, notably the use of face masks by humans in public and when having close contact with camels, and vaccination of both humans and camels. Here, too, the disease-free equilibrium of the resulting extended model (4.6) was shown to be locally asymptotically stable whenever its associated control reproduction number (denoted by RCM) was less than one. The implication of this result is that the implementation of the aforementioned public health interventions could lead to the effective control (or elimination) of MERS-CoV in the Kingdom of Saudi Arabia if the interventions could bring (and maintain) the value of the control reproduction number (RCM) to be below one (provided the initial number of infected individuals and camels lie within the basin of attraction of the DFE).

    The extended model was then simulated to assess the population-level impact of the chosen public health interventions. Specifically, simulations were, first, carried out (for a period of one year) to assess the impact of face mask use as a sole intervention (where the masks are worn by humans when they come in close proximity with camels only) on the transmission dynamics and control of MERS-CoV in the KSA. The extended model was simulated for the special case without vaccination and humans do not wear face masks in public. The simulation results showed that the associated control reproduction number ˜RCMmc decreased with increasing values of the mask efficacy and compliance, and that ˜RCMmc lied in the range ˜RCMmc[0.85,1.1] (suggesting that the use of face masks as a sole intervention can lead to the elimination of the disease in the population, depending on the efficacy of, and compliance in, mask usage in the community. For example, this figure shows that if surgical masks (with an estimated efficacy of about 70%) is prioritized, the disease can be effectively controlled or eliminated (in both the human and camel populations) if at least 45% of the human population consistently wear such masks whenever they are within close proximity to camels. Then, we studied the effects of vaccination of humans as a sole intervention strategy implemented in the Kingdom (i.e., humans do not wear face masks in public or when having close contact with camels, and camels are not vaccinated) under three effectiveness levels of this human vaccination-only strategy (namely, low, moderate, and high). The results obtained, for the case where only humans are vaccinated, showed that although the the burden of MERS-CoV in the human population decreases with an increasing effectiveness level of this (human vaccination-only) intervention (as measured in terms of the reduction in the associated control reproduction number of the extended model or the cumulative number of new cases and cumulative mortality in humans), in comparison to the baseline scenario, this reduction was quite marginal. In other words, these simulations showed that, based on the parameter values in the simulations, vaccinating humans (as a sole intervention) will only offer a marginal impact in reducing the burden of the MERS-CoV pandemic in the Kingdom of Saudi Arabia. On the other hand, the simulations for the case where vaccination of camels was the sole intervention implemented (and humans are not vaccinated and do not wear masks) showed, over the same one-year period, a significant reduction in both the associated reproduction number and the cumulative number of new MERS-COV cases in camels with an increasing effectiveness level of this (camels vaccination only) intervention strategy, in comparison to the baseline scenario. For instance, these simulations show that while the moderate effectiveness level of this strategy could reduce the baseline cumulative new cases at the end of the one-year simulation period by about 24%, the high effectiveness level of this intervention reduces the cumulative number of new cases by about 34%.

    Finally, the extended model was used to simulate the population-level impact of a hybrid strategy, where both humans and camels are vaccinated and humans wear mask in public and when having close contact with camels. Here, too, three effectiveness levels (low, moderate, and high) of the hybrid strategy were considered. The simulation results obtained showed that while this strategy only offers a marginal decrease in the disease burden in humans, with an increasing effectiveness level, in comparison to the baseline, it resulted in a marked decrease in the cumulative number of new MERS-CoV cases in camels (for instance, the moderate and high effectiveness levels of this intervention reduced the baseline cumulative new cases in camels by 32% and 42%, respectively). Thus, this study shows that, while the hybrid strategy only offers a marginal reduction in the burden of the MERS-CoV pandemic in humans, it can result in a significant reduction in the cumulative new cases in camels. The cumulative number of new cases for humans (camels) generated using the hybrid strategy was consistently lower than the cumulative number of new cases obtained using a vaccination-only strategy for humans (camels). For instance, the simulations in this study showed that the high effectiveness level of the hybrid strategy reduced the cumulative number of new cases in humans (camels) by about 4% (15%), in comparison to the corresponding cumulative new cases generated using a vaccination-only strategy (and more of such reductions were recorded in camels than in humans). Thus, this study emphasizes control measures that focus on the camel population (e.g., routine vaccination of camels using an effective vaccine with moderate to high coverage), instead of those that focus on humans. In summary, the analytical and simulation results generated in this study showed that the prospects for the effective control of MERS-CoV in the Kingdom of Saudi Arabia are promising using the public health interventions described in the study, provided their effectiveness levels are high enough. In particular, the implementation of a hybrid strategy that entails the vaccination of both humans and camels and the use of face masks by humans in public and when having close contact with camels can significantly reduce the disease burden in the Kingdom of Saudi Arabia. Greater reductions (in the burden of the disease) were recorded in the camel population than in the human population (the impacts of the interventions in reducing the burden of the pandemic in the latter, from the baseline scenario, were consistently quite marginal). Thus, this study showed that efforts should be focused on prioritizing interventions that effectively control the disease in the camel population, rather than in the human population. In other words, this study suggests that control resources should be prioritized to disease control in camels rather than humans.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    ABG acknowledges the support, in part, of the National Science Foundation (DMS-2052363). AA acknowledges the scholarship received from the University of Tabuk, Saudi Arabia. The authors are grateful to the three anonymous reviewers for their very constructive comments.

    Abba B. Gumel is an editorial board member for Mathematical Biosciences and Engineering and was not involved in the editorial review or the decision to publish this article. All authors declare that there are no competing interests.

    (1) Initial conditions for Figure 8:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (9,900,000, 0, 10,100,100, 10, 1, 3,900,000, 0, 10,977, 1).

    (2) Initial conditions for Figure 10:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (36,617,538, 0, 1, 1, 1, 1, 1, 1,400,000, 0, 1, 7, 1).

    (3) Initial conditions for Figure 12:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (36,617,538, 0, 2196, 1, 1, 1, 1, 1,400,000, 0, 60,977, 0).

    (4) Initial conditions for Figure 13:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (36,617,538, 0, 1, 1, 1, 1, 1, 1,400,000, 0, 1, 7, 1).

    (5) Initial conditions for Figure 14:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (36,617,538, 0, 2196, 1, 1, 1, 1, 1,400,000, 0, 60,977, 0).

    (6) Initial conditions for Figure 15:

    (SH(0), VH(0), EH(0), IHA(0), IHS(0), IHH(0), RH(0), SC(0), VC(0), IC(0), RC(0), ) = (9,000,000, 10, 10,100,100,100, 1, 50,000,100, 10, 9700, 1).



    [1] R. Abraham, J. Marsden, Foundations of Mechanics, 2nd edition, Addison-Wesley, 1978, (with the assistance of Tudor Ratiu and Richard Cushman).
    [2] P. A. Absil, R. Mahony, R. Sepulchre, Optimization Algorithms on Matrix Manifolds, Princeton University Press, Princeton, NJ, 2008. https://doi.org/10.1515/9781400830244
    [3] V. Arnold, On the differential geometry of Lie groups of infinite dimension and its applications to the hydrodynamics of perfect fluids, Ann. Fourier Inst., 16 (1966), 319–361.
    [4] V. Arnold, Mathematical Methods of Classical Mechanics, vol. 60 of Graduate Texts in Mathematics, Springer-Verlag, 1989, Translated from the 1974 Russian original by K. Vogtmann and A. Weinstein.
    [5] M. Barbero Liñán, D. Martín de Diego, Bäcklund transformations in discrete variational principles for Lie–Poisson equations, in Discrete Mechanics, Geometric Integration and Lie–Butcher Series (eds. K. Ebrahimi-Fard and M. Barbero Liñán), Springer International Publishing, Cham, 2018,315–332.
    [6] A. Bloch, L. Colombo, F. Jiménez, The variational discretization of the constrained higher-order Lagrange–Poincaré equations, Discrete Contin. Dyn. Syst., 39 (2019), 309–344. https://doi.org/10.3934/dcds.2019013 doi: 10.3934/dcds.2019013
    [7] A. Bobenko, Y. Suris, Discrete Lagrangian reduction, discrete Euler–Poincaré equations, and semidirect products, Lett. Math. Phys., 49 (1999), 79–93. https://doi.org/10.1023/A:1007654605901 doi: 10.1023/A:1007654605901
    [8] M. I. Caruso, J. Fernández, C. Tori, M. Zuccalli, Discrete mechanical systems in a Dirac setting: a proposal, 2022, URL https://arXiv.org/abs/2203.05600.
    [9] J. Fernández, M. Zuccalli, A geometric approach to discrete connections on principal bundles, J. Geom. Mech., 5 (2013), 433–444. https://doi.org/10.3934/jgm.2013.5.433 doi: 10.3934/jgm.2013.5.433
    [10] E. Hairer, C. Lubich, G. Wanner, Geometric Numerical Integration: Structure-preserving algorithms for ordinary differential equations, vol. 31 of Springer Series in Computational Mathematics, 2nd edition, Springer-Verlag, 2006.
    [11] S. Jalnapurkar, M. Leok, J. Marsden, M. West, Discrete Routh reduction, J. Phys. A: Math. Gen., 39 (2006), 5521. https://doi.org/10.1088/0305-4470/39/19/S12 doi: 10.1088/0305-4470/39/19/S12
    [12] S. Lall, M. West, Discrete variational Hamiltonian mechanics, J. Phys. A: Math. Gen., 39 (2006), 5509. https://doi.org/10.1088/0305-4470/39/19/S11 doi: 10.1088/0305-4470/39/19/S11
    [13] T. Lee, M. Leok, N. Mcclamroch, Lagrangian mechanics and variational integrators on two-spheres, Int. J. Numer. Methods Eng., 79 (2009), 1147–1174. https://doi.org/10.1002/nme.2603 doi: 10.1002/nme.2603
    [14] M. Leok, J. Marsden, A. Weinstein, A discrete theory of connections on principal bundles, 2005, URL https://arXiv.org/abs/math/0508338.
    [15] M. Leok, T. Ohsawa, Discrete Dirac structures and implicit discrete Lagrangian and Hamiltonian systems, AIP Conference Proceedings, 1260 (2010), 91–102. https://doi.org/10.1063/1.3479325 doi: 10.1063/1.3479325
    [16] M. Leok, T. Ohsawa, Variational and geometric structures of discrete Dirac mechanics, Found. Comput. Math., 11 (2011), 529–562. https://doi.org/10.1007/s10208-011-9096-2 doi: 10.1007/s10208-011-9096-2
    [17] M. Leok, J. Zhang, Discrete Hamiltonian variational integrators, IMA J. Numer. Anal., 31 (2011), 1497–1532. https://doi.org/10.1093/imanum/drq027 doi: 10.1093/imanum/drq027
    [18] Z. Ma, C. Rowley, Lie–Poisson integrators: A Hamiltonian, variational approach, Int. J. Numer. Methods Eng., 82 (2010), 1609–1644. https://doi.org/10.1002/nme.2812 doi: 10.1002/nme.2812
    [19] J. Marrero, D. Martín de Diego, E. Martínez, Discrete Lagrangian and Hamiltonian mechanics on Lie groupoids, Nonlinearity, 19 (2006), 1313. https://doi.org/10.1088/0951-7715/19/6/006 doi: 10.1088/0951-7715/19/6/006
    [20] J. Marsden, Lectures on Mechanics, Lecture note series. London Mathematical Society, Cambridge University Press, 1992.
    [21] J. Marsden, S. Pekarsky, S. Shkoller, Discrete Euler–Poincaré and Lie–Poisson equations, Nonlinearity, 12 (1999), 1647–1662. https://doi.org/10.1088/0951-7715/12/6/314 doi: 10.1088/0951-7715/12/6/314
    [22] J. Marsden, S. Pekarsky, S. Shkoller, Symmetry reduction of discrete Lagrangian mechanics on Lie groups, J. Geom. Phys., 36 (2000), 140–151. https://doi.org/10.1016/S0393-0440(00)00018-8 doi: 10.1016/S0393-0440(00)00018-8
    [23] J. Marsden, T. Ratiu, Introduction to Mechanics and Symmetry, vol. 17 of Texts in Applied Mathematics, 2nd edition, Springer-Verlag, 1999.
    [24] J. Marsden, J. Scheurle, Lagrangian reduction and the double spherical pendulum, Z. angew. Math. Phys., 44 (1993), 17–43. https://doi.org/10.1007/BF00914351 doi: 10.1007/BF00914351
    [25] J. E. Marsden, J. Scheurle, The reduced Euler-Lagrange equations, Fields Institute Communications, 1.
    [26] J. Marsden, A. Weinstein, Reduction of symplectic manifolds with symmetry, Rep. Math. Phys., 5 (1974), 121–130. https://doi.org/10.1016/0034-4877(74)90021-4 doi: 10.1016/0034-4877(74)90021-4
    [27] J. Marsden, M. West, Discrete mechanics and variational integrators, Acta Numer., 10 (2001), 317–514. https://doi.org/10.1017/S096249290100006X doi: 10.1017/S096249290100006X
    [28] K. R. Meyer, Symmetries and integrals in mechanics, Dynamical Systems (ed. M. M. Peixoto), Academic Press, 1973,259–272. https://doi.org/10.1016/B978-0-12-550350-1.50025-4
    [29] J. Moser, A. Veselov, Discrete versions of some classical integrable systems and factorization of matrix polynomials, Commun.Math. Phys., 139 (1991), 217–243. https://doi.org/10.1007/BF02352494 doi: 10.1007/BF02352494
    [30] A. Natale, C. Cotter, A variational H(div) finite-element discretization approach for perfect incompressible fluids, IMA J. Numer. Anal., 38 (2017), 1388–1419. https://doi.org/10.1093/imanum/drx075 doi: 10.1093/imanum/drx075
    [31] H. Parks, M. Leok, Variational integrators for interconnected Lagrange–Dirac systems, J. Nonlinear Sci., 27 (2017), 1399–1434. https://doi.org/10.1007/s00332-017-9364-7 doi: 10.1007/s00332-017-9364-7
    [32] D. Pavlov, P. Mullen, Y. Tong, E. Kanso, J. Marsden, M. Desbrun, Structure-preserving discretization of incompressible fluids, Phys. D, 240 (2011), 443–458, https://doi.org/10.1016/j.physd.2010.10.012 doi: 10.1016/j.physd.2010.10.012
    [33] S. Smale, Topology and mechanics. Ⅰ, Inventiones Math., 10 (1970), 305–331. https://doi.org/10.1007/BF01418778 doi: 10.1007/BF01418778
    [34] W. M. Tulczyjew, Geometric Formulations of Physical Theories: Statics and Dynamics of Mechanical Systems, vol. 11 of Monographs and Textbooks in Physical Science Lecture Notes, Bibliopolis, 1989.
    [35] A. J. van der Schaft, Port-Hamiltonian systems: an introductory survey, in International Congress of Mathematicians. Vol. III, Eur. Math. Soc., Zürich, 2006, 1339–1365.
    [36] J. Vankerschaver, Euler–Poincaré reduction for discrete field theories, J. Math. Phys., 48 (2007), 032902. https://doi.org/10.1063/1.2712419 doi: 10.1063/1.2712419
    [37] J. Vankerschaver, F. Cantrijn, Discrete Lagrangian field theories on Lie groupoids, J. Geom. Phys., 57 (2007), 665–689. https://doi.org/10.1016/j.geomphys.2006.05.006 doi: 10.1016/j.geomphys.2006.05.006
    [38] V. S. Varadarajan, Lie Groups, Lie Algebras, and Their Representations, Graduate Texts in Mathematics, Springer New York, 2013.
    [39] A. Weinstein, Lagrangian mechanics and groupoids, Fields Inst. Proc. AMS, 7 (1996), 207–231.
    [40] H. Yoshimura, J. Marsden, Dirac structures in Lagrangian mechanics Part Ⅰ: Implicit Lagrangian systems, J. Geom. Phys., 57 (2006), 133–156. https://doi.org/10.1016/j.geomphys.2006.02.009 doi: 10.1016/j.geomphys.2006.02.009
    [41] H. Yoshimura, J. Marsden, Dirac structures in Lagrangian mechanics Part Ⅱ: Variational structures, J. Geom. Phys., 57 (2006), 209–250. https://doi.org/10.1016/j.geomphys.2006.02.012 doi: 10.1016/j.geomphys.2006.02.012
    [42] H. Yoshimura, J. Marsden, Dirac cotangent bundle reduction, J. Geom. Mech., 1 (2009), 87–158. https://doi.org/10.3934/jgm.2009.1.87 doi: 10.3934/jgm.2009.1.87
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1611) PDF downloads(60) Cited by(2)

Figures and Tables

Figures(4)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog