Citation: Andrea L. Clements, Matthew P. Fraser, Pierre Herckes, Paul A. Solomon. Chemical mass balance source apportionment of fine and PM10 in the Desert Southwest, USA[J]. AIMS Environmental Science, 2016, 3(1): 115-132. doi: 10.3934/environsci.2016.1.115
[1] | Yan-Xia Dang, Zhi-Peng Qiu, Xue-Zhi Li, Maia Martcheva . Global dynamics of a vector-host epidemic model with age of infection. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1159-1186. doi: 10.3934/mbe.2017060 |
[2] | P. van den Driessche, Lin Wang, Xingfu Zou . Modeling diseases with latency and relapse. Mathematical Biosciences and Engineering, 2007, 4(2): 205-219. doi: 10.3934/mbe.2007.4.205 |
[3] | Yicang Zhou, Zhien Ma . Global stability of a class of discrete age-structured SIS models with immigration. Mathematical Biosciences and Engineering, 2009, 6(2): 409-425. doi: 10.3934/mbe.2009.6.409 |
[4] | Simone De Reggi, Francesca Scarabel, Rossana Vermiglio . Approximating reproduction numbers: a general numerical method for age-structured models. Mathematical Biosciences and Engineering, 2024, 21(4): 5360-5393. doi: 10.3934/mbe.2024236 |
[5] | Zongmin Yue, Yuanhua Mu, Kekui Yu . Dynamic analysis of sheep Brucellosis model with environmental infection pathways. Mathematical Biosciences and Engineering, 2023, 20(7): 11688-11712. doi: 10.3934/mbe.2023520 |
[6] | Haonan Zhong, Wendi Wang . Mathematical analysis for COVID-19 resurgence in the contaminated environment. Mathematical Biosciences and Engineering, 2020, 17(6): 6909-6927. doi: 10.3934/mbe.2020357 |
[7] | Curtis L. Wesley, Linda J. S. Allen, Michel Langlais . Models for the spread and persistence of hantavirus infection in rodents with direct and indirect transmission. Mathematical Biosciences and Engineering, 2010, 7(1): 195-211. doi: 10.3934/mbe.2010.7.195 |
[8] | Xia Wang, Yuming Chen . An age-structured vector-borne disease model with horizontal transmission in the host. Mathematical Biosciences and Engineering, 2018, 15(5): 1099-1116. doi: 10.3934/mbe.2018049 |
[9] | Abdelrazig K. Tarboush, Jing Ge, Zhigui Lin . Coexistence of a cross-diffusive West Nile virus model in a heterogenous environment. Mathematical Biosciences and Engineering, 2018, 15(6): 1479-1494. doi: 10.3934/mbe.2018068 |
[10] | Kangbo Bao, Qimin Zhang, Xining Li . A model of HBV infection with intervention strategies: dynamics analysis and numerical simulations. Mathematical Biosciences and Engineering, 2019, 16(4): 2562-2586. doi: 10.3934/mbe.2019129 |
Johne's disease is the clinical manifestation of infection in animals (mostly ruminants) with Mycobacterium avium subspecies paratuberculosis (MAP). It is a chronic inflammatory disease with a long incubation period. It affects the small intestine of ruminants and ultimately leads to the death of the animal. MAP transmits through MAP-contaminated feces, colostrums, milk, water, and soil [6]. Although in utero transmission was also reported [33], the major mode of transmission was reported to be fecal-oral transmission route [6]. Intermittent shedding of MAP into feces occurs at an early stage in MAP-infected animals and continuous higher level of shedding occurs in a portion of the infected animals [20,29]. The outcome of ingestion of MAP-contaminated feces depends on age of the animal and dose of MAP [19,21].
MAP is a slow-growing bacterium and is resistant to extreme temperature. Hence, it can survive outside in the environment for a year or longer [3,32]. But, MAP is a facultative intracellular pathogen and the number of MAP bacilli does not increase in the environment [3].
MAP infection is important to study due to its economic impact on the dairy industry and its possible connection to a human disease, Crohn's disease [26,28,30]. MAP has been estimated to be present in about 68% of US dairy operations [35]. The main problem for the dairy industry is the reduction in milk production and the premature culling of the infected animals. This causes about $250 million annual loss in the productivity of the dairy industry [22].
An important way of controlling MAP infection is through proper management of animals on farms [10]. Although, MAP can be killed by antibiotics, the treatment is not economically feasible to dairy farmers [10]. There are two main ways to diagnosis MAP infection. One is by testing for the bacteria (MAP) and the other is by testing for the antibodies. In the first method of testing, the fecal samples collected from the animals or from the surrounding environment are tested by culturing the bacteria. However, due to the slow growth of MAP these tests could take about two months for the results. Another method of detection of MAP is based on a polymerase chain reaction, which is costly and requires a skilled technician for running the test. In a frequently used test, either blood or milk is tested for the antibodies and the results can be obtained much faster than in MAP detection methods. ELISA (enzyme-linked immunosorbent assays) is an antibody test which is widely used for MAP infection [11]. In 2006, Eda and his colleagues developed a test called EVELISA (ELISA using ethanol extract of MAP) with a much better sensitivity [18,27] than the current ELISA test. Culling the animals with positive test results has been studied to control of MAP infection in dairy farms. In [14], it was shown that with proper farm management, the culling the high-shedding animals can be used to control the disease. Further, it showed that for better results the testing has to be conducted more frequently than once a year to identify and remove high shedders from the farm.
As Johne's disease is a chronic disease with a long incubation period (2-5 years), it is not practical to test different management strategies in actual farms. Mathematical models have therefore been used to predict effectiveness of different control strategies. In a review paper Marcé et al.[15] described eight existing epidemiological models describing the within-herd MAP transmission in dairy farms. Out of those models, four models were discrete-time stochastic models, one was a discrete-time deterministic model and three were continuous-time deterministic models. Some of these models considered horizontal transmission.
But until recently, most models have not considered the MAP transmission through the environment. Humphry et al.[13] considered the ingestion of feces through the environment in their discrete-time stochastic model investigating the dynamics of MAP infection in a suckler-beef herd. In their model, the animals are divided into 4 disease states namely susceptible, subclinical, clinical and resistant. They assumed that the clinically infected animals die within 6 months and further they allow both young and older animals to show clinical signs. The environment was used as a main transmission route in their model. With a 6 months time step in the model, to calculate the probability of infection they considered the bacterial concentration in the environment and the probability of avoiding infection. They assumed that the probability of avoiding infection to be higher for older animals than for younger animals.
In [16], Marcé et al. developed a stochastic model including both the vertical transmission and the horizontal transmission through the ingestion of colostrum, milk or feces in the environment, and they also considered the possibility of calf--to--calf MAP transmission. They studied the dynamics of MAP when no control measures are implemented. For their model they divided the calves into six categories as susceptible, resistant, transiently infectious, latently infected, subclinically infected and clinically infected. They also considered several contaminated environments for each cow category, the environment of the whole farm and the outdoor environment of calves when they graze. They used a model which simulates the population dynamics using specific animal housing arrangements. Further they studied the fadeout and the persistence of the disease starting with small initial conditions for infected animals.
Then, in [17] Marcé et al. studied a closed farm where animals are not purchased from outside to clearly understand the effect of the contact structure of animals within herd on the transmission of MAP. They used the stochastic model developed in [16] and simulated a system without any control in MAP.
In [30], Massaro et al. developed a discrete mathematical model to describe the dynamics of MAP infection in a dairy farm. They have considered the contact structure of the animals in the farm and categorized the animal compartments according to the stage of the disease as well as their age. In their work, they studied the cost effectiveness of EVELISA testing.
Recently, Robins et al. [24] developed an agent-based model to study the dynamics of MAP infection in an dairy farm. They implicitly incorporated the effect of the bacteria in the environment; the direct transmission rate was increased in a location due to the number of high shedders there. Further, they studied the effect of different transmission routes on the prevalence of disease in the system.
Modeling the role of indirect transmission through the environment is an area of current research and modelers are finding it difficult to decide what mechanisms and formats of interaction terms to use. Note that Brehan [2] explained the length of time that a pathogen can survive in the environment would aid in deciding about including such transmission routes and the formats of the interaction terms. Considering the length of time that MAP can survive in the environment, including explicit transmission through MAP in the environment in models seems important, and this is a key component in our paper.
Our goal is to construct a deterministic model to explicitly include compartments for bacteria in different environments and study its effect on transmission of the bacteria among animals in the farm. In studying this effect, we add a novel feature regarding testing. Since testing for MAP usually happens once a year and that is followed by the corresponding culling of infected animals (high shedders) and thenreplacement with susceptible heifers, these culling and replacement actions in the model are implemented as a discrete action at the end of each year. We formulate a system of ordinary differential equations (ODEs) involving progression in this disease and transition in age of the cows on a typical dairy farm. With this system, we are able to show the effects of different transmission routes and management actions on the prevalance. The effects are illustrated with numerical simulations.
The next section gives our model formulation with eleven compartments. Section 3 contains the results of our numerical simulations and discussion of parameters for our model with discrete culling. Our conclusions are presented in the last section.
In our study we develop a continuous deterministic model to describe the dynamics of MAP infection in a dairy farm. We analyze the dynamics of prevalence by taking the environmental persistence of the bacteria into account. In this model we use a system of ODEs to describe the behavior of this infection among dairy cows considering the progression of the disease and the age structure of the cows. We categorize the cows into three age groups as 0 -2 months as calves, 2 -24 months as heifers and above 24 months as adults, similar to previous papers [18,24]. The susceptible cows are in three compartments, Sc, Sh, and Sa, corresponding to calves, heifers and adults. Then, we further divide these age categories according to the stage of the disease with exposed compartments (not able to transmit the disease) Ec, Eh, and Ea. Since calves do not usually progress to shedding of MAP, the low-shedding animals are only in compartments Lh and La; the high shedders are only adults, Ha.
Only the female cows are taken into account and male newborn calves would be removed from the dairy farm. In common dairy farm management, the cows are grouped in different locations according to the age and the milk production. Shortly after the birth the calves are moved into individual hutches and the heifers and adults are kept at separate locations. Usually the adult cows can be in four different locations (housing) as pasture, pregnancy pen, lactation barn and maternity barn. But, in our model for simplicity we consider one location for all the adult cows.
We assume that the calves born to susceptible and exposed adults are not infected with MAP at the birth. But, soon after the birth, some calves born to susceptible or exposed adults can be infected with MAP in the adult's environment (direct contact or from MAP in environment). Also, we assume that some calves born to infected adults can be infected in-utero (vertical transmission) and some will born without the infection. Further, the susceptible calves can be infected by MAP contaminated colostrum and milk. Hence, we subdivided the calves into two categories as susceptible and exposed [9]. We did not consider the calves progressing to the shedding stages. Calves are kept in individual hutches and MAP shed into the environment does not contribute to MAP transmission.
As we assume the heifers can only progress in disease to the level of low-shedding [9,25], the heifers are subdivided into three categories as susceptible, exposed and low-shedding. Since, there are low-shedding animals among heifers, the heifer group housing is considered to have bacteria in that environment (called environment 1 in the model). As a result susceptible heifers not only get infected by direct interaction with low-shedding heifers but also by the bacteria in their housing environment.
In the adult category, infected adults can further progress in disease to the level of high-shedding. Thus, the adults were subdivided into four categories as susceptible, exposed, low-shedding and high-shedding. Also, due to the low-shedding and high-shedding adults, their barns were considered to be contaminated with the bacteria (called environment 2 in the model). Hence, the susceptible adults can be infected due to the direct interaction with low-shedding and high-shedding cows and also due to the interaction with the environment. In a dairy system, testing and culling can be used to control this infection [14]. In common management the cows are tested for MAP infection once or twice a year and if the cows are tested positive as a high shedder they are likely removed from the system.
With all the above assumptions and the total amount of MAP in each of the two environments represented by B1 for heifers and B2 for adults, the flow diagram of our model is given in Figure 1. In Figure 1, the horizontal arrows represent the progression in disease stages and the vertical arrows represent the age transitions. The horizontal dashed arrows represent the sources of bacteria to the heifer environment and the adult environment. Further, the curved dashed arrows (pointing upward) indicate that the bacteria in the environments B1 and B2 interact with the susceptible animals. There is no cross-contamination among calf environment, B1 and B2. Our compartments and parameters are listed in Tables 3 and 4, respectively. We assume that the units of time are per day.
We assume the same birth rates for the susceptible, exposed and infected adults. Also, we assume two different vertical transmission rates for low-shedding and high-shedding animals. Further, different death rates are used for each age category and we assume that the culled animals are replaced by susceptible heifers. Also, we use different transmission rates for the calves due to colostrum and milk from low shedders and high shedders. Further, we use different transmission rates for the susceptible heifers with the low-shedding heifers and for the susceptible adults with the low-shedding adults and high-shedding adults. In the two environments, we assume the same decay rate for the bacteria and two different values for the rates at which the bacteria is added to the corresponding environment from low-shedding cows (heifers and adults) and from high-shedding adults. Further, we assumed different values for the probability of getting infected from the bacteria in the environment for calves, heifers, and adults.
In this study our main goal is to analyze the role of the environments in MAP transmission. Using infrequent testing and diagnostic test sensitivity, later we will add these features as a control into our ODE system.
Since the environmental transmission terms are a key part of our model, we now discuss in detail the choice of this transmission function. The bacteria MAP enters the environment mainly from the feces of the shedding cows.To define the transmission due to the bacteria in the environment we need a function with a threshold level before its transition to its maximum level:
$ f(B) = \frac{1-e^{-\frac{B}{K_1}}}{1+K_2e^{-\frac{B}{K_1}}} \ . $ | (1) |
In (1), B is the total amount of bacteria taken by an animal per unit time (day). Note that K1 changes the slope and shape and K2 changes the threshold value of the function (where it starts to rise sharply). Using a range of effective doses from
In our model, we take B1 and B2 to represent the total amount of MAP in the heifers environment and the adult environment. Now, that B represents the amount of bacteria taken by an animal, we use a scaling factor C to reflect that the animal consumes only a proportion of bacteria in each environment, using
To calculate this proportion of bacteria taken by an animal C, we assume that the bacteria is evenly distributed in the feces and the surrounding environment. Since feces are the main route for MAP to enter and survive in the surrounding environment, the proportion of bacteria taken by an animal per day is approximated by the proportion of the feces taken by an animal per day.
We approximate the amount of new feces entering the environment each day by assuming 800 heifers and 1000 adults in the farm. The low-shedding cows shed less MAP (1-50 MAP/gram of feces)[23] compared to MAP shed by high-shedding cows (100-10,000,000 MAP/gram of feces)[31] We calculate the total amount of new feces in the environment 1 and environment 2 per day as below. On average a cow produces about 29.5 kg of feces per day, which gives, [36]
$
$
|
We assume that a cow will take about 1 mg of feces per day. Hence the proportion of feces taken by an animal per day for the two environments can be calculated as given below. We approximate the proportion of the feces taken by an animal per day, by the proportion of new feces taken by an animal per day:
For environment 1:
$\frac{1}{2.36\times10^{10}} = 4.23\times10^{-11}$ |
For environment 2:
$ \frac{1}{2.95\times10^{10}} = 3.39\times10^{-11}$ |
Hence, for simplicity, we find the average of those 2 values and use
We describe the dynamics of the system with above assumptions using the following system of ODEs:
$ $
|
The birth rates from susceptible, low-shedding and high-shedding adults are represented by b, bLa, and bHa, respectively. The proportions of exposed calves through vertical transmission from low-shedding and high-shedding adults are denoted by
Variable | Defining the variable | Initial value |
Sc | Number of susceptible calves | 130 |
Sh | Number of susceptible heifers | 520 |
Sa | Number of susceptible adults | 650 |
Ec | Number of exposed calves | 70 |
Eh | Number of exposed heifers | 248 |
Ea | Number of exposed adults | 250 |
Lh | Number of low-shedding heifers | 32 |
La | Number of low-shedding adults | 80 |
Ha | Number of high-shedding adults | 20 |
B1 | Amount of bacteria (MAP) in the environment 1(Scaled in 108) | 0.2 |
B2 | Amount of bacteria (MAP) in the environment 2(Scaled in 108) | 590 |
Parameter | Defining the parameter | Parameter value |
b | Birth rate of calves from susceptible and exposed adults | 0.00127 |
bLa | Birth rate of calves from low-shedding adults | 0.00127 |
bHa | Birth rate of calves from high-shedding adults | 0.00127 |
µSc | Death rate of susceptible calves | 0.00028 |
µEc | Death rate of exposed calves | 0.00028 |
µSh | Death rate of susceptible heifers | 0.000063 |
µEh | Death rate of exposed heifers | 0.000063 |
µLh | Death rate of low-shedding heifers | 0.000063 |
µSa | Death rate of susceptible adults | 0.0012 |
µEa | Death rate of exposed adults | 0.0012 |
µLa | Death rate of low-shedding adults | 0.0012 |
µHa | Death rate of high-shedding adults | 0.0012 |
µB1 | Decay rate of bacteria in the heifer environment | 0.0027 |
µB1 | Decay rate of bacteria in the adult environment | 0.0027 |
δ | Culling rate of high-shedding adults | 0.9 |
νL | Probability of getting infected through vertical transmission from low-shedding adults | 0 |
νH | Probability of getting infected through vertical transmission from high-shedding adults | 0.22 |
a1 | Transfer rate from calves to heifers due to age progression | 0.0168 |
a2 | Transfer rate from heifers to adults due to age progression | 0.00151 |
d1 | Transfer rate from exposed heifers to low-shedding heifers | 0.0014 |
d2 | Transfer rate from exposed adults to low-shedding adults | 0.0014 |
d3 | Transfer rate from low-shedding adults to high-shedding adults | 0.00078 |
β1 | Transmission rate for susceptible calves due to the colostrum and milk from low-shedding adults | 0.000021 |
β2 | Transmission rate for susceptible calves due to the colostrum and milk from high-shedding adults | 0.000028 |
γ1 | Transmission rate for susceptible heifers due to direct contact with low-shedding heifers | 0.0000024 |
γ2 | Transmission rate for susceptible adults due to direct contact with low-shedding adults | 0.0000012 |
γ3 | Transmission rate for susceptible adults due to direct contact with high-shedding adults | 0.0000018 |
p | Probability of newborn susceptible calves getting infected by MAP in the adult environment | 0.3 |
r1 | Probability of susceptible heifers getting infected by MAP in the heifer environment | 0.06 |
r2 | Probability of susceptible adults getting infected by MAP in the adult environment | 0.06 |
λ1 | Rate at which the bacteria is added to the heifer environment from the low-shedding heifers | 0.007 |
λ2 | Rate at which the bacteria is added to the adult environment from the low-shedding adults | 0.007 |
λ3 | Rate at which the bacteria is added to the adult environment from the high-shedding adults | 29.5 |
This model has continuous culling of high-shedding adults, represented in the
We use the next generation matrix method [5,4,7,8,34] to analyze the stability of the disease free equilibrium for the model, but an explicit representation for R0 cannot be easily found. We will calculate R0 numerically and present our value after we state the parameter values in the next section.
For the numerical simulations, we have used the initial values and parameter values given in Table 1 and 2. With 2000 cows initially, we chose our initial prevalences (see Table 3), as in [24], which was based on local farms. The birth rates, death rates, vertical transmission rates, and disease progression rates were calculated for this ODE system using the parameter values for the discrete system in [24]. The discrete rates in the agent-based model in [24] were converted to rates for ordinary differential equations (continuous time) with time steps being a day. Since the work in [24] did not have explicit environment compartments, we adjusted the transmission rates (
Susceptible | Exposed | Low-shedding | High-shedding | |
Calves | 65% | 35% | 0% | 0% |
Heifers | 65% | 31% | 4% | 0% |
Adults | 65% | 25% | 8% | 2% |
First to see the dynamics of the disease with no testing or culling we take
The comparison of the dynamics of the total animals in each class with no testing or culling or with annual testing and corresponding culling is shown in Figure 3. With no testing or culling, the number of susceptible cows decreased and at the end of 10 years there were only 26 susceptible calves (20% of initial number), 176 susceptible heifers (34% of initial number) and 5 susceptible adults (0.8% of initial number) in the system. On the other hand, infected animals increased over time and at the end of 10 years there were 301 low-shedding heifers (841% increase), 456 low-shedding adults (470% increase) and 284 high-shedding adults (1320% increase). With the annual testing and culling, the population seemed to approach an equilibrium. Compared to no tasting and culling graphs in Figure 3, with the annual culling, the number of susceptible adults decreased more slowly and about 90 susceptible adults were left in the system at the end of 10 years. Only 5 susceptible adults remained when no culling was applied. Also, with the annual culling, the number of high shedders was much lower than that in no-culling simulation as expected. The increase in susceptible heifers from the annual replacements can also be seen.
The number of total animals in each disease class with no testing or culling and with annual testing and corresponding culling is compared in Figure 4. Over time, with no testing or culling, the number of total susceptible animals decreased significantly as the total number of infected animals (exposed, low-shedding and high-shedding animals) increased in the system. At the end of 10 years there were only 227 susceptible animals left in the system. There were about 723 exposed animals, 757 low-shedding animals and 284 high-shedding animals. Hence, starting with 35% infected animals, 89% of the animals were infected at the end of the 10 year no-testing or culling simulation. But with annual testing and culling, the rate of decrease in total susceptible animals was lower than that in the no-culling simulation. At the end of 10 years, there were about 662 susceptible animals, 664 exposed animals, 632 low-shedding animals and 10 high-shedding animals. Hence, with the annual testing and culling, about 66% of the total animals were infected with MAP at the end of 10 years. This is about 26% decrease in the total infected animals compared to the no-culling simulation.
The comparison of the dynamics of exposed cows from bacteria in each environment when
The effect of culling on the numbers in our compartments is given in the Table 4. Just before culling at the end of 10th year there are about 104 high shedders and after culling at the end of the 10th year there are about 10 high shedders left in the system. Hence, with annual testing and culling, the number of susceptible animals increased while the numbers of shedding animals, especially high-shedding adults decreased.
Compartment | Without culling | With annual testing & culling |
Sc | 26 | 42 |
Ec | 53 | 32 |
Sh | 196 | 530 |
Eh | 349 | 244 |
Lh | 301 | 208 |
Sa | 5 | 90 |
Ea | 321 | 388 |
La | 456 | 424 |
Ha | 284 | 10 |
Figure 6 shows levels of the bacteria in the two environments without and with annual culling. The number of bacteria in environment 2 decreased markedly from about
In Figure 7, the total number of animals getting infected through different transmission routes are given. From the results in Figure 7 with no testing or culling case, most of the cows are infected through environmental bacteria followed by drinking MAP-contaminated milk/colostrum and the direct interactions. In utero transmission contributed the least among the four transmission routes. This order did not change when annual culling is applied.
Further, with the current initial values, if we run the simulation with no testing and culling for 25 years, we get the final values given in column 2 of Table 5, for each infected class and the bacteria in the two environments. Here we can see that the infected classes approach an endemic equilibrium and at the end of 25 years, about 89% of the total animals are infected. Next, if we use these equilibrium values as new initial conditions for our simulation with annual discrete culling (at 90% of high shedders), the new final values for each infected class and the bacteria in the two environments at the end of 10 years are given in Table 5, column 3. Hence, with testing and culling, the system starting from about 89% prevalence decreases to about 69% at the end of 10 years.
Compartment | Equilibrium values after 25 years without culling | Final values with annual testing & culling |
Sc | 25 | 40 |
Ec | 54 | 33 |
Sh | 184 | 498 |
Eh | 349 | 252 |
Lh | 311 | 230 |
Sa | 4 | 67 |
Ea | 308 | 370 |
La | 455 | 439 |
Ha | 296 | 11 |
B1 | 807×108 | 610×108 |
B2 | 3234523×108 | 786526×108 |
Next, using the parameter values in Table 2 and the Disease Free Equilibrium values such that
$R_0 = 6.7$ |
and with continuous culling (
$R_0 = 2.5.$ |
Note that the continuous culling rate of
This work represents the transmission from the environment by using a realistic approximation of bacteria in feces in a farm pasture and barn. The simulations include testing and culling using pulse actions which correspond to actual farm practices.
We calculated the basic reproduction number (R0) for our model using the parameter values. We showed that R0 > 1 and thus the disease free equilibrim is unstable, even with a high level of removal of high-shedding animals from the system. Thus, other management actions besides culling of high shedders are needed to effectively manage this infection. In the future, how to detect and manage low shedders would be important to investigate.
Numerical simulations were run for the model with environmental compartments and annual culling action. Some of the parameters were calculated from the available data and literature and some were determined through knowledge of farm practices. From the numerical simulation results we can conclude that the bacteria in the environment is an important route for MAP transmission. We have analyzed and recorded the effect of discrete (annual) culling on the prevalence in our simulations.
Further modeling and simulation work should be done to see the effect of cleaning of the environment on the prevalence of MAP infection in a farm. Also, we plan to investigate the effect of testing and culling two or three times a year on reducing the prevalence in the system.
[1] | Arizona Department Environmental Quality. Analysis of PM2.5 Exceedences in Pinal County Arizona: Demonstration that PM2.5 Concentration are Driven by Local Sources of PM10 Near the Cowtown Monitor, 2010. Available from: http://www.pinalcountyaz.gov/AirQuality/ Documents/Other%20EPA%20Regulatory%20Actions/CowtownTechnicalWhitePaper.pdf |
[2] | Federal Register, 71 FR 61143 National Ambient Air Quality Standards for PM, 2006. Available from: https://www.gpo.gov/fdsys/granule/FR-2006-10-17/06-8477 |
[3] | Pinal County Air Quality Control District. Pinal County Air Quality Control District Source Apportionment Study, 2003. Available from: http://pinalcountyaz.gov/Departments/AirQuality/ Documents/Monitoring%20Network/pinal_speciation_study.pdf |
[4] | Arizona Department Environmental Quality. Arizona Air Quality Designations Technical Support Document: Boundary Recommendations for the Pinal County 24-hour PM10 Nonattainment Area. Phoenix, AZ: Arizona Department of Environmental Quality, 2010. Available from: http://pinalcountyaz.gov/departments/airquality/documents/other%20epa%20regulatory%20actions/ pinalcountypm10nonattainmtrecommtsd.pdf |
[5] | Clements AL, Fraser MP, Upadhyay N, et al. (2013) Summertime characterization of fine and coarse particulate matter in the desert southwest - Arizona, USA. J Air Waste Manag Assoc 63: 764-772. |
[6] | U.S. Environmental Protection Agency. AirData: Access to Air Pollution Data, 2010. Available from: http://www.epa.gov/airdata/ |
[7] | Clements AL, Fraser MP, Upadhyay N, et al. (2014) Chemical characterization of coarse particulate matter in the desert Southwest - Pinal County Arizona, USA. Atmos Pollut Res 5: 52-61. |
[8] | Upadhyay N, Clements AL, Fraser MP, et al. (2015) Size-differentiated chemical composition of resuspended soil dust from the desert southwest United States. Aerosol Air Qual Res 15: 387-398. |
[9] | Derrick M, Moyers J (1981) Precise and sensitive water-soluble ion extraction method for aerosol samples collected on polytetrafluoroethylene filters. Anal Lett Part A Chem Anal 14: 1637-1652. |
[10] |
Birch ME, Cary RA (1996) Elemental carbon-based method for monitoring occupational exposures to particulate diesel exhaust. Aerosol Sci Technol 25: 221-241. doi: 10.1080/02786829608965393
![]() |
[11] |
Medeiros PM, Simoneit BRT (2007) Analysis of sugars in environmental samples by gas chromatography-mass spectrometry. J Chromatogr A 1141: 271-278. doi: 10.1016/j.chroma.2006.12.017
![]() |
[12] |
Chow JC, Watson JG, Kuhns H, et al. (2004) Source profiles for industrial, mobile, and area sources in the Big Bend Regional Aerosol Visibility and Observational study. Chemosphere 54: 185-208. doi: 10.1016/j.chemosphere.2003.07.004
![]() |
[13] | Schauer JJ, Kleeman MJ, Cass GR, et al. (1999) Measurement of emissions from air pollution sources. 1. C1 through C29 organic compounds from meat charbroiling. Environ Sci Technol 33: 1566-1577. |
[14] | Schauer JJ, Kleeman MJ, Cass GR, et al. (1999) Measurement of emissions from air pollution sources. 2. C-1 through C-30 organic compounds from medium duty diesel trucks. Environ Sci Technol 33: 1578-1587. |
[15] | Schauer JJ, Kleeman MJ, Cass GR, et al. (2002) Measurement of emissions from air pollution sources. 5. C1-C32 organic compounds from gasoline-powered motor vehicles. Environ Sci Technol 36: 1169-1180. |
[16] |
Jia Y, Fraser MP (2011) Characterization of saccharides in size-fractionated ambient particulate matter and aerosol sources: The contribution of Primary Biological Aerosol Particles (PBAPs) and soil to ambient particulate matter. Environ Sci Technol 45: 930-936. doi: 10.1021/es103104e
![]() |
[17] | US Environmental Protection Agency. EPA-CMB8.2 Users Manual, 2004. Available from: https://www3.epa.gov/scram001/models/receptor/EPA-CMB82Manual.pdf |
[18] | Eldred B. Evaluation of the equation for soil composite, 2003. Available from: http://vista.cira.colostate.edu/IMPROVE/Publications/GrayLit/023_SoilEquation/Soil_Eq_Evaluation.pdf |
[19] |
Taylor SR, McLennan SM. (1995) The geochemical evolution of the continental crust. Rev Geophys 33: 241-265. doi: 10.1029/95RG00262
![]() |
[20] | Watson JG, Chow JC (2001) Source characterization of major emission sources in the Imperial and Mexicali Valleys along the US/Mexico border. Sci Total Environ 276: 33-47. |
[21] | Turpin BJ, Lim HJ (2001) Species contributions to PM2.5 mass concentrations: Revisiting common assumptions for estimating organic mass. Aerosol Sci Technol 35: 602-610. |
[22] |
Jia YL, Clements AL, Fraser MP (2010) Saccharide composition in atmospheric particulate matter in the southwest US and estimates of source contributions. J Aerosol Sci 41: 62-73. doi: 10.1016/j.jaerosci.2009.08.005
![]() |
[23] | Zhang KM, Wexler AS (2004) Evolution of Particle Number Distribution near Roadways. Part I: Analysis of Aerosol Dynamics and Its Implications for Engine Emission Measurement. Atmos Environ 38: 6643-6653. |
[24] |
Solomon.PA, Hopke PK, Froines J, et al. (2008) Key Scientific and Policy- and Health-Relevant Findings from EPA’s Particulate Matter Supersites Program and Related Studies: An Integration and Synthesis of Results. J Air Waste Manag Assoc 58: S1-S92. ![]() |
1. | Christina Edholm, Benjamin Levy, Ash Abebe, Theresia Marijani, Scott Le Fevre, Suzanne Lenhart, Abdul-Aziz Yakubu, Farai Nyabadza, 2019, Chapter 5, 978-3-030-22043-3, 109, 10.1007/978-3-030-22044-0_5 | |
2. | Lucie Collineau, Charly Phillips, Brennan Chapman, Agnes Agunos, Carolee Carson, Aamir Fazil, Richard J. Reid-Smith, Ben A. Smith, A within-flock model of Salmonella Heidelberg transmission in broiler chickens, 2020, 174, 01675877, 104823, 10.1016/j.prevetmed.2019.104823 | |
3. | Kaoru Hatate, J. Hunter Rice, Karsten Parker, J. Jayne Wu, Amy Turner, Judith R. Stabel, Shigetoshi Eda, Electrochemical Detection of Serum Antibodies Against Mycobacterium avium Subspecies paratuberculosis, 2021, 8, 2297-1769, 10.3389/fvets.2021.642833 | |
4. | Caroline S. Corbett, Mart C. M. de Jong, Karin Orsel, Jeroen De Buck, Herman W. Barkema, Quantifying transmission of Mycobacterium avium subsp. paratuberculosis among group-housed dairy calves, 2019, 50, 1297-9716, 10.1186/s13567-019-0678-3 | |
5. | IBRAHIM HALIL ASLAN, DAVID BACA-CARRASCO, SUZANNE LENHART, JORGE X. VELASCO-HERNANDEZ, AN AGE STRUCTURE MODEL WITH IMPULSE ACTIONS FOR LEPTOSPIROSIS IN LIVESTOCK CATTLE, 2021, 29, 0218-3390, 75, 10.1142/S0218339021500042 | |
6. | Shi-Yuan Xue, Wei Ma, Meng-Yuan Li, Wei-Kang Meng, Yu-Lin Ding, Bo Yang, Yue-Rong Lv, Rui-Bin Chen, Zhi-Hong Wu, Siqin Tunala, Rong Zhang, Li Zhao, Yong-Hong Liu, The Impact of Mycobacterium avium subsp. paratuberculosis on Intestinal Microbial Community Composition and Diversity in Small-Tail Han Sheep, 2024, 13, 2076-0817, 1118, 10.3390/pathogens13121118 | |
7. | J. Reilly Comper, Karen J. Hand, Zvonimir Poljak, David Kelton, Amy L. Greer, Within-herd mathematical modeling of Mycobacterium avium subspecies paratuberculosis to assess the effectiveness of alternative intervention methods, 2025, 01675877, 106496, 10.1016/j.prevetmed.2025.106496 | |
8. | Rong Zhang, Yue-Rong Lv, Bo Yang, Hao Wang, Jun-Tao Jia, Zhi-Hong Wu, Ming Nie, Lian-Yang Sun, Shi-Yuan Xue, Yu-Lin Ding, Rui-Bin Chen, Siqin Tunala, Li Zhao, Yong-Hong Liu, Prevalence and Genotyping of Mycobacterium avium subsp. paratuberculosis in Sheep from Inner Mongolia, China, 2025, 12, 2306-7381, 326, 10.3390/vetsci12040326 |
Variable | Defining the variable | Initial value |
Sc | Number of susceptible calves | 130 |
Sh | Number of susceptible heifers | 520 |
Sa | Number of susceptible adults | 650 |
Ec | Number of exposed calves | 70 |
Eh | Number of exposed heifers | 248 |
Ea | Number of exposed adults | 250 |
Lh | Number of low-shedding heifers | 32 |
La | Number of low-shedding adults | 80 |
Ha | Number of high-shedding adults | 20 |
B1 | Amount of bacteria (MAP) in the environment 1(Scaled in 108) | 0.2 |
B2 | Amount of bacteria (MAP) in the environment 2(Scaled in 108) | 590 |
Parameter | Defining the parameter | Parameter value |
b | Birth rate of calves from susceptible and exposed adults | 0.00127 |
bLa | Birth rate of calves from low-shedding adults | 0.00127 |
bHa | Birth rate of calves from high-shedding adults | 0.00127 |
µSc | Death rate of susceptible calves | 0.00028 |
µEc | Death rate of exposed calves | 0.00028 |
µSh | Death rate of susceptible heifers | 0.000063 |
µEh | Death rate of exposed heifers | 0.000063 |
µLh | Death rate of low-shedding heifers | 0.000063 |
µSa | Death rate of susceptible adults | 0.0012 |
µEa | Death rate of exposed adults | 0.0012 |
µLa | Death rate of low-shedding adults | 0.0012 |
µHa | Death rate of high-shedding adults | 0.0012 |
µB1 | Decay rate of bacteria in the heifer environment | 0.0027 |
µB1 | Decay rate of bacteria in the adult environment | 0.0027 |
δ | Culling rate of high-shedding adults | 0.9 |
νL | Probability of getting infected through vertical transmission from low-shedding adults | 0 |
νH | Probability of getting infected through vertical transmission from high-shedding adults | 0.22 |
a1 | Transfer rate from calves to heifers due to age progression | 0.0168 |
a2 | Transfer rate from heifers to adults due to age progression | 0.00151 |
d1 | Transfer rate from exposed heifers to low-shedding heifers | 0.0014 |
d2 | Transfer rate from exposed adults to low-shedding adults | 0.0014 |
d3 | Transfer rate from low-shedding adults to high-shedding adults | 0.00078 |
β1 | Transmission rate for susceptible calves due to the colostrum and milk from low-shedding adults | 0.000021 |
β2 | Transmission rate for susceptible calves due to the colostrum and milk from high-shedding adults | 0.000028 |
γ1 | Transmission rate for susceptible heifers due to direct contact with low-shedding heifers | 0.0000024 |
γ2 | Transmission rate for susceptible adults due to direct contact with low-shedding adults | 0.0000012 |
γ3 | Transmission rate for susceptible adults due to direct contact with high-shedding adults | 0.0000018 |
p | Probability of newborn susceptible calves getting infected by MAP in the adult environment | 0.3 |
r1 | Probability of susceptible heifers getting infected by MAP in the heifer environment | 0.06 |
r2 | Probability of susceptible adults getting infected by MAP in the adult environment | 0.06 |
λ1 | Rate at which the bacteria is added to the heifer environment from the low-shedding heifers | 0.007 |
λ2 | Rate at which the bacteria is added to the adult environment from the low-shedding adults | 0.007 |
λ3 | Rate at which the bacteria is added to the adult environment from the high-shedding adults | 29.5 |
Susceptible | Exposed | Low-shedding | High-shedding | |
Calves | 65% | 35% | 0% | 0% |
Heifers | 65% | 31% | 4% | 0% |
Adults | 65% | 25% | 8% | 2% |
Compartment | Without culling | With annual testing & culling |
Sc | 26 | 42 |
Ec | 53 | 32 |
Sh | 196 | 530 |
Eh | 349 | 244 |
Lh | 301 | 208 |
Sa | 5 | 90 |
Ea | 321 | 388 |
La | 456 | 424 |
Ha | 284 | 10 |
Compartment | Equilibrium values after 25 years without culling | Final values with annual testing & culling |
Sc | 25 | 40 |
Ec | 54 | 33 |
Sh | 184 | 498 |
Eh | 349 | 252 |
Lh | 311 | 230 |
Sa | 4 | 67 |
Ea | 308 | 370 |
La | 455 | 439 |
Ha | 296 | 11 |
B1 | 807×108 | 610×108 |
B2 | 3234523×108 | 786526×108 |
Variable | Defining the variable | Initial value |
Sc | Number of susceptible calves | 130 |
Sh | Number of susceptible heifers | 520 |
Sa | Number of susceptible adults | 650 |
Ec | Number of exposed calves | 70 |
Eh | Number of exposed heifers | 248 |
Ea | Number of exposed adults | 250 |
Lh | Number of low-shedding heifers | 32 |
La | Number of low-shedding adults | 80 |
Ha | Number of high-shedding adults | 20 |
B1 | Amount of bacteria (MAP) in the environment 1(Scaled in 108) | 0.2 |
B2 | Amount of bacteria (MAP) in the environment 2(Scaled in 108) | 590 |
Parameter | Defining the parameter | Parameter value |
b | Birth rate of calves from susceptible and exposed adults | 0.00127 |
bLa | Birth rate of calves from low-shedding adults | 0.00127 |
bHa | Birth rate of calves from high-shedding adults | 0.00127 |
µSc | Death rate of susceptible calves | 0.00028 |
µEc | Death rate of exposed calves | 0.00028 |
µSh | Death rate of susceptible heifers | 0.000063 |
µEh | Death rate of exposed heifers | 0.000063 |
µLh | Death rate of low-shedding heifers | 0.000063 |
µSa | Death rate of susceptible adults | 0.0012 |
µEa | Death rate of exposed adults | 0.0012 |
µLa | Death rate of low-shedding adults | 0.0012 |
µHa | Death rate of high-shedding adults | 0.0012 |
µB1 | Decay rate of bacteria in the heifer environment | 0.0027 |
µB1 | Decay rate of bacteria in the adult environment | 0.0027 |
δ | Culling rate of high-shedding adults | 0.9 |
νL | Probability of getting infected through vertical transmission from low-shedding adults | 0 |
νH | Probability of getting infected through vertical transmission from high-shedding adults | 0.22 |
a1 | Transfer rate from calves to heifers due to age progression | 0.0168 |
a2 | Transfer rate from heifers to adults due to age progression | 0.00151 |
d1 | Transfer rate from exposed heifers to low-shedding heifers | 0.0014 |
d2 | Transfer rate from exposed adults to low-shedding adults | 0.0014 |
d3 | Transfer rate from low-shedding adults to high-shedding adults | 0.00078 |
β1 | Transmission rate for susceptible calves due to the colostrum and milk from low-shedding adults | 0.000021 |
β2 | Transmission rate for susceptible calves due to the colostrum and milk from high-shedding adults | 0.000028 |
γ1 | Transmission rate for susceptible heifers due to direct contact with low-shedding heifers | 0.0000024 |
γ2 | Transmission rate for susceptible adults due to direct contact with low-shedding adults | 0.0000012 |
γ3 | Transmission rate for susceptible adults due to direct contact with high-shedding adults | 0.0000018 |
p | Probability of newborn susceptible calves getting infected by MAP in the adult environment | 0.3 |
r1 | Probability of susceptible heifers getting infected by MAP in the heifer environment | 0.06 |
r2 | Probability of susceptible adults getting infected by MAP in the adult environment | 0.06 |
λ1 | Rate at which the bacteria is added to the heifer environment from the low-shedding heifers | 0.007 |
λ2 | Rate at which the bacteria is added to the adult environment from the low-shedding adults | 0.007 |
λ3 | Rate at which the bacteria is added to the adult environment from the high-shedding adults | 29.5 |
Susceptible | Exposed | Low-shedding | High-shedding | |
Calves | 65% | 35% | 0% | 0% |
Heifers | 65% | 31% | 4% | 0% |
Adults | 65% | 25% | 8% | 2% |
Compartment | Without culling | With annual testing & culling |
Sc | 26 | 42 |
Ec | 53 | 32 |
Sh | 196 | 530 |
Eh | 349 | 244 |
Lh | 301 | 208 |
Sa | 5 | 90 |
Ea | 321 | 388 |
La | 456 | 424 |
Ha | 284 | 10 |
Compartment | Equilibrium values after 25 years without culling | Final values with annual testing & culling |
Sc | 25 | 40 |
Ec | 54 | 33 |
Sh | 184 | 498 |
Eh | 349 | 252 |
Lh | 311 | 230 |
Sa | 4 | 67 |
Ea | 308 | 370 |
La | 455 | 439 |
Ha | 296 | 11 |
B1 | 807×108 | 610×108 |
B2 | 3234523×108 | 786526×108 |