Short-term effects of rainfall on childhood hand, foot and mouth disease and related spatial heterogeneity: evidence from 143 cities in mainland China
BMC Public Health volume 20, Article number: 1528 (2020)
Numerous studies have demonstrated the potential association between rainfall and hand, foot and mouth disease (HFMD), but the results are inconsistent. This study aimed to quantify the relationship between rainfall and HFMD based on a multicity study and explore the potential sources of spatial heterogeneity.
We retrieved the daily counts of childhood HFMD and the meteorological variables of the 143 cities in mainland China between 2009 and 2014. A common time series regression model was applied to quantify the association between rainfall and HFMD for each of the 143 cities. Then, we adopted the meta-regression model to pool the city-specific estimates and explore the sources of heterogeneity by incorporating city-specific characteristics.
The overall pooled estimation suggested a nonlinear exposure-response relationship between rainfall and HFMD. Once rainfall exceeded 15 mm, the HFMD risk stopped increasing linearly and began to plateau with the excessive risk ratio (ERR) peaking at 21 mm of rainfall (ERR = 3.46, 95% CI: 2.05, 4.88). We also found significant heterogeneity in the rainfall-HFMD relationships (I2 = 52.75%, P < 0.001). By incorporating the city-specific characteristics into the meta-regression model, temperature and student density can explain a substantial proportion of spatial heterogeneity with I2 statistics that decreased by 5.29 and 6.80% at most, respectively.
Our findings verified the nonlinear association between rainfall and HFMD. The rainfall-HFMD relationship also varies depending on locations. Therefore, the estimation of the rain-HFMD relationship of one location should not be generalized to another location.
Hand, foot and mouth disease (HFMD) is an acute highly contagious viral infection caused by enteroviruses and predominantly affects children under 5 years of age [1, 2]. While symptoms for most children is generally mild, few cases may experience more severe symptoms involving the central nervous system, like acute flaccid paralysis, encephalitis, meningitis, and even fatal cardiopulmonary failure . As a result of continuously large outbreaks and high rates of severe cases, East and Southeast Asia represent the worst hit areas around the world in terms of HFMD epidemics [4,5,6,7]. In mainland China, HFMD is the leading childhood infection among all notifiable diseases and responsible for millions of clinical cases  and greater than 75,000 age-weighted disability adjusted life years (DALYs) annually .
Multiple studies suggested that the water environment plays a crucial role in the survival of enterovirus (i.e., the pathogen of HFMD), thereby potentially affecting the transmission of HFMD. Water environments not only provide the major natural reservoir for enterovirus  but also protect enterovirus from diverse factors resulting in inactivation, such as extreme temperature, ultraviolet radiation, and salinity [11, 12]. In addition, the water environment can facilitate the attachment of enterovirus to suspended particles, which further prolongs enterovirus survival time . As the main source of natural water environments, rainfall was thus considered as an important driving factor of HFMD epidemics.
Although a growing number of epidemiological studies have been recently conducted to quantify the associations between rainfall and HFMD, evidence from previous studies remains inconsistent. Studies in both southern China  and Thailand  reported that the HFMD incidence was positively associated with rainfall. In contrast, studies in Hong Kong, Singapore, Japan, Taiwan , and Korea  suggested a nonsignificant association between rainfall and HFMD. A study in Vietnam  even suggested a negative association. The inconsistency of results among studies may originate from multiple sources. First, the diversity of methodologies and model choices adopted in different studies would inevitably affect the results. In fact, the same data could generate dissimilar conclusions when subject to different analyses . Second, the potential nonlinear exposure-response relationship could be another explanation for the inconsistencies [20, 21] because the estimated associations between rainfall and HFMD may vary depending on the level of rainfall. Finally, the rainfall-HFMD relationship could also be modified by location-specific characteristics [22, 23]. To date, the spatial heterogeneity of rainfall-HFMD relationships across studies remains poorly understood, thus hindering the understanding of potential mechanisms and the usage of rainfall information to predict HFMD epidemics.
To address the above research gap, we proposed a multicity two-stage time series study including 143 cities in mainland China. By estimating the associations between rainfall and HFMD for each of the 143 cities under a common modeling framework, we aimed to assess the spatial heterogeneity of rainfall-HFMD relationships and further explore the potential reasons.
Details of data used in this study have been described previously [22, 23]. Briefly, we retrieved surveillance data of HFMD cases and meteorological variables between from 1 January 2009 to 31 December 2014 in mainland China. We excluded data before 2009 due to the presence of less reliable data , and data after 2014 were excluded to reduce the potential influence of vaccines on the natural history of HFMD . More specifically, surveillance data of HFMD cases were collected from the China Information System for Disease Control and Prevention. Here, HFMD cases were defined as patients with vesicular or papular rashes on hands, feet, mouth or buttocks with or without fever . We further focused our analysis on children aged 0–12 years (accounted for greater than 99% of cases) to allow our results to be more clearly interpretable. Daily surveillance data of meteorological variables were collected from the China Meteorological Data Sharing System based on national ground meteorological stations. Meteorological variables included accumulated rainfall (from 20:00 to 20:00), mean relative humidity, mean temperature, sunshine hours, and mean air pressure.
Given that most of HFMD cases came from urban areas, we aggregated and matched the daily time series of HFMD and meteorological variables of prefectural-level cities (corresponds to medium to large size cities under China’s administrative structure) in mainland China. Confining our analysis to prefectural-level cities also assures better data quality. Finally, we included 143 cities with successful matching in our final analysis. To further explore the sources of spatial heterogeneity, city-specific characteristics were collected from the China city statistical yearbook , including economic characteristics (GDP per person and GDP growth rate), demographic characteristics (population density, population growth rate and student density), health resources (hospital beds and licensed physicians) and traffic (total travel passengers). To represent the climatic differences among cities, the arithmetic means of daily meteorological variables were used for each of the 143 cities as proxy variables.
We adopted a multicity two-stage time series design to robustly estimate rainfall-HFMD relationships and explore potential heterogeneity . In the first stage, we implemented a common time series regression model  for each of the 143 cities and obtained the city-specific estimations of rainfall-HFMD relationships. In the second stage, we pooled the city-specific rainfall-HFMD relationships using meta-analysis and examined whether city-specific characteristics can explain the spatial heterogeneity .
In summary, we applied a common quasi-Poisson time series regression with parametric splines to quantify the associations between rainfall and HFMD for each of the 143 cities. We used the quasi-Poisson distribution was assumed to account for overdispersion of disease counts. The natural cubic spline with 8 degrees of freedom (df) per year was used to control for the seasonality and long-term trends (i.e., time-varying unmeasured confounding) . We also included the indicators of weekdays and holidays and exponentially weighted moving averages of other meteorological variables (temperature, sunshine and air pressure) with the exception of relative humidity to remove residual confounding variables. Humidity was excluded because humidity and rainfall are highly correlated . Given that HFMD is an infectious disease, we also incorporated the autoregressive term of HFMD daily counts at lag 1 and 2 days in logarithmic form to match the transmission mechanism, which was informed by the classical Susceptible-Infectious-Recovery (SIR) model .
In previous studies, rainfall was incorporated into models as either a continuous [32, 33] or categorical variable . For fair comparisons with previous studies, we adopted rainfall as both continuous and categorical variables in this study, whereas the other model choices were identical. For the continuous scale, we applied the natural cubic spline with 5 df to characterize the potential non-linear exposure-response relationship between rainfall and HFMD, which was determined by varying a series of df (see more details in Additional file 3). For the categorical scale, we divided the daily accumulated rainfall into four levels based on the overall percentiles, including the reference group (rainfall ≤10th percentile or 0.07 mm), low rainfall group (rainfall between 10th and 50th percentiles or 0.07 ~ 1.35 mm), high rainfall group (rainfall between 50th and 90th percentiles or 1.35 ~ 7.10 mm) and extreme rainfall group (rainfall >90th percentiles or 7.10 mm) [14, 35]. Then, rainfall was incorporated as a dummy variable with delayed effects in the above time series regression. For both approaches, the same lag range between 0 to 14 lag days and exponentially weighted moving averages were used to characterize the delayed effects based on previous studies [32, 36, 37].
Corresponding to the two approaches used in the first-stage analysis, we applied different methods to obtain the overall rainfall-HFMD relationships depending on the variable types of rainfall. For continuous-scale rainfall, city-specific nonlinear exposure-response relationships between rainfall and HFMD were pooled by the multivariate meta-regression model with random intercepts  to obtain the overall exposure-response curve. For categorical-scale rainfall, city-specific estimates corresponding to different levels of rainfall (i.e., low, high and extreme rainfall versus reference) were pooled by three separate univariate meta-regression models with random intercepts  to obtain the overall estimates of each level of rainfall, separately. We also measured and tested the spatial heterogeneity of city-specific rainfall-HFMD relationships using I2 statistics and the Cochran Q test .
To further explore the reasons for potential heterogeneity, we incorporated the city-specific characteristics as meta-predictors into the above models with random intercepts. We included city-specific characteristics separately rather simultaneously (denoted as single predictor model) given that our main goal was to confirm potential effect modifiers rather than model predictions. Another added benefit of the single predictor model was the avoidance of collinearity issues due to high correlations among city-specific characteristics. Of note, the rainfall distribution was highly right skewed with sparse data in the high value range, which might lead to non-robust estimations of the right tail of exposure-response curves. Therefore, we only reported the results of meta-predictors for the categorical scale rainfall in the main text. Additional details of the corresponding results for continuous-scale rainfall are provided in Additional file 1 and Additional file 2, and these results were consistent with our main findings.
To exclude the possibility that the model choices may substantially affect our conclusions, we implemented additional sensitivity analysis regarding key model choices in the time series regression models. The key model choices included the df of the natural cubic splines used to characterize the rainfall-HFMD relationships, the lag range of the rainfall-HFMD relationships and alternative cut-off points of rainfall groups (5th, 50th, and 95th percentiles). Additional details of our sensitivity analysis results are reported in Additional file 3, Additional file 4 and Additional file 5. In general, our main conclusions were quite robust for numerous model choices.
Statistical analyses were performed using R software (version 3.5.3), specifically dlnm, mvmeta and metaphor packages. Geographic maps were generated using ArcGIS Pro software (version 2.4, authorization number: ESU934769606).
From 2009 to 2014, by including 143 cities in mainland China, we retrieved a total of 3,200,157 HFMD clinical cases in children under 12 years old. Our study area covered the majority of mainland Chinese with great diversity in geographic and climatic conditions. We found that HFMD circulated widely throughout the whole country. In general, HFMD incidences steadily increased from north to south. We particularly observed a high incidence area in the southernmost part of China where the annual HFMD incidences for multiple cities exceeded 1000 cases per 100,000 persons (Fig. 1a). Similarly, the spatial distribution of annual rainfall exhibited an increasing trend from north to south. The whole southern and southeastern coast of China experienced high rainfall levels (Fig. 1c). Regarding temporal trends, we observed an annual major peak in the time series of HFMD incidences, which typically occurred from May to August (Fig. 1b). The annual peak of rainfall roughly echoed that of HFMD incidence (Fig. 1d). Given that similar spatial-temporal patterns were observed between rainfall and HFMD, it was reasonable to assume that rainfall might be an environmental driver of the HFMD epidemic.
Based on continuous-scale rainfall data, we first pooled the city-specific estimates of exposure-response curves between rainfall and HFMD based on the multivariate meta-regression model (Fig. 2). The overall pooled estimates suggested a marginal effect on the rainfall-HFMD relationship with substantial heterogeneity among cities (I2 statistics = 52.75%, P < 0.001). We observed an approximately linear increase in the HFMD risk as rainfall initially increased. Once rainfall exceeded 15 mm, the HFMD risk began to plateau and was even reduced as rainfall further increased. However, we would prefer not to overinterpret the results of extreme rainfalls due to limited data. Overall, we observed a nonlinear association between rainfall and HFMD, and the excessive risk ratio peaked at 21-mm rainfall (ERR = 3.46, 95% CI: 2.05, 4.88).
For categorical-scale rainfall, city-specific estimates for different levels of rainfall were pooled separately by univariate meta-regression model (Fig. 3). The results were consistent with those obtained for continuous-scale rainfall. Compared to the reference group, the overall pooled estimate (measured by ERR) of low, high and extreme rainfall were 1.08% (95% CI: − 0.10, 2.26%), 2.05% (95% CI: 0.52, 3.60%) and 1.25% (95% CI: − 0.74, 3.28%), respectively (Fig. 3d). A steady increase in rainfall-HFMD associations were noted from the low to high rainfall group, whereas no further increase was observed in the extreme rainfall group. Based on the spatial distribution of city-specific estimates, we found that all groups showed similar spatial variations in rainfall-HFMD associations, and the highest values were clustered in south-central China (Fig. 3a-c). Similar to continuous-scale rainfall results, substantial heterogeneity was observed in categorical-scale rainfall results. In addition, an upward trend of heterogeneity was noted, and the I2 statistics increased from 24.98% in the low rainfall group to 45.66% in the extreme rainfall group.
We further examined the potential reasons for spatial heterogeneity in the rainfall-HFMD associations using the univariate meta regression which incorporates city-specific characteristics as explanatory variables. We found that temperature and student density could explain a substantial portion of heterogeneity regardless of rainfall groups. Both variables could reinforce the associations between rainfall and HFMD (Table 1). In the low and high rainfall groups, temperature explained the largest proportion of heterogeneity, and the I2 statistics decreased by 3.23 and 5.29%, respectively (Table 1). Compared to 10th percentile of temperature, the estimated ERR of the 90th percentile of temperature increased from − 0.76% (95% CI: − 2.47, 0.99%) to 3.52% (95% CI: 1.40, 5.69%) and from − 1.21% (95% CI: − 3.41, 1.05%) to 6.09% (95% CI: 3.44, 8.81%) in the low and high rainfall groups, respectively (Fig. 4). However, in the extreme rainfall group, student density explained a slightly larger proportion of heterogeneity compared with temperature, and the I2 statistics decreased by 6.80% (Table 1). Compared to the 10th percentile of student density, the estimated ERR of the 90th percentile of student density increased from − 2.25% (95% CI: − 4.78, 0.35%) to 4.71% (95% CI: 2.08, 7.40%) (Fig. 4).
To our knowledge, this is the first nation-wide study to examine rainfall- HFMD relationships and related heterogeneity. Based on the two-stage time series analysis consisting of 143 cities in mainland China, we observed a significant positive association between rainfall and HFMD overall. This finding was in consistent with the earlier hypothesis that rainfall facilitates enterovirus survival by providing a natural water environment and increasing the humidity of the air, which subsequently increase HFMD risk [39, 40].
However, it is worth noting that the strength of the relationship between rainfall and HFMD was much weaker than that of temperature and humidity reported in similar nation-wide studies [22, 23]. Unlike temperature and humidity, rainfall can be characterized with a spike at zero (i.e., its values at most days were zero), which could further cause the estimation of exposure-response relationship to progress toward a null value . In addition, rainfall also reduces host outdoor activity and effective contacts between susceptible and infectious individuals , which may subsequently hinder the spread of HFMD infection. Both of the above mechanisms could help to explain why many previous studies report a nonsignificant association between rainfall and HFMD and the relatively weaker estimation of an exposure-response relationship reported in this study.
Interestingly, we identified a nonlinear association between rainfall and HFMD. More specifically, we found that the HFMD risk stopped increasingly linearly and began to plateau once rainfall exceeded 15 mm. Our findings were consistent with other studies conducted in Singapore , south China [20, 32, 44] and Hong Kong , which also suggested that these associations would decrease with high levels of rainfall. We hypothesized that two forces drive the nonlinear behavior of the rainfall-HFMD relationship in opposite directions. Given the low level of rainfall, we hypothesized that the impacts of rainfall on host activity would be negligible; thus, the main impact of rainfall is facilitating enterovirus survival and increasing the HFMD risk. Nevertheless, given the high level of rainfall, we hypothesized that the main impact of rainfall would subsequently change to reducing host activity and hindering the spread of HFMD. These findings implied that weather could affect the HFMD risk via multiple pathways and complicated mechanisms.
After implementing a common time series model, we still observed significant spatial heterogeneity in the rainfall-HFMD relationships among different cities. In addition, we found that the heterogeneity of city-specific estimates increased as rainfall increased. As stated earlier, the distribution of rainfall data is highly right skewed with most data clustered in the low rainfall groups. Thus, estimates in the high rainfall lack statistical power and exhibit larger variations compared with the low rainfall group, which may subsequently explain the increased heterogeneity in the high rainfall group .
Furthermore, we found that temperature and student density explain a substantial portion of heterogeneity regardless of the rainfall level. Infectious disease epidemics are mainly determined by a sufficient number of susceptible and infectious individuals as well as effective contacts between susceptible and infectious individuals . Regarding temperature, previous studies suggested that warm temperatures can simultaneously promote the excretion of enterovirus  and host activity , which can increase the concentration of infectious agents in the environment and effective contacts, respectively. Regarding student density, given that children are most vulnerable to the HFMD infections , high student density reflects a high proportion of susceptible individuals in the population, which can directly raise the effective contacts between susceptible and infectious individuals [49, 50]. Both abovementioned mechanisms contribute to the spread of HFMD, thereby reinforcing the associations between rainfall and HFMD. The robust results among different levels of rainfall further verified our hypothesis.
Unlike most previous studies that typically applied a single site analysis, this study included 143 cites in mainland China with great variations in climate, economic development, and health resources. The large sample size and significant diversity allow us to robustly estimate the rainfall-HFMD relationship and thoroughly examine potential heterogeneity. However, several limitations should also be noted and improved in further studies. First, given the time series study design, this study is vulnerable to the issue of ecological fallacies . Although this study may be weak in terms of causal reference, the associations identified in this study can still provide important information on public health practices, like the weather-based early warning system for contagious disease. Second, a large portion of heterogeneity remained unexplained despite the incorporation of a wide range of city-specific characteristics into the meta-regression model. Limited by data accessibility, some other important city-specific characteristics should be taken into consideration in future studies, such as public health interventions, dominant subtypes of enterovirus, vegetation coverage, and soil features. Third, for the sake of simplicity, we did not take the spatial correlation into account when pooling city-specific estimates . While the spatial correlation can be eliminated through incorporating city-specific characteristics into the meta-regression models, we believe that the impacts of spatial correlation on our conclusions should be negligible.
In conclusion, this study verified the nonlinear association between rainfall and HFMD. These findings can help us to elucidate the complicated mechanisms by which rainfall is related to HFMD and provide insights into the establishment of weather-based early warning systems. Especially in the context of global climate change, the distribution of rainfall might experience substantial changes in the near future . In addition, we also found that the rainfall-HFMD relationship is profoundly affected by local conditions. Therefore, the estimation of the rain-HFMD relationship of one location should not be generalized to another location.
Availability of data and materials
The HFMD surveillance data are not publicly available, but access to this data could require permission from the Chinese Center for Disease Control and Prevention or from the corresponding author on reasonable request. The meteorological monitoring data and city-specific characteristics are publicly available, which can be obtained from the China Meteorological Data Sharing System (http://data.cma.cn) and the China city statistical yearbook respectively.
Hand, foot and mouth disease
Age-weighted disability adjusted life years
Excessive risk ratio
Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, Chang Z, Liu F, Fang VJ, Zheng Y, et al. Hand, foot, and mouth disease in China, 2008-12: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18. https://doi.org/10.1016/S1473-3099(13)70342-6.
Solomon T, Lewthwaite P, Perera D, Cardosa MJ, McMinn P, Ooi MH. Virology, epidemiology, pathogenesis, and control of enterovirus 71. Lancet Infect Dis. 2010;10(11):778–90. https://doi.org/10.1016/S1473-3099(10)70194-8.
A Guide to Clinical Management and Public Health Response for Hand, Foot and Mouth Disease (HFMD) [https://www.who.int/vietnam/health-topics/hand-foot-and-mouth-disease-(hfmd)/publications].
Han Z, Zhang Y, Huang K, Wang J, Tian H, Song Y, Yang Q, Yan D, Zhu S, Yao M, et al. Two Coxsackievirus B3 outbreaks associated with hand, foot, and mouth disease in China and the evolutionary history worldwide. BMC Infect Dis. 2019;19(1):466. https://doi.org/10.1186/s12879-019-4107-z.
Liu MY, Liu J, Lai W, Luo J, Liu Y, Vu GP, Yang Z, Trang P, Li H, Wu J. Characterization of enterovirus 71 infection and associated outbreak of hand, foot, and mouth disease in Shawo of China in 2012. Sci Rep. 2016;6:38451. https://doi.org/10.1038/srep38451.
Van Tu P, Thao NTT, Perera D, Truong KH, Tien NTK, Thuong TC, How OM, Cardosa MJ, McMinn PC. Epidemiologic and virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005. Emerg Infect Dis. 2007;13(11):1733–41. https://doi.org/10.3201/eid1311.070632.
Puenpa J, Chieochansin T, Linsuwanon P, Korkong S, Thongkomplew S, Vichaiwattana P, Theamboonlers A, Poovorawan Y. Hand, foot, and mouth disease caused by coxsackievirus A6, Thailand, 2012. Emerg Infect Dis. 2013;19(4):641–3. https://doi.org/10.3201/eid1904.121666.
Overview of the national notifiable communicable disease epidemic from 2016 to 2018 [http://www.nhc.gov.cn/jkj/s2907/new_list.shtml?tdsourcetag=s_pcqq_aiomsg].
Koh WM, Badaruddin H, La H, Chen MI, Cook AR. Severity and burden of hand, foot and mouth disease in Asia: a modelling study. BMJ Glob Health. 2018;3(1):e000442. https://doi.org/10.1136/bmjgh-2017-000442.
Rajtar B, Majek M, Polanski L, Polz-Dacewicz M. Enteroviruses in water environment--a potential threat to public health. Ann Agric Environ Med. 2008;15(2):199–203. https://doi.org/10.1159/000121402.
Pohjanpelto P. Response of enteroviruses to cystine. Virology. 1961;15:225–30. https://doi.org/10.1016/0042-6822(61)90352-x.
Rzezutka A, Cook N. Survival of human enteric viruses in the environment and food. FEMS Microbiol Rev. 2004;28(4):441–53. https://doi.org/10.1016/j.femsre.2004.02.001.
Melnick JL, Gerba CP, Berg G. The ecology of enteroviruses in natural waters. C R C Critical Rev Environ Control. 2009;10(1):65–93. https://doi.org/10.1080/10643388009381677.
Cheng J, Wu J, Xu Z, Zhu R, Wang X, Li K, Wen L, Yang H, Su H. Associations between extreme precipitation and childhood hand, foot and mouth disease in urban and rural areas in Hefei, China. Sci Total Environ. 2014;497:484–90. https://doi.org/10.1016/j.scitotenv.2014.08.006.
Samphutthanon R, Tripathi NK, Ninsawat S, Duboz R. Spatio-temporal distribution and hotspots of hand, foot and mouth disease (HFMD) in northern Thailand. Int J Environ Res Public Health. 2013;11(1):312–36. https://doi.org/10.3390/ijerph110100312.
Lee C-CD, Tang J-H, Hwang J-S, Shigematsu M, Chan T-C. Effect of meteorological and geographical factors on the epidemics of hand, foot, and mouth disease in island-type territory, East Asia. Biomed Res Int. 2015. https://doi.org/10.1155/2015/805039.
Kim BI, Ki H, Park S, Cho E, Chun BC. Effect of Climatic Factors on Hand, Foot, and Mouth Disease in South Korea, 2010-2013. PLoS One. 2016;11(6). https://doi.org/10.1371/journal.pone.0157500.
Phung D, Nguyen HX, Nguyen HLT, Do CM, Tran QD, Chu C. Spatiotemporal variation of hand-foot-mouth disease in relation to socioecological factors: a multiple-province analysis in Vietnam. Sci Total Environ. 2018;610-611:983–91. https://doi.org/10.1016/j.scitotenv.2017.08.158.
Dominici F. Invited commentary: air pollution and health-what can we learn from a hierarchical approach? Am J Epidemiol. 2002;155(1):11–5; discussion 16. https://doi.org/10.1093/aje/155.1.11.
Huang D-C, Wang J-F. Monitoring hand, foot and mouth disease by combining search engine query data and meteorological factors. Sci Total Environ. 2018;612:1293–9. https://doi.org/10.1016/j.scitotenv.2017.09.017.
Liao J, Yu S, Yang F, Yang M, Hu Y, Zhang J. Short-Term Effects of Climatic Variables on Hand, Foot, and Mouth Disease in Mainland China, 2008-2013: A multilevel spatial Poisson regression model accounting for Overdispersion. PLoS One. 2016;11(1). https://doi.org/10.1371/journal.pone.0147054.
Xiao X, Gasparrini A, Huang J, Liao Q, Liu F, Yin F, Yu H, Li X. The exposure-response relationship between temperature and childhood hand, foot and mouth disease: a multicity study from mainland China. Environ Int. 2017;100:102–9. https://doi.org/10.1016/j.envint.2016.11.021.
Bo Z, Ma Y, Chang Z, Zhang T, Liu F, Zhao X, Long L, Yi X, Xiao X, Li Z. The spatial heterogeneity of the associations between relative humidity and pediatric hand, foot and mouth disease: evidence from a nation-wide multicity study from mainland China. Sci Total Environ. 2020;707:136103. https://doi.org/10.1016/j.scitotenv.2019.136103.
Venkatesh S, Memish ZA. Bioterrorism--a new challenge for public health. Int J Antimicrob Agents. 2003;21(2):200–6. https://doi.org/10.1016/s0924-8579(02)00366-7.
Announcement on licensed drugs approved by China Food and Drug Administration (No.4 in 2016). [https://www.nmpa.gov.cn/xxgk/kpzhsh/kpzhshyp/20160414124101246.html].
National guidelines for HFMD prevention and control 2008 [http://www.nhc.gov.cn/jkj/s3577/200805/e73df45b7b1549188b1d4e1efd604da9.shtml].
Chen X: China City Statistical Yearbook [in Chinese]. http://data.cnki.net/yearbook/Single/N2019070173: China Statistics Press, Beijing; 2009-2014.
Peng and Dominici. Statistical Methods for Environmental Epidemiology with R: A Case Study in Air Pollution and Health: Springer Science & Business Media; 2008. https://link.springer.com/book/10.1007/978-0-387-78167-9.
Gasparrini A, Armstrong B, Kenward MG. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat Med. 2012;31(29):3821–39. https://doi.org/10.1002/sim.5471.
Huang Y, Deng T, Yu S, Gu J, Huang C, Xiao G, Hao Y. Effect of meteorological variables on the incidence of hand, foot, and mouth disease in children: a time-series analysis in Guangzhou, China. BMC Infect Dis. 2013;13. https://doi.org/10.1186/1471-2334-13-134.
Imai C, Armstrong B, Chalabi Z, Mangtani P, Hashizume M. Time series regression model for infectious disease and weather. Environ Res. 2015;142:319–27. https://doi.org/10.1016/j.envres.2015.06.040.
Yan S, Wei L, Duan Y, Li H, Liao Y, Lv Q, Zhu F, Wang Z, Lu W, Yin P, et al. Short-Term Effects of Meteorological Factors and Air Pollutants on Hand, Foot and Mouth Disease among Children in Shenzhen, China, 2009–2017. Int J Environ Res Public Health. 2019;16(19). https://doi.org/10.3390/ijerph16193639.
Chen C, Lin H, Li X, Lang L, Xiao X, Ding P, He P, Zhang Y, Wang M, Liu Q. Short-term effects of meteorological factors on children hand, foot and mouth disease in Guangzhou, China. Int J Biometeorol. 2014;58(7):1605–14. https://doi.org/10.1007/s00484-013-0764-6.
Yu G, Li Y, Cai J, Yu D, Tang J, Zhai W, Wei Y, Chen S, Chen Q, Qin J. Short-term effects of meteorological factors and air pollution on childhood hand-foot-mouth disease in Guilin, China. Sci Total Environ. 2019;646:460–70. https://doi.org/10.1016/j.scitotenv.2018.07.329.
Bush KF, O'Neill MS, Li S, Mukherjee B, Hu H, Ghosh S, Balakrishnan K. Associations between extreme precipitation and gastrointestinal-related hospital admissions in Chennai, India. Environ Health Perspect. 2014;122(3):249–54. https://doi.org/10.1289/ehp.1306807.
Zhang Z, Xie X, Chen X, Li Y, Lu Y, Mei S, Liao Y, Lin H. Short-term effects of meteorological factors on hand, foot and mouth disease among children in Shenzhen, China: non-linearity, threshold and interaction. Sci Total Environ. 2016;539:576–82. https://doi.org/10.1016/j.scitotenv.2015.09.027.
Wang P, Zhao H, You F, Zhou H, Goggins WB. Seasonal modeling of hand, foot, and mouth disease as a function of meteorological variations in Chongqing, China. Int J Biometeorol. 2017;61(8):1411–9. https://doi.org/10.1007/s00484-017-1318-0.
Thompson SG, Higgins JP. How should meta-regression analyses be undertaken and interpreted? Stat Med. 2002;21(11):1559–73. https://doi.org/10.1002/sim.1187.
Zhang Q, Zhou M, Yang Y, You E, Wu J, Zhang W, Jin J, Huang F. Short-term effects of extreme meteorological factors on childhood hand, foot, and mouth disease reinfection in Hefei, China: a distributed lag non-linear analysis. Sci Total Environ. 2019;653:839–48. https://doi.org/10.1016/j.scitotenv.2018.10.349.
Jiang FC, Yang F, Chen L, Jia J, Han YL, Hao B, Cao GW. Meteorological factors affect the hand, foot, and mouth disease epidemic in Qingdao, China, 2007-2014. Epidemiol Infect. 2016;144(11):2354–62. https://doi.org/10.1017/s0950268816000601.
Lorenz E, Jenkner C, Sauerbrei W, Becher H. Modeling variables with a spike at zero: examples and practical recommendations. Am J Epidemiol. 2017;185(8):650–60. https://doi.org/10.1093/aje/kww122.
Belanger M, Gray-Donald K, O'Loughlin J, Paradis G, Hanley J. Influence of weather conditions and season on physical activity in adolescents. Ann Epidemiol. 2009;19(3):180–6. https://doi.org/10.1016/j.annepidem.2008.12.008.
Hii YL, Rocklov J, Ng N. Short Term Effects of Weather on Hand, Foot and Mouth Disease. PLoS One. 2011;6(2). https://doi.org/10.1371/journal.pone.0016796.
Zhang W, Du Z, Zhang D, Yu S, Hao Y. Boosted regression tree model-based assessment of the impacts of meteorological drivers of hand, foot and mouth disease in Guangdong, China. Sci Total Environ. 2016;553:366–71. https://doi.org/10.1016/j.scitotenv.2016.02.023.
Wang P, Goggins WB, Chan EY. Hand, foot and mouth disease in Hong Kong: a time-series analysis on its relationship with weather. PLoS One. 2016;11(8):e0161006. https://doi.org/10.1371/journal.pone.0161006.
Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002;21(11):1539–58. https://doi.org/10.1002/sim.1186.
Woolhouse M. Modeling infectious diseases in humans and animals. Nature. 2008;453(7191):34. https://doi.org/10.1038/453034a.
Fong TT, Lipp EK. Enteric viruses of humans and animals in aquatic environments: health risks, detection, and potential water quality assessment tools. Microbiol Mol Biol Rev. 2005;69(2):357–71. https://doi.org/10.1128/MMBR.69.2.357-371.2005.
Wang Y, Feng Z, Yang Y, Self S, Gao Y, Longini IM, Wakefield J, Zhang J, Wang L, Chen X, et al. Hand, foot, and mouth disease in China: patterns of spread and transmissibility. Epidemiology. 2011;22(6):781–92. https://doi.org/10.1097/EDE.0b013e318231d67a.
Bo Y-C, Song C, Wang J-F, Li X-W. Using an autologistic regression model to identify spatial risk factors and spatial risk patterns of hand, foot and mouth disease (HFMD) in mainland China. BMC Public Health. 2014;14. https://doi.org/10.1186/1471-2458-14-358.
Openshaw S. Ecological fallacies and the analysis of areal census data. Environ Plan A. 1984;16(1):17–31. https://doi.org/10.1068/a160017.
Putnam AE, Broecker WS. Human-induced changes in the distribution of rainfall. Sci Adv. 2017;3(5):e1600871. https://doi.org/10.1126/sciadv.1600871.
We specially thank the Chinese Center for Disease Control and Prevention for providing the data.
This study was primarily sponsored by the National Natural Science Foundation of China (grant number 81903415), the Postdoctoral Research Foundation of Sichuan University (grant number 2018SCU12014) and Key R&D Project of Sichuan Province Science and Technology Support Program (grant number 2020YFS0215). Yue Ma was supported by the National Natural Science Foundation of China (grant number 81803332). Xing Zhao was supported by the National Natural Science Foundation of China (grant number 81773548). The funding body has not been involved with the study design nor the conduct of the study.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Multivariate meta-regression models by incorporating city-specific characteristics to explain heterogeneity.
The predicted rainfall-HFMD relationship in ERR for the 10th and 90th percentiles of temperature and student density.
The overall pooled estimates of HFMD-rainfall relationships for different choice of df (the df of natural cubic splines used to characterized the exposure-response relationship).
The overall pooled estimates of HFMD-rainfall relationships for different choice of lag range with the max lag time ranged from 7 to 21 days.
The pooled estimates (measured by ERR) of low, high and extreme rainfall with alternative cut-off points (5th, 50th, 95th percentile as cut-off points).
About this article
Cite this article
Yang, F., Ma, Y., Liu, F. et al. Short-term effects of rainfall on childhood hand, foot and mouth disease and related spatial heterogeneity: evidence from 143 cities in mainland China. BMC Public Health 20, 1528 (2020). https://doi.org/10.1186/s12889-020-09633-1