Time series analysis of mumps and meteorological factors in Beijing, China
BMC Infectious Diseases volume 19, Article number: 435 (2019)
Over the past decades there have been outbreaks of mumps in many countries, even in populations that were vaccinated. Some studies suggest that the incidence of mumps is related to meteorological changes, but the results of these studies vary in different regions. To date there is no reported study on correlations between mumps incidence and meteorological parameters in Beijing, China.
A time series analysis incorporating selected weather factors and the number of mumps cases from 1990 to 2012 in Beijing was performed. First, correlations between meteorological variables and the number of mumps cases were assessed. A seasonal autoregressive integrated moving average model with explanatory variables (SARIMAX) was then constructed to predict mumps cases.
Mean temperature, rainfall, relative humidity, vapor pressure, and wind speed were significantly associated with mumps incidence. After constructing the SARIMAX model, mean temperature at lag 0 (β = 0.016, p < 0.05, 95% confidence interval 0.001 to 0.032) was positively associated with mumps incidence, while vapor pressure at lag 2 (β = ˗0.018, p < 0.05, 95% confidence interval ˗0.038 to ˗0.002) was negatively associated. SARIMAX (1, 1, 1) (0, 1, 1)12 with temperature at lag 0 was the best predictive construct.
The incidence of mumps in Beijing from 1990 to 2012 was significantly correlated with meteorological variables. Combining meteorological variables, a predictive SARIMAX model that could be used to preemptively estimate the incidence of mumps in Beijing was established.
Mumps is a viral respiratory disease that is most likely to occur in children and adolescents. Parotid non-suppurative inflammation and painful swelling of the parotid gland are the main clinical features of mumps. The complications of mumps include meningoencephalitis, meningitis, orchitis, pancreatitis, and ovarian inflammation .
The mumps virus is a single-stranded RNA virus of the paramyxovirus family. It is a moderately to highly contagious virus that only infects humans. The main routes of transmission include direct contact, droplet propagation, and contact with contaminated objects . Administration of the well-established live attenuated vaccine is the primary measure used to prevent mumps.
As at November 2016, of the 192 World Health Organization member states 127 (57%) included mumps vaccine in their national vaccination schedules . Mumps-containing vaccines were introduced in China in 1990. Since 2008 the measles, mumps, and rubella (MMR) vaccine has been included in China’s national immunization programs . There have been outbreaks of mumps in many countries including the United States , Canada , Czech Republic , Belgium , Portugal , and Serbia , even in populations that were vaccinated . A total of 698,092 cases of mumps were reported in mainland China from 2013 to 2015, with an average annual incidence of 17.2 per 100,000 .
Climate change plays an important role in the spread of many infectious diseases, especially vector-borne and water-borne infectious diseases . As early as 2500 years ago, Hippocrates observed the influences of climate change on gastrointestinal infections, tuberculosis, and central nervous system infections . In some recent studies the onset of some infectious diseases including mumps was associated with specific changes in meteorological factors. In a study conducted in Guangzhou in China the incidence of mumps was positively correlated with mean temperature, relative humidity, and wind velocity, and negatively correlated with atmospheric pressure . In a study in Taiwan, which is at a similar latitude to Guangzhou, the occurrence of mumps was significantly correlated with increased temperature and vapor pressure . In a study in Fukuoka Prefecture in Japan the number of pediatric mumps cases increased significantly with increased average temperature and relative humidity . In Jining in China, mumps occurrence was reportedly positively associated with temperature, wind speed, and sunshine duration, and negatively associated with relative humidity . However, large-scale weather changes collectively incorporated into the so called North Atlantic Oscillation phenomenon were reportedly not crucial factors in fluctuations in annual mumps incidence rates in the Czech Republic . The studies described above indicate similarities in different regions that may be related to the locations and climates of the study areas. To date there is no reported study on correlations between mumps incidence and meteorological parameters in Beijing, China.
The purpose of the present study was to assess correlations between weather factors and the incidence of mumps in Beijing, and to establish an accurate model for estimating epidemic trends pertaining to mumps.
Beijing is the capital of China. The city is located in the northern part of the vast North China Plain (39.9° N, 116.3° E), and it is situated in a zone of typical continental monsoonal climate with four clearly distinct seasons. Spring is windy, summer is hot and rainy, autumn is dry, and winter is cold .
Mumps has been a legally notifiable disease in China since September 1989. The monthly mumps data used in this study were obtained from the Chinese Center for Disease Control and Prevention. The number of mumps cases recorded from January 1990 to December 2012 was 197,726, and the diagnoses were based on the criteria used by the National Health and Family Planning Commission of the People’s Republic of China (formerly the Ministry of Health of the People’s Republic of China).
Daily temperature, rainfall, relative humidity, vapor pressure, and wind speed data from 1989 to 2012 were provided by the Beijing Meteorological Bureau. Monthly means of the daily average values of these meteorological characteristics were calculated. The data used in this study are provided as Additional files 1.
A descriptive analysis was conducted to assess the distribution of mumps cases and weather factors in Beijing. Seasonal autoregressive integrated moving average (SARIMA) models were then used to evaluate relationships between monthly numbers of mumps cases and meteorological parameters. SARIMA models were optimal for use in this study because seasonal and non-seasonal trends could be studied . A SARIMA model is described as an autoregressive integrated moving average, p, d, and q, multiplied by P, D, and Q—where the non-seasonal parameters are the number of autoregressive terms (p), the number of differences (d), and the moving average (q), and the seasonal parameters are the number of seasonal autoregressive terms (P), the number of seasonal differences (D), and the seasonal moving average (Q).
The SARIMA model with explanatory variables (SARIMAX) extends the capability of the SARIMA model by integrating external information such as rainfall, wind speed, and other meteorological variables into a time series model . In the present study a SARIMAX model was constructed to investigate mumps cases. The specific method used to generate the model is described below.
The first step was stabilization processing of the sequence. The data were processed by difference and conversion if the sequence was unstable or had a seasonal distribution. The second step was model identification. The order of p, P, q, and Q in the model was determined based on a graph of the autocorrelation function (ACF) and the partial autocorrelation function (PACF). For the pure autoregressive model (p), moving average model (q), and autoregressive moving average model (p, q), the parameter q could not exceed the lag of the ACF, and the parameter p did not exceed the lag of the PACF. The orders of d and D respectively represented the non-seasonal and seasonal difference times . The model parameters were then estimated and verified. The maximum likelihood method and t-tests were used to estimate and test the model parameters. Akaike information criterion (AIC) values were used to measure the model fit. Smaller AIC values indicated a better model fit.
Model diagnostics were then performed. The Ljung-Box Q test was conducted to ascertain whether the residual series were random. A p value less than 0.05 suggested that the residual sequence was not white noise and that it contained information that was inadequately extracted. We then assessed the correlations between pairs of sequences with strong autocorrelations. Pre-whitened data were needed to separate the linear associations from their autocorrelations. Cross-correlation function (CCF) plots were used to evaluate relationships between the number of mumps cases and meteorological factors, and determine which covariates and lags were best for the model . Only covariates that had significant parameter estimates and lowered the AIC value were selected . Lastly, we incorporated the covariates into the model and repeated model parameter estimation/verification and model diagnostics to build the best SARIMAX model. We used the 1990–2010 data to construct the best SARIMA model and SARIMAX models. We then used the 2011–2012 data to assess the predictive capacity of the SARIMAX model. Data analysis was conducted using R 3.3.1 .
There were 197,726 reported mumps cases from 1990 to 2012 in Beijing. As shown in Fig. 1, the peak incidence occurred in 1992 with 38,979 cases, and the lowest incidence occurred in 2003 with 1579 cases. The seasonal variation is clearly depicted in Fig. 2. There was a major peak in the number of mumps cases in the late spring and early summer (April to July) and a minor peak in winter (December to January) during the years included in this study. More descriptive statistics pertaining to meteorological factors and numbers of mumps cases are shown in Table 1.
SARIMA model analysis
A SARIMA model with 252 monthly data-points for mumps cases from 1990 to 2010 without any covariates was developed first. Mumps cases fluctuated within a large range over the time-course (Fig. 1). A logarithmic transformation of the time series of mumps cases was performed to stabilize fluctuations in the data. As the overall logarithmically transformed mumps data exhibited a downward trend and an obvious seasonal distribution, 1-step non-seasonal and 1-step seasonal differences were used separately. The value of both d and D was 1. The ACF and PACF of mumps cases are shown in Fig. 3. The ACF values of lag 2, 12, 14, 26, and 34 exceeded the critical value. The ACF value of lag 14, the neighbor of seasonal lag 12, was caused by the cross effect of the seasonal and non-seasonal autocorrelation. The significant ACF values of lag 26 and 34 may pertain to the presence of a year effect, though 26 months and 34 months are not strictly 2 years or 3 years. Therefore, the respective maximum values of the seasonal parameter Q and the non-seasonal parameter q were 1 and 2. Similarly, the sample PACF values were significant at lag 2, 5, 12, 14, 22, 24, and 36 (Fig. 3), so the respective maximum values of the seasonal parameter P and the non-seasonal parameter p were 3 and 5. We assumed that the maximum values of both P and p were 2 to make the model concise.
We searched all 54 SARIMA models that satisfied the conditions p ≤ 2, P ≤ 2, q ≤ 2, and Q ≤ 1 to find the most suitable model. Only two models yielded statistically significant parameters. Table 2 shows the results for these two models, Model A (SARIMA  [1, 1, 0]12) and Model B (SARIMA  [0, 1, 1]12). The AIC value of Model B (− 63.96) was lower than that of Model A (− 33.22), therefore Model B fit the data better.
SARIMAX model analysis and prediction
The CCF was used to investigate relationships between meteorological factors and mumps cases. The fitted SARIMA model was applied to pre-whiten the data for the monthly averages of the daily mean values of the meteorological factors. Figure 4 shows the cross-correlations between the pre-whitened meteorological variables (temperature, vapor pressure, rainfall, wind speed, and relative humidity) and logarithmically transformed numbers of mumps cases (logmumps) at lags of 1 to 6 months. All of the weather factors except rainfall were significantly associated with logmumps at at least some of the lags. For example, the CCF for vapor pressure and logmumps was significant at lags 1, 2, 4, and 6, and the CCF for wind speed and logmumps was significant at lags 5 and 6. Those significant weather factors were added into the SARIMA model as covariates to establish the SARIMAX model. As shown in Table 3, two SARIMAX models with covariates yielded significant parameters and lowered the AIC value. This result indicated that mean temperature at lag 0 and vapor pressure at lag 2 affected logmumps after fitting the time series regression model. Mean temperature at lag 0 (β = 0.016, p < 0.05, 95% confidence interval 0.001 to 0.032) was positively associated with logmumps, while vapor pressure at lag 2 (β = − 0.018, p < 0.05, 95% confidence interval ˗0.038 to ˗0.002) was negatively associated. SARIMAX (1, 1, 1) (0, 1, 1)12 with temperature at lag 0 was the optimal model with the lowest AIC value (Table 4).
The SARIMAX model described above was used to attempt to retrospectively predict mumps cases from January 2011 to December 2012. The estimated and predicted results are shown in Fig. 5. The predicted monthly numbers of mumps cases all fell within the confidence intervals.
The incidences of many infectious diseases exhibit seasonal variation. In the present study the incidence of mumps cases in Beijing over a time-course exhibited clear seasonal effects. Cases peaked from late spring to early summer (April–July), and in winter (December–January). This is consistent with a study conducted by Li et al.  in Jining, China, in which large peaks were found in May and June, and in winter. In a study of the epidemiology of mumps conducted in China by Cui et al.  most mumps cases occurred between April and July, with a small peak occurring in November and December. Although the mechanisms underlying the seasonality of mumps incidence remain poorly understood, oscillatory changes in infectiousness, contact patterns, pathogen survival, host susceptibility, and population behaviors may contribute to the phenomenon [26, 27]. Seasonal variations in meteorological factors probably also play a role.
The results of the present study are similar to those of several previous studies investigating the effects of weather variables on mumps in Asia [15,16,17,18]. In all of those studies the occurrence of mumps was significantly associated with mean monthly temperature. In several studies there was an approximately linear association between mean temperature and mumps cases, over a certain temperature threshold. For example in Jining, a city in northern China, each 1 °C increase in mean temperature above 4 °C was associated with a 2.72% increase in the risk of mumps . In Taiwan the number of mumps cases started to increase at a temperature of 20 °C, but began to decline when the temperatures exceeded approximately 25 °C . Mumps virus is stable for days at 4 °C . Higher temperatures are more conducive to the survival of mumps virus in the environment, and person-to-person contact . Furthermore, in a study investigating correlations between meteorological conditions and physical activities performed in open-air settings, the number of individuals walking on a public track increased with temperature , which may also be indicative of increased outdoor activities more broadly. Partaking in frequent outdoor activities may increase the risk of mumps infection. In another study conducted in the United States, spring-break college travel was associated with an increase in mumps cases after 01 April .
In the present study vapor pressure at lag 2 was negatively associated with logmumps. In a study conducted in Taiwan the number of mumps cases began to increase at vapor pressures of 5–9 hPa and decreased at vapor pressures > 25–29 hPa . The mechanism by which vapor pressure affects the transmission of mumps virus is poorly understood. Additional studies investigating the underlying mechanisms are warranted.
The survival of viruses outside the host depends partially on relative humidity. Viruses with lipid envelopes survive longer at lower relative humidity (20–30%) . In addition, at higher wind speeds the spread of disease via respiratory droplets is rendered more likely . In the present study, logmumps was significantly correlated with relative humidity and with wind speed at different lags. The varying lag effects associated with weather parameters reported in other studies probably resulted from differences in study locations.
In the present study mean monthly temperature, relative humidity, vapor pressure, and wind speed at different lags were significantly associated with logmumps, with lag effects varying from 0 months to 6 months. After fitting the SARIMAX model, mean temperature at lag 0 was positively associated with logmumps, while vapor pressure at lag 2 was negatively associated. SARIMAX (1, 1, 1) (0, 1, 1)12 with temperature at lag 0 was the optimal model with the highest prediction accuracy, which overcame the hypothesis that the traditional time series model was linearly dependent on the variables included, and improved the accuracy of resulting predictions. The model was established based on mumps incidences and meteorological data in Beijing, China. Accordingly, it is only suitable for use to predict overall trends in Beijing.
The current study had some limitations. More meteorological parameters such as monthly mean, maximum, and minimum temperatures, sunshine duration, and other variables should be included in future studies to comprehensively evaluate relationships between meteorological parameters and mumps. Another limitation was that the mumps incidence data were only available on a per month basis. Weekly or daily incidence data may decrease the accuracy of lagged time estimation. Lastly, potentially confounding variables such as vaccine usage and school and household size may have affected mumps incidence. These data were not available for assessment in the current study.
Various meteorological variables influence the incidence of mumps in Beijing, China. A time series regression model suggested that mean temperature at lag 0 and vapor pressure at lag 2 may influence logmumps. The utilization of a SARIMAX (1, 1, 1) (0, 1, 1)12 model with temperature at lag 0 is recommended for predicting the incidence of mumps in Beijing.
Akaike information criterion
- logmumps :
logarithmically transformed numbers of mumps cases
partial autocorrelation function
seasonal autoregressive integrated moving average
seasonal autoregressive integrated moving average model with explanatory variables
Hviid A, Rubin S, Mühlemann K. Mumps. Lancet. 2008;371(9616):932–44.
MacDonald N, Hatchette T, Elkout L, Sarwal S. Mumps is back: why is mumps eradication not working? Adv Exp Med Biol. 2011;697:197–220.
Vaccine preventable diseases monitoring system. http://www.who.int/immunization_monitoring/en/globalsummary/scheduleselect.cfm.]
Zheng J, Zhou Y, Wang H, Liang X. The role of the China experts advisory committee on immunization program. Vaccine. 2010;28(Suppl 1):A84–7.
Majumder MS, Nguyen CM, Cohn EL, Hswen Y, Mekaru SR, Brownstein JS. Vaccine compliance and the 2016 Arkansas mumps outbreak. Lancet Infect Dis. 2017;17(4):361–2.
Trotz-Williams LA, Mercer NJ, Paphitis K, Walters JM, Wallace D, Kristjanson E, Gubbay J, Mazzulli T. Challenges in interpretation of diagnostic test results in a mumps outbreak in a highly vaccinated population. Clinical and vaccine immunology: CVI. 2017;24:2.
Pazdiora P, Skalova J, Kubatova A, Jezova I, Moravkova I, Podlesna I, Pruchova J, Spacilova M, Svecova M: [Mumps outbreak in the Plzen Region in 2011]. Epidemiologie, mikrobiologie, imunologie : casopis Spolecnosti pro epidemiologii a mikrobiologii Ceske lekarske spolecnosti JE Purkyne 2015, 64(4):242–249.
Braeye T, Linina I, De Roy R, Hutse V, Wauters M, Cox P, Mak R. Mumps increase in Flanders, Belgium, 2012-2013: results from temporary mandatory notification and a cohort study among university students. Vaccine. 2014;32(35):4393–8.
Cordeiro E, Ferreira M, Rodrigues F, Palminha P, Vinagre E, Pimentel JP. Mumps outbreak among highly vaccinated teenagers and children in the central region of Portugal, 2012-2013. Acta Medica Port. 2015;28(4):435–41.
Nedeljkovic J, Kovacevic-Jovanovic V, Milosevic V, Seguljev Z, Petrovic V, Muller CP, Hubschen JM. A mumps outbreak in Vojvodina, Serbia, in 2012 underlines the need for additional vaccination opportunities for young adults. PLoS One. 2015;10(10):e0139815.
Whelan J, van Binnendijk R, Greenland K, Fanoy E, Khargi M, Yap K, Boot H, Veltman N, Swaan C, van der Bij A, et al. Ongoing mumps outbreak in a student population with high vaccination coverage, Netherlands, 2010. Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin. 2010;15:17.
Cui A, Zhu Z, Hu Y, Deng X, Sun Z, Zhang Y, Mao N, Xu S, Fang X, Gao H, et al. Mumps epidemiology and mumps virus genotypes circulating in mainland China during 2013-2015. PLoS One. 2017;12(1):e0169561.
Shuman EK. Global climate change and infectious diseases. The international journal of occupational and environmental medicine. 2011;2(1):11–9.
Falagas ME, Bliziotis IA, Kosmidis J, Daikos GK. Unusual climatic conditions and infectious diseases: observations made by Hippocrates. Enferm Infecc Microbiol Clin. 2010;28(10):716–8.
Yang Q, Yang Z, Ding H, Zhang X, Dong Z, Hu W, Liu X, Wang M, Hu G, Fu C. The relationship between meteorological factors and mumps incidence in Guangzhou, China, 2005-2012. Human vaccines & immunotherapeutics. 2014;10(8):2421–32.
Ho YC, Su BH, Su HJ, Chang HL, Lin CY, Chen H, Chen KT. The association between the incidence of mumps and meteorological parameters in Taiwan. Human vaccines & immunotherapeutics. 2015;11(6):1406–12.
Onozuka D, Hashizume M. Effect of weather variability on the incidence of mumps in children: a time-series analysis. Epidemiol Infect. 2011;139(11):1692–700.
Li R, Lin H, Liang Y, Zhang T, Luo C, Jiang Z, Xu Q, Xue F, Liu Y, Li X. The short-term association between meteorological factors and mumps in Jining, China. Sci Total Environ. 2016;568:1069–75.
Hubálek Z. North Atlantic weather oscillation and human infectious diseases in the Czech Republic, 1951–2003. Eur J Epidemiol. 2005;20(3):263–70.
Zhang DS, Zhang X, Ouyang YH, Zhang L, Ma SL, He J. Incidence of allergic rhinitis and meteorological variables: non-linear correlation and non-linear regression analysis based on Yunqi theory of chinese medicine. Chinese journal of integrative medicine. 2016.
D.Cryer J, Chan K-S: Time Series Analysis with Applications in R, Second edn: Springer; 2008.
Duan Y, Huang XL, Wang YJ, Zhang JQ, Zhang Q, Dang YW, Wang J. Impact of meteorological changes on the incidence of scarlet fever in Hefei City, China. Int J Biometeorol. 2016;60(10):1543–50.
Zhang Y, Bi P, Hiller JE: Weather and the transmission of bacillary dysentery in Jinan, northern China: a time-series analysis. Public Health Rep (Washington, DC : 1974) 2008, 123(1):61–66.
Lin H, Yang L, Liu Q, Wang T, Hossain SR, Ho SC, Tian L. Time series analysis of Japanese encephalitis and weather in Linyi City, China. International journal of public health. 2012;57(2):289–96.
Cui A, Zhu Z, Chen M, Zheng H, Liu L, Wang Y, Ma Y, Wang C, Fang X, Li P, et al. Epidemiologic and genetic characteristics of mumps viruses isolated in China from 1995 to 2010. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases. 2014;21:384–90.
Fisman D. Seasonality of viral infections: mechanisms and unknowns. Clinical microbiology and infection: the official publication of the European Society of Clinical Microbiology and Infectious Diseases. 2012;18(10):946–54.
Fisman DN. Seasonality of infectious diseases. Annu Rev Public Health. 2007;28:127–43.
Suminski RR, Poston WC, Market P, Hyder M, Sara PA. Meteorological conditions are associated with physical activities performed in open-air settings. Int J Biometeorol. 2008;52(3):189–97.
Polgreen PM, Bohnett LC, Yang M, Pentella MA, Cavanaugh JE. A spatial analysis of the spread of mumps: the importance of college students and their spring-break-associated travel. Epidemiol Infect. 2010;138(3):434–41.
Tang JW. The effect of environmental parameters on the survival of airborne infectious agents. J R Soc Interface. 2009;6(Suppl 6):S737–46.
We would like to thank the Chinese Center for Disease Control and Prevention, and the Beijing Meteorological Bureau for providing data used in the study.
The study was supported by the National Natural Science Foundation of China (funding numbers 81574098 and 81503678). The funding bodies played no role in the design of the study, data collection, data analysis, interpretation of data, and writing of the manuscript.
Availability of data and materials
The data used in this study are available from the corresponding author on reasonable request and with permission of the Chinese Center for Disease Control and Prevention and the Beijing Meteorological Bureau. All relevant data are provided as Additional files 2.
Ethics approval and consent to participate
This study was conducted in accordance with the tenets of the Declaration of Helsinki. The research did not involve any direct participation by human subjects. The mumps data were extracted from monthly reports maintained by the Chinese Center for Disease Control. All individual-level data were anonymous. The study was approved by Ethics Committee on Medicine and Laboratory Animals, Beijing University of Chinese Medicine, China (approval number 2017BZHYLL0507).
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.
About this article
Cite this article
Hao, Y., Wang, Rr., Han, L. et al. Time series analysis of mumps and meteorological factors in Beijing, China. BMC Infect Dis 19, 435 (2019). https://doi.org/10.1186/s12879-019-4011-6