Model description.

List of stochastic events, their transitions and rates for a given species e in a given cave c. The total population size is N(e,c) = S(e,c)+E(e,c)+I(e,c)+R(e,c), where S(e,c) is the number of susceptible individuals, E(e,c) is the number of infected but non-infectious individuals, I(e,c) is the number of infectious individuals, and R(e,c) is the number of recovered individuals. The force of infection at time t for the species e in the cave c is given by λ(e,c,t), and α(e,c1,c2) is the rate of exchange of species e from the cave c1 to the cave c2. The rest of the parameters are defined in Table 2.#This event is only defined for migratory species (Miniopterus schreibersii and Myotis capaccinii) during the migratory periods: MSCA→CB (16–30 April), MSCA→CC (16–30 April), MSCB→CA (01–15 September), MSCC→CA (16–30 August), MCCA→CB (16–30 April), MCCA→CC (16–30 April), MCCB→CA (01–15 September), MCCC→CA (16–30 August). For a given simulation, the rates of exchange between caves α(e,c1,c2) and α(e,c2,c1) are such that the same number of individuals migrating from c1 to c2 will later move from c2 to c1 in order to guarantee that population sizes remain constant over time.