Medical and entomological malarial interventions, a comparison and synergy of two control measures using a Ross/Macdonald model variant and openmalaria simulation

Highlights • An adaptation of the classical Ross–Macdonald model for vector disease transmission to incorporate time-dependent medical and entomological control measures.• Modeling both mass drug administration and indoor residual spraying campaigns, the synchronous deployment of both yields a synergy where the impact of a joint intervention exceeds that of isolated campaigns.• Openmalaria simulations, separately run, indicate comparable intervention impacts to the Ross/Macdonald model variant.• The vector reservoir of parasitemia is found to be labile, and this dictates the impacts of the medical and entomological interventions.• A scaling-law level of analysis is performed that estimates the rebound of infections in a community after interventions expire, and not only do higher transmission environments bounce back to prevalent infections faster, communities with stronger interventions are shown to have a slower relapse to parasitemia.


Introduction
Modeling and simulation are eminently appropriate tools for optimizing malaria control strategies [1], helping elucidate their probable mechanisms of action, and optimizing intervention timing and coordination. To compare different methods of suppressing infection in a community with endemic disease, we extend an established Ross/ Macdonald model variant for mosquito-born parasite transmission to simply incorporate time-dependent control interventions. In particular, we focus on mass drug administration (MDA) and indoor residual spraying (IRS), two strongly effective interventions which address, respectively, the human and the mosquito reservoirs of parasitemia. These interventions have direct and complementary effects on the transmission dynamics.
There is context and precedence for this investigation from a long history of both MDA and IRS campaigns, sometimes used in combination. In the post-WWII era, DDT and Chloroquine prompted large scale campaigns with the ensuing Global Malaria Eradication Programme of 1955-1969 [2,3]. Subsequent eradication efforts both succeeded and failed [4,5] and a legacy of this era is acquired resistance [6]. Both Chloroquine and DDT lost efficacy with successive sprayings and the regularly high chemical pressure stemming from frequent treatment with a single pharmaceutical. Later, from 1969 to 1976, the Garki project investigated and recorded many aspects of transmission, immunology, entomology and epidemiology [7] and both IRS and MDA campaigns were carried out as isolated interventions and in tandem. When both IRS and MDA interventions were applied in combination, the regular, high coverage (85%) campaigns were able to achieve a high level of control, though unable to fully interrupt transmission [8]. Such control, though it does not progress directly to elimination, will be relevant for the modeling of these interventions below, as we investigate the mechanism for achieving this pre-elimination level of suppression. Beyond Garki, many MDA [9][10][11] and IRS campaigns [12,13] have been carried out with varying successes, and recently reviewed for Africa [14,15].
We consider the isolated and combined impact of MDA and IRS campaigns in a few different transmission settings. Recent modeling efforts have investigated similar control interventions and combinations of them [16][17][18][19]. Generally, the control measures are carried out as finite programs that are deployed, have a prescribed duration, and promptly end. Post-campaign, in the absence of further interventions or a permanent alteration of the entomological environment, modeled communities eventually rebound to their previous transmission balance where malarial infections are again prevalent. As such, eradication scenarios are not directly considered below but rather the process of large-scale reductions in the host and vector reservoirs is, and this is an important component of any elimination program.
The Ross/Macdonald model with control interventions, below, allows MDA and IRS campaigns to be generalized and understood in an illuminating manner, revealing a significant synergistic effect between them, which, if robust, has potential utility in malaria control policy. The aim of the application of this Ross/Macdonald variant is not so much to emphasize detail or force some precision in its extension, but to look for general trends and to elucidate the dynamics of parasite transfer between the host and vector reservoirs, especially when it is restricted with control interventions. As a result, the modeling is relatively crude and has few parameters. In fact, the IRS intervention below, other than prescribing a coverage and duration, is a singleparameter intervention. This sole parameter is also easily interpreted: it essentially benchmarks the (reduced) reproductive number R I 0 during the intervention. This adapted and extended formalism is relatively easy to use and understand, but correspondingly, has somewhat limited flexibility and muscle. It is best employed, as here, for robust features and trends.
To ensure that the observed synergistic effect and its size are not a quirk of the simple semi-analytical model, we have also run openmalaria simulations. The openmalaria simulations incorporate tremendous detail absent in the simple Ross/Macdonald model, including acquired immunity, demographic information included for transmission intensities and a health system for case management, and so afford a meaningful comparison with the much simpler dynamical model. Interventions are simulated as closely as possible with those described in the text in the semi-analytic Ross/Macdonald model. However this comparison is not perfect. Our intention is not to make an absolutely complete, direct comparison on directly equal footing between the semi-analytic model and the full malarial transmission model of openmalaria, but to see if some simple features are in correspondence, and potentially also find contrast, if it exists.
We also establish a novel scaling result that details a community's collective response, post-intervention, to re-infection. Malaria will invade again when campaigns expire and this timeline of re-establishing parasitemia depends on the transmission intensity and other entomological factors, though it is complicated by many forces wholly absent in a "clean" modeling environment [20]. Estimated from a stability analysis in an accompanying appendix, this rate of re-infection is referenced throughout. The rate enables a general, effective period of intervention efficacy: the waiting period prior to the need for new control measures. To our knowledge, it has not been established in the literature.

Methods: A dynamical variant of the classical model
We begin with a simple population model of local mosquitoes harboring parasites and capable of infecting hosts. Using classical notation, the proportion of this mosquito population is Z, and for simplicity will here consist of a single, short-lived anopheles vector with an average lifespan of − g 1 . This is the sporozoite rate and it grows as more mosquitoes become infectious, acquiring parasites through bloodmeals with infectious hosts, and wanes with their death. Among the very simplest dynamical descriptions of the sporozoite rate is (with others reviewed in [21]) Smith and McKenzie [22], which on the right hand side consists of creation and annihilation terms for these rates. The second is clearly the force of mortality, which is here the only means of diminishing the infectious population. The first creates more infectious mosquitoes through biting events with infectious hosts. If the vector bites hosts periodically at an average rate of a (or the feeding time interval is − a 1 ), and the proportion of infectious humans is X, infectious mosquitoes are created at a rate of acX, where c is a human-to-mosquito transmission efficiency. New infectious mosquitoes are then created at a per-mosquito rate of − acX Z (1 ), essentially the first term in Eq. (1). An extrinsic incubation period for infectiousness n is incorporated here by diminishing the eligible infected vector population so that only −gn exp( ) of newly infected anopheles survive to become infectious, a limitation based on the simple hazard model for mosquito survival. Incorporating the delay in this simple fashion avoids the distinction of infected and infectious populations of mosquitoes, and the need to monitor them separately. Mosquitoes below are thus considered as infectious or not, with the creation rate of the infectious fraction given by − − acX e Z ( ) gn . Mosquitoes are assumed to carry infectiousness, once acquired, to their death.
The evolution of the infectious proportion of humans X in a transmission setting is similarly constructed. First, a host recovers, on average, from a malarial infection at a rate of r, and is only infectious − r 1 days post-initiation. The infectious host population evolves similarly to the above sporozoite rate, and is written [23], The first term creates the infectious populace, again through coupled biting events. The rate that uninfected humans − X (1 )become infected is mabZ, where a is again the rate of biting events, and b is the efficiency coefficient for mosquito-to-human transmission [24]. This term is augmented by the prefactor m, which is the population ratio of mosquitoes to humans (and can be very large in high-transmission environments). This first term, − mabZ X (1 ) is constructed as the creation rate of new, infectious human hosts from the naive host populace and is ignorant of the time lag of parasite development in the hosts. Ignoring this latency, the infectious population of humans is again exactly equal to that of the infected proportion of humans. A summary of the seven classical parameters of transmission is given in Table 1.
These dynamics neglect a large number of complicating influences in both the entomology and epidemiology. First, obviously the spatial component to transmission [25][26][27], and the fact that both mosquito and host populations are moving ones is absent [28][29][30] and is something that will not be discussed here further. Also, host immunity, access to care, and super-infection are wholly neglected, and seem to be culpable, for instance, in one poor performance of Ross/Macdonald models in a high-transmission setting [31]. Simple Ross/Macdonald dynamics also overlook any heterogeneities in host selection, the role of parasite densities in transmission events, in-host dynamics and a more proper inclusion of innoculation latencies. A review of models with these embellishments can be found in reference [32]. Due to these complications, the regime of validity for either one of expressions, Eqs. (1) and (2), is generally not known nor discussed. It is reasonable to suppose, however, that low-to-mid intensity transmission environments, those considered below, might be most amenable to these simple dynamics; infectious transactions may be assumed to be prevalent, and complications such as super-infection play less of a complicating role. Seasonal entomological trends are also a strong dynamical influence and are similarly discarded here in order to focus on the effects of interventions in a very simple setting. Eqs. (1)and (2) describe a dynamical overview of the transmissive elements for infection and provide perhaps the simplest tranmsmission model, for host and vector dynamics. Infectious proportions are boosted with coupled density transmission events and diminished with vector death or the expiry of host infections; these are the only means in which infectiousness is gained or lost.
Analyzing the system Eqs. (1)-(2) for fixed points, there are two, as is well-known [33,34]. The first is the trivial solution, = X 0 and = Z 0, which is an elimination scenario, with neither infectious mosquitoes nor hosts present; there are no parasites transmitted in the host/vector system. The second is a dynamical equilibrium with infectious parasites present in both human hosts and mosquitoes, Two definitions, = γ ac g / and = − P e e gn have been used and the reproductive number R 0 has been introduced, The reproductive number, the number of infected humans created from a single infected host in an otherwise naive population [23,35], is the product of the mismatched ratio of populations m, the mosquito's probability of surviving the extrinsic incubation period = − P gn exp( ), e contributions from the two bites necessary for infection transmission, ac/g and ba/g, and a factor which is the ratio of the amount of time an average mosquito has with an infectious host, g/r. The product of the middle terms is a 2 bc/g, an absolute count of the (lifetime average) successful number of two-fold transmissive events taking place per mosquito. Prefactors diminish or augment this simple count: the mosquito population in m has more mosquitoes (linearly) causing more overall transmission, and their mortality during incubation P e diminishes it. The reproductive number R 0 alone quantifies transmission intensity, a result of the balance of the creation and loss terms of Eqs. (1) and (2).
Rewriting the system above with two transformations, first dividing the infectious populations by their (nontrivial) stable points, = X X X / * and = Z Z Z / *, and second, scaling time by the mosquito lifetime = τ gt. The system of Eqs. (1) and (2), transforms to, with the following redefined parameters: = α r g / is the mosquito lifespan divided by the human infectious period (as mentioned above, this is the ratio of time an average mosquito has with the human infectious period), = β mab g / is the lifetime average number of successful infecting bites per host, and = γ ac g / , as above, is the number of bites in a mosquito's lifetime that infects it. It is clear with Eqs. (5)-(6) that γ mediates human-to-mosquito interactions while β mediates mosquitoto-human ones. The reproductive number is also simply expressed with the scaled parameters, which emphasizes its alternative interpretation which is that it is the ratio of two rates, that of infections invading the community γβP e , divided by the rate they leave, α. The resulting dynamical system in Eqs. (5)-(6) (and thus Eqs. (1) and (2)) is thus fully specified by just four independent parameters: R 0 and γ specify X* via Eq. (3), adding P e specifies Z*, and finally α (β) yields β (α). The complete set is {R 0 , γ, β, P e }, consisting of parameters of the transmission intensity of the setting R 0 , γ and β are scaled interaction parameters for mosquito-to-human mediated transmission and viceversa, and lastly = − P gn exp( ) e is the proportion of mosquitoes that survive the incubation period. The incubation period n, as in Eq. (1) is an obvious alternate specification. Table 1 summarizes the seven classical transmission parameters and the four elemental ones.
The conditions of stability and some other dynamical considerations for the system of Eqs. (5) and (6) are analyzed in Appendix A. In brief, stability about the elimination point, = = X Z 0 is lost when conditions merit R 0 > 1, or when the community has infection rate that outpaces its recovery. All trajectories for reproductive numbers R 0 < 1 attract to the elimination point, = = X Z 0 and all others attract asymptotically to the stable, non-trivial equilibrium point {X*, Z*}, that of Eq. (3). We mention these domains explicitly, though they are widely understood, only because the control interventions considered below will adjust the reproductive number, if temporarily. With good coverage and effective chemoprevention and/or vector control, the reproductive number can fall below unity for the duration of the intervention (high effect size [36]). In this period, dynamics exist on a trajectory that attracts towards elimination. Elimination will however not result from these measures if retraction or intervention expiry causes only a temporary reduction in transmission, and these are the cases explored below. A rebounding R 0 re-introduces transmission and results in a dynamical relaxation to the equilibrium of Eq. (3), re-establishing malaria in the community.

Theory: The addition of interventions
In this section, interventions are modeled as control efforts which change the reproductive number, R 0 , through (say) enhanced mosquito mortality or diminished human infectious periods. These approximated effects perturb the evolution of the infectious populations subject to estimated impacts. For example, periods of intervention activity/duration are approximated, and dynamical trajectories are calculated through the interventions, via Eqs. (5) and (6), but with time periods of suppressed < R R I 0 0 and associated parameters. Upon intervention expiry, the system dynamically relaxes according to the same ambient conditions prior to the intervention. Noted above, and in more detail in Appendix A, globally the system can relax to only one of two fixed points: the trivial equilibrium of elimination, = = X Z * * 0, or the stable transmissive environment given by Eq. (3) where parasites are exchanged freely and there are measurable populations of infectious hosts and mosquitoes at modest R 0 > 1. This relaxation-the slow (or fast) acquisition of malaria in the community, post-intervention-including the characteristic times of regaining equilibrium, is also explored below.
Again, the modeling of MDA and IRS interventions is intentionally simplified. Our interest is to investigate general trends and extract robust, scaling-level effects. For this purpose, three time windows are isolated for the evolution of the infectious proportions X(τ) and Z(τ), which are before, during and after the intervention period. The dynamics evolve according to transmission intensity of each time-period: R R , I 0 0 and R 0 respectively (superscript I denotes an intervention-period value), with the bottlenecked dynamics of the intervention explored in more detail below. Prior to the intervention, the trajectory is static, and that of a community stuck at the stable equilibrium points of Eq. (3). All evolving dynamics are stable with freely exchanged parasites from hosts to mosquitoes and vice-versa. Upon the deployment of either MDA or IRS (or both) campaigns, the second time period, a proportion of the host/vector populations is affected by the intervention, and transmission is immediately reduced. This means, for example, when an MDA is deployed and the coverage is 85% of a local population is treated, ∼ 85% of the infectious pool of humans has their parasite load cleared and their infectious status is exactly voided. This value of the infectious proportion, = X τ τ ( ), 0 with τ 0 a time of intervention deployment, is then used as an initial condition for the subsequent dynamics of the intervention period, with an associated set of parameters appropriate to the intervention: R γ β P { , , , } The time duration of this second period is estimated, which for an MDA is based on the duration of prophylaxis offered by the particular drug administered. Upon expiry of the MDA, the dynamics return to those of the time period prior to the intervention, again with the original parameters, {R 0 , γ, β, P e }. The evolved densities at the expiry of the intervention period with duration τ are used as initial conditions for post-intervention system evolution. This is the starting point for the system's third window, to dynamically relax to one of the two fixed points post-intervention: elimination, or proliferation and the return to an infected populace. In the below, since the pre-elimination parameter values are instantaneously switched back on after = + τ τ τ 0 (for example R 0 is once again suddenly large), elimination can only be achieved with exactly It should be mentioned with these simplifications that decaying intervention effects are not present here. There is no decaying efficacy of an insecticide for example; interventions are taken to be working or wholly ineffective. We choose not to parameterize individual insecticides or medical interventions, but only prescribe an initial impact and effective duration. As will be seen, these generalizations not only enable a simple analysis, they illuminate some of the forces responsible for an intervention's impact, and in particular engender an understanding of the interactions between multiple interventions.
All non-analytic and non-simulation results are full solutions to the coupled nonlinear Eqs. (5) and (6), integrated numerically with a fourth-order Runge-Kutta algorithm as the system itself is not (at least not easily) separably integrable. Parameter values are indicated in the text or figure captions.

Mass drug administration
The effects of an MDA campaign are incorporated through estimates of coverage and duration, as well as adjustments to two internal parameters which account for the altered transmission during the campaign's effective period. First, infections, or here, infectiousness, is squelched at the time of administration and X τ ( ) 0 is set to just 15% of its equilibrium value, a number which corresponds loosely to the campaign's coverage, but also necessarily involves a chosen drug's efficacy, and compliance and/or adherence to its prescriptive treatment. Second, the duration of the MDA also depends on the selection of administered drugs (any parasite resistance is here clearly disregarded), and for simplicity the period of the MDA is taken to be = τ 2, or two mosquito lifetimes (maybe three weeks with ≈ − g 10 1 ). Lastly, systematic effects of the MDA are also ascribed through two internal parameters, r and b, discussed below. Adjustments to these two parameters change the dynamical course of Eq. (2) (or Eq. (5)) through changes in the coefficients, but the prefactors of Eqs. (1) and (6) remain as before. The MDA does not change the entomology-Z(τ 0 ) remains unchanged by the MDA-at the deployment time of the MDA, and the sporozoite rate responds only after this moment. The local infectious proportion of anopheles changes with the MDA only through the coupled dynamics of the infectious populations.
As host infections are cleared for the recipients of the MDA, the average duration of the human infectious period, − r 1 (contained in = α r g / above), decreases as a result. Those infections cleared have shorter durations and, on average, the host infectious period reduces. Transmission mosquito-to-human is also correspondingly squelched, a chemoprophylactic effect for treated individuals, so b (in β above) also diminishes. Through these two modifications, the reproductive number, plummets. To accomplish these reductions, the host healing rate is amplified, r → ξr (or α → ξα in the above) with ξ > 1. This rate amplification diminishes the average duration of a host infection, an obvious effect of administered anti-malarials. The mosquito-to-human transmission for the short period of the campaign is also greatly reduced, a prophylactic effect, and b → b/μ ( = → β mab g β μ / / ) is set to accomplish this. As a consequence, the reproductive number is squelched R 0 → R 0 /μξ for the short duration of the intervention, corresponding to an effect size of 36,37]. This reduction is not entirely independent of the coverage, or vice-versa, in that if coverage is weak, the reduction in R 0 must correspondingly be modest. We focus here on relatively high impact interventions and high coverage in an effort to focus on relatively potent interventions and what may be achieved with them.
Just after the effective time period of the MDA, these effects promptly expire and the parameters return to their pre-intervention values. Succinctly, the MDA campaign has the following settings: an initial impact related to the intervention coverage X(τ 0 ), the duration of the MDA effective period = τ 2, the reduction in the human in- follows from these, again, only during the short time window of the campaign.
With these parameter simplifications, infectious host dynamics during the MDA period < < + − τ τ τ g 2 0 0 1 are determined from, and the intervention values of these prefactors are = α ξα I and = β β μ / I . The sporozoite rate evolves as before, according to Eq. (6) as both γ and P e are unaffected by the intervention: the medical intervention alters only host dynamics directly. The sporozoite rate changes only "downstream," through preempted parasite transmission. It is worth noting that the proportions in Eq. (7) are still normalized by their pre-intervention equilibrium values, X* and Z*.
After the MDA initially causes human infectiousness to plummet with the clearing of the host reservoir, = X τ ( ) 0.15, 0 subsequent host dynamics are determined from the evolution of Eq. (7) with ξ and μ set with two conditions. The first is that these rescalings reduce the reproductive number so ξ and μ must be set so that I 0 as host infectiousness should diminish and attract towards the elimination point = X * 0 during a campaign that is effectively administered with robust anti-malarials. The second condition is to roughly balance the terms on the right hand side of Eq. (7). This selection is made as a sensible estimate: an MDA generally maintains the clearing of the host reservoir during its effective period. That is, if 85% of the population is properly treated with an anti-malarial that has a prophylactic period post-administration, we might expect that host infectiousness X will (at least) not increase, but also hold relatively steady during the period of effectiveness of the drug. If, on the other hand, a particular drug does not have a robust prophylactic effect, μ can be modified (reduced) to accomplish this. Values for ξ and μ, regardless of R 0 , in the cases below are always set with = R 0.5, with R 0 , meaning that here is maintained regardless of the initial, pre-intervention R 0 .
The integrated dynamical trajectories of the MDA campaign are shown in Fig. 1 dynamics for both infectious densities decay towards the elimination stable point (albeit somewhat slowly by construction for X ). The initial impact of an 85% reduction in X is striking, but the sporozoite rate also drops with alacrity, i.e. the vector parasite reservoir depletes quickly, in just a few mosquito lifetimes. The sporozoite rate is clearly very labile with respect to changes in the host infectious density. The MDA campaign targets and mostly clears a very large reservoir of parasites in hosts. The infectious portion of the more ephemeral mosquito population adjusts quickly to reflect these changes. At the end of the MDA campaign, for > + = τ τ τ 12, 0 the intervention has expired and host protection is promptly lost. Mosquitoes continue to bite hosts and individuals previously receiving treatment and prophylaxis are eligible to become infected once again. The stable equilibrium of freely exchanged parasites with X* and Z* in Eq. (3) is regained in time. This decaying approach to X*, Z*, prompted by these uninfected individuals from the MDA becoming re-infected takes much more time in the lowtransmission = R 2.5 0 setting than the high-transmision one, a relatively obvious result, but one which supports the use of MDA in the context of elimination. A stability analysis, described in more detail in Appendix A, estimates the characteristic recovery time(in units of − g 1 ) of the infected host population to be, so that large R 0 clearly bounces back faster, as is seen in Fig. 1. This recovery time, a non-trivial function of transmission parameters, is a decreasing function of γ, as well as β and P e . Boosting the human-tomosquito (or mosquito-to-human) transmission coefficients decreases the time to regain equilibrium. The relation in Eq. (8) is also somewhat unique in that post-campaign, α is conspicuously absent: only those fast dynamics mediated by the entomology and host-transmission, γ, β and P e , are relevant (at first order in R 0 ). As such, while the medical treatment of an MDA campaign clearly affects the recovery time of treated individuals (information contained in α), it is the entomology and parameters of transmission, absent these in-host effects, that causes the community as a whole to lapse once again into a state with prevalent malarial infection. The use of pharmaceuticals for the in-host clearance of infections can establish impressive gains, but the forces for invading parasitemia post-prophylaxis are all first-order dependent on transmission from the vector, despite its depleted reservoir. This foreshadows the needed role for vector control that will be considered in the next sections.
Post-intervention, the rate malarial infections invade the populace may also be written, Here f I (X) is the force of infection, the rate of successful infectious bites received per human. Here, as classically, it is a function of the infected host proportion. A comparable estimate (note expression (8) is approximate) is, The restoration rate, the rate malarial infections return to the community, is essentially the force of infection of the pre-intervention entomological setting. It is the rate of planted, viable infectious bites, that of the entomological setting into which the system relaxes. Before moving on, it is important to remark on Eq. (8) (or Eqs. (9) and (11)) as a scaling relation. It is derived in the linearized approximation of the Eqs. (5) and (6), and as such has relevance in the neighborhood of fixed points such as X*, Z*, while neglecting the nonlinear coupling in these regions. For this matter, the MDA campaign above upon expiry is not a small perturbation from the fixed point of Eq. (3), and the scaling estimate for recovery post-intervention is relatively numerically poor here. A better estimate is found with some further analysis in Appendix A. It is, is the infectious host proportion at intervention expiry. This recovery time is dependent on the suppressed host infectious proportion X f at the end of the intervention: the more parasitemia is suppressed from the campaign, the longer the period of restoration, a sensible result. For this matter, it takes longer for the malarial infections to return in a community where a more successful MDA campaign has taken place. This effect amounts to essentially enhancing the recovery time to a few multiples ≡  Despite these rates and recovery times being valid for a linearized approximation, the response time of Eq. (8) (or (12)) is certainly not without merit given that it illuminates the entomology and transmission properties that are responsible for restoring equilibrium, which is especially true as X(τ) approaches X* ( → X 1). Also, perhaps the most conspicuous attribute of Fig. 1 is the difference in recovery times for the two transmission settings post-intervention. The scaling of this rebound in Eq. (8) with − R 0 1 indicates the relative recoveries of X* differs by an order of magnitude, which matches the recovery trends in the figure extremely well: one is almost exactly ten times longer than the other. One final interesting property explored in Appendix A is that the relapse time of Eq. (8) does not endlessly shorten with R 0 but saturates, ultimately achieving the asymptote ≈ + + τ γ 1 , another consequence of the fast turnover of the vector population.
Of course, in a real setting the force of infection does not instantaneously revert back to pre-intervention levels, as has also been assumed here, but more likely continuously regains momentum as the effect of an intervention deteriorates. This waning influence of an intervention should delay, or slow, reemergence from this estimate. Since immunity is additionally absent in this analysis and it potentially protects some individuals from re-infection, slowing the influx of parasitemia, this estimate is also likely fast from this oversight. Eq. (12) can thus be regarded as a fast estimate, or short threshold, for the characteristic resurgence time.

Indoor residual spraying
An IRS campaign changes the local, ambient vector population and inter alia the reproductive number. Contrasting the MDA campaign, IRS does not change human infectiousness with deployment, so X τ ( ) 0 remains unchanged at this point; IRS affects human infectiousness dynamics only "downstream," as fewer infectious mosquitoes survive to bite humans and transmit infections.
To account for the killing effects of an IRS campaign, a skeletal population model of mosquitoes is required. Killing anopheles with insecticide is an effect contained in m, the ratio of populations, mosquitoes to humans. This ratio m, embedded in the composite parameter β, has thus far been regarded as constant, but for the next few paragraphs alone, we regard it as changing in time to account for the insecticidal killing effect (it will be shown that this time-dependence is spurious). Among the simplest of population models [21], mosquitoes have an emergence/mortality balance, with ϵ a per-host emergence rate and g the force of mortality as above.
In correspondence with the treatment for MDA, all terms are normalized by an unchanging human population of N, i.e. = m M N / , with M the total relevant mosquito population: those proximal, female, and host-seeking. The equilibrium population ratio is = m g * ϵ/ , a balance of birth and death rates. The population dynamics of Eq. (13) has some deficiencies (one emphasized below), but is nonetheless capable of generating insight into the main drivers of transmission reduction. According to this population model, if an IRS campaign diminishes the population to m 0 from the equilibrium value of m* at = τ τ , 0 the mosquito population recovers as, Note that the characteristic time of recovery is (of course) − g , 1 the average mosquito lifetime (again, = τ gt, as above and all time is scaled with the equilibrium mosquito lifecycle). Constant mosquito emergence replaces the lost population in a matter of just a few mosquito lifetimes. This highlights one inherent weakness of the emergence/mortality balance of Eq. (13): the emergence rate is undeterred despite a reduced population of mosquitoes. If perhaps an application of insecticide kills a substantial number of adult mosquitoes, their (likely) fewer offspring do not result. For the same reason, as designed, vector extinction is similarly impossible. Shortcomings such as these will be overlooked in the following in favor of a coarse look at vector control that simply maintains a lower population of mosquitoes when active.
Since IRS kills mosquitoes, with m falling with its deployment, the mosquito mortality rate g must increase. To capture this, for the effective length of time of the IRS campaign, i.e. roughly the viable period of the insecticide, the mortality rate is boosted, g → κg, where κ > 1 here augments the force of mortality and correspondingly diminishes the average lifetime of anopheles. The reproductive number R 0 thus decreases with both fewer mosquitoes m and their shortened lifecycle, diminishing transmission. Augmenting the force of mortality in Eq. (13), with g → κg, the mosquito population re-equilibrates to m*/κ, reducing the equilibration above by a factor of κ. The recovery time, that to regain the new equilibrium m*/κ is short, changing from − g 1 to − κg ( ) , 1 so that not only does the population rebound to a reduced equilibrium, it does so quickly. Because the population recovery time is shorter than an average mosquito lifetime, the second term in Eq. (15) is a fast transient effect, and = m m κ */ alone is sufficient to capture the population reduction for all but the very shortest (irrelevant) time scales. In the following, m no longer contains this time-dependence but during an IRS campaign their ambient population is reduced by the augmented mortality rate m → m/κ.
The rescaling of g during the IRS adjusts the mosquito-to-human transmission prefactor, = β mab g / as well as the fraction of mosquitoes that survive the incubation period, κgn . With the same scale normalization as before, the evolution equations for the IRS transformations are, and note the time = τ gt is still measured in units of − g 1 and not − κg ( ) 1 with this construction. The IRS evolution is modified with a boosted mortality rate g → κg, an altered mosquito-to-human transmission prefactor β, and a diminished survival during the incubation period, P e → P e κ . Fewer survive to pass on parasites with this amplified mortality.
Other than the initial reduction and the campaign duration common to all interventions, it should be noted that the IRS evolution is configured with only a single parameter, κ, which changes these factors. It can be shown the reproductive number is reduced by, a scaling with κ that has been pointed out previously [22]. A last detail is needed for the IRS, which is the connection between a campaign's coverage and how the sporozoite rate adjusts with its deployment, setting Z τ ( ) 0 . This is a rather subtle point, and we might expect that the reduction in the number of ambient mosquitoes m → m/κ be adequate for a description of the insecticide's effect: mosquitoes would be killed indiscriminately and the fraction harboring infections to be unaffected. This is however not the case because the boosted mortality rate g → κg of the campaign modifies the age structure of mosquitoes, and infectious mosquitoes are by requirement older, having survived the latency period. To this end, the hazard model for mosquito survivorship has an age distribution given by and is the fraction affected by a perfect (coverage of unity) IRS campaign. If an IRS campaign affects a fraction of this mosquito population with coverage c 0 , the total reduction at deployment is, the reduction in Z τ ( ) 0 is 71% at deployment for 85% coverage, a substantial reduction.
The effects of this model IRS campaign for two transmission settings are shown in Fig. 2 with 85% coverage, synonymous to the initial reduction in X for the MDA. The intervention duration is much longer for the IRS than the immediate cleansing of an MDA campaign, as insecticide is generally viable for months, and set here to = τ 18 (or ≈ for all transmission environments (which through Eq. (19) sets κ). It is reasonable to expect the sporozoite rate to decline during the effective period of the IRS yielding a net killing effect, diminishing Z and attracting to the stable node of elimination, = Z* 0. Upon retraction, however, prevalent infectious densities again return. Expiry of the IRS here causes the community to be subject to re-infection from the fast dynamics of the recovering vector parasite reservoir. Z quickly bounces back towards the proliferation stable point and host infectiousness follows, scaling with the entomology, the ambient force of infection, as before in Eq. (12).
Another important comparison with an MDA campaign illustrated here is that the host infectiousness X is slow and even comparatively sluggish to respond to changes in the sporozoite rate Z . This is in direct contrast with the MDA campaign above, where fast changes in X result in fast changes in Z . The reverse is not true. Fig. 2 indicates a slow decay in host infectiousness after the mass killing of infectious mosquitoes during the IRS. In fact, the sporozoite rate plummets to roughly a tenth of its equilibrium value for nearly half a year, and host infectiousness decays consistently and comparably slowly. This asymmetry is a result of infections lost only through expiry at a rate of r (see Eq. (2)), a rate directly modified by the MDA but untouched by the IRS. It is only the long duration of the IRS that allows host infections to clear and X to wane. The IRS period has only a trickle of new infections due to a strongly suppressed mosquito population, but with no mechanism or program in place to actually clear existing host infections, waiting out the infectious period is the only means to a reduced presence of parasitemia in the community. The clearing of host infections in an MDA results in a correspondingly fast clearing of the vector reservoir, but the cleansing of the vector reservoir with an IRS does not result in the fast clearing of the host reservoir. Host infections expire only at rate r (contained in α), i.e. slow.
A final note on the model IRS campaign here regards the domain of applicability of the intervention. With the labile nature of the sporozoite rate, low transmission environments R 0 ≲ 10 (not shown in Fig. 2) present a situation where high host infection ≈ X 1, and suppressed Z post-deployment, do not balance the creation and destruction terms of Eq. (1). In these low-transmission environments, anopheles infection can still outpace the relatively (and necessarily) modest κ reduction responsible for the accelerated killing of mosquitoes for a short time. As a consequence, for small pre-intervention R 0 there is a transient period immediately post-IRS deployment where the sporozoite rate grows (responding quickly to the yet-still high X ) but then subsequently falls, as it must attract to the elimation fixed point =

Synchronous IRS and MDA, modeled by the Ross/Macdonald variant
Trajectories for a synchronous deployment of an IRS and MDA are plotted in Fig. 3. Synchronous deployment has the infectious host proportion depleted to just 15% and the sporozoite rate rather comparably reduced from their pre-intervention equilibrium values, in correspondence with the individual campaigns considered above. The interventions are deployed exactly as described in Section 3, with the same durations and reduction parameters. Here, the MDA campaign is finished at = τ 12, and dynamics are then determined by the IRS evolution equations, Eqs. (17) and (18), as it continues beyond the MDA period to = τ 28. Intervention durations are indicated on the figure as before. Comparable to the cases above, at the expiry of the IRS, the parameters instantaneously revert to their environmental settings for transmission prior to the campaigns. Fig. 3 shows this rebound for both host and vector prevalence post-campaign. The lower transmission setting regains equilibrium in the long time limit while the high transmission setting rebounds quickly, again in agreement with the scaling result of Eq. (12). Also plotted are integrable solutions found by exploiting the fact that the sporozoite rate responds very quickly to  changes in the host infectiousness. Appendix A has more details on these trajectories. The most important feature of Fig. 3 is that during the period of the campaign, prevalence is low, nearly in the pre-elimination regime, even for the high-transmission setting. Host infectiousness/infections are first cleared by the MDA to low levels (compare with Fig. 1), and rather than this depletion of the host reservoir being short-lived, this level of prevalence is maintained by the diminished force of infection accomplished by the IRS. In short, the gains of the MDA are sustained by the IRS.
These gains in reduction with combined campaigns are notable because the infection suppression is better than additive. A direct comparison for the high transmission setting, = R 75, 0 is shown in Fig. 4. Host infectiousness trajectories for isolated MDA and IRS campaigns, deployed individually, are compared with that of the combined, synchronous IRS/MDA. It is clear the combined deployment has a synergistic effect, clearing far more infections than either alone.
To measure the robustness of a campaign or campaigns, and to compare effective interventions or intervention sequences, we define their impact, deployed individually or together, as,  and is simply interpreted as the percentage of infections prevented by the campaign times its effective time. As such, a short duration campaign that deeply cleanses the host reservoir, such as a high coverage MDA in a high transmission setting, may have a net impact I comparable to a different, low coverage intervention with a correspondingly long duration. The impact measures the total abilities of the intervention (or a sequence of several) to stem infection.
On this note, the integration limits for Eq. (20) must be large enough to encompass all effects of the intervention(s), especially if their duration is substantial. Provided the time window begins prior to the time of the intervention and is suitably large for long term effects, its duration does not matter. Times τ far in advance or long after the campaign have ≈ X τ X τ ( ) ( ) c 0 and do not contribute; the impact I is constant over a (suitably large) window of integration. There are however a few subtleties for simulations, which are mentioned in Appendix B.
The impacts for an isolated MDA, isolated IRS, and the synchronous MDA+IRS are shown in the inset to Fig. 4, and demonstrate the enhanced suppression of combined, jointly-administered interventions. The impact of the synchronous MDA+IRS is roughly double the impact of an isolated MDA with that of an isolated IRS campaign. This is particularly noteworthy because campaigns that are scheduled and deployed strategically, i.e. jointly administered, vastly outperform two deployed separately, temporally in isolation. This is to say roughly twice the control can be accomplished by jointly deploying the interventions.

Interventions and their impacts with openmalaria simulation
Individual-based simulations are separately run with openmalaria as a comparison for these transmission settings and interventions, and especially to test for the robust synergy in Fig. 4. These simulations include significant detail compared to the semi-analytic model presented above and common ground for comparison is difficult. Several features of openmalaria that are wholly absent in the above model include demographic heterogeneity of transmission, partial immunity of the population, case management in health systems, and variable transmission/infectivity based on a list of factors [40][41][42][43][44][45][46][47]. The  ) and the higher transmission setting bounces back quickly at the expiry of the effective IRS time-span, at = τ 28. A weakly nonlinear solution explored in Appendix A is also drawn. entomology and vector lifecycles are also separably configurable. An appendix below (Appendix B) details several aspects of the simulations run here, and we focus on a comparison that does not necessarily try to strip down openmalaria and compare the base models of transmission and interventions, but rather highlight how reasonably comparable cases yield the same trends. For example, in the below, case management is not entirely absent (though it is minimal), and acquired partial immunity is included in most results, though clearly neither of these are contained in the Ross/Macdonald model above. However, since these additions are integral components of the openmalaria model and to some degree, an elaboration of simpler, underlying transmission models, the intention below is not to force a comparison of exact situations but to see if the more sophisticated simulation approach may also yield a comparable signature of intervention impact and their timings. As will be seen, adhering to general cases and incorporating the most essential features of the interventions enables a reasonable correspondence between the simulation and semi-analytic model. Fig. 5 shows the three situations described above, an isolated MDA, an isolated IRS, and the synchronous deployment of both. Plotted in the tryptic are the Ross/Macdonald semi-analytic theory presented above, and two sets of simulation trajectories for = R 25 0 (corresponding to an annual EIR of = 25 a E bites/host*annum with no seasonal variation, see notes in Appendix B). Here, only the normalized infected host proportion is shown; the sporozoite rate is absent. These trajectories are for hosts with any parasite density in their blood, i.e. those with 0.01 parasites/μL or greater, and not for patent hosts. This standard has been chosen for a congruent comparison with the Ross/Macdonald model, where simply a host either harbors parasites or does not. Intervention coverages of 85% are used and the durations of both interventions are those specified above. No decay is specified for the IRS in the simulation, but it is configured generically as a step function, on during its effective period and then off. Trajectories are normalized by the average prevalence in this community absent interventions. Separate, otherwise equivalent, intervention-free simulations are run to determine these averages, which amount to a determination of X* (and Z* if desired), the prevalence balance points of Eq. (3). Stochastic noise is clearly evident in the eight trajectories plotted. Some of the simulated infected host prevalence trajectories shown in Fig. 5 are for a community absent any immune protection, and others with the openmalaria default for partial immunity [41]. Again, more details of the simulation procedure and the interventions can again be found in Appendix B. A first comparison of all simulations with the Ross/Macdonald variant is that the recovery times, post-intervention are much longer for the simulation. Eq. (12) indicates the community regains half of its infected proportion for τ m ≈ 2 ( = R 25 0 ), agreeing with the Ross/Macdonald trajectories of the plot. In contrast, the simulations here with immunity easily predict three times that. Furthermore, at nearly three months ( ≈ τ 8 post-intervention), the Ross/Macdonald theory indicates the effects of the MDA are essentially gone while it takes perhaps slightly more than a year for the equivalent equilibration in openmalaria. These restoration times, the hangover of these interventions, are substantially different. Malaria invades the community post-intervention with openmalaria at a much slower rate, with and without acquired immunity in the populace. Thus post-MDA, the left panel of the figure, shows one apparent role of acquired immunity in the simulations: it protects the community from the rapid return of widespread infections. A faster return to the proliferation of parasitemia is evident when immune protection is absent. It is, though, not as responsive as the simple Ross/ Macdonald variant model presented above. Which situation better reflects a real (idealized) community suffering from invading malarial infections after a retracted intervention is to us, uncertain. As mentioned above, the Ross/Macdonald rates of Eq. (12) should be regarded as a fast upper limit for recovery due to the abrupt return to the preintervention force of infection (likely an overestimate of transmission at this point in time), and absent immunity. Partial, acquired immunity is instated for all further simulation results below, and the more ponderous return to equilibrium prevalence will be apparent. The IRS deployed in the middle panel of Fig. 5 also presents an interesting comparison. First, IRS simulated by openmalaria is not as effective as IRS modeled by the Ross/Macdonald variant. As mentioned in Section 3.2 above, the rate at which host infections are lost (r in Eq. (2)) essentially determines the diminishing slope, or this effective rate of infection loss in X during the IRS, as only a trickle of new infections amend that rate. Here, in simulation, infections appear to either expire much more slowly, or more than a small number of new infections take place in the IRS period. As discussed in Appendix B, the IRS potency is maximized in this application, with pre-and post-prandial mortality affecting essentially all contacted mosquitoes. The heterogeneity of host selection, an included sophistication in the transmission model built into openmalaria, may also play a role in the weaker overall impact of the IRS.
As a comparison, in both left and center panels, it is intriguing that the rate infections are acquired by a community post-intervention is slow for openmalaria, which is in part due to effects of acquired immunity, but overall the effective force of infection is simply less potent. The combined MDA+IRS campaigns on the right indicates the same trend: the interventions indicate a synergy with joint deployment. Significant gains are again established with synchronous interventions. An inset in Fig. 5  . There is first a comparatively weak IRS, but more significantly, an impact that is better than double for synchronous deployments of a MDA and IRS. The openmalaria simulations demonstrate a synergy that is even larger than that of the Ross/Macdonald model. This is the same simple mechanism at work as seen above: the host reservoir is cleared by the MDA and the vector control prevents it from being refilled. Deployed separately, this coordination is clearly absent and the overall impact is diminished.

Discussion
Both MDA and IRS are potent transmission suppressing interventions, if performed with high population coverage. Given their complementary impacts on the human and vector parasite reservoirs, the uncovered synergy between these interventions is both logical and potentially valuable for malaria control programming. In order to explore the nature and size of this potential synergy, we have extended a simple Ross/Macdonald variant to incorporate these two interventions, MDA and IRS, and have carried out a concise analysis of their overall impacts on population infection.
Each intervention is first modeled in isolation, with intensities of either of the interventions set to the same R I 0 through the duration of the campaigns, an attempt to put them on equal ground considering transmission. The factors of amplification: {ξ, μ, κ} are all introduced as accelerants which manually depress the reproductive number R 0 during the interventions, and are all fairly logical: in an MDA ξ decreases the duration of host infections and μ provides a prophylactic effect, while κ is responsible for killing mosquitoes during an IRS campaign.
The MDA campaign cleanses the host parasites of 85% of the population and having a modest chemoprophylactic effect, protects the recipients for a short duration. After this, individuals are again susceptible and infections are re-established. The dynamics of the Ross/ Macdonald model dictate that these infections restore an equilibrium, or a balance of parasite transactions, post-intervention. The timeline to this re-equilibration, that of relaxing to the stable equilibrium point of Eq. (3), is simply related to the force of infection of the entomological setting the system relapses into. High transmission has a fast relapse to the proliferation of parasitemia, a phenomenon explored in greater detail in Appendix A. This is an approximation, one that is perhaps too fast as an estimate, but one that enables us to discern the forces responsible for restoring the equilibrium, and how this rate scales with R 0 .
The IRS campaign impacts the vector, changing the effective entomological setting of the community. A simple population model for the mosquito population is incorporated which, when diminished, affects the effective R I 0 of the intervention. The IRS, with its much longer duration set by the effective period of the insecticide, inhibits the transmission of parasites by the vector, and leaves hosts with (potentially) no other means of purging infections except to clear on their own slow timescale. Host infections/infectiousness wane very slowly but the populace has protection from new infections through the reduced vector population.
It is sensible that a strategy with combined campaigns would have more impact than isolated campaigns. If the IRS can only offer the protection of reduced biting, while leaving individuals to fight infections off on their own, it is clear that a simultaneously deployed MDA offers the additional, and otherwise absent, therapeutic aspect. The analysis above underlines that a clearing of the host reservoir with an MDA but aided with the added protection of the IRS is a good strategy; the MDA clears infections while IRS effectively protects the cleared individuals from reinfection, adding endurance to these initial gains. The openmalaria simulations confirm this effect and its size are not specific to the Ross/Macdonald model. Lastly, the impact of synergy here is measured in terms of its effect on host-prevalence alone, in terms of those infections (of any parasite load) averted with interventions. An analysis of clinical incidence of disease, instead of the impact on prevalence, would introduce a layer beyond the simple analysis presented here. This synergy, with its extremely simple mechanism, is here revealed through its impact on prevalence alone; a standard metric of malaria in a community.
Simulations were carried out using intentionally simplified settings in order to unmask the fundamental transmission dynamics, and do not necessarily represent a real community or reflect a specific entomological environment. For example, vector biting in our simulations has no seasonal oscillations, and is carried out by a single species. There is also very weak case management. More sophisticated (and real) environments could be simulated but these embellishments might obscure the result (case management, for example, may be represented through a health system parameterization but functions as an additional "background" intervention). The Ross/Macdonald model employs the most basic of transmission dynamics, contains essentially no elaborations, has very few parameters and indicates a strong synergy, together with the full transparency of the forces and effects that enable it. The more complex openmalaria simulation, incorporating a much broader picture of transmission and the many pertinent forces that shape it, confirms this synergy, and in fact predicts a comparable, if not greater impact. Our intent has been to eludicate generalities that are scalinglevel trends and, we hope and expect, are model-independent by nature. We anticipate other transmission models would echo substantively similar results.

Conclusion
We have considered the isolated and combined impacts of two interventions, MDA and IRS, and through two disparate modeling efforts shown the community effects of these interventions is greatly enhanced by their combined application in the setting. This synergy is found to be robust both by the semi-analytic Ross/Macdonald variant and openmalaria. The mechanism of this cooperative impact is illustrated by the semi-analytic model, and also readily understood in terms of basic transmission dynamics.
Additionally, we have also presented a scaling relationship that shows how quickly malaria infections are re-established following an intervention in a simple scenario absent any immunity effects. Not only is resurgence in a strong transmission environment swifter than resurgence in a weak transmission environment, scaling with − R , 0 1 resurgence is shown to be slower for more effective interventions. In the simple scenarios we are modeling, it is plainly beneficial to permanently alter the entomological environment in order to prevent resurgence: this is an obvious result, but one that also highlights the need for integrated vector management.
Eqs. (1)and (2) describe a simple first order, autonomous nonlinear set of equations, and they mimic predator/prey models. Combined with initial conditions, = = X t t X ( ) 0 0 and = = Z t t Z ( ) , 0 0 it constitutes (at least in principle) a solvable system with a unique trajectory. Dynamical stability is assessed from a Jacobian matrix whose eigenvalues categorize linearized trajectory behavior. For the system of Eqs. (5) and (6), it is written, which, evaluated at the trivial, elimination point = = X Z 0 has a larger eigenvalue that determines stability (the other is always negative). The larger is positive for − > γβP α 0, e or R 0 > 1 as is well known. All trajectories for reproductive numbers R 0 < 1 thus attract to the elimination point and all others greater than unity are unstable with respect to elimination. As long as R 0 < 1 is maintained, elimination is the only long-term outcome for system evolution. The analysis below is an extension of work pioneered by Lotka [49] and other, comparable analyses have been performed by other authors [50,51].
The eigenvalues of the Jacobian for the second equilibrium point {X*, Z*} may be written, which shows the discriminant is always positive. Eigenvalues are thus never imaginary and as such there are no cyclic, i.e. periodic solutions common to some predator/prey systems. All trajectories instead exist in attraction basins and terminate at equilibrium points, or potentially run off to infinity. Further inspection indicates that the latter is only possible for a negative determinant, which upon simplification is only the case for R 0 < 1, precisely the region stable for elimination. Thus, the elimination point = X * 0, = Z* 0 and the nontrivial one given in Eq. (3) are both stable fixed points and as such the system asymptotes to either, based on R 0 alone.
The quantities αY and R 0 /Y are nondimensional "bare" rates at which infections are acquired near X*, Z* for humans and mosquitoes. Closer examination of Eq. (A.2) shows that the eigenvalues have a crossover region when R 0 /Y ≈ αY, which corresponds to a curious symmetry point where these parasite acquisition rates are nearly the same in both hosts and mosquitoes. Pre-crossover, expanding for small αY 2 /R 0 , these rates are, and αY . These expansions, the crossover, and the viability of these expansions can all be seen in Fig. A.6. Again, these expressions swap regions of applicability with the crossover, i.e. ↔ + − λ λ , so that the approximation for + λ pre-crossover is an adequate approximation for − λ post-crossover. These eigenvalues are the non-dimensional rates, in units of g, at which the infectious densities regain equilibrium values. To a good approximation (for small α, i.e. momentarily neglecting the small cure rate), host infections are acquired with the elemental solution, The constant term in Eq. (A.4) is dwarfed for modest R 0 and dropped. Eq. (A.6) indicates the sensible result that a larger R 0 has a faster recovery. In contrast, post-crossover with large R 0 , the rate saturates to that of − λ above and asymptotes to → + − λ γ 1 (see Fig. A.6) which is curiously R 0 independent. There is a limit to how fast equilibrium can be recovered-an inertia of the vector reservoir-present even in large R 0 , which is a limiting rate due to the ephemeral lifecycle of the vector. It is noteworthy that pre-crossover this rate grows with R 0 and then post-crossover, saturates. Even for the highest transmission settings, the rate at which equilibirum is recovered never exceeds + γ 1 .

A1. Approximating the time to re-equilibration
In the linearized approximation, post-intervention the system decays to the equilibrium points of Eq. (3) with the general solutions + λ τ exp( ) and − λ τ exp( ), which (somewhat unfortunately) tend to be accurate for small perturbations from the equilibrium. As the heavy suppression of infections with an MDA campaign or the loss of infectious vectors with an IRS are neither small dynamical effects, recovery tracked with this exponential decay is a poor estimate.
One reasonable approximation to estimate recovery trajectories is to exploit the fact that the sporozoite rate Z τ ( ) responds quickly to changes in the host reservoir of X τ ( ), and thus all recovery trajectories in Z τ ( ) in Figs. 1-3 track X τ ( ) very well. Denoting the host infection proportion at intervention expiry as = + X X τ τ ( ) f 0 and simply equating = X τ Z τ ( ) ( ) (which is naturally best when they are reasonably close as in Fig. 3), Eq. (5) simplifies to, species, anopheles gambiae, whose parameters are chosen as previously specified [55] (in an included supplement). The one exception to these parameter assignments is that the extrinsic incubation period is set to = n 10 days in accordance with the average lifetime − g 1 mentioned in the text (which sets = = − − P e e , e gn 1 also, for comparison it is = n 11d in the reference). A correspondence between the EIR E and the reproductive number R 0 is needed to align the simulations with the semi-analytic model. Ref. [48] details different effective R 0 values for openmalaria model simulations, with varying treatment probabilities. For the case of the base model without case management, the reproductive number very nearly corresponds with the annual EIR (total bites per host per annum) a E =365dE . To this end, simulations were run with annual EIR values equal to the reproductive number. A health system is somewhat generically parameterized, as detailed in the openmalaria wiki pages [56]. Defaults for simulations where case management is not a primary interest of the simulation are used, including those for access, adherence, and compliance. Uncomplicated cases of malaria are treated for only 4% of the populace and for 48% of severe episodes. Cure rates are unity for those who access care. This very low level of care allows reasonable correspondence with cases of ref. [48]. Other case fatality rates are as is typical as well as specified sequelae rates.
Finally, interventions are modeled as closely as possible with those described in the text in the semi-analytic model. The effects of the IRS do not slowly decay but are on or off through the use of a step function which activates the intervention from deployment for a duration of 36 timesteps, or 18 − g 1 . The IRS is modeled and deployed as a GVI (generic vector intervention) and maximized in its mosquito killing impact, with preprandial and postprandial killing effects set to 0.99; deterrency is set to zero. The MDA is modeled as clearing the liver for 5 days and blood for 4 timesteps, or 2 − g 1 . Coverage is 85% for both interventions, corresponding to the text. Since about 15 days are required for a full response in the sporozoite rate with vector control in openmalaria, to compare the timing of interventions, the time index of intervention deployment is shifted for an apt comparison. This means, for the Ross/Macdonald variant model, an IRS deployed at τ 0 has the sporozoite rate instantaneously reduced at this moment. An equivalent openmalaria deployment is then begun at τ 0 -(15d/10d), since this accounts for a latent response in the sporozoite rate [53]. We take, for reasons of simplicity, the timing of an IRS campaign as defined to be the very moment that the sporozoite rate plummets. On this note, the MDA timing is also shifted. An MDA deployed within openmalaria at time = t k t Δ 0 has its first effects at timestep + k 1, and thus comparisons are made for a MDA deployed at τ 0 -(5d/10d) with a Ross/Macdonald campaign deployed at τ 0 .
Simulations are performed on a population of = N 10000 individuals with a demographic of that of Ifakara, as detailed in reference [41]. As none of the analytic results in the text above depend on this demography, the distribution is rather irrelevant in this context, but necessary as the force of infection is dependent on an age distribution in openmalaria. The distribution was selected for ease of reproducibility and generality, and it also is representative of the general youthfulness of many sub-Saharan populations. The choice of population size is not selected with a community in mind but rather to simulate a large enough populace to have statistical relevance with small fractional random error. The multiple simulation results plotted demonstrate the stochastic variance in trajectories.
One complication with openmalaria and measuring the impacts of interventions is that in higher transmission settings, high-impact interventions will cause a resurgence in parasitemia that rebounds and can exceed the previous, pre-campaign equilibrium. The impact measure I of Eq. (20) breaks down in these circumstances. As can be seen in Fig. 5, this is not the case here, largely due to the relatively brief interventions. This does not mean that the impact measure I is trouble-free in that post-intervention the long-time limit for integration must be set. With the simulated X τ ( ) equilibrium not being exactly that of pre-intervention, the noisy post-intervention trajectory likely bounces around a nearby equilibrium. To this end, we set the impact integration window to roughly 8 years post-intervention. This means we cut off effects of noise ∼ 8 years (580 timesteps) postintervention, and this is certainly long after these campaigns. Larger windows add greater dispersal about the values of the impact in Fig. 5, as more noise contributes to the sum, but the magnitude of the absolute impact is essentially conserved.
Lastly, cases without immunity are run for 80 years without immune protection and interventions are deployed in the last 5 years of this time window. This is essential to equilibrate the system without the added protection of partial, acquired immunity.