Skip to main content

Sub-national tailoring of malaria interventions in Mainland Tanzania: simulation of the impact of strata-specific intervention combinations using modelling

Abstract

Background

To accelerate progress against malaria in high burden countries, a strategic reorientation of resources at the sub-national level is needed. This paper describes how mathematical modelling was used in mainland Tanzania to support the strategic revision that followed the mid-term review of the 2015–2020 national malaria strategic plan (NMSP) and the epidemiological risk stratification at the council level in 2018.

Methods

Intervention mixes, selected by the National Malaria Control Programme, were simulated for each malaria risk strata per council. Intervention mixes included combinations of insecticide-treated bed nets (ITN), indoor residual spraying, larval source management, and intermittent preventive therapies for school children (IPTsc). Effective case management was either based on estimates from the malaria indicator survey in 2016 or set to a hypothetical target of 85%. A previously calibrated mathematical model in OpenMalaria was used to compare intervention impact predictions for prevalence and incidence between 2016 and 2020, or 2022.

Results

For each malaria risk stratum four to ten intervention mixes were explored. In the low-risk and urban strata, the scenario without a ITN mass campaign in 2019, predicted high increase in prevalence by 2020 and 2022, while in the very-low strata the target prevalence of less than 1% was maintained at low pre-intervention transmission intensity and high case management. In the moderate and high strata, IPTsc in addition to existing vector control was predicted to reduce the incidence by an additional 15% and prevalence by 22%. In the high-risk strata, all interventions together reached a maximum reduction of 76%, with around 70% of that reduction attributable to high case management and ITNs. Overall, the simulated revised NMSP was predicted to achieve a slightly lower prevalence in 2020 compared to the 2015–2020 NMSP (5.3% vs 6.3%).

Conclusion

Modelling supported the choice of intervention per malaria risk strata by providing impact comparisons of various alternative intervention mixes to address specific questions relevant to the country. The use of a council-calibrated model, that reproduces local malaria trends, represents a useful tool for compiling available evidence into a single analytical platform, that complement other evidence, to aid national programmes with decision-making processes.

Background

Since 2000, increased funding towards the universal scale-up of malaria control prevention mainly through insecticide-treated bed nets (ITNs), and treatment with artemisinin-based combination therapy (ACT), substantially reduced the malaria transmission and burden in sub-Saharan Africa [1, 2]. However, in recent years, progress has stalled, and many countries are not on track to achieving national and global targets for 2020 and 2025 as defined in the World Health Organization (WHO) Global Technical Strategy [3]. Mainland Tanzania and nine other countries in Africa contribute together to 66% of the global malaria burden [4]. In 2018 the persistent high burden has instigated the ‘high burden to high impact’ (HBHI) approach for an improved allocation of limited resources in national malaria strategic plans [4]. This approach, encourages national malaria control programmes (NMCPs) in high malaria burden countries to include sub-national stratification of malaria risk with targeted interventions, thereby allowing to intensify control efforts in high transmission areas, while maintaining the gains achieved in low transmission areas. In this context, a rigorous approach for sub-national tailoring of interventions that consist in appropriately selecting intervention mixes for specific risk areas remains a challenge and requires a good understanding of the local context. This can be informed with the use of mathematical modelling to predict the impact that different strategies might have.

Data from epidemiological, clinical, and operational studies along with routine health information systems represent a valuable data source for informing on malaria trends and retrospective intervention impact. However, they can be limited in predicting the impact of interventions over time in specific geographies and for specific combinations of interventions [5]. Mathematical models represent a powerful tool for simulating setting specific malaria transmission dynamics and quantify, with some level of uncertainty, the impact of interventions and their combination in such settings. Mathematical modelling has been used to inform global and national strategies [6,7,8,9], develop target product profiles for new interventions, investigate potential intervention combinations and alternatives (e.g. vector control [10], or vector control and [11]), predict the impact of new interventions such as vaccines [12], and to understand the potential role of surveillance-response [13], among other. The application of mathematical modelling at country level with simulations of suggested intervention mixes for specific geographies has a strong potential to aid decision-making and facilitate a better strategic approach in the selection of interventions.

In mainland Tanzania, the overall malaria prevalence decreased from 18% in 2008 [14] to 7.3% in 2017 [15]. As prevalence declines, the heterogeneity in malaria transmission has increased with 40% of the population living in areas of low or very low malaria risk (prevalence among school children five to sixteen years old PfPR5-16 < 1 and < 5%), 23% living in moderate risk areas (PfPR5-16 5–30%), and 37% in high-risk areas (PfPR5-16 > 30%) that are predominantly located in the north-west and south-east of the country [16]. In the 2015–2020 National Malaria Strategic Plan (NMSP) intervention packages were defined at the regional (admin 1) level with ITNs distributed through mass campaigns or school net programmes (SNP), indoor residual spraying (IRS) in the Lake Zone and larval source management (LSM) in some urban areas with a national target prevalence of less than 1% by 2020 [17]. In 2017, mathematical modelling was applied to assess technical feasibility of reaching the NMSP target and to determine which intervention mixes would maximize impact to meet the target and within the constraint of cost-effectiveness [18]. In alignment with the modelling results, the mid-term strategic review of the NMSP in 2017, concluded that national prevalence targets would not be achieved with the current strategy [19].

In response to the outcomes of the mid-term strategic review, the Tanzanian NMCP developed an innovative approach to stratify the epidemiological risk of the country at the council (admin 2) level as described in Thawer et al. [16]. In the stratification multiple malaria risk indicators from school survey and routine health surveillance data were combined into an overall malaria risk score per council. The resulting risk map grouped the 184 councils of the mainland into four epidemiological strata: very low (PfPR5-16 < 1%), low (PfPR5-16 1–5%), moderate (PfPR5-16 5–30%), and high (PfPR5-16 > 30%), and one operational stratum, urban.

In 2018, recommendations from a consultative meeting [21] with a group of national and international malaria experts, selected a series of intervention mixes to be simulated for each strata including interventions not yet implemented in countries and/or for which WHO guidance was lacking (see complete list in Additional file 1).

The selection of appropriate interventions for each of these strata was performed in partnership with mathematical modellers who provided simulated evidence to compare the impact of various pre-determined intervention mixes.

In this process, three questions were of particular interest: (a) What would be the impact of stopping ITN mass distribution campaigns in the very-low, low, and urban strata? Discontinuing mass campaigns in these areas might allow redistribution of resources for intensified efforts in high-risk areas. (b) What would be the additional benefit of intermittent preventive therapy in school-aged children (IPTsc) when ITNs (and/or IRS) are already deployed in moderate and high-risk strata? School children were found to have high malaria prevalence with high contribution to the transmission in the community [22], and preventive therapy in children might be a valuable addition to the intervention mix deployed in countries [23]. (c) What combination of interventions would be required to substantially reduce malaria in the high risk-stratum? Despite the successes in the past, malaria transmission is still intolerably high in many parts of the country, mainly rural areas [24]. This question aimed to assess the hypothetical maximum impact with the interventions considered for the new strategic plan that could be reachable if funding was not limited. In addition to these questions, it was of interest to compare the predicted impact of interventions defined in the 2015–2020 NMSP with a new candidate strategy for a supplementary strategic plan for the years 2018–2020.

This paper presents the modelling approach conducted in 2018, based on a previous model calibrated to the councils until 2016, used to support answering these specific programmatic questions. The results of the present analysis provided additional information for the NMCP to update their national strategic plan for the period 2018–2020 [20].

Methods

Update of the National Malaria Strategic Plan

Following the mid-term strategic review in 2017, the NMCP decided to update their NMSP and introduce a new malaria risk stratification in mainland Tanzania [16, 20]. The selection of appropriate interventions in each stratum was discussed in consultative meetings held in 2018 [21] with various stakeholders including local researchers, funders, implementers, interregional collaborators, and international partners. A ‘Strength, Weakness, Opportunity and Threat’ (SWOT) analysis was conducted by focus area (integrated malaria vector control, case management, surveillance and monitoring, social behaviour and change communication, and program management), to help in the selection of potential intervention packages for each of the strata [21]. Mathematical modelling was then used to simulate the impact of these suggested intervention mixes per council in order to support the revision of the strategic plan. The revised strategic plan was released in the Supplementary Malaria Midterm Strategic Plan 2018–2020 [20], in the following referred to ‘2018–2020 NMSP’. An overview of the process and methods in shown in Fig. 1 and comparison of the simulated interventions in previous and revised NMSP in Fig. 2.

Fig. 1
figure 1

Overview of modelling support in the process of sub-national tailoring of malaria interventions in Mainland Tanzania. Initial strategic planning meeting [21] that followed the expert consultative and midterm strategic planning meeting (see Runge et al. 2020 [25] for full timeline), that led to the a malaria risk stratification [16] as well as b pre-selection of potential intervention mixes per strata. d The malaria risk stratification and the intervention mixes together with a previously calibrated malaria transmission model [18] were used to generate intervention impact predictions per council. e The objective for the simulations and respective results were 1. the exploration of various intervention mixes at target coverage levels with or without strengthened case management, 2. the comparison of specific intervention scenarios to address some of the questions most relevant to the NMCP at the time of the analysis (questions denoted as ac in results), and 3. the comparison of the simulated current strategy against the potential new strategy. f Modelling results together with NMCP performed analysis and review were discussed. g Formulation of the supplementary malaria midterm strategic plan (2018–2020 NMSP) [20]

Fig. 2
figure 2

Overview of A administrative boundaries, B malaria risk stratification and C, D simulated interventions corresponding to simplified versions of the 2015–2020 and 2018–2020 NMSP per council in Mainland Tanzania. Continuous ITN distribution through schools (ITN-SNP) was operational in fourteen regions (Katavi, Kagera, Mara, Mwanza, Simiyu, Shinyanga, Geita, Lindi, Mtwara, Ruvuma, Morogoro, Tabora, Kigoma, Pwani), and in those areas no ITN-MRC was planned for 2019. ITN-SNP was simulated with 70% coverage in moderate to high transmission strata to reflect school net distributions, or with 40% in the low-malaria risk strata and in urban councils to reflect cumulative coverage of other routine distributions schemes. The final intervention allocation can be obtained from the supplementary malaria midterm strategic plan (2018–2020 NMSP) [20]. ITN, Insecticide Treated Nets; IPTsc, Intermittent preventive therapy in school children; IRS, Indoor Residual Spraying, MRC, Mass Replacement Campaign; SNP, School Net Programme

Model parameterization and calibration

The microsimulation platform OpenMalaria [26] was used to simulate the dynamics of malaria and the impact of interventions in each council. This malaria transmission model represents the dynamics of malaria in humans with an individual-based model and includes a population model simulating the vector dynamics [27,28,29,30,31,32,33,34] as well as the effect size of interventions [10, 29, 35,36,37]. The within-host component of the model had been previously calibrated to historical studies [27, 31]. Transmission seasonality and intensity, historical intervention coverage and vector bionomics are key parameters that characterize a council and were informed by country specific data. Data sources included national-level household surveys and malaria indicator surveys, as well as entomological surveillance reports and intervention distribution information available from the NMCP. Details about the model parameterization and calibration have been described elsewhere [18]. In brief, parameter estimates for historical ITN usage for 2003–2010 were obtained from the Malaria Atlas Project (MAP) [38] and for 2010, 2012 and 2016 from the national-level household survey data [39,40,41]. Model assumptions on treatment seeking behaviour were informed by estimates from the national-level household survey data 2016 combined with a country specific scaling factor for effective treatment given access to health services from Galactionova et al. [42]. Geospatial predictions of Plasmodium falciparum prevalence among children aged 2 to 10 years (PfPR2to10) were derived for each council between 1990 and 2017, based on various community surveys conducted in Tanzania [2]. Using a Bayesian framework, the model was fitted to these predictions to represent historic trends of malaria transmission in each council [18]. Estimated parameters during the fitting process were the pre-intervention annual entomological inoculation rate (EIR) in 2003, and the ITN decay of nets distributed in 2012 for councils that had a mass campaign, or the ITN coverage between 2012 and 2016 for councils with annual school net distributions [43]. Population estimates were obtained from the national census in 2012, assuming linear growth between 2012 and 2016, and constant population afterwards [44]. A constant importation rate of 5 cases per 1000 population per year was used, corresponding to the estimated range of importations with Zanzibar [45]. The model was calibrated for each council and fitted to the model-based prevalence estimates per year between 2003 and 2016. The total simulation time spanned from 2003 to 2022, capturing historical trends until 2016 and projected impact of intervention combinations of interest from 2017 to 2020, with 2 additional years for follow up. In the model, no explicit distinction between urban and rural councils was made. Details on fitting performance per strata are provided in Additional file 2.

Characteristics of simulated interventions 2017–2020

Seven different interventions were simulated either alone or in combinations. The intervention coverage, efficacy and deployments are summarized in Table 1. In the simulations, “effective treatment rates” also referred to as “case management” (CM) coverage were defined as the proportion of uncomplicated malaria cases effectively cured after treatment with an anti-malarial using a parameterization corresponding to ACT, with 98% efficacy in clearing blood stage parasites. As for historical case management, the effective coverage parameter indirectly took into account treatment-seeking and treatment failure rates [42]. The effective treatment coverage for severe malaria cases was fixed at 48% [46].

Table 1 Simulated interventions, deployment, coverage, and efficacy

Coverage of ITNs was defined as the ratio of the population protected by bed nets amongst the total population at risk. The ITN distributions were simulated so that (a) coverage levels (proportion of population protected by a net) increased either at once, representing mass replacement campaigns (ITN-MRC), or (b) coverage were maintained at constant levels, representing continuous distribution mechanisms [43]. In Tanzania, these continuous delivery mechanisms include school net programmes (SNP) and distribution through reproductive and child health clinics (RCH). In the model, the effect of both continuous ITN distribution mechanisms were simulated with a higher and lower cumulative coverage level and both abbreviated using ITN-SNP for simplicity. The ITNs were assumed to have a half-life of 3 years, inferring that only 50% of these would still be used after three years [47] and for conservative purposes, it was assumed that high insecticide resistance to pyrethroids was established throughout the country for all vectors.

The simulations for IRS assumed the insecticide had organophosphate (Actellic 50EC) or carbamate (Bendiocarb) as active ingredients, as these were the used insecticides in recent years [48]. The effect size parameterization was already established previously, with Actellic 50EC parameterized based on an experimental hut study In Côte d’Ivoire [49] and Bendiocarb based on a study in Benin [50] (Briët et al., pers. commun.). Annual rotation between these insecticides was assumed starting with Bendiocarb in September 2017.

Although mass drug administration (MDA) had not been implemented in the country, the intervention was explored in the simulations, assuming an immediate clearance of blood-stage parasites without any lasting prophylactic effect after administration and targeting the whole population.

LSM, specifically with larviciding, was simulated with an effective coverage defined as the proportion of larvae killed (when compared to the number that should have emerged).

For intermittent preventive treatment in school-aged children (IPTsc), targeting children aged 5 to 16 years it was assumed that the drug would lead to immediate parasite clearance and have a prophylactic effect of fourteen days. Seasonal malaria chemoprevention (SMC) was considered for only six eligible councils was not specifically simulated, although roughly approximated with IPTsc in the strategic plan simulations. Intermittent preventative treatment in pregnant women (IPTp) is included in the national strategy but not in the simulations as its impact on transmission is limited [35]. Similarly, behaviour change communication was not explicitly modelled as its direct impact is difficult to quantify. For all simulated interventions it was assumed that the coverage and the effectiveness would be homogeneous within each council.

Simulated intervention mixes and national malaria strategic plans

For the explorative analysis, full factorial combinations of the interventions described above were simulated for all councils. The resulting dataset included model predicted annual outcomes per council per intervention combination for 2016 (baseline) to 2020 (comparison), and 2022. To form a specific intervention mix for the country, the dataset was filtered to match selected intervention mixes to the respective councils. Hence, the simulated national malaria strategic plans were simplified versions, excluding IPTp and IPTi as well as novel insecticides, such as Piperonyl butoxide treated nets (2 councils) or SumiShield. SMC, planned in six councils, was approximated with IPTsc. IRS or LSM did not distinguish between targeted or blanket deployment. Strengthened CM was simulated with a target coverage of 85% effective treatment. The country stratification and simulated interventions are shown in Fig. 2.

Analysis of simulation results

The model outputs were summarized with prevalence rates among children aged between two and ten years and incidence in the whole population, defined as the total number of new cases (i.e. uncomplicated and severe malaria episodes) within a year per 1000 population. Model estimates for council predictions are summarized with median and credible intervals from the posterior distribution from the model calibration. In addition, estimates per strata and at national level were summarized with population-weighted means and 95% confidence intervals (95%CI, shown in squared brackets) of the council median estimates.

Relative reduction between 2016 and 2020 was calculated as (x2016-x2020/ x2016) *100, with x being either prevalence or incidence per council or strata. The impact of discontinuation of ITNs mass campaigns (very low-risk, low-risk, and urban strata) was described by comparing predicted prevalence and percent point change in prevalence between the scenario with strengthened CM only (discontinued vector control interventions starting in 2017) and the scenario with ITNs deployed as a mass campaign (ITN-MRC in 2019) aggregated per strata or per council and pre-intervention EIR within the very low risk strata.

The impact of IPTsc (moderate and high risk-strata) was estimated by calculating the relative reduction in prevalence and incidence compared to no IPTsc, using predictions for 2020. The impact of strengthened CM was assessed in the same way using current CM levels as the counterfactual scenario.

The incremental benefit of adding interventions in the high-risk strata was described using the relative reduction for each intervention mix between 2016 and 2020. The intervention scenarios were: (1) no intervention other than strengthened CM; (2) ITN-MRC; (3) ITN-MRC in combination with IPTsc; (4) ITN-MRC in combination with IRS; (5) ITN-MRC in combination with ITN-SNP; (6) ITN-MRC in combination with ITN-SNP and IPTsc; (7) ITN-MRC in combination with ITN-SNP and IRS; and (8) ITN-MRC in combination with ITN-SNP, IRS and IPTsc. All eight scenarios were simulated assuming strengthened CM.

To compare the impact of the intervention mixes of the two NMSP’s, the difference in predictions was computed for each council (X2015-2020_NMSP—X2018-2020_NMSP, with X being either prevalence or incidence in 2020). The council predictions were aggregated per strata (and nationally) using unweighted and population weighted means for comparison. All analyses were performed using R and RStudio (R Core Team, 2020).

Results

Impact of strata-specific intervention combinations for 2017–2020

Four to eleven scenarios with different intervention mixes and two levels of CM were simulated for an explorative comparison per strata, with fewest scenarios in the moderate risk and most scenarios in the urban stratum (Fig. 3).

Fig. 3
figure 3

Predicted prevalence for intervention mixes per strata for 2020. The error bars and shaded area show the mean and 95% confidence intervals based on heterogeneity among councils. The black and highlighted error bars correspond to current CM level simulations and the grey error bars to strengthened CM. The highlighted error bars show the intervention combinations selected for the 2018–2020 NMSP. The vertical solid line indicates a prevalence of 1% and the dashed line shows the simulated prevalence for 2016

In the very low-risk stratum, the predicted mean prevalence for 2020 ranged from < 0.1% to 1.5% across the simulated interventions and for both CM levels. Maintaining current CM level without additional interventions (counterfactual) was the least effective and implementing MDA was found to be the most effective intervention. Strengthening of case management, the intervention scenario selected for the 2018–2020 NMSP showed a predicted prevalence in 2020 close to the baseline prevalence in 2016.

In the low-risk stratum, the predicted prevalence for 2020 ranged from 0.7% to 4.7% across the simulated interventions. In this stratum, the coverage of ITNs and delivery mechanisms were drivers of impact. Prevalence in 2020 had comparable low levels either when mass campaigns were implemented or if the coverage was at least of 70% when distributed continuously. Amongst the suggested strategies, the ITN mass campaign combined with implementation of LSM had the most impact with predicted prevalence of 0.7% [95% CI 0.5–0.9%] in 2020. On the contrary, a continuous distribution of nets with low coverage (40%) showed the lowest impact, especially with low CM levels (4.6% [95% CI 4.0–5.4%]). All scenarios had on average a lower prevalence compared to the 2016 baseline level (PfPR2to102016: 3.8% [95% CI 2.9–4.9%]), except for ITN continuous distribution with 40% coverage and current CM.

In the moderate-risk stratum, all four simulated scenarios resulted in lower prevalence in 2020 than in 2016 (PfPR2to102016: 15.3% [95% CI 13.2–17.3%]), that was further reduced when assuming strengthened CM, with prevalence ranging from 5.0% to 10.8%. The implementation of both ITN delivery mechanisms simultaneously, as included in the 2018–2020 NMSP, and assuming 70% coverage, was predicted to have similar impact as ITN school campaigns coupled with the implementation of IRS (5.1%, [95% CI 4.6–5.6%] vs 5.0% [95% CI 4.5–5.5%]).

In the high-risk stratum, the predicted prevalence for 2020 ranged between 4.5% and 16.2% across the suggested scenarios, all lower than the baseline prevalence of 25.2% [95% CI 23.6–26.7%]. The simulations showed that the implementation of both ITN delivery mechanisms simultaneously demonstrated a lower prevalence compared to either one alone or in combination with other interventions (7.5% [95% CI 7.0–8.1%] for ITN-MRC + ITN-SNP; 9.2% [95% CI 8.5–9.8%] for ITN-MRC; 10.7% [95% CI 10.1–11.5%] for ITN-SNP). The most impactful intervention mix corresponded, as expected, to the scenario with the most interventions (both ITN distribution mechanism, IPTsc and IRS) and was expected to lower the prevalence to 4.5% [95% CI 4.0–5.0%] in 2020. The intervention mixes simulated for the 2018–2020 NMSP were ITN-SNP, IPTsc and IRS in some councils, the associated predicted prevalence for 2020 was 5.6% [95% CI 5.2–6.2%] with IRS and 9.0% [95% CI 8.4–9.7%] without IRS.

The urban stratum included a broad mix of intervention scenarios, ranging from single interventions such as strengthened CM, LSM, or MDA to combinations of those in addition to either ITN mass or continuous campaigns, depending on the epidemiological strata. The predicted prevalence of the simulated interventions ranged between 0.6% and 7.7%. All intervention combinations were predicted to reduce the prevalence compared to 2016 (PfPR2to102016: 8.1% [95% CI 5.9–10.2%]), except LSM with current CM. In this stratum, four intervention combinations were suggested for the 2018–2020 NMSP depending on the epidemiological risk-strata.

Impact of discontinuation of ITN mass campaigns in very low, low, and urban strata

Without an ITN mass campaign in 2019, an increase in prevalence was predicted in all the three strata with moderate increase by 2020 and considerable increase by 2022. For 2020, the absolute increase was marginal in the very low-risk strata (0.53 [95% CI 0.37–0.70] percent points), slightly higher in the low-risk strata (2.17 [95% CI 1.71–2.62] percent points), and around 5 percent points in the urban strata (5.2 [95% CI 4.3–6.1] percent points) (Fig. 4A). The magnitude in increase differed depending on the pre-intervention transmission intensity. Among the councils in the very low-risk strata, the prevalence remained below the 1% threshold until at least 2022 at pre-intervention EIRs lower than 7 infectious bites per person per annum (ibpa). At pre-intervention EIRs lower than 10 ibpa, the prevalence remained below 1% by 2020, and at pre-intervention transmission levels above 10 ibpa, a prevalence below 1% could not be maintained or reached in 2020 and increased further by 2022 (Fig. 4B).

Fig. 4
figure 4

Predicted impact on prevalence for discontinuation of ITN mass campaigns in the very low, low, and urban strata. A The bar charts show the mean and 95% confidence interval of council prevalence aggregated per strata comparing improved CM without ITN-MRC distribution in 2019 (blue) to improved CM with ITN-MRC distribution in 2019 (orange). B Prevalence per pre-intervention transmission intensity per council in the very low-risk stratum compared for the two scenarios with or without additional ITNs for the years 2016, 2020 and 2022. The points indicate the posterior median, the error bars and the credible intervals resulting from the model calibration. The dashed line highlights a prevalence of one percent

Potential benefit of adding IPTsc in the high and moderate-risk strata combined

IPTsc reduced the predicted prevalence for 2020 on average by 22.3% [95% CI 21.8–22.8%] and the incidence by 15.2% [95% CI 14.7–15.6%], when combined with vector control intervention mixes (ITN-MRC, ITN-SNP, IRS). In comparison to IPTsc, strengthening CM was less effective in reducing incidence in all ages (8.2%, [95% CI 7.8–8.5%], while more effective in reducing prevalence in children between 2 and 10 years (38.9%, [95% CI 38.6–39.4%]) (Fig. 5).

Fig. 5
figure 5

Predicted impact of IPTsc and CM in the moderate and high-risk strata combined. A The impact on incidence in the total population for 2020 and B on PfPR2to10 for 2020. The impact of IPTsc was evaluated when deployed on top of vector control interventions and the impact of strengthened CM was plotted as reference. The boxplots show the distribution among councils in the moderate- and high-risk strata combined. The x-axis presents the different combinations of vector control interventions

Incremental benefit of interventions in the high-risk stratum

The choice of the intervention mixes considered in the high-risk stratum were made to understand the incrementable benefit of each additional interventions. The simulated scenario that corresponded to a maximum number of interventions (both ITN delivery mechanism, strengthened CM, IRS and IPTsc) predicted an average reduction in prevalence by 76.6% [95% CI 72.6–80.7%] between 2016 and 2020. Of that reduction, approximately 70% was attributed to strengthened CM combined with ITN-MRC, which was predicted to reduce the prevalence on average by 52.0%, [95% CI 47.1–57.9%] (Fig. 6A). The other intervention mixes, as shown in Fig. 6A, were predicted to reduce the prevalence by an additional 8% to 25%. The relative reduction in prevalence of additional IRS was comparable to that of additional IPTsc or ITN-SNP (relative reduction in PfPR2to10 63.1% [95%CI 58.4–67.9%] for IRS, 60.9% [95%CI 56.3–65.5%] for IPTsc, 61.1% [95%CI 56.7–65.4%] for ITN-SNP when administered in addition to strengthened CM and ITN-MRC) (Fig. 6A). Strengthened CM and ITN-MRC were predicted to reduce the prevalence by at least 50% in most councils (n = 52 out of 73 districts in high strata including urban high-risk councils). However, in at least 13 councils, even the implementation of all interventions would not be sufficient to reduce prevalence by more than 50% (Fig. 6B). The 13 councils were located in the Southern Zone for which estimated pre-intervention EIR’s were highest, and prevalence predictions started to increase before 2020 (Additional file 2, Fig. A2.5).

Fig. 6
figure 6

Predicted reduction in prevalence in 2020 compared to 2016 per intervention combinations for councils in malaria high-risk strata. A Predicted prevalence reduction between 2016 and 2020 per incremental intervention mix, each point represents a council (n = 73, including four urban councils). The solid line in the boxplot shows the median and the dashed line the mean. B Number of councils per incremental intervention mix grouped by prevalence reduction

Predicted impact of the simulated 2015–2020 NMSP and 2018–2020 NMSP

In the simulation of the 2018–2020 NMSP, the prevalence reductions between 2016 and 2020 were highest in the high and moderate-risk strata (mean PfPR2to10 reduction 63.5%, [95% CI 57.8–69.2%] and 58.9%, [95% CI 55.7–62.1%], respectively), followed by the low- and urban risk strata (57.8%, [95% CI 53.4–62.3%] and 60.1%, [95% CI 52.8–67.3%], respectively) and the lowest in the very low-risk strata (1.5%, [95% CI -12.7–15.8%]) (Additional file 2, Fig A2.4, Table A2.1). Figure 7 shows the ratio of the predicted prevalence and incidence for 2020 of both simulated NMSPs. Compared to the simulated 2015–2020 NMSP scenario, the simulated 2018–2020 NMSP scenario projected lower prevalence and incidence values for most of the councils in the moderate and high-risk strata (with mean ratio of for both strata respectively). In the low-risk strata the simulated 2015–2020 NMSP performed better, while the impact was heterogeneous in the low-risk and in the urban strata. The mean across councils was the same or very similar to the population weighted mean across councils in all strata, except in the urban strata (Fig. 7). In the simulation of the 2018–2020 NMSP, the proportion of the population that would live at a high malaria risk (PfPR2to10 greater than 10%) in 2020 would be around 10% less than in the simulation of the 2015–2020 NMSP (Additional file 2, Table A2.2).

Fig. 7
figure 7

Comparison of predicted impact of the simulated 2015–2020 NMSP and 2018–2020 NMSP per council for 2020. Ratio in predicted prevalence (left panel) and incidence (right panel) for 2020 between both simulated NMSPs. Each horizontal line represents one council. The vertical lines show the mean ratio per strata with solid line for unweighted mean and dot-dashed line for population-weighted mean

Discussion

In this work, mathematical modelling was used to simulate a set of selected intervention mixes between 2017 and 2020 (+ 2022), tailored to the malaria risk strata at council level in mainland Tanzania. The simulated strata-specific intervention mixes were selected to (i) represent the 2015–2020 [17] and the 2018–2020 NMSP [20]; (ii) and to address specific questions relevant to the NMCP. These questions were: (a) what would be the impact of stopping ITN mass distribution campaigns in the very-low, low, and urban strata; (b) what would be the additional benefit of IPTsc when ITNs (and/or IRS) are already deployed in moderate and high-risk strata; and (c) what intervention mix would be required to substantially reduce malaria in the high risk-stratum. The simulations and main analyses were conducted during the strategic planning process in 2018 [25] that accompanied the development of the council-level malaria risk stratification [16]. This work utilizes a previously parameterized model, calibrated to each council in mainland Tanzania, to provide timely model predictions interactively discussed with the NMCP. A table for each of the simulated intervention scenario per strata and the corresponding strategic response was included in the 2018–2020 NMSP [20] (Additional file 3). The process and lessons learned have been published elsewhere [25] and the discussion below focuses on the technical aspects of the modelling analysis.

Overall, both simulated NMSPs had similar impact predictions at national level, that highly varied among and within risk strata. For instance, in the high- and moderate-risk strata the 2018–2020 NMSP was predicted to achieve higher reductions than the 2015–2020 NMSP, which in turn performed slightly better in the very low-risk stratum. Impact predictions for the low-risk and urban strata were mixed and did not show a clear tendency for either strategy. Although, substantial reductions in prevalence for all strategies were predicted at national level between 2017 and 2020, the reductions were not enough to reach the 2015–2020 NMSP target prevalence of less than 1% by 2020. Notably, in both simulated NMSPs the impact on prevalence, with strengthened CM in addition to vector control, was higher than the impact gained through the reallocation of vector control or additional IPTsc for the time frame in consideration. Such a high effect from strengthened CM in combination with vector control was also observed in other countries [8, 51].

In mainland Tanzania, most of the councils in the very low-risk stratum have low receptivity due to unsuitable climate and environmental conditions [52], hence a sub-analysis of the simulations assessed the scenario of discontinuing large scale ITN mass campaigns in the very-low and low risk strata as well as urban strata. The results showed clear increasing trends by 2020 in the low-risk and urban strata for the scenario without an ITN-MRC in 2019 and an even higher increase by 2022, whereas predictions were not homogeneous among councils. In the very low-risk strata, the increase by 2020 was on average low but highly varied depending on the pre-intervention EIR. At low pre-intervention EIR, the prevalence predicted for 2020 was not substantially higher than with a ITN mass campaign in 2019. However, at higher pre-intervention EIR, high CM alone was not enough to keep prevalence low, which aligns with previous modelling studies [53, 54]. While the results could suggest that universal coverage and large-scale deployment of nets might not be required to maintain baseline levels in low-risk settings with persistent low transmission, the results should be interpreted with caution as local contextual factors, especially surveillance and monitoring capacities, were not captured in the model, and as model stochasticity is usually higher in low transmission simulations.

So far, only very few areas in malaria-endemic countries discontinued ITNs. For instance, in Kenya, ITNs that were previously deployed countrywide [55], are deployed in only 50% of counties in the low-risk zone via routine ITN distributions [56]. Zambia is moving from universal ITN to more targeted campaigns [57], whilst in Namibia, ITNs are mostly recommended for personal protection or in transmission foci in malaria risk free or low transmission areas [58], and in Zimbabwe nets are not distributed in the highland areas that are at very low risk of malaria [59]. Analysis on cancellation of ITN mass campaign due to the COVID-19 pandemic in 2020, in low as well as high transmission areas across Sub-Saharan Africa predicted potential substantial increase in cases and deaths, varying across countries [60]. In practice discontinuation of interventions in persistently low transmission areas should only be considered when having a strong surveillance system in place [53], as well as strategies for foci detection and protection of the vulnerable population as included in the 2018–2020 NMSP [20].

In mainland Tanzania as well as other countries, school children were found to be a large reservoir for malaria parasites [61,62,63], with malaria prevalence as high as 76% in some councils [61]. A randomized control trial conducted in Ugandan schoolchildren demonstrated high efficacy of preventive therapies in these groups as well as community-wide effect [64, 65]. Results of a recent systematic review showed beneficial impact of preventive treatment in low and high transmission areas [23] and highlights the importance targeting school children in malaria control. IPTsc, simulated in the moderate and high-risk strata, showed an additional benefit on prevalence and incidence independent of the underlying vector control interventions. These results are lower than estimated effectiveness between 66 to 83% combined across field studies [23]. However, a direct comparison is not possible as the in the present analysis IPTsc was simulated on top of strong vector control interventions, and the effectiveness was calculated compared to counterfactual after three years of deployment. Interestingly, IPTsc showed a higher impact on incidence (all ages) but lower impact on prevalence when compared to the impact of strengthened CM. This could be explained by the differences in target population between the two strategies. Simulated IPTsc targets school children aged 5 to 16-years-old, one of the most vulnerable groups that gets highly infected and symptomatic during the transmission season, whereas CM affects the whole population that might have built disease immunity despite high prevalence. Hence, effectively treating infections might reduce transmission in children and adults but affect the burden, especially for the whole population, less. A simplified IPTsc parameterization was used assuming immediate parasite clearance and prophylactic effect of two weeks with four deployment rounds during the transmission season.

In practice, the post-treatment prophylactic effect was found to range between 14 to 35 days and the deployment schedule differed across trials [23], whereas the optimal deployment of IPTsc across different settings is unclear [66]. In practice, concerns also exist about development of resistance in the parasite against the anti-malarials used for IPTsc [67], and further in-depth analysis with varying anti-malarials and deployment schemes would be useful for more accurate IPTsc impact predictions.

In the third sub-analysis of the presented work, it was assessed what it would take to substantially reduce malaria in the high-risk strata by comparing the incremental impact of additional interventions. In the analysis, a high CM coverage was assumed since strengthening the health care system is a priority independent from vector control and other malaria interventions. The model predictions showed a reduction of around 50% for high coverage levels in CM and ITN-MRCs, and a maximum reduction of around 70% when ITN-SNP, IPTsc and IRS were added. In some councils however, additional interventions of one or more of ITN-SNP, IRS, or IPTisc would be needed to achieve at least a 50% reduction in prevalence. While these findings with greatest impact for CM and ITN align with expectations based on previous published studies [38], the high impact predicted for ITNs is to some extent surprising, as high pyrethroid resistance was assumed for all vectors in all areas. This could be due to the high coverage of at least 80% immediately after the mass campaign simulated in 2019 and the continued protection provided by the physical net barrier [68, 69]. The evidence on reduced net effectiveness at different resistance intensities is limited [70] and varies between modelling studies [29, 71]. In addition, the simulated IRS deployments in 2018 and 2020 used the more short-lived Actellic 50EC [72, 73] instead of Actellic 300CS [74, 75], which likely underestimates predicted impact of IRS. An important additional in-depths analysis in mainland Tanzania would be using vector and location specific resistance parameters and including new insecticides for ITNs, such as pyrethroid piperonyl butoxide [76] and IRS (i.e. SumiShield [77]). Moreover, future reassessment of intervention mixes will need to include improved chemoprevention interventions as well as the malaria vaccine, recently endorsed by the WHO [78].

Utilizing a previously calibrated model to address questions that are relevant to the current situation in the country, enabled delivering timely results to meet the needs of the NMCP during the funding application and strategic planning cycle. In this process, the model results provided an additional layer of information for the selection process of targeted interventions per malaria risk strata. While an impactful approach, the modelling analysis has several limitations that affect the accuracy of the results.

First, the simulated NMSPs are a simplified version of the actual NMSPs [17, 20], as not all interventions were simulated. For instance, IPTp and IPTi were excluded as they were assumed not to have a major impact on malaria transmission in the community [35]. Piperonyl butoxide treated nets were not simulated as, at the time of analysis, planned in only two councils [79]. Similarly, SMC, planned in six councils, was approximated with IPTsc. In the model larviciding did not distinguish between targeted or blanket deployment and due to large uncertainties around feasible larviciding coverages and related effectiveness [80], simulation outputs were excluded from results presented to the NMCP. Strengthened CM was simulated with a target coverage of 85% effective treatment, whereas health system strengthening will likely differ across councils as they depend on baseline performance and strategies to improve CM in respect to the local contexts.

Second, the specific questions addressed would benefit from further in-depth analyses at a more granular level and including more sources of uncertainty, i.e. uncertainty in intervention parameters and model structure, stochasticity and increased population sizes especially in councils with small populations and low transmission intensities. With strata and intervention specific characteristics requiring additional parameters, one large simulation-model becomes inefficient and insufficient to address all programmatic questions. However, to allow for a constant re-evaluation and adjustments in interventions, a fast and flexible modelling approach using a parsimonious model that can be easily updated would be an advantage.

Third, while the results include heterogeneity within strata, the model results are likely less accurate at council level because seasonality, entomology as well as historical intervention parameters were assumed to be homogeneous for councils within a region. Sub-region as well as sub-council heterogeneity are increasingly relevant factors, as risk stratifications are being developed at higher resolution [81] and interventions will be deployed at finer spatial scale. These factors will also be increasingly feasible to account for in future models, as data quantity and quality, especially in routine data, as well as modelling capacities continue to improve. Similarly, urban and rural councils only differed in their transmission intensity and setting specific parameters for seasonality and intervention coverage, whereas socio-demographic factors [82], population density [83], local environment and infrastructure [84], as well as human mobility malaria case importation rates [85, 86] among other factors that are important considerations for intervention implementation were not accounted for.

Finally, the applied model primarily addresses the technical feasibility, not the operational or financial feasibility [87], hence model predictions are likely overestimated and should be interpreted in relative terms. The model predictions obtained in this analysis were designed for comparing the relative impact of the intervention mixes of interest at the time of the analysis rather than making predictions about the future impact of the strategy. This distinction is especially important since long-term temporal effects might bias the intervention impact predictions. Examples for long-term temporal effects include inter-annual variations in climate, that might lead to local malaria epidemics [88,89,90,91], or public health emergencies, such as the COVID-19 pandemic in 2020 [92] that caused service delivery interruptions and delays in several countries and was predicted to potentially substantially increase malaria burden and deaths [60, 93,94,95].

In mainland Tanzania, the IRS campaign in few councils got delayed by around six months due to COVID-19, while the ITN campaign initially planned for 2019, was rescheduled to 2020 due to reasons not related to COVID-19. Recent routine health facility data trends indicate a reduction in outpatient health care seeking and testing whereas test positivity ratio and incidence remained relatively constant. The next round of DHS has not been completed at the time of writing the manuscript and a validation of modelling results was not possible and would have been challenging due to the reasons outlined above. It is however critical to update and recalibrate these models as a dynamic process as new data becomes available to maintain an up-to date country specific model.

The results of this work extend previous work where modelling was applied to assess the technical feasibility of reaching the 2015–2020 NMSP target prevalence of less than 1% by 2020 and to explore alternative intervention allocations at council level that would lead to most impact on prevalence or be most cost-effective in reducing incidence [18]. In contrast to these previous objectives where modelling was used to obtain a new council stratification based on modelled impact, this work used modelling to obtain a comparison for selected interventions for a fixed council stratification based on malaria risk [16, 20]. Both use cases demonstrate the potential value of modelling to support the development of malaria strategic plans, which under the HBHI initiative [4] finds increasing application in other high burden countries [96].

In Tanzania, the shift from almost fully unconstrained mathematical modelling analysis (free combination of interventions and allocation to councils [18]) to a modelling analysis under meaningful constraints set by the NMCP (present analysis) demonstrates two distinct yet related use cases for applying modelling to inform a national malaria control strategy. Although risk stratification should be based on local data, geographical patterns in intervention impact predictions could play a supportive role for sub-national tailoring of interventions to guide intervention allocation for each risk stratum.

The process of subnational intervention impact modeling to support NMCPs will become more reliable and dynamic with the use of increasing high-quality routine data as basis for stratification and modelling. The increased use of mathematical modelling outputs in consultation with NMCPs will result in strengthened strategic and operational planning that will lead towards burden reduction and ultimately elimination.

Conclusion

A modelling approach was presented for predicting the impact of intervention mixes targeted to malaria risk strata in mainland Tanzania as defined by the NMCP during the strategic planning process. By using a previously calibrated model, the model could readily address emerging questions and provided a powerful analytical insight into likely trends of intervention impacts on malaria prevalence and incidence across and within malaria risk strata. The application of modelling for exploring alternative intervention scenarios is likely to increase confidence in the selection of intervention mixes when developing a new national malaria control strategy. Continuous model updates and improvements in the approach will be crucial when scaling up the application of modelling for strategic planning processes in countries.

Availability of data and materials

The simulation datasets generated and analysed during the current study are available in the Zenodo repository (https://doi.org/10.5281/zenodo.5701550), and analysis script available in the GitHub repository, https://github.com/ManuelaRunge/tza_nmsp2018_modeling. All other data are available in the paper and additional information files.

Abbreviations

ACT:

Artemisinin-based combination therapy

CM:

Case management

EIR:

Entomological Inoculation Rate

HBHI:

High Burden to High Impact

ibppa:

Infectious bites per person per annum

iCCM:

Integrated Community Case Management

IPTp:

Intermittent preventive therapy in pregnancy

IPTsc:

Intermittent preventive therapy in school children

IRS:

Indoor residual spraying

ITN:

Insecticide-treated bed nets

LSM:

Larval source management

MoHCDGEC:

Ministry of Health, Community Development, Gender, Elderly and Children

MDA:

Mass drug administration

MRC:

Mass replacement campaigns

NMCP:

National Malaria Control Programme

NMSP:

National Malaria Strategic Plan

PfPR:

Plasmodium falciparum Parasite rate

SNP:

School net programmes

SMC:

Seasonal Malaria Chemoprevention

References

  1. WHO. World Malaria Report 2019. Geneva: World Health Organization; 2019. https://www.who.int/malaria/world_malaria_report_2019/en/.

  2. Snow RW, Sartorius B, Kyalo D, Maina J, Amratia P, Mundia CW, et al. The prevalence of Plasmodium falciparum in sub Saharan Africa since 1900. Nature. 2017;550:515–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. WHO. Global Technical Strategy for Malaria 2016–2030. Geneva, World Health Organization, 2015. https://www.who.int/docs/default-source/documents/global-technical-strategy-for-malaria-2016-2030.pdf.

  4. WHO, RBM Partnership to End Malaria. High burden to high impact: a targeted malaria response. Geneva, World Health Organization, 2019. Report No.: WHO/CDS/GMP/2018.25.

  5. Chubb MC, Jacobsen KH. Mathematical modeling and the epidemiological research process. Eur J Epidemiol. 2010;25:13–9.

    Article  PubMed  Google Scholar 

  6. Walker PGT, Griffin JT, Ferguson NM, Ghani AC. Estimating the most efficient allocation of interventions to achieve reductions in Plasmodium falciparum malaria burden and transmission in Africa: a modelling study. Lancet Glob Health. 2016;4:e474–84.

    Article  PubMed  Google Scholar 

  7. Scott N, Hussain SA, Martin-Hughes R, Fowkes FJI, Kerr CC, Pearson R, et al. Maximizing the impact of malaria funding through allocative efficiency: using the right interventions in the right locations. Malar J. 2017;16:368.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Hamilton M, Mahiane G, Werst E, Sanders R, Briët O, Smith T, et al. Spectrum-Malaria: a user-friendly projection tool for health impact assessment and strategic planning by malaria control programmes in sub-Saharan Africa. Malar J. 2017;16:68.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Otieno G, Koske JK, Mutiso JM. Cost effectiveness analysis of optimal malaria control strategies in Kenya. Mathematics. 2016;4:14.

    Article  Google Scholar 

  10. Chitnis N, Schapira A, Smith T, Steketee R. Comparing the effectiveness of malaria vector-control interventions through a mathematical model. Am J Trop Med Hyg. 2010;83:230–40.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Selvaraj P, Wenger EA, Gerardin J. Seasonality and heterogeneity of malaria transmission determine success of interventions in high-endemic settings: a modeling study. BMC Infect Dis. 2018;18:413.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Penny MA, Verity R, Bever CA, Sauboin C, Galactionova K, Flasche S, et al. Public health impact and cost-effectiveness of the RTS, S/AS01 malaria vaccine: a systematic comparison of predictions from four mathematical models. Lancet. 2016;387:367–75.

    Article  PubMed  PubMed Central  Google Scholar 

  13. RBM. Mathematical Modelling to Support Malaria Control and Elimination. Roll Back Malaria Partnership, UNDP; 2010. http://www.rollbackmalaria.org/microsites/wmd2014/report5.html.

  14. Tanzania Commission for AIDS, Zanzibar AIDS Commission, National Bureau of Statistics, MEASURE DHS, Macro International Inc. Tanzania HIV/AIDS and Malaria Indicator Survey 2007–08. Dar es Salaam, Tanzania, 2008.

  15. Ministry of Health, Community Development, Gender, Elderly and Children; Ministry of Health Zanzibar; National Bureau of Statistics; Office of Chief Government Statistician; The DHS Program. Tanzania Malaria Indicator Survey (TMIS) Malaria Atlas 2017. Dar es Salaam, Tanzania and Rockville, USA. 2018. https://dhsprogram.com/pubs/pdf/ATR20/ATR20.pdf.

  16. Thawer SG, Chacky F, Runge M, Reaves E, Mandike R, Lazaro S, et al. Sub-national stratification of malaria risk in mainland Tanzania: a simplified assembly of survey and routine data. Malar J. 2020;19:177.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Ministry of Health, Community Development, Gender, Elderly and Children; National Malaria Control Programme (NMCP). National malaria strategic plan 2014–2020. Dar es Salaam, 2014. https://www.out.ac.tz/wp-content/uploads/2019/10/Malaria-Strategic-Plan-2015-2020-1.pdf.

  18. Runge M, Snow RW, Molteni F, Thawer S, Mohamed A, Mandike R, et al. Simulating the council-specific impact of anti-malaria interventions: a tool to support malaria strategic planning in Tanzania. PLoS ONE. 2020;15:e0228469.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ministry of Health, Community Development, Gender, Elderly and Children, National Malaria Control Programme, World Health Organization (WHO Tanzania). Report of the mid-term review of the National Malaria Control Strategic Plan 2015–2020. Dar es Salaam, Tanzania: 2017.

  20. Ministry of Health Community Development Gender Elderly & Children; National Malaria Control Programme. Supplementary malaria midterm strategic plan 2018–2020. Dar es Salaam, 2019.

  21. Ministry of Health, Community Development, Gender, Elderly and Children; National Malaria Control Programme. Consultative malaria expert meeting report 2018. Dar es Salaam, Tanzania: 2018.

  22. Nankabirwa J, Brooker SJ, Clarke SE, Fernando D, Gitonga CW, Schellenberg D, et al. Malaria in school-age children in Africa: an increasingly important challenge. Trop Med Int Health. 2014;19:1294–309.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Cohee LM, Opondo C, Clarke SE, Halliday KE, Cano J, Shipper AG, et al. Preventive malaria treatment among school-aged children in sub-Saharan Africa: a systematic review and meta-analyses. Lancet Glob Health. 2020;8:e1499–511.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Ministry of Health, Community Development, Gender, Elderly and Children; Zanzibar Ministry of Health; Tanzania National Bureau of Statistics; Office of Chief Government Statistician; ICF International. Tanzania Malaria Indicator Survey 2017. Dar es Salaam, Tanzania, and Rockville, USA; 2017. https://dhsprogram.com/pubs/pdf/MIS31/MIS31.pdf.

  25. Runge M, Molteni F, Mandike R, Snow RW, Lengeler C, Mohamed A, et al. Applied mathematical modelling to inform national malaria policies, strategies and operations in Tanzania. Malar J. 2020;19:101.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Swiss TPH. OpenMalaria. A simulator of malaria epidemiology and control. Swiss Tropical and Public Health Institute (Swiss TPH); 2020. https://github.com/SwissTPH/openmalaria. Accessed 31 Mar 2020.

  27. Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, et al. Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008;135:1507–16.

    Article  CAS  PubMed  Google Scholar 

  28. Chitnis N, Hardy D, Smith T. A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012;74:1098–124.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Briët OJ, Penny MA, Hardy D, Awolola TS, Van Bortel W, Corbel V, et al. Effects of pyrethroid resistance on the cost effectiveness of a mass distribution of long-lasting insecticidal nets: a modelling study. Malar J. 2013;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Stuckey EM, Smith TA, Chitnis N. Estimating malaria transmission through mathematical models. Trends Parasitol. 2013;29:477–82.

    Article  PubMed  Google Scholar 

  31. Chitnis N, Hyman JM, Cushing JM. Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull Math Biol. 2008;70:1272.

    Article  PubMed  Google Scholar 

  32. Stuckey EM, Stevenson JC, Cooke MK, Owaga C, Marube E, Oando G, et al. Simulation of malaria epidemiology and control in the highlands of western Kenya. Malar J. 2012;11:357.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Penny MA, Maire N, Bever CA, Pemberton-Ross P, Briët OJT, Smith DL, et al. Distribution of malaria exposure in endemic countries in Africa considering country levels of effective treatment. Malar J. 2015;14:384.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, et al. Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012;9:e1001157.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Ross A, Penny M, Maire N, Studer A, Carneiro I, Schellenberg D, et al. Modelling the epidemiological impact of intermittent preventive treatment against malaria in infants. PLoS ONE. 2008;3:e2661.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Camponovo F, Ockenhouse CF, Lee C, Penny MA. Mass campaigns combining antimalarial drugs and anti-infective vaccines as seasonal interventions for malaria control, elimination and prevention of resurgence: a modelling study. BMC Infect Dis. 2019;19:920.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Reiker T, Chitnis N, Smith T. Modelling reactive case detection strategies for interrupting transmission of Plasmodium falciparum malaria. Malar J. 2019;18:259.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Tanzania National Bureau of Statistics (NBS), ICF Macro. Tanzania Demographic and Health Survey 2010. Dar es Salaam, Tanzania, and Rockville, USA; 2011. http://dhsprogram.com/publications/publication-FR243-DHS-Final-Reports.cfm.

  40. Tanzania Commission for AIDS, Zanzibar AIDS Commission, National Bureau of Statistics, Office of the Chief Government Statistician, ICF International 2013. Tanzania HIV/AIDS and Malaria Indicator Survey 2011–12. Dar es Salaam, Tanzania, and Rockville, USA; 2013.

  41. Ministry of Health, Community Development, Gender, Elderly and Children, Zanzibar Ministry of Health, Tanzania National Bureau of Statistics, Office of Chief Government Statistician, ICF International. Tanzania Demographic and Health Survey and Malaria Indicator Survey 2015–2016. Dar es Salaam, Tanzania and Rockville, Maryland, USA, 2016.

  42. Galactionova K, Tediosi F, de Savigny D, Smith T, Tanner M. Effective coverage and systems effectiveness for malaria case management in sub-Saharan African countries. PLoS ONE. 2015;10:e0127818.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Lalji S, Ngondi JM, Thawer NG, Tembo A, Mandike R, Mohamed A, et al. School distribution as keep-up strategy to maintain universal coverage of long-lasting insecticidal nets: implementation and results of a program in Southern Tanzania. Glob Health Sci Pract. 2016;4:251–63.

    Article  PubMed  PubMed Central  Google Scholar 

  44. National Bureau of Statistics, Tanzania, Office of Chief Government Statistician, Zanzibar. 2012 Population and Housing Census. Dar es Salaam, Tanzania, 2013.

  45. Tatem AJ, Qiu Y, Smith DL, Sabot O, Ali AS, Moonen B. The use of mobile phone data for the estimation of the travel patterns and imported Plasmodium falciparum rates among Zanzibar residents. Malar J. 2009;8:287.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, et al. An approach to model the costs and effects of case management of Plasmodium falciparum malaria in Sub-Saharan Africa. Am J Trop Med Hyg. 2006;75:90–103.

    Article  PubMed  Google Scholar 

  47. Kilian A, Byamukama W, Pigeon O, Atieli F, Duchon S, Phan C. Long-term field performance of a polyester-based long-lasting insecticidal mosquito net in rural Uganda. Malar J. 2008;7:49.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Stelmach R, Millikan E, Lalji S, Colaço R. Tanzania Vector Control Scale-up Project (TVCSP): cost of indoor residual spraying 2008–2015. RTI International, PMI, 2016.

  49. Tchicaya ES, Nsanzabana C, Smith TA, Donzé J, de Hipsl ML, Tano Y, et al. Micro-encapsulated pirimiphos-methyl shows high insecticidal efficacy and long residual activity against pyrethroid-resistant malaria vectors in central Côte d’Ivoire. Malar J. 2014;13:332.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Agossa FR, Aïkpon R, Azondékon R, Govoetchan R, Padonnou GG, Oussou O, et al. Efficacy of various insecticides recommended for indoor residual spraying: pirimiphos methyl, potential alternative to bendiocarb for pyrethroid resistance management in Benin, West Africa. Trans R Soc Trop Med Hyg. 2014;108:84–91.

    Article  CAS  PubMed  Google Scholar 

  51. Farooqui HH, Hussain MA, Zodpey S. Malaria control in India: has sub-optimal rationing of effective interventions compromised programme efficiency? WHO South-East Asia J Public Health. 2012;1:128.

    Article  PubMed  Google Scholar 

  52. National Malaria Control Programme; World Health Organization Tanzania; Ifakara Health Institute; KEMRI-Wellcome Trust; Report funded by Roll Back Malaria and Department for International Development-UK. An epidemiological profile of malaria and its control in mainland Tanzania. Dar es Salaam, Tanzania, 2013. http://www.inform-malaria.org/wp-content/uploads/2014/05/Tanzania-Epi-Report-060214.pdf.

  53. Yukich JO, Chitnis N. Modelling the implications of stopping vector control for malaria control and elimination. Malar J. 2017;16:411.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Crowell V, Hardy D, Briët O, Chitnis N, Maire N, Smith T. Can we depend on case management to prevent re-establishment of P. falciparum malaria, after local interruption of transmission? Epidemics. 2012;4:1–8.

    Article  PubMed  Google Scholar 

  55. Noor AM, Alegana VA, Patil AP, Snow RW. Predicting the unmet need for viologically targeted coverage of insecticide-treated nets in Kenya. Am J Trop Med Hyg. 2010;83:854–60.

    Article  PubMed  PubMed Central  Google Scholar 

  56. President’s Malaria Initiative (PMI). Malaria Operational Plan 2019 - Kenya. https://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy19/fy-2019-kenya-malaria-operational-plan.pdf?sfvrsn=3.

  57. President’s Malaria Initiative (PMI). Malaria Operational Plan FY 2019 - Zambia. https://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy19/fy-2019-zambia-malaria-operational-plan.pdf?sfvrsn=3.

  58. National Vector-borne Diseases Control Programme; Ministry of Health and Social Services. Malaria Strategic Plan (2010–2016). Republic of Namibia: 2010. https://endmalaria.org/sites/default/files/namibia2010-2016.pdf.

  59. President’s Malaria Initiative (PMI). Malaria Operational Plan FY 2018 – Zimbabwe. https://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy-2018/fy-2018-zimbabwe-malaria-operational-plan.pdf?sfvrsn=7

  60. Weiss DJ, Bertozzi-Villa A, Rumisha SF, Amratia P, Arambepola R, Battle KE, et al. Indirect effects of the COVID-19 pandemic on malaria intervention coverage, morbidity, and mortality in Africa: a geospatial modelling analysis. Lancet Infect Dis. 2021;21:59–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Chacky F, Runge M, Rumisha SF, Machafuko P, Chaki P, Massaga JJ, et al. Nationwide school malaria parasitaemia survey in public primary schools, the United Republic of Tanzania. Malar J. 2018;17:452.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Gitonga CW, Edwards T, Karanja PN, Noor AM, Snow RW, Brooker SJ. Plasmodium infection, anaemia and mosquito net use among school children across different settings in Kenya. Trop Med Int Health. 2012;17:858–70.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Coalson JE, Walldorf JA, Cohee LM, Ismail MD, Mathanga D, Cordy RJ, et al. High prevalence of Plasmodium falciparum gametocyte infections in school-age children using molecular detection: patterns and predictors of risk from a cross-sectional study in southern Malawi. Malar J. 2016;15:527.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Staedke SG, Maiteki-Sebuguzi C, Rehman AM, Kigozi SP, Gonahasa S, Okiring J, et al. Assessment of community-level effects of intermittent preventive treatment for malaria in schoolchildren in Jinja, Uganda (START-IPT trial): a cluster-randomised trial. Lancet Glob Health. 2018;6:e668–79.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Nankabirwa JI, Wandera B, Amuge P, Kiwanuka N, Dorsey G, Rosenthal PJ, et al. Impact of intermittent preventive treatment with dihydroartemisinin-piperaquine on malaria in Ugandan schoolchildren: a randomized, placebo-controlled trial. Clin Infect Dis. 2014;58:1404–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Matangila JR, Doua JY, Mitashi P, da Luz RI, Lutumba P, Van Geertruyden JP. Efficacy and safety of intermittent preventive treatment in schoolchildren with sulfadoxine/pyrimethamine (SP) and SP plus piperaquine in Democratic Republic of the Congo: a randomised controlled trial. Int J Antimicrob Agents. 2017;49:339–47.

    Article  CAS  PubMed  Google Scholar 

  67. O’Meara WP, Breman JG, McKenzie FE. The promise and potential challenges of intermittent preventive treatment for malaria in infants (IPTi). Malar J. 2005;4:33.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Lindblade KA, Mwandama D, Mzilahowa T, Steinhardt L, Gimnig J, Shah M, et al. A cohort study of the effectiveness of insecticide-treated bed nets to prevent malaria in an area of moderate pyrethroid resistance. Malawi Malar J. 2015;14:31.

    Article  PubMed  Google Scholar 

  69. Ochomo EO, Bayoh NM, Walker ED, Abongo BO, Ombok MO, Ouma C, et al. The efficacy of long-lasting nets with declining physical integrity may be compromised in areas with high levels of pyrethroid resistance. Malar J. 2013;12:368.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Pryce J, Richardson M, Lengeler C. Insecticide-treated nets for preventing malaria. Cochrane Database Syst Rev. 2018;11:CD000363.

    PubMed  Google Scholar 

  71. Churcher TS, Lissenden N, Griffin JT, Worrall E, Ranson H. The impact of pyrethroid resistance on the efficacy and effectiveness of bednets for malaria control in Africa. Elife. 2016;5:e16090.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Fuseini G, Ebsworth P, Jones S, Knight D. The efficacy of ACTELLIC 50 EC, Pirimiphos Methyl, for indoor residual spraying in Ahafo, Ghana: Area of high vector resistance to pyrethroids and organochlorines. J Med Entomol. 2011;48:437–40.

    Article  CAS  PubMed  Google Scholar 

  73. Aïkpon R, Sèzonlin M, Tokponon F, Okè M, Oussou O, Oké-Agbo F, et al. Good performances but short lasting efficacy of Actellic 50 EC Indoor Residual Spraying (IRS) on malaria transmission in Benin. West Africa Parasit Vectors. 2014;7:256.

    Article  PubMed  Google Scholar 

  74. Haji KA, Thawer NG, Khatib BO, Mcha JH, Rashid A, Ali AS, et al. Efficacy, persistence and vector susceptibility to pirimiphos-methyl (Actellic® 300CS) insecticide for indoor residual spraying in Zanzibar. Parasit Vectors. 2015;8:628.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Mashauri FM, Manjurano A, Kinunghi S, Martine J, Lyimo E, Kishamawe C, et al. Indoor residual spraying with micro-encapsulated pirimiphos-methyl (Actellic® 300CS) against malaria vectors in the Lake Victoria basin. Tanzania PLoS One. 2017;12:e0176982.

    Article  PubMed  Google Scholar 

  76. Protopopoff N, Mosha JF, Lukole E, Charlwood JD, Wright A, Mwalimu CD, et al. Effectiveness of a long-lasting piperonyl butoxide-treated insecticidal net and indoor residual spray interventions, separately and together, against malaria transmitted by pyrethroid-resistant mosquitoes: a cluster, randomised controlled, two-by-two factorial design trial. Lancet. 2018;391:1577–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Agossa FR, Padonou GG, Koukpo CZ, Zola-Sahossi J, Azondekon R, Akuoko OK, et al. Efficacy of a novel mode of action of an indoor residual spraying product, SumiShield® 50WG against susceptible and resistant populations of Anopheles gambiae (s.l.) in Benin. West Africa Parasit Vectors. 2018;11:293.

    Article  PubMed  Google Scholar 

  78. WHO recommends groundbreaking malaria vaccine for children at risk. Geneva, World Health Organization, 2021. https://www.who.int/news/item/06-10-2021-who-recommends-groundbreaking-malaria-vaccine-for-children-at-risk. Accessed 6 Oct 2021.

  79. President’s Malaria Initiative (PMI). Malaria Operational Plan FY 2019 - Tanzania. https://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy19/fy-2019-tanzania-malaria-operational-plan.pdf?sfvrsn=3.

  80. Runge M, Mapua S, Nambunga I, Smith TA, Chitnis N, Okumu F, et al. Evaluation of different deployment strategies for larviciding to control malaria: a simulation study. Malar J. 2021;20:324.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Alegana VA, Suiyanka L, Macharia PM, Ikahu-Muchangi G, Snow RW. Malaria micro-stratification using routine surveillance data in Western Kenya. Malar J. 2021;20:22.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Smith JL, Mumbengegwi D, Haindongo E, Cueto C, Roberts KW, Gosling R, et al. Malaria risk factors in northern Namibia: The importance of occupation, age and mobility in characterizing high-risk populations. PLoS ONE. 2021;16:e0252690.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Kabaria CW, Gilbert M, Noor AM, Snow RW, Linard C. The impact of urbanization and population density on childhood Plasmodium falciparum parasite prevalence rates in Africa. Malar J. 2017;16:49.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Kabaria CW, Molteni F, Mandike R, Chacky F, Noor AM, Snow RW, et al. Mapping intra-urban malaria risk using high resolution satellite imagery: a case study of Dar es Salaam. Int J Health Geogr. 2016;15:26.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Lowa M, Sitali L, Siame M, Musonda P. Human mobility and factors associated with malaria importation in Lusaka district, Zambia: a descriptive cross sectional study. Malar J. 2018;17:404.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, et al. Quantifying the impact of human mobility on malaria. Science. 2012;338:267–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. The malERA Consultative Group on Modeling. A research agenda for malaria eradication: modelling. PLoS Med. 2011;8:e1000403.

    Article  PubMed Central  Google Scholar 

  88. Lindsay SW, Bødker R, Malima R, Msangeni HA, Kisinza W. Effect of 1997–98 El Niño on highland malaria in Tanzania. Lancet. 2000;355:989–90.

    Article  CAS  PubMed  Google Scholar 

  89. Kilian AHD, Langi P, Talisuna A, Kabagambe G. Rainfall pattern, El Niño and malaria in Uganda. Trans R Soc Trop Med Hyg. 1999;93:22–3.

    Article  CAS  PubMed  Google Scholar 

  90. Bangs MJ, Subianto DB. El Niño and associated outbreaks of severe malaria in highland populations in Irian Jaya, Indonesia: a review and epidemiological perspective. Southeast Asian J Trop Med Public Health. 1999;30:608–19.

    CAS  PubMed  Google Scholar 

  91. Mabaso MLH, Kleinschmidt I, Sharp B, Smith T. El Niño Southern Oscillation (ENSO) and annual malaria incidence in Southern Africa. Trans R Soc Trop Med Hyg. 2007;101:326–30.

    Article  PubMed  Google Scholar 

  92. Dawood FS, Ricks P, Njie GJ, Daugherty M, Davis W, Fuller JA, et al. Observations of the global epidemiology of COVID-19 from the prepandemic period using web-based surveillance: a cross-sectional analysis. Lancet Infect Dis. 2020;20:1255–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Aborode AT, David KB, Uwishema O, Nathaniel AL, Imisioluwa JO, Onigbinde SB, et al. Fighting COVID-19 at the expense of malaria in Africa: the consequences and policy options. Am J Trop Med Hyg. 2021;104:26–9.

    Article  PubMed  Google Scholar 

  94. Heuschen A-K, Lu G, Razum O, Abdul-Mumin A, Sankoh O, von Seidlein L, et al. Public health-relevant consequences of the COVID-19 pandemic on malaria in sub-Saharan Africa: a scoping review. Malar J. 2021;20:339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Sherrard-Smith E, Hogan AB, Hamlet A, Watson OJ, Whittaker C, Winskill P, et al. The potential public health consequences of COVID-19 on malaria in Africa. Nat Med. 2020;26:1411–6.

    Article  CAS  PubMed  Google Scholar 

  96. WHO. World malaria report 2020. Geneva: World Health Organization; 2020. https://www.who.int/publications/i/item/9789240015791.

Download references

Acknowledgements

The authors are grateful to all the participants who contributed to the expert discussion during the strategic planning workshop in Bagamoyo in May 2018, especially the members of the National Malaria Control Programme, of the Ministry of Health, Community Development, Gender, Elderly and Children of mainland Tanzania, President’s Office Regional Administration and Local Government offices, WHO country office and development partners. The authors thank Munir Winkel for helpful comments on earlier versions of the manuscript. The simulations in OpenMalaria were performed at sciCORE (http://scicore.unibas.ch/), the scientific computing core facility at University of Basel. MR acknowledges Northwestern University for providing resources to finalize the manuscript.

Funding

Funding was received for EP and MR from the Global Fund to Fight Aids, Tuberculosis and Malaria (GFATM), the NETCELL Project funded by the Swiss Agency for Development and Cooperation, and the Swiss Tropical and Public Health Institute. RWS is supported as a Wellcome Trust Principal Fellow (# 079080, 103602 and 212176). RWS also acknowledges support from the UK’s Department for International Development under a project entitled Strengthening the Use of Data for Malaria Decision Making in Africa (DFID Programme Code # 203155).

Author information

Authors and Affiliations

Authors

Contributions

MR performed the analysis and wrote the draft manuscript. ST, FM, FC, RM, AM, SM supported interpretation of modelling results during the strategic planning workshop and contributed to the discussions of the selection of scenarios to simulate. RWS, CL and EP provided overall guidance in the project and in the preparation of the draft. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Manuela Runge or Emilie Pothin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Strategic planning questions for modelling in mainland Tanzania.

Additional file 2.

Historical trends per strata and fitting performance.

Additional file 3.

Results table as included in the revised strategic plan.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Runge, M., Thawer, S.G., Molteni, F. et al. Sub-national tailoring of malaria interventions in Mainland Tanzania: simulation of the impact of strata-specific intervention combinations using modelling. Malar J 21, 92 (2022). https://doi.org/10.1186/s12936-022-04099-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-022-04099-5

Keywords