
==== Front
BMJ Open
BMJ Open
bmjopen
bmjopen
BMJ Open
2044-6055
BMJ Publishing Group BMA House, Tavistock Square, London, WC1H 9JR

39160100
10.1136/bmjopen-2023-082503
bmjopen-2023-082503
Original Research
Public Health
1724
1506
Spatial-temporal analysis of climate and socioeconomic conditions on cholera incidence in Mozambique from 2000 to 2018: an ecological longitudinal retrospective study
https://twitter.com/jose_chaibo
http://orcid.org/0000-0002-7884-3302
Armando Chaibo Jose 12cjarmando.jose@gmail.com

Rocklöv Joacim 13joacim.rockloev@uni-heidelberg.de

Sidat Mohsin 2mmsidat@gmail.com

Tozan Yesim 4tozan@nyu.edu

Mavume Alberto Francisco 2amavume@gmail.com

Bunker Aditi 5aditi.bunker@uni-heidelberg.de

Sewe Maquins Odhiambo 16maquines.sewe@umu.se

1 Department of Public Health and Clinical Medicine, Sustainable Health Section, Umeå University, Umea, Sweden
2 Eduardo Mondlane University, Maputo, Mozambique
3 Heidelberg Institute of Global Health and Interdisciplinary Centre for Scientific Computing, Heidelberg University, Heidelberg, Germany
4 School of Global Public Health, NYU, New York, New York, USA
5 Heidelberg Institute of Global Health, Heidelberg University, Heidelberg, Germany
6 Universitatsklinikum Heidelberg Heidelberg Institute of Global Health, Heidelberg, Germany
Supplemental material: Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.

None declared.

Mr; cjarmando.jose@gmail.com
2024
19 8 2024
14 8 e08250311 1 2024
30 7 2024
Copyright © Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY. Published by BMJ.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Abstract

Objectives

This study aims to assess both socioeconomic and climatic factors of cholera morbidity in Mozambique considering both spatial and temporal dimensions.

Design

An ecological longitudinal retrospective study using monthly provincial cholera cases from Mozambican Ministry of Health between 2000 and 2018. The cholera cases were linked to socioeconomic data from Mozambique Demographic and Health Surveys conducted in the period 2000–2018 and climatic data; relative humidity (RH), mean temperature, precipitation and Normalised Difference Vegetation Index (NDVI). A negative binomial regression model in a Bayesian framework was used to model cholera incidence while adjusting for the spatiotemporal covariance, lagged effect of environmental factors and the socioeconomic indicators.

Setting

Eleven provinces in Mozambique.

Results

Over the 19-year period, a total of 153 941 cholera cases were notified to the surveillance system in Mozambique. Risk of cholera increased with higher monthly mean temperatures above 24°C in comparison to the reference mean temperature of 23°C. At mean temperature of 19°C, cholera risk was higher at a lag of 5–6 months. At a shorter lag of 1 month, precipitation of 223.3 mm resulted in an 57% increase in cholera risk (relative risk, RR 1.57 (95% CI 1.06 to 2.31)). Cholera risk was greatest at 3 lag months with monthly NDVI of 0.137 (RR 1.220 (95% CI 1.042 to 1.430)), compared with the reference value of 0.2. At an RH of 54%, cholera RR was increased by 62% (RR 1.620 (95% CI 1.124 to 2.342)) at a lag of 4 months. We found that ownership of radio RR 0.29, (95% CI 0.109 to 0.776) and mobile phones RR 0.262 (95% CI 0.097 to 0.711) were significantly associated with low cholera risk.

Conclusion

The derived lagged patterns can provide appropriate lead times in a climate-driven cholera early warning system that could contribute to the prevention and management of outbreaks.

Epidemiology
Infection control
Public health
Risk Factors
Swedish International Development Agency (SIDA). Not applicable for award/grant number.
==== Body
pmcSTRENGTHS AND LIMITATIONS OF THIS STUDY

Employs an ecological longitudinal retrospective design to analyse temporal and spatial trends in cholera incidence, yielding insights into the dynamics.

The method comprehensively elucidates cholera incidence in Mozambique by exploring the influence of climate and socioeconomic conditions employing advanced statistical techniques to explore the complex relationships.

Understanding the lagged effects of climate and socioeconomic factors is essential for discerning the temporal delay between risk exposure and cholera incidence, offering insights into long-term dynamics for effective management.

The use of cross-sectional Demographic and Health Survey data, collected every 5 years and aggregation to provincial, may have masked the relationship with cholera.

Introduction

Cholera is a waterborne disease caused by the bacteria Vibrio cholerae with the whole population at risk. Humans infection of cholera mostly occurs through faecal-oral route by the ingestion of contaminated food or water.13 Improved access to water, sanitation and hygiene (WASH) has been shown to decrease cholera occurrence by interrupting transmission routes.47 However, cholera continues to affect vulnerable populations in low-resource settings with inadequate water and sanitation infrastructure.810 Cholera imposes a substantial global burden, with an annual estimated incidence of 1.3–4 million cases worldwide and results in high mortality rates, ranging from 21 000 to 143 000 deaths each year.41113 Globally, the number of cholera cases showed a significant decline from 2014 to 2017. In 2015, there was a 9.49% reduction in reported cases compared with the 190 549 cases in 2014. This downward trend continued with a further 30.66% decrease in 2016, followed by an additional 37.02% decline in 2017. These figures indicate a consistent and substantial reduction in cholera cases worldwide over this period.11 14 Cholera remains endemic in sub-Saharan Africa (SSA)15 and more than 87 million people live in areas with high risk of cholera transmission.16 Between 1970 and 2011, there were a total of 3 221 050 cholera cases, and 9.79% (315 295) of the cases were reported in Mozambique.17 Between 2000 and 2015, 83% of the cholera mortalities reported by the WHO occurred in SSA.18 In 2017, there were 3220 deaths among the reported 1 79 835 cases, a case fatality rate of 1.8%.19 SSA including Mozambique experiences year-round transmission of cholera.19 20 Mozambique registered the first cholera cases in early 1973.21 22 Notified cholera cases from Mozambique account for about one-fifth to one-third of all cases reported in Africa across the years.21 23 In 2019, Sofala Province in Mozambique reported 6766 cholera cases with 8 associated fatalities.24 According to the United Nations Office for the Coordination of Humanitarian Affairs report, Mozambique recorded approximately 26 841 cholera cases and 123 fatalities in April 2023.25

The ongoing cholera burden in Mozambique is attributed to inadequate sanitation and extreme weather events, such as tropical cyclones and droughts2629 though the proportion of the population who have access to adequate sanitation in Mozambique has increased from 28% in 201527 to 38% in 2022.30 Studies have investigated how WASH improvements can reduce the risk of cholera transmission in Mozambique,31 32 Congo,33 Niger34 and Kenya35 and yet the effectiveness of these interventions has been not established.36 The key association between cholera risk and WASH factors can help to identify and prioritise areas with the highest need for intervention.3739 In addition to the WASH factors, cholera transmission is affected by various climatic factors such as temperature and precipitation.4042 Temperature is considered as an important factor for V. cholerae occurrence and propagation in the water reserves such as river, lake, oceans, stagnated water and groundwater.43 44 Furthermore, precipitation influences cholera transmission by increasing the risk of water contamination and the concentration of V. cholerae through runoff.45 46 These climate factors affect the temporal and spatial range of the V. cholerae and influence exposure pathways.47 In Mozambique, cholera transmission generally peaks between December and April during the hot and rainy season.48 49

Mozambique is a low-income country with weak WASH infrastructure and limited resources to scale up interventions against cholera outbreaks. Cholera has been linked to the households with low socioeconomic conditions, poverty, inequalities, education level and poor sanitation.31 50 The incubation period of cholera is relatively short ranging from just 2 hours to 5 days, thus the cases can increase quite substantially in a relatively short time.51 52 Identification of vulnerable areas for multisectoral interventions can be a great first step for cholera control and mitigation.53 Several interventions are used for cholera control and prevention including distribution of oral cholera vaccine, improved drinking water and sanitation facilities for communities at risk areas, as well as provisions of adequate healthcare.54 55 In this study, we aim to assess how climate and socioeconomic conditions affect cholera incidence in Mozambique considering the lagged patterns over space and time. Through analysis of temporal lag patterns, we can better understand and respond to the dynamics of cholera transmission, leading to more effective prevention and control strategies.

Materials and methods

Study area

Mozambique is one of the most vulnerable African countries to climate change along its coasts,56 with the majority of people living in rural areas along the coast.57 Mozambique is bounded by Tanzania on the north, Malawi and Zambia on the north-west, Zimbabwe on the west, South Africa and Swaziland on the south-west and the Indian Ocean on the east (see online supplemental figure S1).

Data

Weekly cholera cases at the province level were extracted from the Mozambique Ministry of Health disease surveillance system for the years 2000–2018 and aggregated to monthly totals. Data on the climatic factors, including daily precipitation, relative humidity (RH), Normalised Difference Vegetation Index (NDVI), and minimum and maximum temperature (Tmin, Tmax), were sourced from National Center for Environmental Prediction (NCEP).58

The NCEP employs advanced computer models that integrate observational data to simulate and predict weather patterns, including temperature, precipitation, NDVI and RH. These diverse data collection methods enable NCEP to generate reliable and detailed climate data which is available globally at different resolutions.59 Monthly summaries of the daily climatic data were computed and summarised at the provincial level based on shapefiles.60 In addition, the population data used in the analysis was sourced from Woldpop.61

Socioeconomic data from the Demographic Health Surveys (DHS) for the years 2003, 2009 and 2015 were included in the analysis.62 The socioeconomic variables included in this study followed similar variables used in Armando et al63 and were aggregated to the provincial level by applying the method used in Armando et al63 (see online supplemental tables S1 and S2).

Patient and public involvement

No patients or the public were involved in the design, conduct, reporting or dissemination of this study.

Statistical analysis

We employed an ecological longitudinal retrospective study design to analyse the temporal and spatial pattern of cholera incidence rates in Mozambique. We applied distributed lag non-linear models (DLNM)64 in a Bayesian framework with integrated nested Laplace approximation (INLA)65 to model the delayed and non-linear relationship between minimum temperature (Tmin), mean temperature (Tmean), maximum temperature (Tmax), total precipitation, RH and NDVI and cholera incidence adjusting for the DHS-derived socioeconomic indicators. In DLNM methodology, a bidimensional cross-basis function is created in order to simultaneously capture the lag and the variable dimension spaces.64 A natural cubic spline with 3 df was used for both the lag and the variable dimension. Lags of 0–6 months for climatic variables were assessed. The overdispersed cholera cases were modelled using a negative binomial distribution. Backward elimination based on significance at a 95% credible interval was used to select the socioeconomic variables in the final model.

A significance level of 5% means that there is a 5% chance of incorrectly rejecting the null hypothesis. The final model consisted of the cross-basis functions of the climatic variables and the significant socioeconomic factors.

We used the following reference values 18°C, 23°C, 28°C, 77 mm, 73%, 0.2 for Tmin, Tmean, Tmax, precipitation, RH and NDVI, respectively, when defining cross-basis functions. Interpretations are made in reference to these values.

The final model selected is represented in equation (1) below:

Yit∼Negbin

log⁡(Yit)=βo+ui+νi+ωt+fxj,vardf,lagdf+log⁡(Et)1

Yit represents cholera incidence rate for province i in month t, βo is the intercept, ui is the provincial exchangeable area effect while νi is the spatially structured random area effect. ωt is the monthly random effect modelled with random walk of order 1. fxj,vardf,lagdf represents the defined cross-basis function for climate and socioeconomic variable xj with vardf and lagdf representing the df for the variable and lag basis function. Et is the population offset added to the model with a coefficient of 1.

We used default INLA prior specifications for the provincial and monthly random effects.

All the analyses were done using R V.4.2.066 while DLNM64 and INLA65 packages in R were used for the analysis.

Results

Cholera cases and environmental variables

There were 153 941 cholera cases reported between 2000 and 2018 in the whole of Mozambique. The annual cholera incidence rates ranged from a high of 181.5 in the year 2002 to a low of 1.8 cases per 100 000 population in the year 2014 (see online supplemental tables S3 and S4). In the periods 2005–2009, 2010–2014 and 2015–2018, Mozambique experienced huge reductions in cholera incidence rates. For example, between 2005 and 2009, cholera incidence rate was 531.2 per 100 000 population with a decrease by 68.91% when considering the period 2000–2004 when it was 167.3 per 100 000 population. This downward trend continued in 2010–2014, with a further reduction of 92.54% to a rate of 40.1 per 100 000. By 2015–2018, cholera incidence rate had decreased by 92.74% compared with the 2000–2004 period to 39.1 per 100 000 (table 1).

Table 1 Summary of monthly cholera cases and environmental factors in Mozambique from 2000 to 2018

	2000–2018	2000–2004	2005–2009	2010–2014	2015–2018	
Cholera cases	
 Min	0	11	0	0	0	
 Max	8349	8349	5579	1529	1750	
 Mean (SD)	675.18 (1408.68)	1612.95 (2203.69)	613.42 (1099.62)	160.68 (305.78)	223.29 (422.78)	
Cholera incidence rate per 100 000	
 Min	0	0.61	0	0	0	
 Max	494.70	494.70	270.73	72.17	71.85	
 Mean (SD)	37.00 (81.12)	95.38 (130.24)	30.67 (53.65)	7.37 (14.17)	8.96 (17.26)	
Minimum temperature	
 Min	13.32	14.09	13.96	13.32	14.50	
 Max	23.0	22.75	23.0	22.69	22.76	
 Mean (SD)	18.98 (2.94)	19.14 (2.92)	19.11 (3.00)	18.70 (3.03)	18.97 (2.83)	
Mean temperature	
 Min	18.65	18.98	19.20	18.65	19.59	
 Max	27.34	26.84	27.34	26.99	27.34	
 Mean (SD)	23.73 (2.44)	23.64 (2.46)	23.98 (2.45)	23.47 (2.50)	23.85 (2.37)	
Maximum temperature	
 Min	23.65	23.66	24.36	23.65	24.59	
 Max	31.92	31.30	31.74	31.56	31.92	
 Mean (SD)	28.48 (2.06)	28.14 (2.11)	28.85 (2.00)	28.24 (2.06)	28.74 (2.01)	
Normalised Different Vegetation Index	
 Min	0.121	0.121	0.171	0.165	0.129	
 Max	0.365	0.365	0.343	0.349	0.331	
 Mean (SD)	0.240 (0.050)	0.240 (0.050)	0.250 (0.060)	0.240 (0.050)	0.230 (0.050)	
Relative humidity	
 Min	59.17	60.41	59.79	59.38	59.17	
 Max	91.48	91.48	89.79	89.52	88.62	
 Mean (SD)	73.88 (9.19)	74.86 (9.26)	72.75 (9.18)	74.15 (9.35)	73.72 (9.03)	
Precipitation	
 Min	1.35	2.42	2.48	1.35	2.89	
 Max	345.01	345.01	262.25	344.84	260.31	
 Mean (SD)	77.60 (80.85)	90.28 (90.52)	74.68 (79.34)	73.79 (78.48)	70.18 (73.20)	

The seasonal variation in cholera incidence rates followed the seasonal variation in environmental factors (figure 1, see online supplemental figure S2). During the period 2000–2018, cholera incidence rate decreased, and cholera burden was geographically concentrated in the central and northern parts of Mozambique (figure 1, see online supplemental figure S3-C and S4). The year 2014 had the lowest national average of 40 cholera cases per month while the year 2002 had the highest mean of about 2809 cases, as shown in online supplemental table S4. Cholera outbreaks in Mozambique show a strong seasonal pattern with an increased burden from November to April (figure 1). Cholera seasonality in Mozambique varies geographically, with peak transmission occurring earlier in the year in the northern regions, such as Cabo Delgado and Nampula provinces, and central regions, including Zambezia, Sofala and Manica provinces (see online supplemental figure S4). Over the study period, the northern and central regions of Mozambique received more precipitation compared with other parts of the country. The driest provinces were Maputo, Gaza and Maputo City (see online supplemental figure S5). The highest temperatures (minimum, mean and maximum) are observed along the central and northern coast while the lowest temperatures occur in Gaza and Maputo province (see online supplemental figure S6). Mozambique exhibits a simple seasonal temperature profile, with the lowest temperatures in July and the highest in December. Between 2000 and 2018, the annual mean RH in Mozambique ranged from 63.6% to 79.8% (see online supplemental figure S9), indicating significant fluctuations in atmospheric moisture levels over the years. Lower NDVI values are seen in Maputo City, Gaza, Maputo, Tete and Niassa provinces (see online supplemental figure S10) while the higher NDVI values are found in the central provinces such as Sofala, Manica and Zambezia (see online supplemental figures S10 and S11).

Figure 1 Seasonal variation in the monthly cholera incidence rates (red) and monthly precipitation (light blue), minimum temperature (orange), mean temperature (black), maximum temperature (dark red), relative humidity (light green) and NDVI (green) in Mozambique from 2000 to 2018. NDVI, Normalised Difference Vegetation Index.

Socioeconomic indicators

Summaries of the DHS-derived socioeconomic and WASH factors aggregated at the provincial level are shown in online supplemental table S1.

Specifically considering access to clean water and adequate sanitation, two key promoters of overall health and well-being, the majority of households (56.8%) in the provinces reported no access to treated drinking water (see online supplemental table S1). The proportion of households that shared a toilet facility with other households differed across the provinces ranging from 4% to 17% with an average of 10.1%. The average proportion of households with radio and mobile phone ownership was 39%, 67% with a range of 27%–55% and 45%–97%, respectively. As for educational attainment, on average, 25.1% of households in the provinces reported no education while the highest was 38.6% (see online supplemental table S1).

Model results

Temperature

Table 2 displays the relationship between mean temperature and cholera risk at different lags and percentiles. The relative risk (RR) of cholera gradually increased as minimum, mean and maximum temperature increased (figure 2A, see online supplemental figure S12). We observed a lower risk of cholera for mean temperatures between 19°C and 22°C at lags 0–4 months, compared with the reference of 23°C. The highest RR of cholera was at 28°C with a lag of 0 months (table 2). At 19°C, cholera risk was also high compared with that the reference but with a much longer lag of 5–6 months. Lower risk of cholera was, however, observed at much shorter lags. We observed an increasing RR of cholera for a mean temperature above 24°C (table 2). We observed an increased risk of cholera at maximum temperatures between 24°C and 33°C with lags of 2–3 months, compared with the reference value of 28°C, although this finding was not statistically significant (see online supplemental table S5).

Table 2 The association between mean temperature, precipitation, relative humidity, NDVI and cholera risk at lags 0–6 months

Lag	Mean temperature (Ref=23°C)	
Percentile	
5th	10th	25th	50th	75th	90th	99th	
19	20	22	24	26	27	28	
0	0.72	0.70	0.85	1.15	1.28	1.24	1.54	
1	0.79	0.78	0.90	1.11	1.22	1.21	1.19	
2	0.87	0.87	0.93	1.10	1.17	1.20	1.22	
3	0.94	0.94	0.97	1.04	1.14	1.19	1.25	
4	1.00	1.00	0.99	1.02	1.12	1.20	1.29	
5	1.07	1.05	1.00	1.01	1.11	1.21	1.33	
6	1.13	1.10	1.02	1.00	1.12	1.23	1.37	
	Precipitation (Ref=77 mm)	
1.2	2.4	7.4	33.7	114.5	223.3	395.1	
0	0.75	0.76	0.77	0.85	1.14	1.57*	1.80	
1	0.78	0.78	0.79	0.87	1.22*	1.46*	1.53	
2	0.86	0.86	0.87	0.92	1.07	1.21	1.12	
3	0.99	0.99	0.99	0.99	1.00	0.95	0.81	
4	1.05	1.05	1.05	1.03	0.97	0.88	0.77	
5	0.96	0.96	0.96	0.98	1.02	1.02	0.91	
6	0.76	0.77	0.78	0.86	1.12	1.37	1.15	
	Relative humidity (Ref=73%)	
54	61	67	74	82	88	92	
0	0.56	0.61	0.78	1.03	1.13	1.00	0.89	
1	0.85	0.87	0.92	1.01	1.08	1.09	1.09	
2	1.22	1.16	1.07	0.99	1.04	1.18	1.32	
3	1.53*	1.40*	1.17	0.98	1.03	1.26	1.50	
4	1.62*	1.45*	1.19*	0.98	1.04	1.31	1.62	
5	1.51*	1.36	1.15	0.98	1.07	1.35	1.65	
6	1.32	1.20	1.08	0.99	1.11	1.37	1.64	
	NDVI (Ref=0.2)	
0.137	0.157	0.189	0.238	0.295	0.337	0.400	
0	1.02	1.00	1.00	0.99	1.00	1.03	1.09	
1	1.12	1.09	1.03	0.91	0.81	0.79	0.83	
2	1.20*	1.15*	1.04*	0.85*	0.69	0.65*	0.67	
3	1.22*	1.17*	1.05*	0.83*	0.65*	0.59*	0.59*	
4	1.16*	1.13*	1.03*	0.87*	0.70*	0.63*	0.59*	
5	1.05	1.04	1.01	0.94	0.84	0.76	0.65	
6	0.93	0.94	0.98	1.05	1.06	0.97	0.74	
* Significant.

NDVINormalised Difference Vegetation Index

Figure 2 Overall effect of (A), precipitation (B), relative humidity (C) and NDVI (D) on cholera risk in Mozambique, 2000–2018. The reference values for, precipitation, relative humidity and NDVI were 23°C, 77 mm, 73% and 0.2, respectively. NDVI, Normalised Difference Vegetation Index.

Precipitation

Table 2 shows the estimated relationship between precipitation and relative cholera risk at lags 0–6 months. There were delayed effects on cholera risk with precipitation, with the maximum effects of heavy precipitation reached at shorter lags of less than 1 month. For example, total precipitation of 223.3 mm with a lag of 1 month, a resulted in 57% increase in cholera risk (RR 1.57 (95% CI 1.06 to 2.31)). We also observed that the effect of heavy precipitation on cholera risk was attenuated at lags of 3–5 months. The overall association between cumulative precipitation and cholera risk is displayed in figure 2B. Increasing amounts of precipitation were associated with increased cholera risk at lags 0–2 months (table 2). The RR of cholera decreased as precipitation increased at a lag of 4 months (table 2). At precipitation of above 500 mm, cholera risk was lower (figure 2B).

Relative humidity

Table 2 shows exposure-lag response surface for RH and cholera risk. For RH between 54% and 67%, the risk of cholera increased at longer lags for RH of 54% (RR 1.620 (95% CI 1.124 to 2.342)) at a lag of 4 months. This association is more intense at a lag of 3–5 months. Higher RH values were found to be significantly associated with cholera risk at longer lags. At an RH of 67%, cholera risk increased by 19% (RR 1.190 (95% CI 1.017 to 1.392) at a lag of 4 months. The non-linear relationship between RH and cholera risk is shown in figure 2C. For RH greater than 80% and below 70%, cholera risk consistently increased though the association was not statistically significant (figure 2C).

Normalised Difference Vegetation Index

Table 2 displays the lag–response association between NDVI and cholera risk. For NDVI of 0.137, cholera risk was highest at a lag of 3 months (RR 1.220 (95% CI 1.042 to 1.430)). We observed a decreased RR for NDVI values from 0.238 to 0.400 compared with the reference value of 0.2, for example, at an NDVI value of 0.337 cholera risk was 37% lower (RR 0.630 (95% CI 0.468 to 0.857)), compared with the reference. The overall relationship between NDVI and cholera risk is displayed in figure 2D. In comparison to the reference value of 0.2, cholera risk was significantly lower for NDVI values above 0.22 (figure 2D). Specifically, at a monthly mean NDVI of 0.38, cholera risk was 88.9% lower, (RR 0.110 (95% CI 0.012 to 0.984)), and at 0.27, it was 66.78% lower (RR 0.332 (95% CI 0.118 to 0.928)) (figure 2D). For NDVI values below the reference value of 0.2, the risk of cholera slightly increased though it did not show a significant relationship with cholera risk (figure 2D).

Socioeconomic factors for cholera

Table 3 shows the association between household socioeconomic indicators and cholera risk. A high proportion of households with a radio (used for delivering outbreak messages) was significantly associated with very low risk, 70.1% lower (RR 0.29 (95% CI 0.109 to 0.776)), compared with having a low percentage of radio ownership. Similarly, high proportion of mobile phones ownership was associated with a 73.76% decreased risk of cholera compared households with low ownership (RR 0.262 (95% CI 0.097 to 0.711)). Conversely, we observed an increased risk of cholera in households sharing a toilet, though this was not statistically significant.

Table 3 The association between socioeconomic variables and cholera risk

	Radio ownership (Ref=26.68%)	
Percentile	
5th	10th	25th	50th	75th	90th	
27	37	39	40	41	50	
RR	1.01	0.56	0.39	0.33*	0.28*	0.29*	
	Mobile phone ownership (Ref=45%)	
45.4	45.7	46.9	69.9	86.1	87.9	
RR	0.997	0.994	0.986	0.697	0.262*	0.292*	
	Toilet sharing (Ref=7%)	
4.02	5.02	7.76	9.23	13.13	14.2	
RR	1.096	1.034	1.043	1.212	1.302	1.108	
The reference value from households with a radio, mobile phones and those who share toilet facilities is 26.68%, 45% and 7%, respectively

* Significant.

RRrelative risk

Discussion

We analysed the non-linear relationship between delayed climatic conditions and cholera risk in Mozambique while adjusting for socioeconomic conditions, space and time dependencies in a Bayesian framework making inferences using the computationally efficient INLA methodology. We show that temperature, precipitation, RH and NDVI influence the spatiotemporal distribution of cholera incidences in Mozambique.

Meteorological factors play a crucial role in cholera transmission pathways, for example, high temperatures in a warmer season may provide a suitable condition for V. cholerae to proliferate.424467 In this study, we found that mean temperature above 24°C increased cholera risk at lags 0–6 months. At mean temperatures between 19°C and 22°C, the risk was lower at shorter lags. A study looking at countries in SSA linked elevated cholera risk with mean temperature at lags above 2 months1 which is consistent with our study. In Zanzibar, cholera risk increased with 1°C rise in temperature with a delay of 4 months.68 In Tanzania, they also found a significant relationship between temperature and cholera risk, however, they did not consider the lagged effect of temperature.69

Similar to temperature, precipitation also modulates cholera risk acting at different spatial and temporal scales. In this study, we found that increased precipitation was associated with an increased risk of cholera occurrence at a 1-month lag. With monthly precipitation of between 150 and 500 mm, the risk of cholera was high at a shorter lag of 1 month and reduced at higher lags of 3–5 months. This resonates with our finding as the incubation period for V. cholerae causing cholera, is usually less than a month with most symptoms appearing within 1–3 days after ingestion of contaminated food or water. In South Sudan, precipitation was found to be the most important driver for cholera, however, this study also did not explore the delayed effect of precipitation.45 In Bangladesh, a significant 1-month delay, similar to our finding, was found between precipitation and cholera risk where a 1 mm/day increase in mean precipitation was linked with a 6.5% increase in cholera transmission,70 while in this study, we found increase of 57% at precipitation of 223.3. In Zambia, higher cholera risk was found at a lag less than 1 month for precipitation above 50 mm.71 Though in South Africa,72 precipitation had the strongest association with cholera risk at 2-month lag, while in Haiti,73 shorter lags less than 1 month were found, very similar to our findings.

We found a significant association between NDVI values and cholera risk, with lower NDVI values linked to higher cholera risk and higher NDVI values associated with reduced risk. Lower NDVI values are often associated with dry or drought conditions, which can exacerbate cholera transmission. During droughts, water scarcity forces populations to rely on contaminated water sources, such as stagnant ponds, rivers and lakes for drinking and household purposes, thereby increasing the risk of cholera transmission. NDVI is the most used vegetation index in different areas including epidemiology and climate variability.74 NDVI increases with the amount of green biomass and precipitation. In this study, we found NDVI to be positively associated with precipitation.

Our findings indicated that RH levels between 54% and 67% were linked to an increased risk of cholera, with effects observed after delays ranging from 0 to 6 months. Specifically, the risk of cholera was notably higher when longer lag periods were considered within this RH range.

This is consistent with the study done in India.75 Similarly, in Zanzibar, cholera outbreaks were found to be significantly associated with RH at a lag of 5 months.68 RH critically influences cholera transmission dynamics by modulating the survival of V. cholerae bacteria in environmental reservoirs and impacting human susceptibility to infection. RH directly influences the persistence of V. cholerae in various water sources, including rivers, ponds and coastal areas, potentially extending the bacterium’s survival duration in water and thereby augmenting the risk of exposure through contaminated water sources.

Household’s socioeconomic disparities can help elucidate the fluctuations in cholera transmission between and within provinces in Mozambique during the outbreak in addition to the climatic factors. Studies done in Zambia76 77 and Kenya78 79 showed that households sharing toilet facilities were at higher risk of cholera incidences, a finding consistent with our study. Sharing toilet facilities increases cholera risk due to factors such as inadequate sanitation, increased cross-contamination and limited access to hygiene resources. Shared toilet facilities often lack proper sanitation and maintenance, resulting in the accumulation of faecal matter and increased likelihood of faecal-oral transmission of V. cholerae.

A study conducted in Haiti found that household with radio had lower cholera risk80 which is similar to this study. Additionally, mobile phone ownership has been shown to play an important role in cholera interventions and prevention strategies81 and our results corroborate similar findings that a higher prevalence of mobile phone ownership results in reduced cholera risk as it enhances communication and information dissemination, which can lead to better hygiene practices and faster response to cholera outbreaks. An example of this was shown in Mozambique where a short message service was used to gather information to support cholera response activities and vaccination campaigns.81

We have shown that temperature precipitation, NDVI and RH contributed to increased cholera risk at different lag periods in Mozambique. The results of the study highlight the need for identifying vulnerable populations to further support cholera control efforts and the utility of combining climate, environmental conditions, regional spatial stratification, socioeconomic factors and public health interventions related to cholera risk. Appropriate interventions for control and elimination of cholera in Mozambique require a multidisciplinary innovative approach, to prevention and sustained political commitment at national, provincial and district levels, as well as continued investment in improving the availability of WASH services and infrastructure with appropriate quality and quantity across urban and rural areas of the country.

Results from our current study support the uptake of measures that can also help achieve the Sustainable Development Goals, especially those related to (1) good health and well-being, (2) clean water and sanitation and (3) climate action.

This study had some limitations that need to be acknowledged. The cross-sectional nature of the Demographic and Health Survey (DHS) data limits our ability to draw precise conclusions about their influence on cholera transmission in time. Since DHS data are collected every few years, it may not accurately capture the temporal dynamics of cholera transmission and outbreak patterns within the study area. This temporal mismatch means that the data often do not coincide with the actual periods of cholera occurrence.

Spatial variability is another limitation, as the cholera, climate and DHS data were aggregated to the broader provincial level in our study. This aggregation may have masked the important local-level variations in cholera risk and transmission dynamics. Consequently, our analysis may have overlooked important microlevel heterogeneities that are critical for understanding the nuanced patterns of cholera spread and identifying targeted intervention strategies. Additionally, the study’s reliance on secondary sources of socioeconomic data such as the DHS limits the inclusion of certain key indicators such as infrastructural factors, for example, water quality, sanitation infrastructure and healthcare access.

Our study offers important contributions to understanding cholera transmission, the limitations related to the cross-sectional nature of the DHS data, potential recall bias, spatial aggregation, lack of certain key variables and inability to assess long-term trends must be carefully considered when interpreting the findings.

Conclusions

This study explored lag patterns of monthly climate variables and cholera morbidity in Mozambique. Results show how the use of climate variables could help in the early warning system and control of cholera. The public health measures for the prevention and control of cholera occurrence must factor in climate considerations if it is to reduce vulnerability and increase the adaptive capacity of the population in Mozambique. Evidence presented here can support the planning, monitoring and evaluation of cholera control efforts.

supplementary material

10.1136/bmjopen-2023-082503 online supplemental file 1

Acknowledgements

The first author is thankful to the Ministry of Health in Mozambique for providing the data. We are deeply grateful to all those who played a role in the success of this research paper.

Data availability statement

Data are available on reasonable request.

Review Process File
19 08 2024

Funding: The study is funded by Swedish International Development Agency (SIDA). Not applicable for award/grant number.

Prepublication history and additional supplemental material for this paper are available online. To view these files, please visit the journal online (https://doi.org/10.1136/bmjopen-2023-082503).

Provenance and peer review: Not commissioned; externally peer reviewed.

Patient consent for publication: Not applicable.

Data availability free text: The cholera datasets analysed during the current study are available from the corresponding author on reasonable request.

Map disclaimer: The inclusion of any map (including the depiction of any boundaries therein), or of any geographic or locational reference, does not imply the expression of any opinion whatsoever on the part of BMJ concerning the legal status of any country, territory, jurisdiction or area or of its authorities. Any such expression remains solely that of the relevant source and is not endorsed by BMJ. Maps are provided without any warranty of any kind, either express or implied.

Patient and public involvement: Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
==== Refs
References

1 Perez-Saez J Lessler J Lee EC et al The seasonality of cholera in sub-Saharan Africa: a statistical modelling study Lancet Glob Health 2022 10 e831 9 10.1016/S2214-109X(22)00007-9 35461521
2 Cliff JL Zinkin P Martelli A A hospital outbreak of cholera in Maputo, Mozambique Trans R Soc Trop Med Hyg 1986 80 473 6 10.1016/0035-9203(86)90349-4 3798545
3 Goh KT Teo SH Lam S et al Person-to-person transmission of cholera in a psychiatric hospital J Infect 1990 20 193 200 10.1016/0163-4453(90)90994-j 2341728
4 Jones N Bouzid M Few R et al Water, sanitation and hygiene risk factors for the transmission of cholera in a changing climate: using a systematic review to develop a causal process diagram J Water Health 2020 18 145 58 10.2166/wh.2020.088 32300088
5 Mukandavire Z Liao S Wang J et al Estimating the reproductive numbers for the 2008-2009 cholera outbreaks in Zimbabwe Proc Natl Acad Sci U S A 2011 108 8767 72 10.1073/pnas.1019712108 21518855
6 Travis Yates TY Jelena Allen Vujcic JAV Myriam Leandre Joseph MLJ et al Water, sanitation, and hygiene interventions in outbreak response: a synthesis of evidence Waterlines 2018 37 5 30 10.3362/1756-3488.17-00015
7 Clemens JD Nair GB Ahmed T et al Cholera The Lancet 2017 390 1539 49 10.1016/S0140-6736(17)30559-7
8 Musoke D Namata C Lubega GB et al The role of environmental health in preventing antimicrobial resistance in low- and middle-income countries Environ Health Prev Med 2021 26 100 10.1186/s12199-021-01023-2 34610785
9 Bwire G Munier A Ouedraogo I et al Epidemiology of cholera outbreaks and socio-economic characteristics of the communities in the fishing villages of Uganda: 2011-2015 PLoS Negl Trop Dis 2017 11 e0005407 10.1371/journal.pntd.0005407 28288154
10 Sodjinou VD Sodjinou VD Sodjinou VD et al The 2021 cholera outbreak in West Africa: epidemiology and public health implications Arch Clin Biomed Res 2022 6 296 307
11 World Health Organization World Health Statistics 2018: Monitoring Health for the SDGs, Sustainable Development Goals World Health Organization 2018
12 Ali M Nelson AR Lopez AL et al Updated global burden of cholera in endemic countries PLoS Negl Trop Dis 2015 9 e0003832 10.1371/journal.pntd.0003832 26043000
13 Ali M Lopez AL You YA et al The global burden of cholera Bull World Health Organ 2012 90 209 218A 10.2471/BLT.11.093427 22461716
14 mondiale de la Santé, O. and W.H. Organization Weekly epidemiological record= relevé épidémiologique hebdomadaire, 2016 Wkly Epidemiol Rec 2016 91 432 40
15 Zerbo A Castro Delgado R González PA A review of the risk of cholera outbreaks and urbanization in sub-Saharan Africa J Biosaf Biosec 2020 2 71 6 10.1016/j.jobb.2020.11.004
16 Lessler J Moore SM Luquero FJ et al Mapping the burden of cholera in sub-Saharan Africa and implications for control: an analysis of data across geographical scales The Lancet 2018 391 1908 15 10.1016/S0140-6736(17)33050-7
17 Mengel MA et al Cholera Outbreaks in Africa Cholera Outbreaks 2014 117 44
18 WHO Number of Reported Deaths - Data by Country WHO 2023
19 World Health Organization Cholera, 2017 Weekly Epidemiological Record Geveva 2017
20 Legros D Partners of the Global Task Force on Cholera Control Global cholera epidemiology: opportunities to reduce the burden of cholera by 2030 J Infect Dis 2018 218 S137 40 10.1093/infdis/jiy486 30184102
21 W.R.H Organization Global Task Force on Cholera Control Cholera Country Profile: Mozambique WHO 2013
22 WHO Cholera Country Profile: Mozambique. Global Task on Cholera Contro World Health Organization 2009
23 Langa JP Sema C De Deus N et al Epidemic waves of cholera in the last two decades in Mozambique J Infect Dev Ctries 2015 9 635 41 10.3855/jidc.6943 26142674
24 Mongo E Cambaza E Nhambire R Outbreak of cholera due to cyclone Idai in central mozambique IntechOpen 2019 10.3390/ijerph16162925
25 OCHA Mozambique: Tropical Cyclone Freddy, Floods and Cholera Situation Report No.2 2023 2
26 Lequechane JD Mahumane A Chale F et al Mozambique’s response to cyclone Idai: how collaboration and surveillance with water, sanitation and hygiene (WASH) interventions were used to control a cholera epidemic Infect Dis Poverty 2020 9 68 10.1186/s40249-020-00692-5 32546268
27 Borja-Vega C et al Findings of the Mozambique Water Supply, Sanitation, and Hygiene Poverty Diagnostic The World Bank 2018
28 OCHA Mozambique Response Plan: Cyclone Freddy, Floods & Cholera (March - September 2023) 2023a
29 Bié AJ de Camargo R Mavume AF et al Numerical modeling of storm surges in the coast of Mozambique: the cases of tropical cyclones Bonita (1996) and Lisette (1997) Ocean Dyn 2017 67 1443 59 10.1007/s10236-017-1095-7
30 UNICEF and WHO Progress on Household Drinking Water, Sanitation and Hygiene 2000–2022: Special Focus on Gender 2023
31 Collins AE Lucas ME Islam MS et al Socio‐economic and environmental origins of cholera epidemics in Mozambique: guidelines for tackling uncertainty in infectious disease prevention and control Int J Environ Stud 2006 63 537 49 10.1080/00207230600963122
32 WHO Cholera - Mozambique 2023
33 Gallandat K Jeandron A Ross I et al The impact of improved water supply on cholera and diarrhoeal diseases in Uvira, Democratic republic of the Congo: a protocol for a pragmatic stepped-wedge cluster randomised trial and economic evaluation Trials 2021 22 408 10.1186/s13063-021-05249-x 34154636
34 Graveleau J Reserva ME Keita A et al Influence of community-led total sanitation and water coverages in the control of cholera in Madarounfa, Niger (2018) Front Public Health 2021 9 643079 10.3389/fpubh.2021.643079 33996720
35 Kiama C Okunga E Muange A et al Mapping of cholera hotspots in Kenya using epidemiologic and water, sanitation, and hygiene (WASH) indicators as part of Kenya’s new 2022-2030 cholera elimination plan PLOS Negl Trop Dis 2023 17 e0011166 10.1371/journal.pntd.0011166 36930650
36 Taylor DL Kahawita TM Cairncross S et al The impact of water, sanitation and hygiene interventions to control cholera: a systematic review PLoS One 2015 10 e0135676 10.1371/journal.pone.0135676 26284367
37 Ratnayake R Finger F Azman AS et al Highly targeted spatiotemporal interventions against cholera epidemics, 2000-19: a scoping review Lancet Infect Dis 2021 21 e37 48 10.1016/S1473-3099(20)30479-5 33096017
38 Burrowes V Perin J Monira S et al Risk factors for household transmission of Vibrio cholerae in Dhaka, Bangladesh (CHoBI7 trial) Am J Trop Med Hyg 2017 96 1382 7 1382 10.4269/ajtmh.16-0871 28719281
39 Wolfe M Kaur M Yates T et al A systematic review and meta-analysis of the association between water, sanitation, and hygiene exposures and cholera in case-control studies Am J Trop Med Hyg 2018 99 534 45 534 10.4269/ajtmh.17-0897 29968551
40 Change, I.P.O.C., Ipcc Climate Change 2014
41 Islam MS Zaman MH Islam MS et al Environmental reservoirs of vibrio cholerae Vaccine (Auckl) 2020 38 A52 62 10.1016/j.vaccine.2019.06.033
42 Usmani M Brumfield KD Magers BM et al Combating cholera by building predictive capabilities for pathogenic vibrio cholerae in Yemen Sci Rep 2023 13 2255 10.1038/s41598-022-22946-y 36755108
43 Kruger SE Lorah PA Okamoto KW Mapping climate change’s impact on cholera infection risk in Bangladesh PLOS Glob Public Health 2022 2 e0000711 10.1371/journal.pgph.0000711 36962590
44 Asadgol Z Mohammadi H Kermani M et al The effect of climate change on cholera disease: the road ahead using artificial neural network PLoS One 2019 14 e0224813 10.1371/journal.pone.0224813 31693708
45 Lemaitre J Pasetto D Perez-Saez J et al Rainfall as a driver of epidemic cholera: comparative model assessments of the effect of intra-seasonal precipitation events Acta Trop 2019 190 235 43 10.1016/j.actatropica.2018.11.013 30465744
46 Christaki E Dimitriou P Pantavou K et al The impact of climate change on cholera: a review on the global status and future challenges Atmosphere (Basel) 2020 11 449 10.3390/atmos11050449
47 UNICEF Cholera: A Global Call to Action 2023
48 Baltazar CS Pezzoli L Baloi LD et al Conditions to eliminate cholera in Mozambique - the pathway for the development of the national cholera plan Pan Afr Med J 2022 42 279 10.11604/pamj.2022.42.279.36368 36405663
49 WHO Weekly Regional Cholera Bulletin WHO: WHO - World Health Organisation 2023
50 Das R Nasrin S Palit P et al Vibrio choleraein rural and urban Bangladesh, findings from hospital-based surveillance, 2000–2021 Sci Rep 2023 13 6411 10.1038/s41598-023-33576-3 37076586
51 World Health Organization Ending Cholera a Global Roadmap to 2030, in Ending Cholera a Global Roadmap to 2030 2017 32
52 Azman AS Rudolph KE Cummings DAT et al The incubation period of cholera: a systematic review J Infect 2013 66 432 8 10.1016/j.jinf.2012.11.013 23201968
53 Organization WH Cholera vaccines: WHO position paper WER 2010 85 117 28 20349546
54 Zuckerman JN Rombo L Fisch A The true burden and risk of cholera: implications for prevention and control Lancet Infect Dis 2007 7 521 30 10.1016/S1473-3099(07)70138-X 17584531
55 Anetor GO Abraham F Knowledge of cholera and its prevention amongst urban residents of a district in Abuja: the pivotal role of health education REJHS 2020 8 102 12 10.4314/rejhs.v8i2.6
56 INGC Main Report: INGC Climate Change. Study on the Impact of Climate Change on Disaster Risk in Mozambique Synthesis Report, in N.I.f.D. Management Maputo-Mozambique INGC, Mozambique 2009
57 UN-HABITAT Mozambique Country Brief 2023: To Promote Sustainable Urbanization in Mozambique as Driver of Socioeconomic Development, Climate Resilience and Peace 2023
58 NCEP NCEP-reanalysis 2022 Available https://psl.noaa.gov/data/gridded/help.html#FTP
59 Campos RM Analysis of environmental prediction data of temperature and relative humidity from NCEP/NOAA IJESNR 2021 27 IJESNR 10.19080/IJESNR.2021.27.556217
60 ROSEA O Mozambique Administrative Boundary Common Operational Database (COD-AB) OCHA 2019
61 Tatem AJ WorldPop, open data for spatial demography Sci Data 2017 4 170004 10.1038/sdata.2017.4 28140397
62 DHS Survey search 2015 Available http://dhsprogram.com/What-We-Do/survey-search.cfm
63 Armando CJ Rocklöv J Sidat M et al Climate variability, socio-economic conditions and vulnerability to malaria infections in Mozambique 2016-2018: a spatial temporal analysis Front Public Health 2023 11 1162535 10.3389/fpubh.2023.1162535 37325319
64 Gasparrini A Armstrong B Kenward MG Distributed lag non-linear models Stat Med 2010 29 2224 34 10.1002/sim.3940 20812303
65 Rue H Martino S Chopin N Approximate bayesian inference for latent gaussian models by using integrated nested laplace approximations J R Stat Soc Series B Stat Methodol 2009 71 319 92 10.1111/j.1467-9868.2008.00700.x
66 R Core Team, A. and R.C. Team R: A Language and Environment for Statistical Computing Vienna, Austria R Foundation for Statistical Computing 2022
67 Shackleton D Memon FA Nichols G et al Mechanisms of cholera transmission via environment in India and Bangladesh: state of the science review Rev Environ Health 2024 39 313 29 10.1515/reveh-2022-0201 36639850
68 Reyburn R Kim DR Emch M et al Climate variability and the outbreaks of cholera in Zanzibar, East Africa: a time series analysis Am J Trop Med Hyg 2011 84 862 9 10.4269/ajtmh.2011.10-0277 21633020
69 Trærup SLM Ortiz RA Markandya A The costs of climate change: a study of cholera in Tanzania Int J Environ Res Public Health 2011 8 4386 405 10.3390/ijerph8124386 22408580
70 Constantin de Magny G Murtugudde R Sapiano MRP et al Environmental signatures associated with cholera epidemics Proc Natl Acad Sci U S A 2008 105 17676 81 10.1073/pnas.0809654105 19001267
71 Luque Fernández MA Bauernfeind A Jiménez JD et al Influence of temperature and rainfall on the evolution of cholera epidemics in Lusaka, Zambia, 2003-2006: analysis of a time series Trans R Soc Trop Med Hyg 2009 103 137 43 10.1016/j.trstmh.2008.07.017 18783808
72 Mendelsohn J Dawson T Climate and cholera in KwaZulu-Natal, South Africa: the role of environmental factors and implications for epidemic preparedness Int J Hyg Environ Health 2008 211 156 62 10.1016/j.ijheh.2006.12.002 17383231
73 Eisenberg MC Kujbida G Tuite AR et al Examining rainfall and cholera dynamics in Haiti using statistical and dynamic modeling approaches Epidemics 2013 5 197 207 10.1016/j.epidem.2013.09.004 24267876
74 Garai S Khatun M Singh R et al Assessing correlation between rainfall, Normalized Difference Vegetation Index (NDVI) and Land Surface Temperature (LST) in Eastern India Saf Extreme Environ 2022 4 119 27 10.1007/s42797-022-00056-2
75 Rajendran K Sumi A Bhattachariya MK et al Influence of relative humidity in vibrio cholerae infection: a time series model Indian J Med Res 2011 133 138 45 138 21415487
76 DuBois AE Sinkala M Kalluri P et al Epidemic cholera in urban Zambia: hand soap and dried fish as protective factors Epidemiol Infect 2006 134 1226 30 10.1017/S0950268806006273 16623992
77 Sasaki S Suzuki H Igarashi K et al Spatial analysis of risk factor of cholera outbreak for 2003-2004 in a peri-urban area of Lusaka, Zambia Am J Trop Med Hyg 2008 79 414 21 18784235
78 Shultz A Omollo JO Burke H et al Cholera outbreak in Kenyan refugee camp: risk factors for illness and importance of sanitation Am J Trop Med Hyg 2009 80 640 5 19346392
79 Mahamud AS Ahmed JA Nyoka R et al Epidemic cholera in Kakuma Refugee Camp, Kenya, 2009: the importance of sanitation and soap J Infect Dev Ctries 2012 6 234 41 10.3855/jidc.1966 22421604
80 Grandesso F Allan M Jean-Simon PSJ et al Risk factors for cholera transmission in Haiti during inter-peak periods: insights to improve current control strategies from two case-control studies Epidemiol Infect 2014 142 1625 35 10.1017/S0950268813002562 24112364
81 Almeida S et al Enhancing Community Engagement Through Data Collection: Controlling the Cholera Epidemic in Mozambique 2020
