Baptism of Fire: Modeling the Effects of Prescribed Fire on Lyme Disease

Recently, tick-borne illnesses have been trending upward and are an increasing source of risk to people's health in the United States. This is due to range expansion in tick habitats as a result of climate change. Thus, it is imperative to find a practical and cost-efficient way of managing tick populations. Prescribed burns are a common form of land management that can be cost-efficient if properly managed and can be applied across large amounts of land. In this study, we present a compartmental model for ticks carrying Lyme disease and uniquely incorporate the effects of prescribed fire using an impulsive system to investigate the effects of prescribed fire intensity (high and low) and the duration between burns. Our study found that fire intensity has a larger impact in reducing tick population than the frequency between burns. Furthermore, burning at high intensity is preferable to burning at low intensity whenever possible, although high-intensity burns may be unrealistic due to environmental factors. Annual burns resulted in the most significant reduction in infectious nymphs, which are the primary carriers of Lyme disease.


Introduction
Many ticks are disease vectors that significantly impact public health. Reports of overall tick-borne diseases doubled from 2006 to 2018 [1] while the incidence of Lyme disease in the United States has been steadily increasing, from a little less than four cases per 100,000 people in the 1990s to close to 10 cases per 100,000 people in the early 2000s [2]. New pathogens continue to emerge, including heartland virus, Bourbon virus, Borrelia miyamotoi, Borrelia mayonii, and Ehrlichia muris eauclairensis [3]. Climate change has expanded the northern borders of tick habitats and increased winter tick activity, contributing to the prevalence of tickborne diseases [4]. erefore, finding a practical and costefficient way to manage tick populations has become extremely important. e majority of ticks that carry Lyme disease are infected through mice or other small rodents [5], so most methods that have looked at tick reduction are focused on either host reduction or tick elimination [6]. Ticks can have up to 2-to 3-year long life cycles, with four general stages that each requires one blood meal after the eggs hatch: egg, larvae, nymph, and adult [6]. e majority of human infections come from tick nymphs, which are much smaller than adult ticks (less than 2 mm long or about the size of a poppy seed), making them more difficult to spot on the human body and therefore more likely to remain undetected [6]. Nymphs are also more numerous than adult ticks and are most active during the spring and summer months, when the number of people who spend time outside is substantially larger than those during other months. e black-legged tick, also known as Ixodes scapularis, has a life cycle that generally lasts two years, while the life cycle of the lone star tick (Amblyomma americanum) is around three years long. Most ticks hatch from their eggs in the spring and have the ability to live for three to five months between each blood meal [7].
Prescribed fires, or controlled burns, are a common and necessary form of land management in many different environments that are also effective in controlling tick populations. is is through both directly killing ticks along with destroying their leaf litter habitat [7]. Larvae, nymphs, and adults spend the vast majority of their time in leaf litter other than the few days that they are feeding on their hosts. Controlled burns are appealing due to their time and cost efficiency along with their ability to be applied across a large amount of land. ey are generally most effective in the late spring and early summer, as that time coincides with when nymph ticks are questing for hosts (although this is heavily dependent on the type of land that is being burned) [8]. Primary concerns around prescribed fire include air quality (due to smoke) and the potential for the fire to burn out of control; however, these can be prevented when proper precautions are taken.
Many studies have looked at the impact of prescribed fires on tick populations, with conflicting results. e majority of these studies agree that tick populations decrease immediately after a burn but recover to pre-burn abundance after around one year [6]. Other studies have found that although the nymph population decreased, the risk of encountering infectious nymphs remained the same [9] or that the tick population even increased [10]. However, these studies often fail to account for the logistics of true prescribed burning (long term and over lots of lands on a regular basis) or other predictors of tick abundance such as host abundance, climate, or vegetation structure [6]. ere is also the possibility that postburn recolonization rates vary based on tick species, habitat type, climate, and burn intensity [11,12].
A study done by Allan [11] in the oak-hickory ecosystems of the Missouri Ozarks looked at the relationship between lone star tick larval populations and deer abundance under long-term burn management. e sites were burned in the spring at low intensity every 3-5 years. e ticks were depleted but then rapidly grew starting two years postburn, coming back down to pre-burn abundance around five years postburn [11]. e researchers attributed this increase to the high host populations postburn, as freshly burned areas are better for deer to forage in. ese issues could be countered by more frequent, longer, and larger-scale burns, which correlate with other studies that also believe that burns at higher intensity are most effective in countering ticks than those at low intensity [13]. Gleim et al. [3] found that long-term prescribed fire (regular burning for 10+ years) significantly reduced tick abundance, regardless of burn interval, host abundance, or vegetation structure. is is primarily due to the change in vegetation structure, creating a hotter and drier environment that is less appealing for ticks [6]. ese burnings decreased the encounter rate with infectious ticks by 98% in plots in southwestern Georgia and northwestern Florida. However, more research in a variety of environments needs to be done regarding realistic prescribed burning as a tick management technique. e goal of this study was to develop a compartmental model for ticks carrying Lyme disease to see how they are affected by prescribed burns. Compartmental epidemic models are a modeling technique that segments individuals in the populations into different classes according to individual's disease status [14]. We look at both fire intensities (high and low) and the duration between fires to understand how this common land practice affects tick populations and the prevalence of Lyme disease among them. We also investigate whether intensity or duration plays a more significant role in tick population reduction overall. To the best of our knowledge, this is the first study to use a mathematical model for Lyme disease to examine the effects of prescribed fire. A preprint of this study has already been published [15]. e remainder of the work in this study is organized as follows. In Section 2, we formulate our baseline tick/Lyme disease model, compute the model basic reproduction number, and carry out basic stability analysis including sensitivity analysis to determine the parameter with the most impact on the basic reproduction number. In Section 2.3, we describe the tick model with the effect of prescribed fire using an impulsive system of ordinary differential equations and present some stability analysis results of the impulsive system. In Section 2.3.1, we discuss the estimation of parameters related to prescribed fire from literature. In Section 3, we present some simulation results, and in Section 4, we discuss our findings and close with conclusions.

Materials and Methods
is model was created by incorporating two subgroups: mice and ticks. e mice population is divided into susceptible (S M (t)) and infected mice (I M (t)). e tick population is divided by life stage (eggs, larvae, nymph, and adult) and further divided into susceptible and infected groups for larvae (S L (t) and I L (t)), nymphs (S N (t) and I N (t)), and adults (S A (t) and I A (t)). Since ticks must take a blood meal before they become infected and there is no vertical transmission for the disease, all eggs remain susceptible (S E (t)). Individuals move between compartments according to their life stage and disease status. A detailed description of the model formulation is given in Appendix A. e related model schematic is given in Figure 1, and the description of the model variables and parameters is stated in Table 1.
e basic qualitative properties of the tick model shown in Figure 1 and described in Appendix A, its positivity, and the boundedness of solutions are given in Appendix B.

Reproduction Number
e associated reproduction number, denoted by R 0 , obtained using the next-generation matrix method [19,20] for the Lyme disease model (A-5) shown in Figure 1 and described in Appendix A, is given as follows: However, we made a simplifying assumption that μ L � μ N � μ A � μ T . en, the reproduction number in (1) becomes See details in Appendix C. e basic reproduction number, R 0 , is defined as the expected number of new infections that result from one infectious individual in a population that is fully susceptible [19,20].
is value is extremely significant because if the reproduction number is less than unity (R 0 ≤ 1) then the disease cannot invade the population and it will die out in the community. Conversely, if R 0 > 1 then the disease will continue to persist in the population. is determines whether there is a possibility of disease elimination or whether the goal should be to manage transmission within the community.

Sensitivity
Analysis. Sensitivity analysis was conducted to determine the contribution of each of the model parameters to the reproduction number R 0 . e results of this help identify which parameters are the best to target regarding interventions and future data collection. A normalized forward sensitivity index was used to determine the ratio of the relative change in R 0 based on a relative change in a parameter. e sensitivity indices of R 0 are derived as follows: where Υ R 0 p is the forward sensitivity index of R 0 with respect to parameter p. Parameter p is a parameter within R 0 . e results of this sensitivity analysis are shown in Figure 2. e larval development rate, the carrying capacity of the environment, and the disease transmission probability for both mice and ticks all positively affect R 0 . As these values increase, R 0 will also increase. We are more interested in the mice birth rate and the tick death rate, which negatively affect R 0 . As these parameters decrease, R 0 will also decrease. e most influential parameter on R 0 is the tick death rate. is suggests that control strategies that effectively target the spread of Lyme disease will focus on increasing the tick death rate and, to a lesser extent, also decreasing the mice birth rate. e sensitivity analysis aligns with our results from the model, as prescribed fire is a mechanism that most significantly affects the tick death rate.

Tick Model with Prescribed
Fire. In this section, we consider the effect of fire on ticks and the small mammal population. We do not explicitly incorporate fire into the model shown in Figure 1 and described in Appendix A; rather, we consider the effect of fire on population size Figure 1: Flow diagram of the Lyme disease model for ticks and mice. e mice population is represented by the red boxes and is divided into susceptible (S M (t)) and infected (I M (t)) compartments. e tick population is represented by the blue boxes and is divided by life stage and infection status. It consists of susceptible eggs (S E (t)), susceptible larvae (S L (t)), infected larvae (I L (t)), susceptible nymphs (S N (t)), infected nymphs (I N (t)), susceptible adults (S A (t)), and infected adults (I A (t)).
after the burns. To introduce the effect of prescribed fire into the Lyme disease model, we have the following system of nonautonomous impulsive differential equations.
Subject to the prescribed fire impulsive condition: where t n is the times that prescribed fire is implemented, which may be fixed or non-fixed; within this study, we will consider the case with fixed times. e parameters ] j , where j � E, L, N, A, M are the proportion of the tick and mice population that is reduced by the fire. In Section 2.3.1, we discuss how these parameters are estimated using data from low-and high-intensity fires. e existence and stability of the impulsive model (3) are given in Appendix D.

Prescribed Fire Parameter Estimation.
To estimate the parameters (] L , ] N , ] A , ] M ) that quantify the reduction in each tick life stage and mice population after the different burn intensities (low and high burns), we use data from [13,21]. e parameters for each life stage relating to the low-intensity burn were estimated using data from [13], while the parameters relating to the high-intensity burn were estimated using data from [21]. Each group of parameters was separated into the larval, nymph, and adult life stages, but only data in [21] provided data to estimate the parameter for mice population. Below, we give a summary description of the study sites in each study, the amount of ticks and mice collected, and how these parameters are estimated from the data collected.
(1) High-Intensity Fire. e study in [21] was conducted in chaparral habitat at the University of California, Hopland Research and Extension Center, in Mendocino County, CA. e study took advantage of two prescribed fires ignited on June 1, 1995, that were intended to reduce fire load in two chaparral plots, Maude's Glade (MG) and Don's Brush Plot (DBP). e fires were ignited by hand crews using drip torches in a strip headfire configuration to produce relatively uniform fire behavior that left no live branches in the shrub line, which we assume as high-intensity fire. For a period of 13 months beginning a month before the burn, control and treatment areas were monitored for the presence of ticks by flagging the vegetation or ground between 0800 and 1000 hours. Each flag sample involved sweeping a standard 1 m 2 white flannel tick flag in a 180°arc from a particular bush. e number of flag samples varied between 24 and 184 samples per day per plot. Ticks collected from these flag samples had their species and sex recorded and then that were released at the site of capture. Carbon dioxide (CO 2 ) trapping was done monthly, starting one month before the burn, and ending one year after the burn date. Ten enamelware pans were buried at 15 m intervals as pitfall traps. 5-10 g of dry ice was added to the center of the pan between 1000 and 1200 hours at control sites and between 1200 and 1400 hours at treatment sites. Ticks were collected 2 hours after the   Table 1: Description of the variables and parameters for the Lyme disease model in Figure 1.
Number Nymph-to-adult development rate of (adult development rate) 0.00491 [16] addition of dry ice and had their numbers, species, sex, and life stage recorded and then that were released at the site of capture. To assess the abundance of rodents and the associated ticks on them, live traps were set to catch the rodents. Non-folding Sherman traps were set at 7.5 m intervals along the ecotones of the area, starting one week before the burn and continuing for two simultaneous nights each month throughout the duration of the study. At the two study sites (MG and DBP), six tick species (namely Ixodes pacificus, Ixodes jellisoni, Ixodes spinipalpis, Ixodes woodi, Dermacentor occidentalis, and Dermacentor parumapertus) were removed from the six different rodents species caught; these include California kangaroo rat (Dipodomys californicus californicus), brush mouse (Peromyscus boylii), pinyon mouse (P. truei sequoiensis), deer mouse (P. maniculatus gambelii), dusky-footed woodrat (Neotoma fuscipes), and western harvest mouse (Reithrodontomys megalotis longicaudus). After the fire treatment, about half as many rodents were trapped at the treated sites compared with control sites.
All the ixodid tick species (Ixodes pacificus, Ixodes jellisoni, Ixodes spinipalpis, Ixodes woodi) are competent hosts that are able to transmit Lyme disease, but we use only the data of Ixodes pacificus to estimate the parameter used to quantify the reduction in the tick population as a result of fire. Table 2 gives a breakdown of Ixodes pacificus collected on the rodents in the control and treated areas at both MG and DBP sites, as well as the number of rodents collected at these sites pre-and postburn.
(2) Computing Parameters. ] L , ] N , ] A , ] M for High-Intensity Fire. To compute these parameters, we assume that equal numbers of ticks and rodents are in the sites (control and treated sites) pre-burn, and the difference in number is due to the burn, since we have no way to measure exactly the number of ticks in all the study sites. First, we take the difference between the total number of ticks and mice in the control and treatment sites and divide it by the total number of ticks and mice in the control sites. en, we subtract these proportions from 1 to give the proportions reduced as a result of the burn: adult: (3) Low-Intensity Fire. e study in [13] was conducted in an open oak woodland barren complex at Western Illinois University's Alice L. Kibbe Field Station located in Warsaw, in Hancock County, IL, USA. ese areas are comprised of multiple habitat types including oak-hickory woodlands, early successional woodlands, oak barrens, floodplain forests, restored tallgrass prairies, and hill prairies. e entire study site was last burned in 2004 (B04), and two additional burns were carried out in spring of 2014 (B14) and 2015 (B15). e burns were considered low intensity because most flame heights were less than 1 meter and plant mortality was limited to the understory vegetative community [22].
Ticks were collected through a flagging method similar to that of the high-intensity fire, where a 1 m 2 flannel cloth was attached to a dowel (although they used bamboo instead of wood) and swept through the vegetation. is was done every two weeks when the vegetation was dry between 1200 and 1800 hours during two consecutive years (May 9, 2015, to October 30, 2015, and April 22, 2016, to November 4, 2016. e cloth flag was then placed into a sealed bag and frozen for three days to ensure tick mortality. e ticks were removed and identified using taxonomic keys; ticks collected in 2015 and 2016 also had their DNA extracted. A total of 2788 Amblyomma americanum, 54 Ixodes scapularis, and 23 Dermacentor variabilis ticks were collected in 2015 and 2016. Amblyomma americanum ticks collected in B04 made up 51% of the collection (n � 1433), while those collected in B14 made up 37% (n � 1045) of the collections, and those collected in B15 constituted 11% (n � 307) of the collection. Of these ticks, 2% (n � 67) were adults, 4% (n � 107) were nymphs, and 93% (n � 2614) were larvae. Of the 23 D. variabilis collected, 74% (n � 17) were adults, 9% (n � 2) were nymphs, and 17% (n � 4) were larvae, while 4% (n � 2) of the 54 I. scapularis collected were adults, 22% (n � 12) were nymphs, and 74% (n � 40) were larvae. Here, we use the data collected for Ixodes scapularis to estimate the parameter quantifying the reduction in the tick population due to fire since I. scapularis can transmit Lyme disease. e study did not indicate whether these were data for the pre-burn or postburn number of I. scapularis collected, nor did it provide data for the number of mice that were caught. To estimate these parameters, we assume that the 54 I. scapularis collected are the ticks left after the burn. We then divide the numbers collected in each age group by the total number collected and subtract the proportion obtained from 1 to obtain the proportion reduced by fire.
Although these were the best sources of data we could find (that accounted for tick life stages pre-and postburn along with fire intensity), unfortunately, they were not ideal.
e burns were conducted in vastly different environments-the high-intensity fire took place in California chaparral, while the low-intensity fire took place in Illinois oak woodland. e high-intensity burn was only conducted once in the summer, while the low-intensity burn was conducted in the spring for two consecutive years. We were unable to find any sources of data from similar geographic locations or number of burns that were detailed enough in terms of tick data and fire intensity data to be used. ese variances between the two data sets have the potential to drastically affect the burn results.

Results
To address our research goals, we started by looking at how different burn frequencies and intensities affected tick populations, focusing specifically on the nymphs, since the primary mode of transmission for Lyme disease from ticks to humans is through infectious nymphs. High-intensity fires substantially change the aboveground structure, with no live branches and few shrub skeletons leftover [12,22]. ese burns are much more uniform than those at lower intensities, which are patchy and have vegetation cover and woody debris leftover that ticks are able to survive the burn-in [13]. Our first simulation in Figure 3 shows the substantial difference between high-intensity and low-intensity burns regarding how effective they are at reducing the infectious nymph population.
Although both burns start out at around the same effectiveness regardless of intensity, the high-intensity burn proves to reduce the infectious nymph population much more efficiently than the low-intensity burn. e infectious nymph population affected by the low-intensity fire continues to increase despite the burns, while the infectious nymph population affected by the high-intensity fire remains consistently low. In our second simulation shown in Figure 4, we consider the duration between burns for both high-and lowintensity fires. We consider the effects of burning for a period of six years once every six years, once every three years, once every other year, and annually. We found in general that as the duration between burns decreases the infectious nymph population also decreases regardless of the burn intensities. Furthermore, we found that the duration between burns has a more significant effect on ticks with higher-intensity fires than with lower-intensity fires. Fire intensity appears to have a larger influence on tick reduction than the duration of the burns, as burning fewer times at a higher intensity is more effective than burning more times at a lower intensity. For example, high-intensity burns once every three years reduce the infectious nymph population more than low-intensity burns once every two years. However, high-intensity burns might be unrealistic due to environmental factors. In that case, annual burns at low intensity result in the most significant reduction in infectious nymphs in its category.

Discussion.
e most effective simulation in minimizing infectious nymph populations is having annual high-intensity burns, although this is not always practical due to environmental factors such as weather, burn location, and fuel loads. e geographic location of the burn is also  [21] under assumed high-intensity fire. C � control and B � burn. DBP  79  101  3  1  14  5  56  27  MG  80  17  0  2  8  7  54  25  Total  159  118  3  3  22  12 Table 1.

Low intensity
High intensity  Canadian Journal of Infectious Diseases and Medical Microbiology important to consider, as it is riskier to have higher-intensity fires closer to human settlements. Most prescribed fires are at low to moderate intensity and are repeated every 1-5 years, with their primary objective being to reduce fuel loads and act as a wildfire prevention tool [21]. ere have been a few negative outcomes associated with intense annual burning such as oak tree mortality, soil compaction, and an increased number of tree cankers and tree colonization by root fungus [23], but these depend heavily on the type of environment that is being burned. Overall, infectious tick populations substantially decrease as the time between burns also decreases, regardless of the intensity that the burn is conducted at. ese findings were confirmed in a private conversation with Gallagher of the U.S. Forest Service Research and Development team. He recently finished working on a study that looked at how long-term prescribed burns affected tick populations in New Jersey oak woodlands. ey found preliminary data that strongly suggest that different tick species are impacted differently by prescribed burns due to their varying moisture sensitivities [12]. It appears that blacklegged ticks, the primary vector of Lyme disease in humans, are the most sensitive to moisture loss, while Lone Star ticks and Gulf Coast ticks seem to be more resilient. Gallagher believes that one of the reasons for this is the thickness of the scutum, which is a shield-like plate on the back of hardback ticks [12]. In general, initial data from this study agree with the results of our model. eir high-intensity burns, occurring once every 20 years with data collection for four consecutive years after the burn, saw the greatest reduction in tick populations. Annual low-intensity burns for over 25 years had the second largest tick reduction, while the single low-intensity burn showed the smallest amount of tick reduction.
Few research studies have examined the effects of prescribed fire on tick populations and disease prevalence [11]. e studies do often choose to focus on one aspect, usually that of how the burns affect the ticks, with the fire itself being of secondary interest [7]. Future research that records detailed information about fires, such as flame height, difference in vegetation presence before and after the burn, how the fire was ignited, weather conditions during the burns, and overall fire behavior, would be extremely helpful in ensuring the accuracy of the parameters used in the model. Knowledge about the exact dates of the burns and the specifics of the duration between burns is also useful, as seasonality has the potential to play a large role in the effectiveness of burns on disease prevalence. is is especially true with diseases that are primarily found to be transmitted to humans in a certain life stage of the tick, such as with nymphs and Lyme disease. Knowing the exact dates of the burns along with the time in between is useful when comparing the burn results to tick life cycles. For example, spring burns will mostly kill nymphs, although this can vary by species. Noting the geographic context is also important, as the environment plays a role in how these burn impact tick populations. Finally, increasing the time frame of studies that look at how prescribed fire impacts ticks will also be useful. Many studies that claim to study how ticks are affected by prescribed fires only involve a single burn, which is vastly different from true prescribed burning, which involves regular fires over many years.
Land geography seems to have a role in the effectiveness of prescribed burns on tick population reduction regarding how fast the environment is able to regrow back to its original state before the burn, although to the best of our knowledge there have not been any studies that specifically look at this. Certain environmental factors might make certain geographic locations better candidates for prescribed fires when the main goal is to reduce tick populations, and this is a subject that is worth looking into in the future. It also appears that different tick species have different levels of heat and moisture resistance, causing them to be affected differently by the fires [12,13]. Only very preliminary research has been done on this so far, although it is a topic that will substantially influence the choice of whether or not to use prescribed fires for tick reduction.

Conclusion.
To conclude, in this study we develop a simple model for Lyme disease transmission and used it to investigate the impact of prescribed burn frequency and fire intensity on the spread of Lyme disease. To the best of our knowledge, this is the first model to incorporate the effect of fire into a mathematical model of ticks and Lyme disease.
e key findings from this study are summarized as follows. e simulations of the Lyme disease model (A-5) with prescribed burns show that the following: (i) e most influential parameters impacting the reproduction number, R 0 , from the sensitivity analysis are π M , σ T , K, μ T , β T , β M . (ii) Intensity appears to be a larger influence on tick reduction than duration is-it is better to burn fewer times at a high intensity than to burn more often at a lower intensity. (iii) Burning at high intensity is preferable to burning at low intensity whenever possible, although the high intensity may be unrealistic due to environmental factors. (iv) For low intensity, annual burns resulted in the most significant reduction in infectious nymphs, which are the primary carriers of Lyme disease.
is study has shown the effect of prescribed fire on ticks and the spread of Lyme disease. In a future study, we will consider the impact of seasonality on the effectiveness of prescribed burns and the impact of the timing of the burns on disease spread. Geographic landscape and tick habitat are important factors affecting the tick population. In another future work, we will consider different types of geographic landscapes and determine which landscape is most effective in using prescribed burns for tick population control and reduction.

A. Model Formulation
is model was created by incorporating two subgroups: mice and ticks. e mice population is divided into susceptible (S M (t)) and infected mice (I M (t)). e tick population is divided by life stage (eggs, larvae, nymph, and adult) and further divided into susceptible and infected groups for larvae (S L (t) and I L (t)), nymphs (S N (t) and I N (t)), and adults (S A (t) and I A (t)). Since ticks must take a blood meal before they become infected and there is no vertical transmission for the disease, all eggs remain susceptible (S E (t)). Individuals move between compartments according to their life stage and disease status. We assume that all transition rates are of the current population and remain steady, with no migration into or out of the overall population. e force of infection in the mice (or the rate that susceptible mice become infected) is given as follows: where the parameter β M is the probability that infection will occur if a mouse is bitten by an infectious tick, multiplied by the number of all infectious ticks-the sum of the larvae, nymphs, and adults-and then divided by the total number of mice where N M � S M + I M . For simplicity, we assume that there is a homogenous mixing of both mice and tick populations. New mice are born at a rate of π M , and the susceptible mice move into the infected compartment at a rate of λ M . Since the disease does not affect the mice, they remain in the infected compartment for the rest of their lives and death related to the disease is not incorporated into the model. erefore, natural death, given as the rate μ M , is the only factor decreasing the population of both susceptible and infected mice. e equations for the susceptible mouse population and infected mouse population are shown as follows: e force of infection in the ticks (or the rate that ticks become infectious) is given as follows: where the parameter β T is the probability that infection will occur if a tick bites an infected mouse, multiplied by the total number of infected mice, and divided by the overall mice population (both susceptible and infected). Since there is no vertical transmission of the disease between adult ticks and their eggs, we assume that there is no infected egg compartment. We also assume that all adult ticks are capable of reproduction, regardless of whether they are susceptible or infected. e eggs mature into the larvae category at a rate of σ T , with a certain percentage dying naturally at the rate μ E .
e susceptible larvae then have the possibility of moving into the infected larvae category at a rate of λ T or remaining in the susceptible larval compartment. Both susceptible and infected larval populations are affected by the natural death rate of μ L . ey also both move into their respective nymph compartments at a rate of τ T , regardless of whether they are susceptible or infected. is leads to the following system of equations for larvae: (A.5) After taking another blood meal, susceptible nymphs move into the infected nymph compartment at the rate λ T . Both susceptible and infectious nymph populations are reduced by the natural death rate of μ N and continue to mature into adults at the rate of c T . e equations for the development rate of nymphs are given as follows: Although it is less common than at other life stages, adult ticks are still able to become infectious. Susceptible adults move into the infected compartment at the rate λ T , after they take a blood meal at that life stage. Both susceptible and infected adult populations are removed by natural death rate of μ A . e equations for both susceptible and infected adult ticks are given as follows: Incorporating all the assumptions and equations above, we have the following system of differential equations: Canadian Journal of Infectious Diseases and Medical Microbiology (A.8)

B. Analysis of the Tick Model (A-5)
Basic Qualitative Properties. Positivity and boundedness of solutions: for the tick model (A-5) to be epidemiologically meaningful, it is essential to prove that all its state variables are nonnegative for all time. In other words, solutions of the model system (A-5) with nonnegative initial data will remain nonnegative for all time t > 0. where λ M (t) � (β M (I L (t) + I N (t) + I A (t))/N M (t)).
Using the integrating factor method, we can rewrite equation (B.9) as follows: Hence, so that 10 Canadian Journal of Infectious Diseases and Medical Microbiology Similarly, it can be shown that F > 0 for all t > 0. Invariant Regions e Lyme disease model (A-5) will be analyzed in a biologically feasible region as follows. Consider the feasible region where where N T � S L + I L + S N + I N + S A + I A .
□ Lemma B.2. e region Ω ⊂ R 9 + is positively invariant for the model (A-5) with nonnegative initial conditions in R 9 + .
Proof. Summing the first two equations of model (A-5), we have the following: Next, the last seven equations of model (A-5) give the following after summing the equations representing the larval, nymph, and adult stages: where μ T � min μ L , μ N , μ A . Since K is the carrying capacity, it follows that S E ≤ K. Hence, equation (B.18) becomes the following: us, where (C. 22) e stability of E 0 can be established by calculating the reproduction number R 0 using the next-generation operator method on system (A-5). Taking I L , I N , I A , and I M as the infected compartments and then using the notation in [20], the Jacobian F and V matrices for new infectious terms and the remaining transfer terms, respectively, are defined as follows: where ρ is the spectral radius. We made a simplifying assumption that en, the reproduction number in (C.24) becomes e expression R 0 is the number of secondary infections in completely susceptible population due to infections from one introduced tick or mouse with Lyme disease. Further, using eorem 2 in [20], the following result is established.

D. Existence and Stability of Disease-Free Periodic Solution
and in case t � nT, S M (t) � S * M � S M ((n + 1)T), so (D.33) is periodic with period T.
us, the solution of the first equation of (D.28) is a solution not only in the time interval [0, T) but also for all t ≥ 0. Hence, the solution of (D.33) in the time interval [0, T) is as follows: (D.34) Following the above argument, the solution of the second equation of (D.28) for all t ≥ 0 is given as follows:  [24], it can be shown that (S M (t), S T (t)) is globally asymptotically stable using the stroboscopic map. Hence, we summarize the results as follows.

Lemma D.4.
e model (D.28) (D.29) has a unique diseasefree periodic solution given as follows:  [25] and define the following matrices: (D.37) Let A be a n × n matrix, Φ A(·) (t) be the fundamental solution matrix of the linear ordinary differential system x ′ � Ax, and ρ(Φ A(·) (w)) be the spectral radius of Φ A(·) (w). Let S M � s m (t) + S M (t), S T � s t (t) + S T (t), I M � i m (t), I T � i t (t). en, system (D.28) can be written as follows: x ′ (t) � Qx(t), t ≠ nT, n ∈ N, can see that the eigenvalues of P 1 e UT are less than 1. Furthermore, if spectral radius ρ(P 2 Φ (F− V) (T)) < 1, then DFPS is stable. is leads to the following theorem. We denote the spectral radius as R p � ρ(P 2 Φ (F− V) (T)). Note that R p does not produce the number of individuals infected by a single infected carrier or infectious individual. In particular, it does not produce the average number of secondary infections [25]. However, it works as a threshold such that the disease persists as R p > 1 [25].

Data Availability
All data used in this study came from published cited sources and are included in the text.

Conflicts of Interest
e authors declare that there are no conflicts of interest.