A regional-scale, high resolution dynamical malaria model that accounts for population density, climate and surface hydrology

Influence of environmental, geographic, socio-demographic, and epidemiological factors on presence of malaria at the community level in two continents

The interactions of environmental, geographic, socio-demographic, and epidemiological factors in shaping mosquito-borne disease transmission dynamics are complex and changeable, influencing the abundance and distribution of vectors and the pathogens they transmit. In this study, 27 years of cross-sectional malaria survey data (1990-2017) were used to examine the effects of these factors on Plasmodium falciparum and Plasmodium vivax malaria presence at the community level in Africa and Asia. Monthly long-term, open-source data for each factor were compiled and analyzed using generalized linear models and classification and regression trees. Both temperature and precipitation exhibited unimodal relationships with malaria, with a positive effect up to a point after which a negative effect was observed as temperature and precipitation increased. Overall decline in malaria from 2000 to 2012 was well captured by the models, as was the resurgence after that. The models also indicated higher malaria in regions with lower economic and development indicators. Malaria is driven by a combination of environmental, geographic, socioeconomic, and epidemiological factors, and in this study, we demonstrated two approaches to capturing this complexity of drivers within models. Identifying these key drivers, and describing their associations with malaria, provides key information to inform planning and prevention strategies and interventions to reduce malaria burden.

The Impact of Climatic Factors on Temporal Mosquito Distribution and Population Dynamics in an Area Targeted for Sterile Insect Technique Pilot Trials

It is widely accepted that climate affects the mosquito life history traits; however, its precise role in determining mosquito distribution and population dynamics is not fully understood. This study aimed to investigate the influence of various climatic factors on the temporal distribution of Anopheles arabiensis populations in Mamfene, South Africa between 2014 and 2019. Time series analysis, wavelet analysis, cross-correlation analysis, and regression model combined with the autoregressive integrated moving average (ARIMA) model were utilized to assess the relationship between climatic factors and An. arabiensis population density. In total 3826 adult An. arabiensis collected was used for the analysis. ARIMA (0, 1, 2) (0, 0, 1)12 models closely described the trends observed in An. arabiensis population density and distribution. The wavelet coherence and time-lagged correlation analysis showed positive correlations between An. arabiensis population density and temperature (r = 0.537 ), humidity (r = 0.495) and rainfall (r = 0.298) whilst wind showed negative correlations (r = -0.466). The regression model showed that temperature (p = 0.00119), rainfall (p = 0.0436), and humidity (p = 0.0441) as significant predictors for forecasting An. arabiensis abundance. The extended ARIMA model (AIC = 102.08) was a better fit for predicting An. arabiensis abundance compared to the basic model. Anopheles arabiensis still remains the predominant malaria vector in the study area and climate variables were found to have varying effects on the distribution and abundance of An. arabiensis. This necessitates other complementary vector control strategies such as the Sterile Insect Technique (SIT) which involves releasing sterile males into the environment to reduce mosquito populations. This requires timely mosquito and climate information to precisely target releases and enhance the effectiveness of the program, consequently reducing the malaria risk.

Estimating the effects of temperature on transmission of the human malaria parasite, Plasmodium falciparum

Despite concern that climate change could increase the human risk to malaria in certain areas, the temperature dependency of malaria transmission is poorly characterized. Here, we use a mechanistic model fitted to experimental data to describe how Plasmodium falciparum infection of the African malaria vector, Anopheles gambiae, is modulated by temperature, including its influences on parasite establishment, conversion efficiency through parasite developmental stages, parasite development rate, and overall vector competence. We use these data, together with estimates of the survival of infected blood-fed mosquitoes, to explore the theoretical influence of temperature on transmission in four locations in Kenya, considering recent conditions and future climate change. Results provide insights into factors limiting transmission in cooler environments and indicate that increases in malaria transmission due to climate warming in areas like the Kenyan Highlands, might be less than previously predicted.

The effect of explicit convection on simulated malaria transmission across Africa

Malaria transmission across sub-Saharan Africa is sensitive to rainfall and temperature. Whilst different malaria modelling techniques and climate simulations have been used to predict malaria transmission risk, most of these studies use coarse-resolution climate models. In these models convection, atmospheric vertical motion driven by instability gradients and responsible for heavy rainfall, is parameterised. Over the past decade enhanced computational capabilities have enabled the simulation of high-resolution continental-scale climates with an explicit representation of convection. In this study we use two malaria models, the Liverpool Malaria Model (LMM) and Vector-Borne Disease Community Model of the International Centre for Theoretical Physics (VECTRI), to investigate the effect of explicitly representing convection on simulated malaria transmission. The concluded impact of explicitly representing convection on simulated malaria transmission depends on the chosen malaria model and local climatic conditions. For instance, in the East African highlands, cooler temperatures when explicitly representing convection decreases LMM-predicted malaria transmission risk by approximately 55%, but has a negligible effect in VECTRI simulations. Even though explicitly representing convection improves rainfall characteristics, concluding that explicit convection improves simulated malaria transmission depends on the chosen metric and malaria model. For example, whilst we conclude improvements of 45% and 23% in root mean squared differences of the annual-mean reproduction number and entomological inoculation rate for VECTRI and the LMM respectively, bias-correcting mean climate conditions minimises these improvements. The projected impact of anthropogenic climate change on malaria incidence is also sensitive to the chosen malaria model and representation of convection. The LMM is relatively insensitive to future changes in precipitation intensity, whilst VECTRI predicts increased risk across the Sahel due to enhanced rainfall. We postulate that VECTRI's enhanced sensitivity to precipitation changes compared to the LMM is due to the inclusion of surface hydrology. Future research should continue assessing the effect of high-resolution climate modelling in impact-based forecasting.

Relative importance of VECTRI model parameters in the malaria disease transmission and prevalence

In this study, a sensitivity analysis on a VECTRI dynamical model of malaria transmission is investigated to determine the relative importance of model parameters to disease transmission and prevalence. Apart from being most climatic prone, Odisha is a highly endemic state for malaria in India. The lack in sufficient modeling studies severely impacts the malarial process studies which further hinder the possibility of malaria early warning systems and preventive measures to be undertaken beforehand. Therefore, modeling studies and investigating the relationship between malaria transmission process studies and associated climatic factors are the need of the hour. Environmental conditions have pronounced effects on the malaria transmission dynamics and abundance of the poikilothermic vectors, but the exact relationship of sensitivity for these parameters is not well established. Sensitivity analysis is a useful tool for ascertaining model responses to different input variables. Therefore, in order to perform the requisite study, a dynamical model, VECTRI, is utilized. The study period ranges from 2000 to 2013, where several sensitivity tests are performed using different model parameters such as infiltration and evaporation rate loss of ponds, degree-days for parasite development, threshold temperature for parasite development, threshold temperature for egg development in the vector, and maximum and minimum temperature for larvae survival. The experiments suggest that the lower value of minimum temperature for larvae survival (rlarv_tmin), i.e., 16 °C, provides higher vector density and entomological inoculation rate (EIR) values. EIR reaches its maximum, when the threshold temperature for parasite development (rtsporo) is 22 °C and degree-days for parasite development (dsporo) is 8 degree-days. No change is observed in the vector density; even when rtsporo is 30 °C, values of EIR are close to 0. A successive increment of infiltration and evaporation rate loss of ponds (rwaterfrac evap126) values from 130 to 200 mm/day result in approximately 5% consistent decline in vector density and EIR. The study concludes that the most sensitive parameters are dsporo, rlarv_tmin, and rwaterfrac evap126. The VECTRI model is rather insensitive to maximum temperature for larvae survival (rlarv_tmin) for vector density and EIR variables. Further certain modifications and improvements are required in VECTRI to predict out variables like vector density and EIR more accurately in highly endemic region.

Effects of Climate Variability on Malaria Transmission in Southern Côte d'Ivoire, West Africa

Malaria continues to be a major public health concern with a substantial burden in Africa. Even though it has been widely demonstrated that malaria transmission is climate-driven, there have been very few studies assessing the relationship between climate variables and malaria transmission in Côte d'Ivoire. We used the VECTRI model to predict malaria transmission in southern Côte d'Ivoire. First, we tested the suitability of VECTRI in modeling malaria transmission using ERA5 temperature data and ARC2 rainfall data. We then used the projected climatic data pertaining to 2030, 2050, and 2080 from a set of 14 simulations from the CORDEX-Africa database to compute VECTRI outputs. The entomological inoculation rate (EIR) from the VECTRI model was well correlated with the observed malaria cases from 2010 to 2019, including the peaks of malaria cases and the EIR. However, the correlation between the two parameters was not statistically significant. The VECTRI model predicted an increase in malaria transmissions in both scenarios (RCP8.5 and RCP4.5) for the time period 2030 to 2080. The monthly EIR for RCP8.5 was very high (1.74 to 1131.71 bites/person) compared to RCP4.5 (0.48 to 908 bites/person). These findings call for greater efforts to control malaria that take into account the impact of climatic factors.

Anopheles vector distribution and malaria transmission dynamics in Gbêkê region, central Côte d'Ivoire

BACKGROUND

A better understanding of vector distribution and malaria transmission dynamics at a local scale is essential for implementing and evaluating effectiveness of vector control strategies. Through the data gathered in the framework of a cluster randomized controlled trial (CRT) evaluating the In2Care (Wageningen, Netherlands) Eave Tubes strategy, the distribution of the Anopheles vector, their biting behaviour and malaria transmission dynamics were investigated in Gbêkê region, central Côte d'Ivoire.

METHODS

From May 2017 to April 2019, adult mosquitoes were collected monthly using human landing catches (HLC) in twenty villages in Gbêkê region. Mosquito species wereidentified morphologically. Monthly entomological inoculation rates (EIR) were estimated by combining the HLC data with mosquito sporozoite infection rates measured in a subset of Anopheles vectors using PCR. Finally, biting rate and EIR fluctuations were fit to local rainfall data to investigate the seasonal determinants of mosquito abundance and malaria transmission in this region.

RESULTS

Overall, Anopheles gambiae, Anopheles funestus, and Anopheles nili were the three vector complexes found infected in the Gbêkê region, but there was a variation in Anopheles vector composition between villages. Anopheles gambiae was the predominant malaria vector responsible for 84.8% of Plasmodium parasite transmission in the area. An unprotected individual living in Gbêkê region received an average of 260 [222-298], 43.5 [35.8-51.29] and 3.02 [1.96-4] infected bites per year from An. gambiae, An. funestus and An. nili, respectively. Vector abundance and malaria transmission dynamics varied significantly between seasons and the highest biting rate and EIRs occurred in the months of heavy rainfall. However, mosquitoes infected with malaria parasites remained present in the dry season, despite the low density of mosquito populations.

CONCLUSION

These results demonstrate that the intensity of malaria transmission is extremely high in Gbêkê region, especially during the rainy season. The study highlights the risk factors of transmission that could negatively impact current interventions that target indoor control, as well as the urgent need for additional vector control tools to target the population of malaria vectors in Gbêkê region and reduce the burden of the disease.

Bias-Corrected CMIP5 Projections for Climate Change and Assessments of Impact on Malaria in Senegal under the VECTRI Model

On the climate-health issue, studies have already attempted to understand the influence of climate change on the transmission of malaria. Extreme weather events such as floods, droughts, or heat waves can alter the course and distribution of malaria. This study aims to understand the impact of future climate change on malaria transmission using, for the first time in Senegal, the ICTP's community-based vector-borne disease model, TRIeste (VECTRI). This biological model is a dynamic mathematical model for the study of malaria transmission that considers the impact of climate and population variability. A new approach for VECTRI input parameters was also used. A bias correction technique, the cumulative distribution function transform (CDF-t) method, was applied to climate simulations to remove systematic biases in the Coupled Model Intercomparison Project Phase 5 (CMIP5) global climate models (GCMs) that could alter impact predictions. Beforehand, we use reference data for validation such as CPC global unified gauge-based analysis of daily precipitation (CPC for Climate Prediction Center), ERA5-land reanalysis, Climate Hazards InfraRed Precipitation with Station data (CHIRPS), and African Rainfall Climatology 2.0 (ARC2). The results were analyzed for two CMIP5 scenarios for the different time periods: assessment: 1983-2005; near future: 2006-2028; medium term: 2030-2052; and far future: 2077-2099). The validation results show that the models reproduce the annual cycle well. Except for the IPSL-CM5B model, which gives a peak in August, all the other models (ACCESS1-3, CanESM2, CSIRO, CMCC-CM, CMCC-CMS, CNRM-CM5, GFDL-CM3, GFDL-ESM2G, GFDL-ESM2M, inmcm4, and IPSL-CM5B) agree with the validation data on a maximum peak in September with a period of strong transmission in August-October. With spatial variation, the CMIP5 model simulations show more of a difference in the number of malaria cases between the south and the north. Malaria transmission is much higher in the south than in the north. However, the results predicted by the models on the occurrence of malaria by 2100 show differences between the RCP8.5 scenario, considered a high emission scenario, and the RCP4.5 scenario, considered an intermediate mitigation scenario. The CanESM2, CMCC-CM, CMCC-CMS, inmcm4, and IPSL-CM5B models predict decreases with the RCP4.5 scenario. However, ACCESS1-3, CSIRO, NRCM-CM5, GFDL-CM3, GFDL-ESM2G, and GFDL-ESM2M predict increases in malaria under all scenarios (RCP4.5 and RCP8.5). The projected decrease in malaria in the future with these models is much more visible in the RCP8.5 scenario. The results of this study are of paramount importance in the climate-health field. These results will assist in decision-making and will allow for the establishment of preventive surveillance systems for local climate-sensitive diseases, including malaria, in the targeted regions of Senegal.

The current landscape of software tools for the climate-sensitive infectious disease modelling community

Climate-sensitive infectious disease modelling is crucial for public health planning and is underpinned by a complex network of software tools. We identified only 37 tools that incorporated both climate inputs and epidemiological information to produce an output of disease risk in one package, were transparently described and validated, were named (for future searching and versioning), and were accessible (ie, the code was published during the past 10 years or was available on a repository, web platform, or other user interface). We noted disproportionate representation of developers based at North American and European institutions. Most tools (n=30 [81%]) focused on vector-borne diseases, and more than half (n=16 [53%]) of these tools focused on malaria. Few tools (n=4 [11%]) focused on food-borne, respiratory, or water-borne diseases. The under-representation of tools for estimating outbreaks of directly transmitted diseases represents a major knowledge gap. Just over half (n=20 [54%]) of the tools assessed were described as operationalised, with many freely available online.

Understanding the spatial non-stationarity in the relationships between malaria incidence and environmental risk factors using Geographically Weighted Random Forest: A case study in Rwanda

As found in the health studies literature, the levels of climate association between epidemiological diseases have been found to vary across regions. Therefore, it seems reasonable to allow for the possibility that relationships might vary spatially within regions. We implemented the geographically weighted random forest (GWRF) machine learning method to analyze ecological disease patterns caused by spatially non-stationary processes using a malaria incidence dataset for Rwanda. We first compared the geographically weighted regression (WGR), the global random forest (GRF), and the geographically weighted random forest (GWRF) to examine the spatial non-stationarity in the non-linear relationships between malaria incidence and their risk factors. We used the Gaussian areal kriging model to disaggregate the malaria incidence at the local administrative cell level to understand the relationships at a fine scale since the model goodness of fit was not satisfactory to explain malaria incidence due to the limited number of sample values. Our results show that in terms of the coefficients of determination and prediction accuracy, the geographical random forest model performs better than the GWR and the global random forest model. The coefficients of determination of the geographically weighted regression (R2), the global RF (R2), and the GWRF (R2) were 4.74, 0.76, and 0.79, respectively. The GWRF algorithm achieves the best result and reveals that risk factors (rainfall, land surface temperature, elevation, and air temperature) have a strong non-linear relationship with the spatial distribution of malaria incidence rates, which could have implications for supporting local initiatives for malaria elimination in Rwanda.

Environmental and socio-economic determinants of the occurrence of malaria clusters in Colombia

This study identifies the environmental and socio-economic determinants of clusters of high malaria incidence in Colombia during the period of 2008-2019. The malaria cases were obtained from the National System of Surveillance in Public Health, with 798,897 cases reported in the 986 Colombian municipalities evaluated during the study period. Spatial autocorrelation of incidence was examined with global and local indices. Clusters were identified in the Amazon, Pacific, and Uraba-Bajo Cauca-Alto Sinú regions. The factors associated with a municipality belonging to a high-incidence cluster were identified using a logistic regression model with mixed effects and showed a positive association for the variables (forest coverage and minimum multi-year average rainfall). An inverse relationship was observed for aqueduct coverage and the odds of belonging to a cluster. A 1% increase in forest coverage was associated with a 4.2% increase in the odds of belonging to a malaria cluster. The association with minimum multi-year average rainfall was positive (OR = 1.0011; 95% CI 1.0005-1.0027). A 1% increase in aqueduct coverage was associated with a 4.3% decrease in the odds of belonging to malaria cluster. The identification of malaria cluster determinants in Colombia could help guide surveillance and disease control policies.

Climate Drivers of Malaria Transmission Seasonality and Their Relative Importance in Sub-Saharan Africa

A new database of the Entomological Inoculation Rate (EIR) was used to directly link the risk of infectious mosquito bites to climate in Sub-Saharan Africa. Applying a statistical mixed model framework to high-quality monthly EIR measurements collected from field campaigns in Sub-Saharan Africa, we analyzed the impact of rainfall and temperature seasonality on EIR seasonality and determined important climate drivers of malaria seasonality across varied climate settings in the region. We observed that seasonal malaria transmission was within a temperature window of 15°C-40°C and was sustained if average temperature was well above 15°C or below 40°C. Monthly maximum rainfall for seasonal malaria transmission did not exceed 600 in west Central Africa, and 400 mm in the Sahel, Guinea Savannah, and East Africa. Based on a multi-regression model approach, rainfall and temperature seasonality were found to be significantly associated with malaria seasonality in all parts of Sub-Saharan Africa except in west Central Africa. Topography was found to have significant influence on which climate variable is an important determinant of malaria seasonality in East Africa. Seasonal malaria transmission onset lags behind rainfall only at markedly seasonal rainfall areas such as Sahel and East Africa; elsewhere, malaria transmission is year-round. High-quality EIR measurements can usefully supplement established metrics for seasonal malaria. The study's outcome is important for the improvement and validation of weather-driven dynamical mathematical malaria models that directly simulate EIR. Our results can contribute to the development of fit-for-purpose weather-driven malaria models to support health decision-making in the fight to control or eliminate malaria in Sub-Saharan Africa.

Predicting Plasmodium knowlesi transmission risk across Peninsular Malaysia using machine learning-based ecological niche modeling approaches

The emergence of potentially life-threatening zoonotic malaria caused by Plasmodium knowlesi nearly two decades ago has continued to challenge Malaysia healthcare. With a total of 376 P. knowlesi infections notified in 2008, the number increased to 2,609 cases in 2020 nationwide. Numerous studies have been conducted in Malaysian Borneo to determine the association between environmental factors and knowlesi malaria transmission. However, there is still a lack of understanding of the environmental influence on knowlesi malaria transmission in Peninsular Malaysia. Therefore, our study aimed to investigate the ecological distribution of human P. knowlesi malaria in relation to environmental factors in Peninsular Malaysia. A total of 2,873 records of human P. knowlesi infections in Peninsular Malaysia from 1st January 2011 to 31st December 2019 were collated from the Ministry of Health Malaysia and geolocated. Three machine learning-based models, maximum entropy (MaxEnt), extreme gradient boosting (XGBoost), and ensemble modeling approach, were applied to predict the spatial variation of P. knowlesi disease risk. Multiple environmental parameters including climate factors, landscape characteristics, and anthropogenic factors were included as predictors in both predictive models. Subsequently, an ensemble model was developed based on the output of both MaxEnt and XGBoost. Comparison between models indicated that the XGBoost has higher performance as compared to MaxEnt and ensemble model, with AUC_ROC values of 0.933 ± 0.002 and 0.854 ± 0.007 for train and test datasets, respectively. Key environmental covariates affecting human P. knowlesi occurrence were distance to the coastline, elevation, tree cover, annual precipitation, tree loss, and distance to the forest. Our models indicated that the disease risk areas were mainly distributed in low elevation (75-345 m above mean sea level) areas along the Titiwangsa mountain range and inland central-northern region of Peninsular Malaysia. The high-resolution risk map of human knowlesi malaria constructed in this study can be further utilized for multi-pronged interventions targeting community at-risk, macaque populations, and mosquito vectors.

Malaria metrics distribution under global warming: assessment of the VECTRI malaria model over Cameroon

Malaria is a critical health issue across the world and especially in Africa. Studies based on dynamical models helped to understand inter-linkages between this illness and climate. In this study, we evaluated the ability of the VECTRI community vector malaria model to simulate the spread of malaria in Cameroon using rainfall and temperature data from FEWS-ARC2 and ERA-interim, respectively. In addition, we simulated the model using five results of the dynamical downscaling of the regional climate model RCA4 within two time frames named near future (2035-2065) and far future (2071-2100), aiming to explore the potential effects of global warming on the malaria propagation over Cameroon. The evaluated metrics include the risk maps of the entomological inoculation rate (EIR) and the parasite ratio (PR). During the historical period (1985-2005), the model satisfactorily reproduces the observed PR and EIR. Results of projections reveal that under global warming, heterogeneous changes feature the study area, with localized increases or decreases in PR and EIR. As the level of radiative forcing increases (from 2.6 to 8.5 W.m^-2), the magnitude of change in PR and EIR also gradually intensifies. The occurrence of transmission peaks is projected in the temperature range of 26-28 °C. Moreover, PR and EIR vary depending on the three agro-climatic regions of the study area. VECTRI still needs to integrate other aspects of disease transmission, such as population mobility and intervention strategies, in order to be more relevant to support actions of decision-makers and policy makers.

Distribution and Risk Factors of Malaria in the Greater Accra Region in Ghana

Malaria remains a serious public health challenge in Ghana including the Greater Accra Region. This study aimed to quantify the spatial, temporal and spatio-temporal patterns of malaria in the Greater Accra Region to inform targeted allocation of health resources. Malaria cases data from 2015 to 2019 were obtained from the Ghanaian District Health Information and Management System and aggregated at a district and monthly level. Spatial analysis was conducted using the Global Moran's I, Getis-Ord Gi*, and local indicators of spatial autocorrelation. Kulldorff's space-time scan statistics were used to investigate space-time clustering. A negative binomial regression was used to find correlations between climatic factors and sociodemographic characteristics and the incidence of malaria. A total of 1,105,370 malaria cases were reported between 2015 and 2019. Significant seasonal variation was observed, with June and July being the peak months of reported malaria cases. The hotspots districts were Kpone-Katamanso Municipal District, Ashaiman Municipal Districts, Tema Municipal District, and La-Nkwantanang-Madina Municipal District. While La-Nkwantanang-Madina Municipal District was high-high cluster. The Spatio-temporal clusters occurred between February 2015 and July 2017 in the districts of Ningo-Prampram, Shai-Osudoku, Ashaiman Municipal, and Kpone-Katamanso Municipal with a radius of 26.63 km and an relative risk of 4.66 (p < 0.001). Malaria cases were positively associated with monthly rainfall (adjusted odds ratio [AOR] = 1.01; 95% confidence interval [CI] = 1.005, 1.016) and the previous month's cases (AOR = 1.064; 95% CI 1.062, 1.065) and negatively correlated with minimum temperature (AOR = 0.86, 95% CI = 0.823, 0.899) and population density (AOR = 0.996, 95% CI = 0.994, 0.998). Malaria control and prevention should be strengthened in hotspot districts in the appropriate months to improve program effectiveness.

Spatial analysis to evaluate risk of malaria in Northern Sumatera, Indonesia

Background

As Indonesia aims for malaria elimination by 2030, provisional malaria epidemiology and risk factors evaluation are important in pursue of this national goal. Therefore, this study aimed to understand the risk factor of malaria in Northern Sumatera.

Methods

Malaria cases from 2019 to 2020 were obtained from the Indonesian Ministry of Health Electronic Database. Climatic variables were provided by the Center for Meteorology and Geophysics Medan branch office. Multivariable logistic regression was undertaken to understand the risk factors of imported malaria. A zero-inflated Poisson multivariable regression model was used to study the climatic drivers of indigenous malaria.

Results

A total of 2208 (indigenous: 76.0% [1679] and imported: 17.8% [392]) were reported during the study period. Risk factors of imported malaria were: ages 19–30 (adjusted odds ratio [AOR] = 3.31; 95% confidence interval [CI] 1.67, 2.56), 31–45 (AOR = 5.69; 95% CI 2.65, 12.20), and > 45 years (AOR = 5.11; 95% CI 2.41, 10.84). Military personnel and forest workers and miners were 1,154 times (AOR = 197.03; 95% CI 145.93, 9,131.56) and 44 times (AOR = 44.16; 95% CI 4.08, 477,93) more likely to be imported cases as compared to those working as employees and traders. Indigenous Plasmodium falciparum increased by 12.1% (95% CrI 5.1%, 20.1%) for 1% increase in relative humidity and by 21.0% (95% CrI 9.0%, 36.2%) for 1 °C increase in maximum temperature. Plasmodium vivax decreased by 0.8% (95% CrI 0.2%, 1.3%) and 16.7% (95% CrI 13.7%, 19.9%) for one meter and 1 °C increase of altitude and minimum temperature. Indigenous hotspot was reported by Kota Tanjung Balai city and Asahan regency, respectively. Imported malaria hotspots were reported in Batu Bara, Kota Tebing Tinggi, Serdang Bedagai and Simalungun.

Conclusion

Both indigenous and imported malaria is limited to a few regencies and cities in Northern Sumatera. The control measures should focus on these risk factors to achieve elimination in Indonesia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04262-y.

Spatial patterns and climate drivers of malaria in three border areas of Brazil, Venezuela and Guyana, 2016-2018

In 2020, 77% of malaria cases in the Americas were concentrated in Venezuela, Brazil, and Colombia. These countries are characterized by a heterogeneous malaria landscape and malaria hotspots. Furthermore, the political unrest in Venezuela has led to significant cross-border population movement. Hence, the aim of this study was to describe spatial patterns and identify significant climatic drivers of malaria transmission along the Venezuela-Brazil-Guyana border, focusing on Bolivar state, Venezuela and Roraima state, Brazil. Malaria case data, stratified by species from 2016 to 2018, were obtained from the Brazilian Malaria Epidemiology Surveillance Information System, the Guyana Vector Borne Diseases Program, the Venezuelan Ministry of Health, and civil society organizations. Spatial autocorrelation in malaria incidence was explored using Getis-Ord (Gi*) statistics. A Poisson regression model was developed with a conditional autoregressive prior structure and posterior parameters were estimated using the Bayesian Markov chain Monte Carlo simulation with Gibbs sampling. There were 685,498 malaria cases during the study period. Plasmodium vivax was the predominant species (71.7%, 490,861). Malaria hotspots were located in eight municipalities along the Venezuela and Guyana international borders with Brazil. Plasmodium falciparum increased by 2.6% (95% credible interval [CrI] 2.1%, 2.8%) for one meter increase in altitude, decreased by 1.6% (95% CrI 1.5%, 2.3%) and 0.9% (95% CrI 0.7%, 2.4%) per 1 cm increase in 6-month lagged precipitation and each 1 °C increase of minimum temperature without lag. Each 1 °C increase of 1-month lagged maximum temperature increased P. falciparum by 0.6% (95% CrI 0.4%, 1.9%). P. vivax cases increased by 1.5% (95% CrI 1.3%, 1.6%) for one meter increase in altitude and decreased by  1.1% (95% CrI 1.0%, 1.2%) and 7.3% (95% CrI 6.7%, 9.7%) for each 1 cm increase of precipitation lagged at 6-months and 1 °C increase in minimum temperature lagged at 6-months. Each 1°C increase of two-month lagged maximum temperature increased P. vivax by 1.5% (95% CrI 0.6%, 7.1%). There was no significant residual spatial clustering after accounting for climatic covariates. Malaria hotspots were located along the Venezuela and Guyana international border with Roraima state, Brazil. In addition to population movement, climatic variables were important drivers of malaria transmission in these areas.

Inference and dynamic simulation of malaria using a simple climate-driven entomological model of malaria transmission

Given the crucial role of climate in malaria transmission, many mechanistic models of malaria represent vector biology and the parasite lifecycle as functions of climate variables in order to accurately capture malaria transmission dynamics. Lower dimension mechanistic models that utilize implicit vector dynamics have relied on indirect climate modulation of transmission processes, which compromises investigation of the ecological role played by climate in malaria transmission. In this study, we develop an implicit process-based malaria model with direct climate-mediated modulation of transmission pressure borne through the Entomological Inoculation Rate (EIR). The EIR, a measure of the number of infectious bites per person per unit time, includes the effects of vector dynamics, resulting from mosquito development, survivorship, feeding activity and parasite development, all of which are moderated by climate. We combine this EIR-model framework, which is driven by rainfall and temperature, with Bayesian inference methods, and evaluate the model’s ability to simulate local transmission across 42 regions in Rwanda over four years. Our findings indicate that the biologically-motivated, EIR-model framework is capable of accurately simulating seasonal malaria dynamics and capturing of some of the inter-annual variation in malaria incidence. However, the model unsurprisingly failed to reproduce large declines in malaria transmission during 2018 and 2019 due to elevated anti-malaria measures, which were not accounted for in the model structure. The climate-driven transmission model also captured regional variation in malaria incidence across Rwanda’s diverse climate, while identifying key entomological and epidemiological parameters important to seasonal malaria dynamics. In general, this new model construct advances the capabilities of implicitly-forced lower dimension dynamical malaria models by leveraging climate drivers of malaria ecology and transmission.

Climate plays a fundamental and complex role in malaria transmission, by acting on multiple aspects of mosquito ecology and parasite transmissibility. However, to express malaria transmission pressure, malaria models with implicit vector dynamics have relied on indirect predictors of vector ecology, such as temporal seasonality or interpolations of rainfall/temperature, instead of entomological processes directly informed by ambient conditions. This approach obscures the specific influence of environmental conditions on relevant vector and parasite ecology, as well as meaningful interpretation of climate variability within these models. Here, we demonstrate that both interpretability and ecological effect from climate can be instantiated in lower dimension dynamical models through representation of transmission pressures via a climate-driven Entomological Inoculation Rate (EIR). This process-based model framework is driven by local rainfall and temperature, which regulate multiple aspects of the EIR, namely mosquito density, host-seeking activity, and parasite infectivity. Our results indicate that the climate-driven model construct is able to reproduce regional and local malaria transmission at seasonal and inter-annual time scales, while enabling identification of key entomological determinants of transmission.

Potential future malaria transmission in Odisha due to climate change

Future projections of malaria transmission is made for Odisha, a highly endemic region of India, through numerical simulations using the VECTRI dynamical model. The model is forced with bias-corrected temperature and rainfall from a global climate model (CCSM4) for the baseline period 1975–2005 and for the projection periods 2020s, 2050s, and 2080s under RCP8.5 emission scenario. The temperature, rainfall, mosquito density and entomological inoculation rate (EIR), generated from the VECTRI model are evaluated with the observation and analyzed further to estimate the future malaria transmission over Odisha on a spatio-temporal scale owing to climate change. Our results reveal that the malaria transmission in Odisha as a whole during summer and winter monsoon seasons may decrease in future due to the climate change except in few districts with the high elevations and dense forest regions such as Kandhamal, Koraput, Raygada and Mayurbhanj districts where an increase in malaria transmission is found. Compared to the baseline period, mosquito density shows decrease in most districts of the south, southwest, central, north and northwest regions of Odisha in 2030s, 2050s and 2080s. An overall decrease in malaria transmission of 20–40% (reduction in EIR) is seen during the monsoon season (June-Sept) over Odisha with the increased surface temperature of 3.5–4 °C and with the increased rainfall of 20–35% by the end of the century with respect to the baseline period. Furthermore, malaria transmission is likely to reduce in future over most of the Odisha regions with the increase in future warm and cold nights temperatures.

Assessment of Climate-Driven Variations in Malaria Transmission in Senegal Using the VECTRI Model

Several vector-borne diseases, such as malaria, are sensitive to climate and weather conditions. When unusual conditions prevail, for example, during periods of heavy rainfall, mosquito populations can multiply and trigger epidemics. This study, which consists of better understanding the link between malaria transmission and climate factors at a national level, aims to validate the VECTRI model (VECtor borne disease community model of ICTP, TRIeste) in Senegal. The VECTRI model is a grid-distributed dynamical model that couples a biological model for the vector and parasite life cycles to a simple compartmental Susceptible-Exposed-Infectious-Recovered (SEIR) representation of the disease progression in the human host. In this study, a VECTRI model driven by reanalysis data (ERA-5) was used to simulate malaria parameters, such as the entomological inoculation rate (EIR) in Senegal. Observed malaria data from the National Malaria Control Program in Senegal (PNLP/Programme National de Lutte contre le Paludisme au Senegal) and outputs from the climate data used in this study were compared. The findings highlight the unimodal shape of temporal malaria occurrence, and the seasonal malaria transmission contrast is closely linked to the latitudinal variation of the rainfall, showing a south–north gradient over Senegal. This study showed that the peak of malaria takes place from September to October, with a lag of about one month from the peak of rainfall in Senegal. There is an agreement between observations and simulations about decreasing malaria cases on time. These results indicate that the southern area of Senegal is at the highest risk of malaria spread outbreaks. The findings in the paper are expected to guide community-based early-warning systems and adaptation strategies in Senegal, which will feed into the national malaria prevention, response, and care strategies adapted to the needs of local communities.

Larval flushing alters malaria endemicity patterns in regions with similar habitat abundance

A model of Anopheles gambiae populations dynamics coupled with Plasmodium falciparum transmission dynamics is extended to include mechanisms of larval flushing which are known to occur. Flushing dynamics are modeled using a simulation that incorporates seasonal, autocorrelated, and random components based on 30 years of rainfall data for the Kakamega District of the western Kenya highlands. The model demonstrates that flushing phenomena can account for differences between regions with the same annual larval habitat pattern, changing the World Health Organization endemicity classification from either hyperendemic or holoendemic to hypoendemic disease patterns. Mesoendemic patterns of infection occur at the boundary of the holoendemic to hypoendemic transition. For some levels of flushing the entomological inoculation rate drops to an insignificant amount and disease disappears, while the annual indoor resting density remains well above zero. In these scenarios, the disease is hypoendemic, yet the model shows that outbreaks can occur when disease is introduced at particular time points.

Understanding the effect of climate change in the distribution and intensity of malaria transmission over India using a dynamical malaria model

Efforts have been made to quantify the spatio-temporal malaria transmission intensity over India using the dynamical malaria model, namely, Vector-borne Disease Community Model of International Centre for Theoretical Physics Trieste (VECTRI). The likely effect of climate change in the variability of malaria transmission intensity over different parts of India is also investigated. The Historical data and future projection scenarios of the rainfall and temperature derived from the Coupled Model Intercomparison Project Phase 5 (CMIP5) model output are used for this purpose. The Entomological Inoculation Rate (EIR) and Vector are taken as quantifiers of malaria transmission intensity. It is shown that the maximum number of malaria cases over India occur during the Sept-Oct months, whereas the minimum during the Feb-Apr months. The malaria transmission intensity as well as length of transmission season over India is likely to increase in the future climate as a result of global warming.

Projecting the risk of mosquito-borne diseases in a warmer and more populated world: a multi-model, multi-scenario intercomparison modelling study

BACKGROUND

Mosquito-borne diseases are expanding their range, and re-emerging in areas where they had subsided for decades. The extent to which climate change influences the transmission suitability and population at risk of mosquito-borne diseases across different altitudes and population densities has not been investigated. The aim of this study was to quantify the extent to which climate change will influence the length of the transmission season and estimate the population at risk of mosquito-borne diseases in the future, given different population densities across an altitudinal gradient.

METHODS

Using a multi-model multi-scenario framework, we estimated changes in the length of the transmission season and global population at risk of malaria and dengue for different altitudes and population densities for the period 1951-99. We generated projections from six mosquito-borne disease models, driven by four global circulation models, using four representative concentration pathways, and three shared socioeconomic pathways.

FINDINGS

We show that malaria suitability will increase by 1·6 additional months (mean 0·5, SE 0·03) in tropical highlands in the African region, the Eastern Mediterranean region, and the region of the Americas. Dengue suitability will increase in lowlands in the Western Pacific region and the Eastern Mediterranean region by 4·0 additional months (mean 1·7, SE 0·2). Increases in the climatic suitability of both diseases will be greater in rural areas than in urban areas. The epidemic belt for both diseases will expand towards temperate areas. The population at risk of both diseases might increase by up to 4·7 additional billion people by 2070 relative to 1970-99, particularly in lowlands and urban areas.

INTERPRETATION

Rising global mean temperature will increase the climatic suitability of both diseases particularly in already endemic areas. The predicted expansion towards higher altitudes and temperate regions suggests that outbreaks can occur in areas where people might be immunologically naive and public health systems unprepared. The population at risk of malaria and dengue will be higher in densely populated urban areas in the WHO African region, South-East Asia region, and the region of the Americas, although we did not account for urban-heat island effects, which can further alter the risk of disease transmission.

FUNDING

UK Space Agency, Royal Society, UK National Institute for Health Research, and Swedish Research Council.

Impact of an accelerated melting of Greenland on malaria distribution over Africa

Studies about the impact of future climate change on diseases have mostly focused on standard Representative Concentration Pathway climate change scenarios. These scenarios do not account for the non-linear dynamics of the climate system. A rapid ice-sheet melting could occur, impacting climate and consequently societies. Here, we investigate the additional impact of a rapid ice-sheet melting of Greenland on climate and malaria transmission in Africa using several malaria models driven by Institute Pierre Simon Laplace climate simulations. Results reveal that our melting scenario could moderate the simulated increase in malaria risk over East Africa, due to cooling and drying effects, cause a largest decrease in malaria transmission risk over West Africa and drive malaria emergence in southern Africa associated with a significant southward shift of the African rain-belt. We argue that the effect of such ice-sheet melting should be investigated further in future public health and agriculture climate change risk assessments.

Climate-proofing a malaria eradication strategy

Two recent initiatives, the World Health Organization (WHO) Strategic Advisory Group on Malaria Eradication and the Lancet Commission on Malaria Eradication, have assessed the feasibility of achieving global malaria eradication and proposed strategies to achieve it. Both reports rely on a climate-driven model of malaria transmission to conclude that long-term trends in climate will assist eradication efforts overall and, consequently, neither prioritize strategies to manage the effects of climate variability and change on malaria programming. This review discusses the pathways via which climate affects malaria and reviews the suitability of climate-driven models of malaria transmission to inform long-term strategies such as an eradication programme. Climate can influence malaria directly, through transmission dynamics, or indirectly, through myriad pathways including the many socioeconomic factors that underpin malaria risk. These indirect effects are largely unpredictable and so are not included in climate-driven disease models. Such models have been effective at predicting transmission from weeks to months ahead. However, due to several well-documented limitations, climate projections cannot accurately predict the medium- or long-term effects of climate change on malaria, especially on local scales. Long-term climate trends are shifting disease patterns, but climate shocks (extreme weather and climate events) and variability from sub-seasonal to decadal timeframes have a much greater influence than trends and are also more easily integrated into control programmes. In light of these conclusions, a pragmatic approach is proposed to assessing and managing the effects of climate variability and change on long-term malaria risk and on programmes to control, eliminate and ultimately eradicate the disease. A range of practical measures are proposed to climate-proof a malaria eradication strategy, which can be implemented today and will ensure that climate variability and change do not derail progress towards eradication.

Past eight-year malaria data in Gedeo zone, southern Ethiopia: trend, reporting-quality, spatiotemporal distribution, and association with socio-demographic and meteorological variables

Background

Informed decision making is underlined by all tiers in the health system. Poor data record system coupled with under- (over)-reporting of malaria cases affects the country’s malaria elimination activities. Thus, malaria data at health facilities and health offices are important particularly to monitor and evaluate the elimination progresses. This study was intended to assess overall reported malaria cases, reporting quality, spatiotemporal trends and factors associated in Gedeo zone, South Ethiopia.

Methods

Past 8 years retrospective data stored in 17 health centers and 5 district health offices in Gedeo Zone, South Ethiopia were extracted. Malaria cases data at each health center with sociodemographic information, between January 2012 and December 2019, were included. Meteorological data were obtained from the national meteorology agency of Ethiopia. The data were analyzed using Stata 13.

Results

A total of 485,414 suspected cases were examined for malaria during the previous 8 years at health centers. Of these suspects, 57,228 (11.79%) were confirmed malaria cases with an overall decline during the 8-year period. We noted that 3758 suspected cases and 467 confirmed malaria cases were not captured at the health offices. Based on the health centers records, the proportions of Plasmodium falciparum (49.74%) and P. vivax (47.59%) infection were nearly equivalent (p = 0.795). The former was higher at low altitudes while the latter was higher at higher altitudes. The over 15 years of age group accounted for 11.47% of confirmed malaria cases (p < 0.001). There was high spatiotemporal variation: the highest case record was during Belg (12.52%) and in Dilla town (18,150, 13.17%, p < 0.001) which is located at low altitude. Monthly rainfall and minimum temperature exhibited strong associations with confirmed malaria cases.

Conclusion

A notable overall decline in malaria cases was observed during the eight-year period. Both P. falciparum and P. vivax were found at equivalent endemicity level; hence control measures should continue targeting both species. The noticed under reporting, the high malaria burden in urban settings, low altitudes and Belg season need spatiotemporal consideration by the elimination program.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-05783-8.

Estimating the malaria transmission over the Indian subcontinent in a warming environment using a dynamical malaria model

Malaria is a major public health problem in India. The malaria transmission is sensitive to climatic parameters. The regional population-related factors also influence malaria transmission. To take into account temperature and rainfall variability and associated population-related effects (in a changing climate) on the malaria transmission over India, a regional dynamical malaria model, namely VECTRI (vector-borne disease community model) is used. The daily temperature and rainfall data derived from the historical (years 1961-2005) and representative concentration pathway (years 2006-2050) runs of the Coupled Model Intercomparison Project Phase 5 models have been used for the analysis. The model results of the historical run are compared with the observational data. The spatio-temporal changes (region-specific as well as seasonal changes) in the malaria transmission as a result of climate change are quantified over the India. The parameters related to the breeding cycle of malaria as well as those which estimate the malaria cases are analyzed in the global warming scenario.

Climatological, virological and sociological drivers of current and projected dengue fever outbreak dynamics in Sri Lanka

The largest ever Sri Lankan dengue outbreak of 2017 provides an opportunity for investigating the relative contributions of climatological, epidemiological and sociological drivers on the epidemic patterns of this clinically important vector-borne disease. To do so, we develop a climatologically driven disease transmission framework for dengue virus using spatially resolved temperature and precipitation data as well as the time-series susceptible-infected-recovered (SIR) model. From this framework, we first demonstrate that the distinct climatological patterns encountered across the island play an important role in establishing the typical yearly temporal dynamics of dengue, but alone are unable to account for the epidemic case numbers observed in Sri Lanka during 2017. Using a simplified two-strain SIR model, we demonstrate that the re-introduction of a dengue virus serotype that had been largely absent from the island in previous years may have played an important role in driving the epidemic, and provide a discussion of the possible roles for extreme weather events and human mobility patterns on the outbreak dynamics. Lastly, we provide estimates for the future burden of dengue across Sri Lanka using the Coupled Model Intercomparison Phase 5 climate projections. Critically, we demonstrate that climatological and serological factors can act synergistically to yield greater projected case numbers than would be expected from the presence of a single driver alone. Altogether, this work provides a holistic framework for teasing apart and analysing the various complex drivers of vector-borne disease outbreak dynamics.

A spatio-temporal analysis to identify the drivers of malaria transmission in Bhutan

At a time when Bhutan is on the verge of malaria elimination, the aim of this study was to identify malaria clusters at high geographical resolution and to determine its association with local environmental characteristics. Malaria cases from 2006–2014 were obtained from the Vector-borne Disease Control Program under the Ministry of Health, Bhutan. A Zero-Inflated Poisson multivariable regression model with a conditional autoregressive (CAR) prior structure was developed. Bayesian Markov chain Monte Carlo (MCMC) simulation with Gibbs sampling was used to estimate posterior parameters. A total of 2,062 Plasmodium falciparum and 2,284 Plasmodium vivax cases were reported during the study period. Both species of malaria showed seasonal peaks with decreasing trend. Gender and age were not associated with the transmission of either species of malaria. P. falciparum increased by 0.7% (95% CrI: 0.3%, 0.9%) for a one mm increase in rainfall, while climatic variables (temperature and rainfall) were not associated with P. vivax. Insecticide treated bed net use and residual indoor insecticide coverage were unaccounted for in this study. Hot spots and clusters of both species were isolated in the central southern part of Bhutan bordering India. There was significant residual spatial clustering after accounting for climate and demographic variables.

Modelled and observed mean and seasonal relationships between climate, population density and malaria indicators in Cameroon

Background

A major health burden in Cameroon is malaria, a disease that is sensitive to climate, environment and socio-economic conditions, but whose precise relationship with these drivers is still uncertain. An improved understanding of the relationship between the disease and its drivers, and the ability to represent these relationships in dynamic disease models, would allow such models to contribute to health mitigation and adaptation planning. This work collects surveys of malaria parasite ratio and entomological inoculation rate and examines their relationship with temperature, rainfall, population density in Cameroon and uses this analysis to evaluate a climate sensitive mathematical model of malaria transmission.

Methods

Co-located, climate and population data is compared to the results of 103 surveys of parasite ratio (PR) covering 18,011 people in Cameroon. A limited set of campaigns which collected year-long field-surveys of the entomological inoculation rate (EIR) are examined to determine the seasonality of disease transmission, three of the study locations are close to the Sanaga and Mefou rivers while others are not close to any permanent water feature. Climate-driven simulations of the VECTRI malaria model are evaluated with this analysis.

Results

The analysis of the model results shows the PR peaking at temperatures of approximately 22 °C to 26 °C, in line with recent work that has suggested a cooler peak temperature relative to the established literature, and at precipitation rates at 7 mm day⁻¹, somewhat higher than earlier estimates. The malaria model is able to reproduce this broad behaviour, although the peak occurs at slightly higher temperatures than observed, while the PR peaks at a much lower rainfall rate of 2 mm day⁻¹. Transmission tends to be high in rural and peri-urban relative to urban centres in both model and observations, although the model is oversensitive to population which could be due to the neglect of population movements, and differences in hydrological conditions, housing quality and access to healthcare. The EIR follows the seasonal rainfall with a lag of 1 to 2 months, and is well reproduced by the model, while in three locations near permanent rivers the annual cycle of malaria transmission is out of phase with rainfall and the model fails.

Conclusion

Malaria prevalence is maximum at temperatures of 24 to 26 °C in Cameroon and rainfall rates of approximately 4 to 6 mm day⁻¹. The broad relationships are reproduced in a malaria model although prevalence is highest at a lower rainfall maximum of 2 mm day⁻¹. In locations far from water bodies malaria transmission seasonality closely follows that of rainfall with a lag of 1 to 2 months, also reproduced by the model, but in locations close to a seasonal river the seasonality of malaria transmission is reversed due to pooling in the transmission to the dry season, which the model fails to capture.

Biological Adaptations Associated with Dehydration in Mosquitoes

Diseases that are transmitted by mosquitoes are a tremendous health and socioeconomic burden with hundreds of millions of people being impacted by mosquito-borne illnesses annually. Many factors have been implicated and extensively studied in disease transmission dynamics, but knowledge regarding how dehydration impacts mosquito physiology, behavior, and resulting mosquito-borne disease transmission remain underdeveloped. The lapse in understanding on how mosquitoes respond to dehydration stress likely obscures our ability to effectively study mosquito physiology, behavior, and vectorial capabilities. The goal of this review is to develop a profile of factors underlying mosquito biology that are altered by dehydration and the implications that are related to disease transmission.

Amazon deforestation drives malaria transmission, and malaria burden reduces forest clearing

Deforestation and land use change are among the most pressing anthropogenic environmental impacts. In Brazil, a resurgence of malaria in recent decades paralleled rapid deforestation and settlement in the Amazon basin, yet evidence of a deforestation-driven increase in malaria remains equivocal. We hypothesize an underlying cause of this ambiguity is that deforestation and malaria influence each other in bidirectional causal relationships-deforestation increases malaria through ecological mechanisms and malaria reduces deforestation through socioeconomic mechanisms-and that the strength of these relationships depends on the stage of land use transformation. We test these hypotheses with a large geospatial dataset encompassing 795 municipalities across 13 y (2003 to 2015) and show deforestation has a strong positive effect on malaria incidence. Our results suggest a 10% increase in deforestation leads to a 3.3% increase in malaria incidence (∼9,980 additional cases associated with 1,567 additional km² lost in 2008, the study midpoint, Amazon-wide). The effect is larger in the interior and absent in outer Amazonian states where little forest remains. However, this strong effect is only detectable after controlling for a feedback of malaria burden on forest loss, whereby increased malaria burden significantly reduces forest clearing, possibly mediated by human behavior or economic development. We estimate a 1% increase in malaria incidence results in a 1.4% decrease in forest area cleared (∼219 fewer km² cleared associated with 3,024 additional cases in 2008). This bidirectional socioecological feedback between deforestation and malaria, which attenuates as land use intensifies, illustrates the intimate ties between environmental change and human health.

Thermal biology of mosquito-borne disease

Mosquito-borne diseases cause a major burden of disease worldwide. The vital rates of these ectothermic vectors and parasites respond strongly and nonlinearly to temperature and therefore to climate change. Here, we review how trait-based approaches can synthesise and mechanistically predict the temperature dependence of transmission across vectors, pathogens, and environments. We present 11 pathogens transmitted by 15 different mosquito species - including globally important diseases like malaria, dengue, and Zika - synthesised from previously published studies. Transmission varied strongly and unimodally with temperature, peaking at 23-29ºC and declining to zero below 9-23ºC and above 32-38ºC. Different traits restricted transmission at low versus high temperatures, and temperature effects on transmission varied by both mosquito and parasite species. Temperate pathogens exhibit broader thermal ranges and cooler thermal minima and optima than tropical pathogens. Among tropical pathogens, malaria and Ross River virus had lower thermal optima (25-26ºC) while dengue and Zika viruses had the highest (29ºC) thermal optima. We expect warming to increase transmission below thermal optima but decrease transmission above optima. Key directions for future work include linking mechanistic models to field transmission, combining temperature effects with control measures, incorporating trait variation and temperature variation, and investigating climate adaptation and migration.

Numerical Modeling of the Dynamics of Malaria Transmission in a Highly Endemic Region of India

Using a dynamical model (VECTRI) for malaria transmission that accounts for the influence of population and climatic conditions, malaria transmission dynamics is investigated for a highly endemic region (state of Odisha) in India. The model is first calibrated over the region, and subsequently numerical simulations are carried out for the period 2000–2013. Using both model and observations we find that temperature, adult mosquito population, and infective biting rates have increased over this period, and the malaria vector abundance is higher during the summer monsoon season. Regionally, the intensity of malaria transmission is found to be higher in the north, central and southern districts of Odisha where the mosquito populations and the number of infective bites are more and mainly in the forest or mountainous ecotypes. We also find that the peak of the malaria transmission occurs when the monthly mean temperature is in the range of ~28–29 °C, and monthly rainfall accumulation in the range of ~200–360 mm.

A Dynamical and Zero-Inflated Negative Binomial Regression Modelling of Malaria Incidence in Limpopo Province, South Africa

Recent studies have considered the connections between malaria incidence and climate variables using mathematical and statistical models. Some of the statistical models focused on time series approach based on Box-Jenkins methodology or on dynamic model. The latter approach allows for covariates different from its original lagged values, while the Box-Jenkins does not. In real situations, malaria incidence counts may turn up with many zero terms in the time series. Fitting time series model based on the Box-Jenkins approach and ARIMA may be spurious. In this study, a zero-inflated negative binomial regression model was formulated for fitting malaria incidence in Mopani and Vhembe-two of the epidemic district municipalities in Limpopo, South Africa. In particular, a zero-inflated negative binomial regression model was formulated for daily malaria counts as a function of some climate variables, with the aim of identifying the model that best predicts reported malaria cases. Results from this study show that daily rainfall amount and the average temperature at various lags have a significant influence on malaria incidence in the study areas. The significance of zero inflation on the malaria count was examined using the Vuong test and the result shows that zero-inflated negative binomial regression model fits the data better. A dynamical climate-based model was further used to investigate the population dynamics of mosquitoes over the two regions. Findings highlight the significant roles of Anopheles arabiensis on malaria transmission over the regions and suggest that vector control activities should be intense to eradicate malaria in Mopani and Vhembe districts. Although An. arabiensis has been identified as the major vector over these regions, our findings further suggest the presence of additional vectors transmitting malaria in the study regions. The findings from this study offer insight into climate-malaria incidence linkages over Limpopo province of South Africa.

Forecasting Zoonotic Infectious Disease Response to Climate Change: Mosquito Vectors and a Changing Environment

Infectious diseases are changing due to the environment and altered interactions among hosts, reservoirs, vectors, and pathogens. This is particularly true for zoonotic diseases that infect humans, agricultural animals, and wildlife. Within the subset of zoonoses, vector-borne pathogens are changing more rapidly with climate change, and have a complex epidemiology, which may allow them to take advantage of a changing environment. Most mosquito-borne infectious diseases are transmitted by mosquitoes in three genera: Aedes, Anopheles, and Culex, and the expansion of these genera is well documented. There is an urgent need to study vector-borne diseases in response to climate change and to produce a generalizable approach capable of generating risk maps and forecasting outbreaks. Here, we provide a strategy for coupling climate and epidemiological models for zoonotic infectious diseases. We discuss the complexity and challenges of data and model fusion, baseline requirements for data, and animal and human population movement. Disease forecasting needs significant investment to build the infrastructure necessary to collect data about the environment, vectors, and hosts at all spatial and temporal resolutions. These investments can contribute to building a modeling community around the globe to support public health officials so as to reduce disease burden through forecasts with quantified uncertainty.

Dynamical Malaria Forecasts Are Skillful at Regional and Local Scales in Uganda up to 4 Months Ahead

Malaria forecasts from dynamical systems have never been attempted at the health district or local clinic catchment scale, and so their usefulness for public health preparedness and response at the local level is fundamentally unknown. A pilot preoperational forecasting system is introduced in which the European Centre for Medium Range Weather Forecasts ensemble prediction system and seasonal climate forecasts of temperature and rainfall are used to drive the uncalibrated dynamical malaria model VECTRI to predict anomalies in transmission intensity 4 months ahead. It is demonstrated that the system has statistically significant skill at a number of sentinel sites in Uganda with high-quality data. Skill is also found at approximately 50% of the Ugandan health districts despite inherent uncertainties of unconfirmed health reports. A cost-loss economic analysis at three example sentinel sites indicates that the forecast system can have a positive economic benefit across a broad range of intermediate cost-loss ratios and frequency of transmission anomalies. We argue that such an analysis is a necessary first step in the attempt to translate climate-driven malaria information to policy-relevant decisions.

Exploring the Influence of Daily Climate Variables on Malaria Transmission and Abundance of Anopheles arabiensis over Nkomazi Local Municipality, Mpumalanga Province, South Africa

The recent resurgence of malaria incidence across epidemic regions in South Africa has been linked to climatic and environmental factors. An in-depth investigation of the impact of climate variability and mosquito abundance on malaria parasite incidence may therefore offer useful insight towards the control of this life-threatening disease. In this study, we investigate the influence of climatic factors on malaria transmission over Nkomazi Municipality. The variability and interconnectedness between the variables were analyzed using wavelet coherence analysis. Time-series analyses revealed that malaria cases significantly declined after the outbreak in early 2000, but with a slight increase from 2015. Furthermore, the wavelet coherence and time-lagged correlation analyses identified rainfall and abundance of Anopheles arabiensis as the major variables responsible for malaria transmission over the study region. The analysis further highlights a high malaria intensity with the variables from 1998-2002, 2004-2006, and 2010-2013 and a noticeable periodicity value of 256-512 days. Also, malaria transmission shows a time lag between one month and three months with respect to mosquito abundance and the different climatic variables. The findings from this study offer a better understanding of the importance of climatic factors on the transmission of malaria. The study further highlights the significant roles of An. arabiensis on malaria occurrence over Nkomazi. Implementing the mosquito model to predict mosquito abundance could provide more insight into malaria elimination or control in Africa.

Uncertainty in malaria simulations in the highlands of Kenya: Relative contributions of model parameter setting, driving climate and initial condition errors

In this study, experiments are conducted to gauge the relative importance of model, initial condition, and driving climate uncertainty for simulations of malaria transmission at a highland plantation in Kericho, Kenya. A genetic algorithm calibrates each of these three factors within their assessed prior uncertainty in turn to see which allows the best fit to a timeseries of confirmed cases. It is shown that for high altitude locations close to the threshold for transmission, the spatial representativeness uncertainty for climate, in particular temperature, dominates the uncertainty due to model parameter settings. Initial condition uncertainty plays little role after the first two years, and is thus important in the early warning system context, but negligible for decadal and climate change investigations. Thus, while reducing uncertainty in the model parameters would improve the quality of the simulations, the uncertainty in the temperature driving data is critical. It is emphasized that this result is a function of the mean climate of the location itself, and it is shown that model uncertainty would be relatively more important at warmer, lower altitude locations.

Mathematical modeling of climate change and malaria transmission dynamics: a historical review

Malaria, one of the greatest historical killers of mankind, continues to claim around half a million lives annually, with almost all deaths occurring in children under the age of five living in tropical Africa. The range of this disease is limited by climate to the warmer regions of the globe, and so anthropogenic global warming (and climate change more broadly) now threatens to alter the geographic area for potential malaria transmission, as both the Plasmodium malaria parasite and Anopheles mosquito vector have highly temperature-dependent lifecycles, while the aquatic immature Anopheles habitats are also strongly dependent upon rainfall and local hydrodynamics. A wide variety of process-based (or mechanistic) mathematical models have thus been proposed for the complex, highly nonlinear weather-driven Anopheles lifecycle and malaria transmission dynamics, but have reached somewhat disparate conclusions as to optimum temperatures for transmission, and the possible effect of increasing temperatures upon (potential) malaria distribution, with some projecting a large increase in the area at risk for malaria, but others predicting primarily a shift in the disease's geographic range. More generally, both global and local environmental changes drove the initial emergence of P. falciparum as a major human pathogen in tropical Africa some 10,000 years ago, and the disease has a long and deep history through the present. It is the goal of this paper to review major aspects of malaria biology, methods for formalizing these into mathematical forms, uncertainties and controversies in proper modeling methodology, and to provide a timeline of some major modeling efforts from the classical works of Sir Ronald Ross and George Macdonald through recent climate-focused modeling studies. Finally, we attempt to place such mathematical work within a broader historical context for the "million-murdering Death" of malaria.

Mosquito-Borne Diseases: Advances in Modelling Climate-Change Impacts

Vector-borne diseases are on the rise globally. As the consequences of climate change are becoming evident, climate-based models of disease risk are of growing importance. Here, we review the current state-of-the-art in both mechanistic and correlative disease modelling, the data driving these models, the vectors and diseases covered, and climate models applied to assess future risk. We find that modelling techniques have advanced considerably, especially in terms of using ensembles of climate models and scenarios. Effects of extreme events, precipitation regimes, and seasonality on diseases are still poorly studied. Thorough validation of models is still a challenge and is complicated by a lack of field and laboratory data. On a larger scale, the main challenges today lie in cross-disciplinary and cross-sectoral transfer of data and methods.

Behavioral, climatic, and environmental risk factors for Zika and Chikungunya virus infections in Rio de Janeiro, Brazil, 2015-16

The burden of arboviruses in the Americas is high and may result in long-term sequelae with infants disabled by Zika virus infection (ZIKV) and arthritis caused by infection with Chikungunya virus (CHIKV). We aimed to identify environmental drivers of arbovirus epidemics to predict where the next epidemics will occur and prioritize municipalities for vector control and eventual vaccination. We screened sera and urine samples (n = 10,459) from residents of 48 municipalities in the state of Rio de Janeiro for CHIKV, dengue virus (DENV), and ZIKV by molecular PCR diagnostics. Further, we assessed the spatial pattern of arbovirus incidence at the municipal and neighborhood scales and the timing of epidemics and major rainfall events. Lab-confirmed cases included 1,717 infections with ZIKV (43.8%) and 2,170 with CHIKV (55.4%) and only 29 (<1%) with DENV. ZIKV incidence was greater in neighborhoods with little access to municipal water infrastructure (r = -0.47, p = 1.2x10^-8). CHIKV incidence was weakly correlated with urbanization (r = 0.2, p = 0.02). Rains began in October 2015 and were followed one month later by the largest wave of ZIKV epidemic. ZIKV cases markedly declined in February 2016, which coincided with the start of a CHIKV outbreak. Rainfall predicted ZIKV and CHIKV with a lead time of 3 weeks each time. The association between rainfall and epidemics reflects vector ecology as the larval stages of Aedes aegypti require pools of water to develop. The temporal dynamics of ZIKV and CHIKV may be explained by the shorter incubation period of the viruses in the mosquito vector; 2 days for CHIKV versus 10 days for ZIKV.

Using Rainfall and Temperature Data in the Evaluation of National Malaria Control Programs in Africa

Since 2010, the Roll Back Malaria (RBM) Partnership, including National Malaria Control Programs, donor agencies (e.g., President's Malaria Initiative and Global Fund), and other stakeholders have been evaluating the impact of scaling up malaria control interventions on all-cause under-five mortality in several countries in sub-Saharan Africa. The evaluation framework assesses whether the deployed interventions have had an impact on malaria morbidity and mortality and requires consideration of potential nonintervention influencers of transmission, such as drought/floods or higher temperatures. Herein, we assess the likely effect of climate on the assessment of the impact malaria interventions in 10 priority countries/regions in eastern, western, and southern Africa for the President's Malaria Initiative. We used newly available quality controlled Enhanced National Climate Services rainfall and temperature products as well as global climate products to investigate likely impacts of climate on malaria evaluations and test the assumption that changing the baseline period can significantly impact on the influence of climate in the assessment of interventions. Based on current baseline periods used in national malaria impact assessments, we identify three countries/regions where current evaluations may overestimate the impact of interventions (Tanzania, Zanzibar, Uganda) and three countries where current malaria evaluations may underestimate the impact of interventions (Mali, Senegal and Ethiopia). In four countries (Rwanda, Malawi, Mozambique, and Angola) there was no strong difference in climate suitability for malaria in the pre- and post-intervention period. In part, this may be due to data quality and analysis issues.

Seasonally lagged effects of climatic factors on malaria incidence in South Africa

Globally, malaria cases have drastically dropped in recent years. However, a high incidence of malaria remains in some sub-Saharan African countries. South Africa is mostly malaria-free, but northeastern provinces continue to experience seasonal outbreaks. Here we investigate the association between malaria incidence and spatio-temporal climate variations in Limpopo. First, dominant spatial patterns in malaria incidence anomalies were identified using self-organizing maps. Composite analysis found significant associations among incidence anomalies and climate patterns. A high incidence of malaria during the pre-peak season (Sep-Nov) was associated with the climate phenomenon La Niña and cool air temperatures over southern Africa. There was also high precipitation over neighbouring countries two to six months prior to malaria incidence. During the peak season (Dec-Feb), high incidence was associated with positive phase of Indian Ocean Subtropical Dipole. Warm temperatures and high precipitation in neighbouring countries were also observed two months prior to increased malaria incidence. This lagged association between regional climate and malaria incidence suggests that in areas at high risk for malaria, such as Limpopo, management plans should consider not only local climate patterns but those of neighbouring countries as well. These findings highlight the need to strengthen cross-border control of malaria to minimize its spread.

Modeling and analyzing the impact of temperature and rainfall on mosquito population dynamics over Kwazulu-Natal, South Africa

Malaria parasites are strongly dependent on Anopheles mosquitoes for transmission; for this reason, mosquito population dynamics are a crucial determinant of malaria risk. However, temperature and rainfall play a significant role in both aquatic and adult stages of the Anopheles. Consequently, it is important to understand the biology of malaria vector mosquitoes in the study of malaria transmission. In this study, we develop a climate-based, ordinary-differential-equation model to analyze how rainfall and temperature determine mosquito population size. In the model, we consider in detail the influence of ambient temperature on gonotrophic and sporogonic cycles over Amajuba District, Kwazulu-Natal Province, South Africa. In particular, we further use the model to simulate the spatial distribution of the mosquito biting rate over the study region. Our results reflect high seasonality of the population of An. gambiae over the region and also demonstrate the influence of climatic factors on the mosquito population dynamics.

Assessing the Risk Factors Associated with Malaria in the Highlands of Ethiopia: What Do We Need to Know?

Malaria has been Ethiopia's predominant communicable disease for decades. Following the catastrophic malaria outbreak in 2003⁻2004, the Federal Ministry of Health (FMoH) took drastic public health actions to lower the burden of malaria. The FMoH achieved significant declines in malaria mortality and incidence, and recently declared its objective to achieve malaria elimination in low malaria transmission areas of Ethiopia by 2020. However, while the overall malaria prevalence has decreased, unpredictable outbreaks increasingly occur irregularly in regions previously considered "malaria-free". Such outbreaks have disastrous consequences on populations of these regions as they have no immunity against malaria. The Amhara Region accounts for 31% of Ethiopia's malaria burden and is targeted for malaria elimination by the FMoH. Amhara's epidemiological surveillance system faces many challenges to detect in a timely manner the unpredictable and irregular malaria outbreaks that occur in areas of otherwise low transmission. Despite the evidence of a shift in malaria transmission patterns, Amhara's malaria control interventions remain constrained to areas that are historically known to have stable malaria transmission. This paper discusses the influence of temperature and precipitation variability, entomological parameters, and human population mobility on malaria transmission patterns across the Amhara Region, and in particular, in areas of unstable transmission. We argue that malaria epidemiological surveillance systems can be improved by accounting for population movements in addition to environmental and entomological factors. However, to date, no study has statistically analyzed the interplay of population dynamics on environmental and entomological drivers of malaria transmission.

Modelling the influence of temperature and rainfall on the population dynamics of Anopheles arabiensis

BACKGROUND

Malaria continues to be one of the most devastating diseases in the world, killing more humans than any other infectious disease. Malaria parasites are entirely dependent on Anopheles mosquitoes for transmission. For this reason, vector population dynamics is a crucial determinant of malaria risk. Consequently, it is important to understand the biology of malaria vector mosquitoes in the study of malaria transmission. Temperature and precipitation also play a significant role in both aquatic and adult stages of the Anopheles.

METHODS

In this study, a climate-based, ordinary-differential-equation model is developed to analyse how temperature and the availability of water affect mosquito population size. In the model, the influence of ambient temperature on the development and the mortality rate of Anopheles arabiensis is considered over a region in KwaZulu-Natal Province, South Africa. In particular, the model is used to examine the impact of climatic factors on the gonotrophic cycle and the dynamics of mosquito population over the study region.

RESULTS

The results fairly accurately quantify the seasonality of the population of An. arabiensis over the region and also demonstrate the influence of climatic factors on the vector population dynamics. The model simulates the population dynamics of both immature and adult An. arabiensis. The simulated larval density produces a curve which is similar to observed data obtained from another study.

CONCLUSION

The model is efficiently developed to predict An. arabiensis population dynamics, and to assess the efficiency of various control strategies. In addition, the model framework is built to accommodate human population dynamics with the ability to predict malaria incidence in future.

A Regional Model for Malaria Vector Developmental Habitats Evaluated Using Explicit, Pond-Resolving Surface Hydrology Simulations

Dynamical malaria models can relate precipitation to the availability of vector breeding sites using simple models of surface hydrology. Here, a revised scheme is developed for the VECTRI malaria model, which is evaluated alongside the default scheme using a two year simulation by HYDREMATS, a 10 metre resolution, village-scale model that explicitly simulates individual ponds. Despite the simplicity of the two VECTRI surface hydrology parametrization schemes, they can reproduce the sub-seasonal evolution of fractional water coverage. Calibration of the model parameters is required to simulate the mean pond fraction correctly. The default VECTRI model tended to overestimate water fraction in periods subject to light rainfall events and underestimate it during periods of intense rainfall. This systematic error was improved in the revised scheme by including the a parametrization for surface run-off, such that light rainfall below the initial abstraction threshold does not contribute to ponds. After calibration of the pond model, the VECTRI model was able to simulate vector densities that compared well to the detailed agent based model contained in HYDREMATS without further parameter adjustment. Substituting local rain-gauge data with satellite-retrieved precipitation gave a reasonable approximation, raising the prospects for regional malaria simulations even in data sparse regions. However, further improvements could be made if a method can be derived to calibrate the key hydrology parameters of the pond model in each grid cell location, possibly also incorporating slope and soil texture.