Seasonality Impact on the Transmission Dynamics of Tuberculosis

The statistical data of monthly pulmonary tuberculosis (TB) incidence cases from January 2004 to December 2012 show the seasonality fluctuations in Shaanxi of China. A seasonality TB epidemic model with periodic varying contact rate, reactivation rate, and disease-induced death rate is proposed to explore the impact of seasonality on the transmission dynamics of TB. Simulations show that the basic reproduction number of time-averaged autonomous systems may underestimate or overestimate infection risks in some cases, which may be up to the value of period. The basic reproduction number of the seasonality model is appropriately given, which determines the extinction and uniform persistence of TB disease. If it is less than one, then the disease-free equilibrium is globally asymptotically stable; if it is greater than one, the system at least has a positive periodic solution and the disease will persist. Moreover, numerical simulations demonstrate these theorem results.


Introduction
Tuberculosis (TB) remains one of the world's deadliest communicable diseases. In 2013, it was estimated that 9.0 million people developed TB and 1.5 million died from the disease, and TB is slowly declining each year [1]. According to the online global TB data collection system, China alone accounted for 11% of the total cases, which is the second country of 22 TB high-burden countries, only after India. Shaanxi is one of the more serious TB provinces in China, and its reported new cases each year reach about 25,000, which is the second in the number of cases of infectious diseases in Shaanxi, only next to the hepatitis B virus (HBV).
Some researchers have investigated the influence of seasonal variations on the transmission dynamics of infectious diseases [2][3][4]. And seasonal variation in TB incidence has been described in many countries and cities, such as India, United States, Russia, New York city, and Hong Kong [5][6][7][8][9]. TB is a seasonal disease in China [10,11], but it remains unknown in Shaanxi.  (Table 1) are available on the data-center of China public health science [12]. We apply seasonal filters in MATLAB program, to deseasonalize the time series; then the original time series is decomposed into three components: trend curve and season and irregular noise. The trend curve is the long-term and medium-to-long term movement of the series; it also contains consequential turning points; the seasonal component is within one-year (12 months) fluctuations about the trend that recur in a similar way in the same month or quarter every year; and the irregular component is the residual component that still remains after trend curve and seasonal component are removed from the original series. of this multiplicative decomposition, the original series are compared to a series reconstructed using the component estimates in Figure 3(a) and giving irregular noise component in Figure 3(b). It has shown that the multiplicative decomposition is fit for the TB data of Shaanxi. In Figure 2(a), trend component shows that there is a fast upward trend after 2004 and in 2005 (the data in 2004 is low just because this is the first year for reported data, and some staff may just begin to the report system, so some patients may be omitted), then a downward trend from 2005 to 2007, and then a slowly upward in 2008, a steadily decreasing trend from 2009 to 2012. In Figure 2(b), firstly, seasonal component shows that seasonal amplitude decreases each year from 2004 to 2012; secondly, it illustrates there exists a seasonal period , 12 months, and along with peak and trough months: the first-and second-peak month of TB notification in March and January, respectively, and between them, there exists a trough, February, between these two months (February may be a peak month, but for Chinese lunar new year, Spring Festival, there exists a special reason that patients may not choose diagnosis for regarding illness as an unlucky thing in this traditional festival); the trough month is in December. TB is a seasonal disease in Shaanxi. In addition, the peak and trough of TB transmission actually are in winter and in autumn, respectively, due to the delay which tends to last 4-8 weeks. Understanding TB seasonality may help TB programs better to plan and allocate resources for TB control activities [8]. Motivated by this, we formulate a seasonality epidemic model with periodic coefficients for TB and study its global dynamics in this paper. This paper is organized as follows. In Section 2, a seasonality TB model is formulated. In Section 3, a unique diseasefree equilibrium is obtained, and the basic reproduction number for the periodic model is given in detail. Furthermore, some numerical simulations are used to compare the average basic reproduction number [ ] to in different cases. Threshold dynamics of the TB model is analyzed in Section 4. Meanwhile, some numerical simulations are provided to validate analytical results. Finally, conclusions are given in Section 5.
From Figure 2(b), new TB cases have shown the periodic monthly trend and the possible causes of the seasonal pattern. Shaanxi has a continental monsoon climate and has four distinct seasons. Seasons in Shaanxi are defined as spring (February-April), summer (May-July), autumn (August-October), and winter (November-January). TB is usually acquired through airborne infection from active TB cases; its transmission and progress tend to be effected by the climate within one year (12 months). In particular, the indoor activities are much more in winter than in a warm climate, which improves the probability of susceptible individuals exposed to Mycobacterium tuberculosis (Mtb) from the infectious individuals in a room with windows closed for a longer period of time [13]; thus infection rate may have the periodic   influence. In addition, during these months near or in the highest peak month of TB cases, cold weather and lack of sunshine, which may reduce human immunity, cause a higher disease-induced rate. Individuals with lower Vitamin D level may be more reactivated for TB [14]. Thus, diseaseinduced rate and reactivation rate may also have the periodic influence. To describe and study the TB transmission in Shaanxi, three periodic coefficients are selected: (i) infection rate coefficient ( ) and the bilinear incidence ( ) which are applied in this model; (ii) reactivation rate coefficient ( ), at which an individual leaves the latent compartment for becoming infectious; and (iii) disease-induced rate coefficient ( ), which is the disease-induced death rate coefficients for individuals in compartment . In view of the periodic trend of monthly, ( ), ( ), and ( ) are assumed to be positive periodic continuous function of with period . In some TB models, the fast and slow progression has been considered [15][16][17]. Based on those, a seasonality TB model with fast and slow progression and periodic coefficients is formulated in this section. The transfer among compartments is schematically depicted in Figure 4. It leads to the following model of ordinary differential equations: with initial condition ( (0), (0), (0), (0)) = ( 0 , 0 , 0 , 0 ) ∈ R 4 + , and all parameters are positive. Here, Λ is the recruitment rate and parameter is the natural death rate coefficient; is the rate coefficient at which an infective individual leaves the infectious compartment to the recovered; (0 ≤ < 1) is the fraction of infected individuals who are fast developing into infected cases and enter the infectious compartment directly, while 1 − is the fraction of infected individuals who are slowly developing into infected cases and transferred to the latent compartment.
In the rest of this paper system (2) will be studied in region Ω.

Disease-Free Equilibrium and the Basic Reproduction Number
To study system (2), some notations are introduced.
Let ( ) be a continuous, cooperative, irreducible, and × matrix function with period > 0, and let Φ ( ) be the fundamental solution matrix of the linear ordinary differential equation: And let (Φ ( )) be the spectral radius of Φ ( ). By Perron-Frobenius theorem, (Φ ( )) is the principle eigenvalue of Φ ( ), in the sense that it is simple and admits an eigenvector There is a unique disease-free steady state 0 , that is (Λ/ , 0, 0), for system (2).
In the following, the basic reproduction number will be introduced for system (2) according to the general procedure presented in [2].
With fl ( , , ), system (2) becomes where , Furthermore, here comes Then ( ) is nonnegative, and − ( ) is cooperative in the sense that the off-diagonal elements of − ( ) are nonnegative. Thus, it is easy to verify that system (2) satisfies the assumptions (A1)-(A7) in [2]. Define ( , ), ≥ , which is a 2 × 2 matrix, and is the evolution operator of the linear -periodic system That is, for each ∈ R, ( , ) satisfies where is the 2 × 2 identity matrix. Thus, the monodromy matrix Φ − ( ) of (9) equals ( , 0), ≥ 0. Assume that the population is near the disease-free periodic state 0 . And suppose that ( ), -periodic in , is the initial distribution of infectious individuals. Then ( ) ( ) is the rate of new infections produced by the infected individuals who were introduced at time . Given ≥ , ( , ) ( ) ( ) gives the distribution of those infected individuals who were newly infected at time and remain in infected compartments at . It follows that is the distribution of accumulative new infections at time produced by all those infected individuals ( ) introduced at time ( ≤ ). Let be the ordered Banach space of allperiodic functions from R to R , which is equipped with the maximum norm ‖ ⋅ ‖ ∞ and the positive cone + = { ∈ : ( ) ≥ 0, ∈ R}. Define a linear operator : → by Then, according to Wang and Zhao [2], the basic reproduction number is defined as for the periodic epidemic system (2), where ( ) denotes the spectral radius of the matrix .
Computational and Mathematical Methods in Medicine   5 In the constant case, that is, By van den Driessche and Watmough [18], here comes In the periodic case, in order to characterize , consider the following linear -periodic equation: with parameter ∈ (0, ∞). Let ( , , ) be the evolution operator of system (16) (2), the following statements are valid: (i) If ( ( , 0, )) = 1 has a positive solution 0 , then 0 is an eigenvalue of , and hence > 0.
And the threshold behaviors occur.  Thus, the disease-free periodic solution 0 for system (2) is locally asymptotically stable if < 1 and unstable if > 1.
Define  An example is given to show that the basic reproduction number of the time-averaged autonomous systems may underestimate, estimate, or overestimate the infection risk.
Now Lemma 2 is applied to calculate the basic reproduction number of system (2).     The trends of R T and [R T ] as g 0 is changing superior, inferior, or equal to the basic reproduction number , which is up to value of period . Furthermore, in the next section, we will prove some theoretical results of system (2), in which serves as a threshold parameter: if < 1, then there exists a globally asymptotically stable disease-free periodic state 0 (Λ/ , 0, 0); if > 1, then the disease is persistent in the population and there exists at least one positive periodic solution.

Extinction and Uniform Persistence
The following lemma is useful for our discussion in this section.

Theorem 7.
If > 1, system (2) is uniformly persistent, and there exists at least one positive periodic solution.
It is easy to see that Ω and Ω 0 are positively invariant, Ω 0 is a relatively closed set in Ω, and is point dissipative from Theorem 1.
We claim that for any > 0.

Discussion
Monthly pulmonary TB cases, from January 2004 to December 2012 in Shaanxi province, have been analyzed by the seasonal adjustment program. It has been found that TB cases show seasonal variation in Shaanxi: the peak months (January and March) compare to the lowest trough month (December). It is necessary to study the seasonality TB epidemic model according to the seasonal component of Shaanxi's data. Considering the regularity of peak TB seasonality may help allocate resources for the prevention and treatment of TB activities, a periodic TB epidemic model has been formulated and studied.
The basic reproduction number for the periodic model is given. By numerical simulation, has been compared to the average basic reproduction number [ ] in different parameter values. It has shown that [ ] may overestimate or underestimate or be equal to in different cases. It is also up to the value of periodic . Furthermore, in theory, the threshold dynamics has been studied for the periodic TB model: if < 1, then the disease-free equilibrium is globally asymptotically stable; that is, the TB disease will disappear eventually; if > 1, then there exists at least one positive periodic solution and the disease will be uniformly persistent.
Our numerical simulations are in good accordance with these theoretical results. It should be noted that we have   Figure 18; then = 2.251. These figures show that the disease will be asymptotic to a periodic solution, which is the same as Theorem 7.
not fit Shaanxi's data in these simulations since we cannot accurately estimate some parameters' values in the periodic model according to available data or references. Despite lack of comparing the model results with the Shaanxi's data, our theoretical results have shown that the basic reproduction number with periodic plays a crucial role in determining dynamics of the seasonality TB disease and could be used for controlling the spread of TB epidemic in reality strategies.