The interrelationship between meteorological parameters and leptospirosis incidence in Hambantota district, Sri Lanka 2008-2017 and practical implications

Spatio-temporal modeling of human leptospirosis prevalence using the maximum entropy model

BACKGROUND

Leptospirosis, a zoonotic disease, stands as one of the prevailing health issues in some tropical areas of Iran. Over a decade, its incidence rate has been estimated at approximately 2.33 cases per 10,000 individuals. Our research focused on analyzing the spatiotemporal clustering of Leptospirosis and developing a disease prevalence model as an essential focal point for public health policymakers, urging targeted interventions and strategies.

METHODS

The SaTScan and Maximum Entropy (MaxEnt) modeling methods were used to find the spatiotemporal clusters of the Leptospirosis and model the disease prevalence in Iran. We incorporated nine environmental covariates by employing a spatial resolution of 1 km x 1 km, the finest resolution ever implemented for modeling Human Leptospirosis in Iran. These covariates encompassed the Digital Elevation Model (DEM), slope, displacement areas, water bodies, and land cover, monthly recorded Normalized Difference Vegetation Index (NDVI), monthly recorded precipitation, monthly recorded mean and maximum temperature, contributing significantly to our disease modeling approach. The analysis using MaxEnt yielded the Area Under the Receiver Operating Characteristic Curve (AUC) metrics for the training and test data, to evaluate the accuracy of the implemented model.

RESULTS

The findings reveal a highly significant primary cluster (p-value < 0.05) located in the western regions of the Gilan province, spanning from July 2013 to July 2015 (p-value < 0.05). Moreover, there were four more clusters (p-value < 0.05) identified near Someh Sara, Neka, Gorgan and Rudbar. Furthermore, the risk mapping effectively illustrates the potential expansion of the disease into the western and northwestern regions. The AUC metrics of 0.956 and 0.952 for the training and test data, respectively, underscoring the robust accuracy of the implemented model. Interestingly, among the variables considered, the influence of slope and distance from water bodies appears to be minimal. However, altitude and precipitation stand out as the primary determinants that significantly contribute to the prevalence of the disease.

CONCLUSIONS

The risk map generated through this study carries significant potential to enhance public awareness and inform the formulation of impactful policies to combat Leptospirosis. These maps also play a crucial role in tracking disease incidents and strategically directing interventions toward the regions most susceptible.

SARIMA and ARDL models for predicting leptospirosis in Anuradhapura district Sri Lanka

Leptospirosis is considered a neglected tropical disease despite its considerable mortality and morbidity. Lack of prediction remains a major reason for underestimating the disease. Although many models have been developed, most of them focused on the districts situated in the wet zone due to higher case numbers in that region. However, leptospirosis remains a major disease even in the dry zone of Sri Lanka. The objective of this study is to develop a time series model to predict leptospirosis in the Anuradhapura district situated in the dry zone of Sri Lanka. Time series data on monthly leptospirosis incidences from January 2008 to December 2018 and monthly rainfall, rainy days, temperature, and relative humidity were considered in model fitting. The first 72 months (55%) were used to fit the model, and the subsequent 60 months(45%) were used to validate the model. The log-transformed dependent variable was employed for fitting the Univariate seasonal ARIMA model. Based on the stationarity of the mean of the five variables, the ARDL model was selected as the multivariate time series technique. Residuals analysis was performed on normality, heteroskedasticity, and serial correlation to validate the model. The lowest AIC and MAPE were used to select the best model. Univariate models could not be fitted without adjusting the outliers. Adjusting seasonal outliers yielded better results than the models without adjustments. Best fitted Univariate model was ARIMA(1,0,0)(0,1,1)₁₂,(AIC-1.08, MAPE-19.8). Best fitted ARDL model was ARDL(1, 3, 2, 1, 0),(AIC-2.04,MAPE-30.4). The number of patients reported in the previous month, rainfall, rainy days, and temperature showed a positive association, while relative humidity was negatively associated with leptospirosis. Multivariate models fitted better than univariate models for the original data. Best-fitted models indicate the necessity of including other explanatory variables such as patient, host, and epidemiological factors to yield better results.

Spatial-temporal patterns and risk factors for human leptospirosis in Thailand, 2012-2018

Leptospirosis is a globally important zoonotic disease. The disease is particularly important in tropical and subtropical countries. Infections in humans can be caused by exposure to infected animals or contaminated soil or water, which are suitable for Leptospira. To explore the cluster area, the Global Moran's I index was calculated for incidences per 100,000 population at the province level during 2012-2018, using the monthly and annual data. The high-risk and low-risk provinces were identified using the local indicators of spatial association (LISA). The risk factors for leptospirosis were evaluated using a generalized linear mixed model (GLMM) with zero-inflation. We also added spatial and temporal correlation terms to take into account the spatial and temporal structures. The Global Moran's I index showed significant positive values. It did not demonstrate a random distribution throughout the period of study. The high-risk provinces were almost all in the lower north-east and south parts of Thailand. For yearly reported cases, the significant risk factors from the final best-fitted model were population density, elevation, and primary rice crop arable areas. Interestingly, our study showed that leptospirosis cases were associated with large areas of rice production but were less prevalent in areas of high rice productivity. For monthly reported cases, the model using temperature range was found to be a better fit than using percentage of flooded area. The significant risk factors from the model using temperature range were temporal correlation, average soil moisture, normalized difference vegetation index, and temperature range. Temperature range, which has strongly negative correlation to percentage of flooded area was a significant risk factor for monthly data. Flood exposure controls should be used to reduce the risk of leptospirosis infection. These results could be used to develop a leptospirosis warning system to support public health organizations in Thailand.

The correlation between three teleconnections and leptospirosis incidence in the Kandy District, Sri Lanka, 2004-2019

Background

Leptospirosis is a bacterial zoonosis. Leptospirosis incidence (LI) in Sri Lanka is high. Infected animals excrete leptospires into the environment via their urine. Survival of leptospires in the environment until they enter into a person and several other factors that influence leptospirosis transmission are dependent upon local weather. Past studies show that rainfall and other weather parameters are correlated with the LI in the Kandy district, Sri Lanka. El Niño Southern Oscillation (ENSO), ENSO Modoki, and the Indian Ocean Dipole (IOD) are teleconnections known to be modulating rainfall in Sri Lanka. There is a severe dearth of published studies on the correlations between indices of these teleconnections and LI.

Methods

We acquired the counts of leptospirosis cases notified and midyear estimated population data of the Kandy district from 2004 to 2019, respectively, from weekly epidemiology reports of the Ministry of Health and Department of Census and Statistics of Sri Lanka. We estimated weekly and monthly LI of Kandy. We obtained weekly and monthly teleconnection indices data for the same period from the National Oceanic and Atmospheric Administration (NOAA) of the USA and Japan Agency for Marine-Earth Science and Technology (JAMSTEC). We performed wavelet time series analysis to determine correlations with lag periods between teleconnection indices and LI time series. Then, we did time-lagged detrended cross-correlation analysis (DCCA) to verify wavelet analysis results and to find the magnitudes of the correlations detected.

Results

Wavelet analysis displayed indices of ENSO, IOD, and ENSO Modoki were correlated with the LI of Kandy with 1.9–11.5-month lags. Indices of ENSO showed two correlation patterns with Kandy LI. Time-lagged DCCA results show all indices of the three teleconnections studied were significantly correlated with the LI of Kandy with 2–5-month lag periods.

Conclusions

Results of the two analysis methods generally agree indicating that ENSO and IOD modulate LI in Kandy by modulating local rainfall and probably other weather parameters. We recommend further studies about the ENSO Modoki and LI correlation in Sri Lanka. Monitoring for extreme teleconnection events and enhancing preventive measures during lag periods can blunt LI peaks that may follow.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41182-021-00325-z.

Time series models for prediction of leptospirosis in different climate zones in Sri Lanka

In tropical countries such as Sri Lanka, where leptospirosis—a deadly disease with a high mortality rate—is endemic, prediction is required for public health planning and resource allocation. Routinely collected meteorological data may offer an effective means of making such predictions. This study included monthly leptospirosis and meteorological data from January 2007 to April 2019 from Sri Lanka. Factor analysis was first used with rainfall data to classify districts into meteorological zones. We used a seasonal autoregressive integrated moving average (SARIMA) model for univariate predictions and an autoregressive distributed lag (ARDL) model for multivariable analysis of leptospirosis with monthly average rainfall, temperature, relative humidity (RH), solar radiation (SR), and the number of rainy days/month (RD). Districts were classified into wet (WZ) and dry (DZ) zones, and highlands (HL) based on the factor analysis of rainfall data. The WZ had the highest leptospirosis incidence; there was no difference in the incidence between the DZ and HL. Leptospirosis was fluctuated positively with rainfall, RH and RD, whereas temperature and SR were fluctuated negatively. The best-fitted SARIMA models in the three zones were different from each other. Despite its known association, rainfall was positively significant in the WZ only at lag 5 (P = 0.03) but was negatively associated at lag 2 and 3 (P = 0.04). RD was positively associated for all three zones. Temperature was positively associated at lag 0 for the WZ and HL (P < 0.009) and was negatively associated at lag 1 for the WZ (P = 0.01). There was no association with RH in contrast to previous studies. Based on altitude and rainfall data, meteorological variables could effectively predict the incidence of leptospirosis with different models for different climatic zones. These predictive models could be effectively used in public health planning purposes.