Abstract
Background Malaria transmission is primarily limited to tropical regions where environmental conditions are conducive for the survival of Plasmodium parasites and Anopheles mosquitoes. Adequate rainfall provides breeding sites, while suitable temperatures facilitate vector mosquito life-cycles and parasite development. Evaluating the efficacy of vector control interventions is crucial to determine their effectiveness in reducing malaria transmission. The aim of this study was to explore how these factors affect transmission dynamics at varying levels of vector control efficacy.
Methods We developed a vector-host compartmental mathematical model to compare three published approaches to incorporating weather influences on malaria transmission. The first approach examines mosquito biting behaviour and mortality rates in larval and adult stages. The second focuses on temperature effects on mosquito life-cycle characteristics throughout the aquatic and adult stages. The third considers how temperature and rainfall influence adult mosquito behaviour, environmental carrying capacity, and survival during the aquatic stages. Model simulations were conducted at different annual vector control coverage levels, to identify variations in transmission patterns and seasonal variability in daily and annual incidence across three climate regions.
Results The first approach indicates sustained seasonal transmission, with lower cases per 1,000 in tropical regions compared to semi-arid and sub-tropical regions, even with enhanced vector control reducing cases. The second approach predicts sharp, brief annual peaks, with zero transmission during winter in tropical and semi-arid regions, while sub-tropical regions experience ongoing transmission. In contrast, the third approach projects multiple irregular peaks, with transmission ceasing in winter across all regions.
Conclusions Simulations indicate that climatic events like heatwaves or flooding, can trigger mosquito population surges and malaria outbreaks, even in areas previously free of malaria, despite strong vector control efforts. However, the results demonstrate that sustained and effective vector control, particularly in regions with moderate temperatures, can substantially reduce malaria incidence. Effective malaria control requires incorporating weather predictions into intervention plans, and enhancing current vector control strategies with supplementary measures like larval source management. Accurate timing and targeting of these interventions, based on transmission season projections, are crucial for maintaining robust control as weather conditions evolve and to prepare for future challenges.
Background
Malaria is an infectious disease caused by Plasmodium parasites transmitted to humans through the bites of infectious female Anopheles mosquitoes [1]. This disease poses a significant public health concern, impacting the health and livelihoods of individuals in tropical regions worldwide [2, 3]. Malaria remains a significant cause of mortality, with around half a million deaths attributed to the disease annually [4–6]. The World Health Organization (WHO) categorises the prevalence of malaria in a country based on the number of reported cases per population at risk. Countries with fewer than 100 cases per 1000 population per year are classified as low transmission, while those reporting more than 450 cases per 1000 population per year are considered high transmission settings [7]. Several intervention strategies have been implemented to combat malaria transmission, these include case management through effective treatment with artemisinin-based combination therapies (ACTs) together with vector control measures such as indoor residual spraying (IRS) of insecticides and widespread distribution of long-lasting insecticide-treated nets (LLINs) [8–15]. These interventions play a crucial role in reducing the incidence of malaria and improving public health outcomes in affected areas. Despite this, seasonal outbreaks frequently occur due to favourable weather conditions and reduced efficacy of interventions. The influence of temperature and rainfall on factors such as mosquito abundance or changes in mosquito breeding, survival and behaviour contribute to these challenges [2, 16–18].
Malaria transmission is primarily restricted to tropical regions where temperatures are conducive for both parasite and mosquito development, and adequate rainfall to facilitate the availability of breeding sites for mosquitoes [8–10]. Mosquito breeding involves the laying of eggs in different water bodies, from temporary puddles to permanent lakes and dams with varying levels of salinity, turbidiy and pollution. After hatching, the larvae undergo four developmental stages, known as instars, before they pupate. These life stages are aquatic, while adult mosquitoes are the only terrestrial phase, having emerged from the pupae. The rate of maturation from one stage to the next is primarily dependent on temperature [19]. Higher temperatures can accelerate mosquito development within their upper thermal limits, leading to a higher likelihood of survival during the aquatic stages and an increase in the adult mosquito population. Similarly, the malaria parasite matures more quickly within mosquitoes at warmer temperatures, with optimal sporogony occurring around 25°C [20]. However, if temperatures exceed the upper thermal limit of 30°C, sporogony is halted, and neither the mosquitoes nor the malaria parasite are likely to survive [2, 5, 6, 20]. Additionally, adequate rainfall generates sufficient habitats for mosquito larvae to grow, leading to abundant adult mosquitoes capable of spreading the disease. However, limited rainfall results in fewer breeding sites and heavy rainfall leads to flushing out of mosquito larvae, both of which result in reduced mosquito populations available to spread the disease [11, 13–15]. It is important to consider the impact of weather factors like temperature and rainfall on mosquito development, survival, and behaviour when implementing vector control measures to manage malaria outbreaks. Accurate seasonal weather forecasts of these factors can be used in malaria models as early warning systems in endemic regions. These models can also assess potential changes in malaria prevalence due to seasonal weather variations.
Several mathematical models incorporating empirical data and statistical approaches have shown associations between variations in meteorological factors with malaria incidence [4, 12, 16, 18, 21–23]. While extensive research into the influence of temperature variability on malaria transmission has been conducted, rainfall has often been relegated to a secondary factor, despite efforts to develop suitable models that account for the combined impact of temperature and rainfall on vector and parasite development [2, 5, 24, 25]. The significance of these models has grown because while statistical models have been valuable in revealing relationships between environmental variables and transmission intensity, process-based mathematical models offer a more explanatory insight into the balance between internal factors (resulting from biological processes) and external factors (such as changes in environmental variables) that drive transmission. These dynamic models are essential because they account for the biological processes driving malaria transmission within an environment that changes over different time scales [22, 26–29]. Furthermore, a dynamic model is crucial for capturing invasion dynamics and effectively forecasting the emergence of new outbreaks. This includes scenarios where changes in temperature and rainfall patterns may render previously unsuitable areas conducive to transmission, as well as instances of human migration (and thereby Plasmodium parasite migration) or mosquito spread into previously unaffected regions [30–33].
A mechanistic model calibrated with weather data from multiple regions in Africa found that malaria infection tends to rise within temperatures ranging from 16 to 25°C, and decreases within the 25 to 28°C range [29]. In contrast, a study using a climate-based vector-host model, which considers the aquatic stages in the mosquito lifecycle, indicated that transmission is optimized in the temperature range of 21 to 25°C, accompanied by 95 to 125 mm of rainfall [25]. A vector-host modeling approach was employed to investigate the dynamics of malaria transmission, specifically focusing on the aquatic stages of mosquito vectors [11]. The model categorized adult mosquitoes based on their behaviours such as biting, resting, and host-seeking, while also considering the impact of ambient air temperature and water body temperature on mosquito breeding success. By quantifying the seasonality of Anopheles arabiensis population densities, the model accurately simulated the observed trends in larval density [1, 2, 11]. This alignment with laboratory experiments highlight the robustness of the model in capturing real-world conditions. In contrast, models that neglect the aquatic stages of mosquito vectors often fail to account for these crucial aspects of malaria transmission dynamics [4, 24, 33, 34].
While mechanistic models provide a valuable framework for understanding malaria dynamics, various approaches have been developed to incorporate the effects of temperature and rainfall on malaria transmission dynamics, focusing on different critical stages of transmission. [24] developed a weather-based model that examines malaria transmission by focusing on an age-structured vector population with periodic birth rates. Their model integrates factors such as temperature-dependent egg production, biting rates and mortality rates of both aquatic and adult mosquitoes. Similarly, [11] employed a weather-based mathematical model to investigate the impact of temperature on mosquito development and behaviours. This model incorporates temperature-dependent factors such as egg production, as well as the development and mortality rates of eggs, larvae and pupae, alongside adult mosquito activities such as resting, mating and host-seeking behaviours. Furthermore, [2] utilized a weather-based mathematical model to explore the combined effects of temperature and rainfall on malaria transmission dynamics. Their model includes mosquito recruitment dependent on temperature and rainfall, infection dynamics between humans and mosquitoes, and explores the impact of interventions such as insecticide spraying. Although these studies have explored the impact of temperature and rainfall on malaria transmission, none do so completely, and [24] do not explicitly consider aquatic stages. Instead they treat the combined aquatic population as a single unit and assume that its mortality rate is equivalent to that of the larvae only. Additionally, [11] explicitly account for these stages but focus solely on mosquito population dynamics, neglecting transmission. Moreover, [2] include the effects of rainfall but also neglect the aquatic population. Because these studies each focus on different aspects of malaria dynamics, it is essential to identify the most appropriate approach for investigating mosquito population dynamics and guiding resource allocation and decision-making for malaria control and prevention under varying weather conditions. This study aims to develop a compartmental mathematical model to evaluate the efficacy of vector control interventions by exploring how different approaches to incorporating weather influences affect malaria transmission dynamics. The research compares these approaches to identify the most effective framework for modeling the impacts of weather on malaria transmission and mosquito population dynamics across diverse settings. The study seeks to inform evidence-based strategies for optimizing vector control interventions, ensuring their timely and targeted implementation to sustain control and address emerging malaria transmission challenges under changing climatic conditions.
Methods
We constructed a compartmental model by grouping individuals with similar characteristics related to developmental stage, risk of infection, infectiousness, treatment seeking behaviour and recovery from disease into compartments within mosquito and human populations. The model diagram (see Figure 1) illustrates the transmission of malaria between mosquito vectors and human hosts, accounting for the aquatic stages of the mosquito life-cycle which includes eggs (Ea) that hatch into larva (La) that develop into pupae (Pa) that emerge as adult mosquitoes. We consider a collective population of larvae without explicitly accounting for the instar stages of larval development. The population of adult mosquitoes is grouped into mosquitoes susceptible to malaria infection (Sm), infected mosquitoes that cannot transmit infection to humans (Em), and infectious mosquitoes that can transmit the infection to humans (Im). Those mosquitoes that cannot transmit infections can become infectious if Plasmodium sporozoites (following sporogony in the midgut oocyts) have migrated to their salivary glands. The host population is grouped into individuals at risk of infection (S), infected humans that cannot transmit the infection to mosquitoes (E), and infected people that can transmit the infection to mosquitoes but either do not show symptoms (A), show uncomplicated symptoms (Iu) or are experiencing severe symptoms (Is). Furthermore, individuals treated for uncomplicated symptoms (Tu) or severe symptoms (Ts) and individuals that have recovered from the disease (R) are also considered.
Transmission flow diagram illustrating the development of mosquitoes from egg to larvae then pupae and the transmission of malaria between a human and mosquito population.
We adopt the technique discussed in [35], to investigate the role of vector control in effectively preventing transmission between mosquitoes and humans. Vector control coverage is considered as a compartment accounting for annual intervention deployment rate and loss of operational efficacy of insecticides used. The mosquito development and disease pathways illustrated in the model diagram as well as the intervention compartment are represented by the following system of ordinary differential equations (ODEs) describing malaria transmission between vector mosquitoes and human hosts. Full descriptions of the parameters and compartments are provided in Tables 1 and 2.
Model Compartments and Parameters
The transmission model is parameterized with values provided and described in Table 2. The total number of mosquito eggs laid by a population of female mosquitoes is represented by ne, with these eggs being oviposited into water bodies by individual female mosquitoes at a rate of θ and cleared at a rate of µe. Mosquito eggs hatch into larvae at a rate of κe, and competition for nutrients and predation leads to larvae dying at a rate of µl, whereas those that survive develop into pupae at a rate of κl. Clearance of pupae from water bodies occurs at a rate of µp and those that survives emerge as adult mosquitoes at a rate of κp. Adult mosquitoes die at a rate of µmand susceptible mosquitoes become infected at a rate of λm. The malaria parasite develops from gametocytes to sporozoites within infected mosquitoes at a rate of γm.
We consider that recruitment of individuals into the host population occurs at a birth rate of µb, with the mortality rates in humans given by µh. Susceptible humans become infected through infectious mosquito bites at a rate of λh and malaria parasite development within individuals exposed to infection occurs at a rate of γh. The proportion of infectious individuals that do not show symptoms is given by pa and the loss of uncomplicated symptoms occurs at a rate of ω. Infectious people with uncomplicated symptoms receive treatment at a rate of τu, whereas those that do not receive treatment recover from disease at a rate of δr. Infected people experience severe malaria symptoms at a rate of ν and are lost at a rate of α. People with severe symptoms receive treatment at a rate of τs but suffer malaria induced death at rate of µs. People that have been treated for uncomplicated or severe symptoms recover from disease at rates δu and δs, respectively. Finally, recovered people are considered to possess temporary immunity acquired from infection and lose this immunity to become susceptible to infection again at a rate of ρ.
Comparison of Modelling Approaches
Three approaches were adopted to investigate the effects of temperature and rainfall on malaria dynamics while accounting for different aspects of malaria transmission and the mosquito life-cycle. Approach A adopted in [24], focuses on aspects of transmission that are critical to the effectiveness of vector control that include mosquito biting behaviour, as well as mortality of adult mosquitoes and larvae. Vector control through IRS or ITNs limits contact between susceptible individuals and mosquitoes, which protects against infectious mosquito bites. Temperature also limit exposure to infected mosquitoes and the number of mosquitoes available to transmit infection by influencing the number of larvae that die and the lifespan of adult mosquitoes. These factors are also considered in Approach B used in [11], which extends the role of temperature to the development and mortality of mosquito eggs, larvae, and pupae. The effects of temperature on the total number of eggs laid and the rate of ovipositing by a female mosquito during her lifespan are also considered. Despite overlap with both these approaches, [2] adopted Approach C to account for the developmental period during aquatic stages and the survival probability of eggs, larvae and pupae in habitats with changing temperatures. Furthermore, temperature dependence of the extrinsic incubation period for parasites to develop within mosquitoes and rainfall dependence of the environmental carrying capacity are considered.
Approach A: Temperature-Regulated Mosquito Mortality
We adopt the approach of [24] and [1] by modelling aspects of the mosquito life-cycle that are crucial to mosquito survival and malaria transmission. How many eggs survive the larval stage of the mosquito life-cycle to become pupae influences the growth of the adult mosquito population which is regulated through fluctuating mortality at different temperatures. Mosquitoes that survive into adulthood and come into contact with humans drive transmission through biting and blood feeding. Seasonal changes in temperature will affect this behaviour, especially when outside the optimal thermal limits of the mosquito population. Evidence from Mpumalanga province, South Africa, suggests that the adults of the malaria vector species An. arabiensis and An. parensis do not undergo a period of suspended development during the colder and drier winter months. Instead, adult female mosquitoes continue to seek for blood meals in winter but their population densities are significantly reduced by the limited number of breeding sites available. Furthermore, colder temperatures likely reduce their blood feeding propensity and therefore their fertility [37]. It is assumed that mosquito larvae and adult mortality rates are at least partially temperature dependent, that is µl(T) and µm(T) for larvae and adult mosquitoes, respectively, and that the biting rate, a(T) (see Figure 2 in Appendix 4 for plots of these functions).
Projected aquatic mosquito populations across regions.
Approach B: Temperature-Sensitive Mosquito Life-Cycle
Here the approach of [11] is extended to account for the temperature dependence of mosquito characteristics related to the development of eggs, larvae and pupae are explicitly considered. We consider the total number of eggs laid by a population of female mosquitoes ne(T), the daily rate of egg oviposition θ(T), the development rates of eggs κe(T), larvae κl(T) and pupae κp(T), the mortality rates of eggs µe(T), larvae µl(T) and pupae µp(T), as well as the biting rate of adult mosquitoes a(T) (see Figure 3 in Appendix 4 for plots of these functions).
Projected adult mosquito populations across regions.
Approach C: Temperature and Rainfall Dependent Aquatic Survival
This last approach is presented in [2], where the environmental carrying capacity of mosquito eggs Ke(R) and the survival probabilities of eggs pe(R), larvae pl(R, T) and pupae pp(R) depend on rainfall. The rainfall threshold for flushing out of aquatic mosquitoes due to heavy rainfall is denoted by Rl. Furthermore, the biting rate a(T), mortality rate of aquatic and adult mosquitoes, µl(T) and µm(T) are included. The development period of larvae tl(T), and the incubation period of malaria parasites within mosquitoes γm(T) are all considered temperature dependent (see Figure 4 in Appendix 4 for plots of these functions).
Projected host populations across regions.
This study uses historical mean monthly surface air temperature and precipitation datasets for the 1950 to 2022 period in three regions with semi-arid, tropical and subtropical climate conditions. Temperatures in the semi-arid region that had the highest peaks around 30 degrees Celsius, followed by tropical then the sub-tropical region which also experienced the lowest temperatures around 15 degrees Celsius during the 2016 to 2022 period. The tropical region accumulated the largest amount of rainfall peaking around 250 to 300 mm during months of heavy rainfall, followed by the semi-arid region with slightly less rainfall peaking around 150 to 250 mm despite experiencing the longest periods without rainfall. The sub-tropical region had the least amount of rainfall peaking below 150 mm during periods of heavy rainfall between 2016 and 2022. The temperature and rainfall trends for the period investigated in the study are presented in Figure 1 in Appendix 4. The model is simulated with the same initial conditions and varying parameters across three modelling approaches A, B and C. Furthermore, we explore vector control scenarios with annual coverage levels ranging from 25 % to 90 % in three climate regions. Simulations are conducted for the entire period of seventy two years from 1950 to 2022 to observe annual trends in malaria incidence and variations in seasonal patterns are analysed for the last six years from 2016 to 2022. The impact of temperature and rainfall on malaria transmission in three regions are compared through projected aquatic and adult mosquitoes populations as well as the number of new cases, populations at-risk of infection, asymptomatic populations and annual incidence rates. We conducted our comparative study using formulations of the temperature and rainfall dependent model parameters derived in studies by [2, 11] and [24].
Results
Aquatic Mosquito Population Dynamics
Three approaches were adopted to integrate the influence of weather conditions on the development and survival of Anopheles mosquito eggs, larvae and pupae provided vector control intervention in these regions. Approach A projects relatively low but persistent aquatic mosquito populations with eggs and larvae populations peaking around 10 million in semi-arid and tropical regions, and 12 million in the sub-tropical region, whereas pupa populations remain substantially lower in all three regions. This approach projects a larger egg population across all regions due to more sustained adult populations, leading to higher exposure to infection than Approaches B and C. In addition, Approach A suggests that cooler temperatures in the sub-tropical region lead to moderate adult mosquito mortality, sustaining a more stable adult population capable of continuous egg production. However, the warmer temperatures in the semi-arid and tropical regions increase mosquito mortality, thereby reducing adult populations. Approach B suggests seasonal growth in aquatic mosquitoes populations with egg populations peaking around 40 million in tropical and semi-arid regions despite the short lived abundance of aquatic mosquitoes in the semi-arid region. Although the egg populations in the sub-tropical region peak around 25 million, these populations remain consistent throughout the year while experiencing momentary periods of low aquatic populations compared to semi-arid and tropical regions. Despite the higher mortality of adult mosquitoes in warmer climates, these temperatures accelerate egg-laying and the number of eggs laid, supporting growth in egg populations. On the other hand, Approach C highlights that sufficient rainfall creates ideal aquatic environments, allowing mosquito eggs, larvae, and pupae to thrive across a range of temperatures. Approach C projects very low aquatic populations throughout while experiencing seasonal spikes in egg populations peaking around 60, 75 and 80 million in semi-arid, sub-tropical and tropical regions, respectively. Additionally, this approach suggests that variations in rainfall help temporarily maintain egg populations, especially at moderate temperatures in the tropical region.
Adult Mosquito Population Trends
We observe different population trends in adult mosquito populations across different regions, with susceptible mosquitoes accounting for the majority of the population for all approaches. Approach A projects persistent aquatic mosquitoes populations and Approach B suggests seasonal peaks in these populations, whereas short periods with abundant in aquatic mosquitoes observed in Approach C. Approach B suggests that while warmer temperatures in the tropical region increase mosquito mortality, they ultimately reduce the risk of infection compared to the sub-tropical region, where adult populations remain relatively stable. Approach C, however, shows that lower mosquito mortality rates in the sub-tropical region result in a sustained mosquito population, whereas warmer conditions in the semi-arid region typically lead to higher mortality rates. Approach A suggests that susceptible mosquito populations range between 1 to 3 million in the tropical region whereas this population may grow to 4 million mosquitoes in semi-arid and sub-tropical regions. This approach further indicates that cooler temperatures in the sub-tropical region limit mosquito biting, while moderate temperatures in the tropical region increase the risk of mosquito infection. Exposed and infectious mosquito populations remain relatively low while peaking around 1 million mosquitoes across regions. These populations persist throughout the year contributing to malaria transmission in these areas, particularly in semi-arid and tropical regions where these populations fluctuate rapidly. Approach B projects transitory adult mosquito populations in the semi-arid and tropical regions peaking between 1 and 2 million susceptible mosquitoes followed by exposed and infectious mosquitoes at relatively lower populations. This approach projects stable adult mosquito populations in the sub-tropical region, with populations of susceptible, exposed and infectious mosquitoes remaining relatively similar across the years below a 1 million mosquitoes in each population. Approach C suggests that adult mosquito populations experience sharp, short-lived spikes throughout the years, which are characterized by irregular surges in infectious mosquito populations that remain very low across all regions during this period. These populations peak below 1 million mosquitoes in semi-arid and tropical regions, whereas susceptible mosquito populations may grow up to 4 million in the sub-tropical region.
Host Population Dynamics
All approaches project seasonal peaks in susceptible populations that alternate with the asymptomatic populations in semi-arid, sub-tropical and tropical regions. These populations contribute to the majority of the total host population, followed by exposed individuals and those that have recovered from malaria. Approach A suggests that susceptible populations peak around 400 thousand in semi-arid and tropical regions but may rise to 600 thousand in the sub-tropical region. In contrast, Approach B projects that susceptible populations peak around 600 thousand in semi-arid and tropical regions and may drop to 400 thousand in the sub-tropical region, whereas Approach C suggests that populations at risk of infection remain around 600 thousand. Asymptomatic populations in all regions stabilise around 250 thousand despite differing transmission seasons across Approaches. Approach A indicates that warmer temperatures in the semi-arid and tropical regions promote increased mosquito biting, which increases exposure risk compared to the cooler sub-tropical region. In contrast, Approach B indicates that lower mortality rates in the sub-tropical region enhance mosquito biting, whereas biting is restricted in the semi-arid and tropical regions due to higher mosquito mortality at warmer temperatures. Approach C suggests that adequate rainfall may sustain mosquito populations which may continue to transmit malaria infection even in region with less favourable temperatures. Approach A projects annual peaks in asymptomatic populations which persist for longer in tropical compared to semi-arid and sub-tropical regions. However, Approach B suggests that these populations are sustained for longer in the sub-tropical rather than the tropical region. Approach C predicts several irregular peaks in infectious mosquito populations throughout the year across all regions, which are driven by more frequent periods of optimal weather conditions as a result of ideal rainfall and temperature aligning to support mosquito development. The sporadic surges in infectious mosquito populations when Approach C is adopted, lead to sharp, short-lived peaks in new malaria cases in humans. Nonetheless, Approach B indicates that seasonal growth rather than irregular peaks in infectious mosquitoes, results in greater exposure and higher peaks in new malaria cases even compared to the persistent transmission in Approach A.
Daily Malaria Incidence Projections
Figure 5 shows that all three approaches project different transmission seasons where Approach A projects sustained seasonal transmission peaking around 20 new cases per 1000 in the tropical region and 25 new cases per 1000 in semi-arid and sub-tropical regions. Improved vector control intervention from 25 % to 90 % annual coverage leads to reductions in new malaria cases, particularly in the sub-tropical region. Approach B projects sharp annual peaks in new malaria cases peaking around 25, 60 and 75 cases per 1000 in sub-tropical, tropical and semi-arid regions, respectively. These peaks are short lived in the semi-arid region with cases dropping to zero in winter as a result of low temperatures which limit contact between mosquitoes and humans in populations at risk which reduces the number of successful infectious bites during this period. Approach B also projects no transmission occurs in winter when temperatures decline in the tropical region, malaria transmission persists slightly longer after the surge in new cases. In contrast, transmission in the sub-tropical regions continues throughout winter as temperatures do not drop low enough to stop mosquito biting and as a result the spread of malaria infection. Approach C also projects short lived peaks in new malaria cases and that transmission stops during winter across all regions, however this approach projects multiple peaks throughout the year. Provided there is irregularity in cases across regions, cases peak around 60 to 75 new cases per 1000 in tropical and semi-arid regions, whereas in the sub-tropical region where cases peak below 75 cases per 1000, cases may rise up to 100 new cases per 1000 as experienced in 2018.
Daily malaria incidence rates across regions.
Annual Malaria Transmission Trends
Figure 6 shows plots of annual trends in malaria transmission in semi-arid, sub-tropical and tropical regions during the 1950 to 2022 period. Approach A suggests that malaria incidence ranges between 40 to 60 new cases per 100 000 per year in semi-arid and tropical regions, with a slight increase throughout the simulation period. Whereas cases remain around 10 to 40 new cases per 100 000 per year in the sub-tropical region despite fluctuations and reduction in annual cases as vector control coverage improves. Approach B suggests that annual cases in the sub-tropical region remain fairly consistent and there is a decline in annual cases from 80 to around 40 cases per 100 000 per year in the tropical and around 60 to below 20 cases per 100 000 per year in the semi-arid region. Approach C projects the least incidence with relatively stable transmission leading to a gradual rise in new cases peaking below 40 cases per 100 000 in all regions. Cases remain around 20 new cases per 100 000 per year in semi-arid and tropical regions, while cases may rise up to 40 new cases per 100 000 per year. Although malaria transmission is halted around 1960 and again around 1980 in the sub-tropical region when vector control is sustained at 90 % annual coverage, malaria transmission persisted until 2022.
Annual malaria incidence rates across regions.
Discussion
This study used monthly mean rainfall and temperature datasets obtained from the Climate Change Knowledge Portal (CCKP) for the 1950 to 2022 period in three areas representing semi-arid, tropical and sub-tropical regions. Model simulations were conducted for this period and annual outcomes of malaria incidence are observed at this scale, however the seasonal patterns are explored at a shorter period from 2016 to 2022. Seasonal increases in asymptomatic populations play a key role in sustaining malaria transmission, particularly in Approach C, which demonstrates irregular peaks in adult mosquito populations. This pattern is influenced by temperature and rainfall, as higher temperatures accelerate mosquito development and parasite maturation, leading to an increased risk of malaria transmission. As observed in Approach A, moderate temperatures in the tropical region regulate mosquito survival, creating a larger mosquito population at risk compared to the sub-tropical region, which promotes more frequent biting in comparison to the semi-arid region. However, in the semi-arid region, as shown in Approach B, warmer temperatures accelerate mosquito biting, but the impact is less pronounced than in the tropical region. Approach C suggests that while rainfall plays a critical role in mosquito population dynamics, heavy rainfall periods in the semi-arid and tropical regions can inhibit mosquito growth by flushing out mosquito larvae, reduce mosquito populations and potentially limit the spread of malaria. On the other hand, insufficient rainfall reduces available breeding sites, which further lowers mosquito abundance and, consequently, malaria transmission. On the other hand, although Approach C projects sharp peaks in daily new cases, these are more sporadic and irregular. This pattern reflects the impact of rainfall on mosquito populations, as highlighted in previous studies: periods of heavy rainfall can flush out mosquito larvae, while insufficient rainfall limits breeding sites, both leading to fluctuating mosquito populations. Therefore, while Approach B’s seasonal growth in malaria cases indicates a more stable transmission cycle, the sporadic surges predicted by Approach C highlight the role of irregular environmental factors in driving short-term increases in malaria incidence.
Approach A projects the lowest daily incidence rates in the sub-tropical region due to cooler temperatures that limit mosquito biting, a finding consistent with past studies highlighting the role of temperature in regulating mosquito survival and development. Cooler temperatures reduce both mosquito and parasite maturation, limiting transmission, as supported by research indicating that mosquito development and malaria parasite maturation are hindered outside optimal temperature ranges [8, 11, 12, 14]. On the other hand, Approach B projects higher annual malaria cases in the semi-arid and tropical regions due to extended peak transmission periods. This pattern aligns with findings from past studies that associate warmer temperatures with accelerated mosquito development, which in turn prolongs the periods of high mosquito activity and malaria transmission [2, 11, 13, 14]. This explains the higher cumulative cases under Approach B, where extended warm periods foster more consistent transmission. Approach C projects sharp peaks in daily new cases which are more sporadic and irregular, reflecting fluctuating mosquito populations driven by rainfall variability. Heavy rainfall can wash away mosquito larvae, reducing mosquito numbers, while insufficient rainfall results in fewer breeding sites, both contributing to the irregular surges seen in Approach C. This pattern highlights the importance of process-based models that consider both temperature and rainfall impacts on mosquito and parasite development. Dynamic models of malaria transmission, which account for these interactions provide key insights into the duration of transmission periods as well as transmission intensity and frequency. These explain how prolonged warm conditions can sustain malaria transmission over time, as observed in Approach A and B, and how environmental fluctuations contribute to the sporadic surges in Approach C. Although rainfall has often been considered a secondary factor in malaria transmission, these results highlight its critical role in shaping mosquito populations, especially in regions with irregular rainfall patterns, as seen in Approach C.
Effective vector control targeting adult mosquitoes reduces annual incidence rates across all regions, with Approach A showing the most substantial reduction. A mechanistic model calibrated with weather data from various African regions found that malaria infection rates rise when temperatures range from 16 to 25°C but decrease when temperatures exceed 25°C. In the semi-arid region, temperatures peaked around 30°C, creating a challenging environment for mosquito survival. This suggests that adult mosquito populations, which are crucial for transmission, are most abundant and active within the lower range of these temperatures, supporting the substantial reductions in transmission observed under Approach A. Studies have shown that temperatures above 28°C tend to reduce transmission due to the stress on both mosquito populations and the malaria parasite [19, 25, 29]. However, the combination of high temperatures and intermittent rainfall results in a variable environment, allowing mosquito populations to thrive during wetter periods, leading to seasonal peaks in transmission, as seen in Approach C. These irregular surges in transmission reflect the fluctuating nature of the environment in this region.In contrast, the tropical region, with temperatures ranging from 25°C to 30°C and the highest rainfall (250 to 300 mm), creates optimal conditions for mosquito breeding. The malaria parasite also matures best at temperatures between 21 and 25°C, suggesting that the tropical region likely experiences a consistent transmission cycle during heavy rainfall months. This consistency is reflected in Approach B, where seasonal variability in annual cases corresponds to favourable environmental conditions. Approach B also shows a notable reduction in incidence, although seasonal variability influences annual cases. The high rainfall and suitable temperatures provide ideal conditions for both the aquatic stages of the mosquito life cycle and the adult mosquito populations, resulting in higher malaria incidence in this region. This is consistent with findings from climate-based vector-host models reveal that transmission is optimised when temperatures are between 21 and 25°C, and rainfall ranges from 95 to 125 mm, conditions that support mosquito breeding [2, 5, 6, 20]. Seasonal transmission patterns observed in Approach B likely reflect these environmental factors, with transmission peaking during favourable periods and waning during less optimal times. These findings emphasise the importance of sustained vector control in regions where malaria persists during cooler periods to prevent resurgence when conditions become favourable.
Increasing vector control coverage to 90% results in periods of low to no transmission, particularly in the tropical region when Approach C is adopted. This aligns with the results from vector-host models, which indicate that high levels of control can significantly reduce mosquito populations, leading to low transmission. However, as seen in Approach C, resurgences in malaria cases occur after periods of low transmission, driven by sustained asymptomatic populations and irregular spikes in infectious mosquitoes. These resurgences reflect findings from studies that show the seasonal fluctuations in larval density and mosquito breeding success, with cooler temperatures hindering mosquito activity, particularly biting, and slow parasite development [6, 11, 24, 29, 30]. The sub-tropical region, with temperatures around 15°C, is less conducive to year-round mosquito survival and transmission. Additionally, with the least rainfall, peaking below 150 mm, this region has fewer breeding sites, further contributing to lower transmission rates. However, despite these lower environmental risks, malaria transmission can still occur, particularly with asymptomatic individuals or seasonal surges in mosquito populations, which is reflected in the sporadic surges predicted in Approach C. All three approaches indicate low malaria infection rates within mosquito populations across all regions, suggesting that vector control measures are effective at reducing the parasite load in mosquitoes. However, Approach A projects a consistent adult mosquito population throughout the year, while Approach B suggests a seasonal pattern with localized peaks, and Approach C shows sporadic, localized surges driven by mismatches in temperature and rainfall seasonal patterns. Despite these differences, all approaches indicate large, sustained populations at risk of infection, emphasizing the ongoing need for comprehensive vector control strategies that account for environmental factors to address both aquatic and adult mosquito populations. These modelling approaches have shown that changes in environmental factors like temperature and rainfall can shift the suitability of an area for malaria transmission, sometimes causing outbreaks in previously unaffected regions. The seasonal growth in malaria cases under Approach B contrasts with the sporadic surges under Approach C, where the influence of rainfall and environmental fluctuations is more pronounced, illustrating the complex and dynamic nature of malaria transmission in different regions.
Although this study provides important insights into how weather conditions influence malaria transmission, it has several limitations. The modelling approaches account for specific assumptions regarding the relationship between mosquito life cycles and environmental factors. Additionally, relying on monthly mean weather data without accounting for daily fluctuations may restrict the applicability of the model to localized regions and limit its utility for short-term forecasting of transmission periods. Moreover, the models primarily emphasize weather influences while neglecting critical factors such as vector and parasite species diversity, healthcare access, and human movement, which significantly impact malaria dynamics. Future research should aim to develop more comprehensive malaria control strategies by examining the effectiveness of treatments, larval habitat management, and the impact of case importation due to travel. Incorporating these elements into the models and validating them with vector surveillance data on abundance and distribution would enhance their practical relevance for malaria prevention and control.
Conclusions
This comprehensive analysis underscores the importance of context-specific malaria control strategies and enhances our understanding of how temperature and rainfall influence various factors in malaria transmission dynamics, especially in regions with variable weather conditions. The three modeling approaches focus on distinct aspects of malaria transmission, highlighting the importance of tailoring interventions to specific local conditions, such as rainfall patterns, vector behaviour, and the impacts of climate change. The simulation results indicate that climatic fluctuations such as heat shocks or flooding can trigger surges in mosquito populations and subsequent malaria outbreaks, even with high levels of vector control. This has potentially far-reaching consequences because extreme weather events may enable the introduction of malaria into malaria-free areas or areas recently cleared of the disease. However, the findings also show that consistent and effective vector control significantly reduces malaria incidence, particularly in regions with moderate temperatures. Effective malaria control requires the integration of weather predictions into intervention plans to maintain efficacy in the face of changing environmental conditions. To accurately predict malaria transmission dynamics and inform intervention strategies, malaria control programs should prioritize selecting climate-vector modeling frameworks that are tailored to local climatic and transmission conditions. In light of these considerations, the following recommendations are proposed to address challenges and ensure sustained malaria control amidst evolving weather patterns:
The selection of models should align with the unique temperature and rainfall patterns, and the transmission dynamics of a region.
Malaria control programs should use model outputs to inform the timing, targeting, and scaling of primary malaria control measures, such as IRS, ITN distribution, and supplementary larval source management.
Combining multiple modelling approaches may provide more robust predictions in complex environments with overlapping environmental drivers.
Future models should integrate additional factors such as vector and parasite species diversity, human movement, and healthcare access to improve accuracy and utility in decision-making.
By selecting the appropriate climate-vector model, malaria control programs can ensure effective and targeted interventions, mitigating the impacts of environmental fluctuations on malaria transmission dynamics.
Declarations
Funding
This work was supported, in whole or in part, by the Bill & Melinda Gates Foundation [INV047-048]. Under the grant conditions of the Foundation, a Creative Commons Attribution 4.0 Generic License has already been assigned to the Author Accepted Manuscript version that might arise from this submission.
Competing Interests
The authors declare that they have no competing interests.
Ethics Approval
The research protocol for this study was reviewed and approved by the University of Cape Town’s Faculty of Health Sciences Human Research Ethics Committee (UCT HREC) [058/2024] and the Limpopo Provincial Department of Health [LP 2024- 05-026].
Code Availability
The code developed to generate simulation results for the model is publicly available at the Malaria Weather Models Github repository. https://github.com/maditogladstone/regionalmalariaweathermodels.git
Author Contributions
GT developed the mathematical model, implemented the code, conducted the analysis, and prepared the manuscript draft. SP contributed to the conceptual design of the study, provided guidance throughout the research process, and reviewed and provided feedback on the manuscript. BD reviewed and revised the manuscript. All authors read and approved the final manuscript.
Acknowledgements
This research was financially supported by the Bill and Melinda Gates Foundation (BMGF). We express our gratitude for the insights and resources provided by the Malaria Modelling and Analytics, Leaders in Africa (MMALA) programme. Additionally, we extend our thanks to the Limpopo Malaria Control Programme, Clinton Health Access Initiative (CHAI) South Africa and the SADC Elimination 8 (E8) Initiative Technical Working Group (TWG) for their invaluable support and collaboration throughout this study.
Footnotes
The manuscript was revised to investigate malaria incidence different climate regions rather than disease dynamics in various transmission settings. Basil D. Brooke also provided feedback and contributed in the review and revision of the manuscript. The methodology remains the same, however, analysis focuses on malaria vector population dynamics and the impact of climate variability on aquatic and adult population abundance and seasonal growth. Further analyses on the role of climate in driving transmission on host populations was conducted, with incidence and prevalence trends investigated across different climate regions.
List of Abbreviations
- WHO
- World Health Organization
- ACTs
- Artemisinin-based Combination Therapies
- IRS
- Insecticide Residual Spraying
- LLINs
- Long-Lasting Insecticide Nets
- ODEs
- Ordinary Differential Equations
- CCKP
- Climate Change Knowledge Portal