Spatio-temporal dynamics of three diseases caused by *Aedes*-borne arboviruses in Mexico

Background The intensity of transmission of Aedes-borne viruses is heterogeneous, and multiple factors can contribute to variation at small spatial scales. Illuminating drivers of heterogeneity in prevalence over time and space would provide information for public health authorities. The objective of this study is to detect the spatiotemporal clusters and determine the risk factors of three major Aedes-borne diseases, Chikungunya virus (CHIKV), Dengue virus (DENV), and Zika virus (ZIKV) clusters in Mexico. Methods We present an integrated analysis of Aedes-borne diseases (ABDs), the local climate, and the socio-demographic profiles of 2469 municipalities in Mexico. We used SaTScan to detect spatial clusters and utilize the Pearson correlation coefficient, Randomized Dependence Coefficient, and SHapley Additive exPlanations to analyze the influence of socio-demographic and climatic factors on the prevalence of ABDs. We also compare six machine learning techniques, including XGBoost, decision tree, Support Vector Machine with Radial Basis Function kernel, K nearest neighbors, random forest, and neural network to predict risk factors of ABDs clusters. Results DENV is the most prevalent of the three diseases throughout Mexico, with nearly 60.6% of the municipalities reported having DENV cases. For some spatiotemporal clusters, the influence of socio-economic attributes is larger than the influence of climate attributes for predicting the prevalence of ABDs. XGBoost performs the best in terms of precision-measure for ABDs prevalence. Conclusions Both socio-demographic and climatic factors influence ABDs transmission in different regions of Mexico. Future studies should build predictive models supporting early warning systems to anticipate the time and location of ABDs outbreaks and determine the stand-alone influence of individual risk factors and establish causal mechanisms.