- Research article
- Open Access
Trends of reported human brucellosis cases in mainland China from 2007 to 2017: an exponential smoothing time series analysis
Environmental Health and Preventive Medicine volume 23, Article number: 23 (2018)
The main objective of this study was to describe the temporal distribution of monthly reported human brucellosis cases in mainland China and develop an appropriate time series model for short-term extrapolation forecast.
Surveillance data of the monthly reported human brucellosis cases occurring from April 1, 2007, to March 31, 2017, in mainland China were obtained. The spectrum analysis was first adopted to find the cyclic and seasonal features, the existence of the seasonality and trend was determined by exponential smoothing method and the seasonal-trend decomposition. The candidate models of exponential smoothing included the additive model and multiplicative model; R2 was selected as the indicator for the selection of candidate model, and the stability of the model was verified by adjusting the training data and test data set. Finally, the extrapolations of monthly incident human brucellosis cases in 2017 were made.
From April 1, 2007, to March 31, 2017, a total of 435,108 cases of Brucellosis occurred in mainland China were reported, with an average of 3626 cases per month and a standard deviation of 1834 cases. The R2 of the exponential smoothing method that based on additive model increased steadily from 0.927 to 0.949 with the increase of the data volume. Ten of 12 actual values fell in the confidence interval of predicted value.
Human brucellosis cases peaked during the months from March to August in mainland China, with clear seasonality. The exponential smoothing based on the additive model method could be effectively used in the time series analysis of human brucellosis in China. Control methods, such as vaccination, quarantine, elimination of infected animals, and good hygiene within the production cycle, should be strengthened with paying more attention to the seasonality. Further research is warranted to explore the drivers behind the seasonality.
The transmission of infectious diseases between domestic and wild animals is posing great challenges to public health [1,2,3]. Among the numerous emerging infectious diseases, including zoonoses, human brucellosis is one of the most common zoonotic diseases. Although there has been a noticeable decrease in the incidence and prevalence of brucellosis in many countries due to their domestic animal brucellosis control or eradication programs, human brucellosis is a growing public health concern in China with the increasing number of reported cases and more widespread natural foci . In mainland China, human brucellosis was first reported in 1905, the incidence of human brucellosis in China was quite severe before the 1980s and declined later [4, 5]. Comprehensive prevention and control measures have been gradually implemented in mainland China since 1950 . Nationwide reporting of human brucellosis was established during 1950–1963, during 1964–1976, vaccination for animals and humans was applied in the high-risk areas, and then during 1977–1988, a national brucellosis control program was administered with the introduction of clinical features, case definition, laboratory examination criteria, treatment options, and control measures; the vaccination of livestock was selected as the main control method. Since 1990, sentinel surveillance of the seroprevalence of brucellosis in humans and animals has been added to the national brucellosis control program . However, the incidence of human brucellosis has been increasing greatly since 1990, against the background that the incidence of animal brucellosis remained unchanged . Also, the geographic distribution of affected regions gradually expanded, from traditional pasturing areas to agricultural areas, from the north part to the south part . These kinds of variation could be attributed to various kinds of risk factors. The number of livestock increased significantly in the past decade to meet the growing demands for meat consumption in mainland China , which would give rise to the increase of the total population of infected domestic animals, even with relatively low-level seroprevalence in livestock. Lack of quarantine or pasteurization in the livestock product supply also increased the risk of infections in non-occupational populations and urban settings [9, 10].
Multiple publications have been dealing with the epidemiological characteristics and spatio-temporal distribution of human brucellosis [11, 12], while the studies focusing on the forecasting models were limited. The epidemic of human brucellosis is associated with the variation features of periodicity and seasonality; the major part of cases can be seen during the peak period, while sporadic cases can be observed during the non-peak period. Early recognition of the epidemics may provide strong support for the prevention and control of human brucellosis. The short-term forecast can evaluate the prevention or control measures; meanwhile, it can adopt timely and effective countermeasures for the epidemic peak that may occur. Among the available prediction models, time series analysis  was the most commonly adopted method, and some characteristics could be found from the annual data.
In the present study, the reported number of human brucellosis cases on a monthly basis in mainland China from 2007 to 2017 was collected and the short-term forecast was made with the seasonal-trend decomposition (STL) and exponential smoothing model.
Design and ethical considerations
The current study utilized a longitudinal, non-experimental design. This study was conducted on the basis of the nationwide surveillance of human brucellosis in mainland China and was exempt from institutional review board assessment. All the data in the present study were publicly open and were supplied and analyzed in an anonymous format, without access to any personal identifying information.
The monthly reported numbers of human brucellosis cases (clinically diagnosed and laboratory confirmed) occurring from April 1, 2007, to March 31, 2017, in mainland China were obtained from the National Health Commission of the People’s Republic of China and the CDC of the People’s Republic of China. In China, the Law on the Prevention and Control of Infectious Diseases  requires health-care staff to report any of the 39 infectious diseases to the Center for Disease Control and Prevention (CDC) through the National Noticeable Infectious Disease Reporting system (NIDR). Human brucellosis is listed as a Class II legally mandated notifiable disease and brucellosis is listed as a Class II key disease by the Disease Prevention and Control of Livestock and Poultry in China .
The clinical diagnosis of human brucellosis cases was made according to the clinical manifestations (undulant fever, sweats, fatigue, myalgia, arthralgia, arthritis, hepatomegaly, splenomegaly, epididymo-orchitis, etc.) and confirmed by serologic test or isolation of the organism, in accordance with the case definition of the World Health Organization .
The monthly reported number of human brucellosis cases was exported to Microsoft Excel 2016 for the distribution description. The spectrum analysis was first adopted to find the cyclic and seasonal features; the existence of the seasonality and trend was determined by exponential smoothing method and the seasonal-trend decomposition. The decomposition model assumes that the data has the following form: Times series = Pattern + Error = Trend cycle + Seasonality + Error. The candidate models of exponential smoothing included additive model and multiplicative model; R2 was selected as the indicator for the selection of candidate model, and the stability of the model was verified by adjusting the training data and test data set. Finally, the extrapolations of monthly incident human brucellosis cases in 2017 were made.
All statistical analyses were conducted by using the software IBM SPSS Statistics 21.0 for Windows (IBM, Asian Analytics Shanghai).
Reported cases of human brucellosis cases in China
From April 1, 2007, to March 31, 2017, a total of 435,108 cases of brucellosis occurred in mainland China were reported, with an average of 3626 cases per month and a standard deviation of 1834 cases. A heat map of the monthly cases was created to visualize the long term of reported cases during the 10 years (Fig. 1). The annual increase of the reported cases is 7.8%. The peak of monthly reported cases was 8102 cases, which appeared in July 2014; the bottom was 728 cases, which appeared in January 2009.
Seasonal features of reported brucellosis cases
Brucellosis cases peaked during the months from March to August, as observed in nearly all of the years of the study, while the most noticeable annual increase was observed from 2014 to 2015. There was a clear seasonality of monthly human brucellosis cases during the study period. During the peak season, a total of 302,128 cases were found, accounting for 69.44% of the total cases, with an average of 5035 cases and a standard deviation of 1440 cases. Overall, June was the month with the highest number of reported cases; there were 435,108 cases reported across the study period, where the peak of reported cases (60,782 cases) was in 2015. During the peak season, the peak of monthly reported cases was 8102 cases, which appeared in July 2014; the bottom was 2293 cases, which appeared in March 2008. By contrast, January had the lowest number of reported cases (17,965 cases) across the period, with the lowest monthly number (728 cases) in 2009.
According to spectrum analysis, the frequency ◊ period = 1; it is decomposed into many sine waves by the Fourier series. From Fig. 2, it can be seen that among these waves, the amplitude of the wave with the frequency of 0.0833 is the strongest, and its reciprocal value (1/0.0833) is the period of 12 months. Taking into account the peak time of reported cases, the seasonality of human brucellosis cases was thus determined.
Trend and periodicity of reported brucellosis cases
The seasonal decomposition of additive model indicated that there was a clear trend of increase in the monthly number of human brucellosis cases and there was an obvious upward trend from April 2007 to July 2015 (Fig. 3).
From Fig. 3, it can be seen that in the additive model, the first cycle is from April 2008 to April 2010, the second cycle is from April 2010 to July 2012, the third cycle is from July 2012 to April 2015, and the last cycle has been ongoing since April 2015. The error between observed and fitted values in additive model ranged from − 628 to 811 (Fig. 4).
Additive model sensitivity and forecasting performance
A good model requires stable predictions, thus 6 to 9 years of data was adopted to predict the next year’s incidence, respectively, to test the validity of the additive model of exponential smoothing method. From Table 1, it can be seen that the R2 of the model increased steadily from 0.927 to 0.949 with the increase of the data volume. From Fig. 5 and Table 2, it can be seen that all the actual values except for two in April and May 2017 fall in the confidence interval of predicted value.
In the present study, data accumulated in mainland China was used to provide an overview of the temporal distribution of human brucellosis cases and the feasibility of exponential smoothing method in the forecasting of incident human brucellosis cases was explored.
The study adopted the number of reported human brucellosis cases to present the temporal distribution in mainland China and the numbers that provided could give some hints to the allocation of brucellosis-related healthcare resources. The clear seasonality of reported human brucellosis was found in the present study, indicating more contact between human and infected livestock from March to August which is the cultivation time in the northern agricultural areas and the lambing time in the pastoral areas in China, also the warmer temperature might be suitable for the zoonotic transmission of pathogens [17, 18]. Temperature, monthly hours of sunshine and rainfall have shown the highest probabilities of influencing brucellosis transmission in China according to the Granger causality tests .
Models dealing with time series data have been traditionally used in econometrics, and recently their use is increasing in medical fields under the background of increasing data quality of disease surveillance database [19, 20]. Generally, the time series of disease surveillance sourced data can be decomposed into four components: the trend, the cyclical, the seasonal (systematic, calendar related movements), and the irregular pattern. This study targeted the application of exponential smoothing method to predict human cases of brucellosis based on historical time series. For exponential smoothing method, also known as exponential weighted average method, the major advantage of this model is that it chooses the weighted average of each period as a decreasing exponential sequence, giving greater weight to the historical data closer to the forecast period. It might be suitable for the forecast of time-series data with trends and seasonality. During the extrapolation prediction, all the actual values except for two (April and May 2017) fell in the confidence interval of predicted value (Fig. 5). If in the real world an actual incident value exceeds the upper confidence limit of the forecast, it will alert us the high emerging risk of the outbreak or epidemic which is necessary to put more special effort into the coping strategy and solutions of brucellosis, not only in medical and veterinary fields but also in economic areas [21,22,23,24].
The limitations of this study should be acknowledged. First, as brucellosis is a zoonotic disease caused by Brucella spp. and the transmission routes to humans are direct with infected animals, consumption of meat and dairy products, or intentional or accidental aerosol exposure [25, 26]. Monthly data of the abovementioned factors could not be obtained in the present study; we could only collect these data on an annual basis. The risk of the disease is mainly affected by the amount of livestock herds [27,28,29]. We have found that the Pearson correlation coefficient between the annual number of human brucellosis cases and cattle slaughter amount, sheep slaughter amount, beef yield, mutton yield, and wool yield is 0.803, 0.781, 0.642, 0.708, and 0.564, respectively, indicating the existence of strong correlations. Second, the analysis and forecast were based on the data that were collected from a passive surveillance system of infectious diseases. During the study period, the modifications of the laboratory diagnostic test, the diagnosis criteria of cases, and the availability of health-care providers and hospital resources might have an impact on the data quality. Third, data about the pathogens and their variations were not available in the present study. Fourth, we are not able to utilize data at the individual province level. However, the monthly data that we used were the most nationally official data of human brucellosis in mainland China.
Human brucellosis cases peaked during the months from March to August in mainland China, with clear seasonality. The exponential smoothing that based on additive model method could be effectively used in the time series analysis of human brucellosis in China. Control methods, such as vaccination, quarantine, elimination of infected animals and good hygiene within the production cycle should be strengthened with paying more attention to the seasonality. Further research is warranted to explore the drivers behind the seasonality.
Bayesian Information Criterion
Center for Disease Control and Prevention
Mean absolute error
Mean absolute percentage error
Maximum absolute deviation
Maximum absolute percentage error
National noticeable infectious disease reporting system
Root mean square error
Peiris JS, Cowling BJ, Wu JT, et al. Interventions to reduce zoonotic and pandemic risks from avian influenza in Asia. Lancet Infect Dis. 2016;16(2):252–8.
Coker R, Rushton J, Mounier-Jack S, et al. Towards a conceptual framework to support one-health research for policy on emerging zoonoses. Lancet Infect Dis. 2011;11(4):326–31.
Franco MP, Mulder M, Gilman RH, et al. Human brucellosis. Lancet Infect Dis. 2007;7(12):775–86.
Lai S, Zhou H, Xiong W, et al. Changing epidemiology of human brucellosis, China, 1955-2014. Emerg Infect Dis. 2017;23(2):184–94.
Deqiu S, Donglou X, Jiming Y. Epidemiology and control of brucellosis in China. Vet Microbiol. 2002;90(1–4):165–82.
Gang S, Zhao Y, Zhang S. Analysis of 1990 to 2001 surveillance effect of national major surveillance place for brucellosis [in Chinese]. Chinese Journal of Control of Endemic. Diseases. 2002;17(5):285–8.
Zhong Z, Yu S, Wang X, et al. Human brucellosis in the People’s Republic of China during 2005-2010. Int J Infect Dis. 2013;17(5):e289–92.
National Development and Reform Commission of China. National beef and mutton production development plan (2013–2020) [cited 27 April 2018]. http://www.ndrc.gov.cn/zcfb/zcfbghwb/201402/P020140221362074574291.pdf.
Tan Z, Huang Y, Liu G, Zhou W, Xu X, Zhang Z, et al. A familial cluster of human brucellosis attributable to contact with imported infected goats in Shuyang, Jiangsu Province, China, 2013. Am J Trop Med Hyg. 2015;93:757–60.
Chen Q, Lai S, Yin W, et al. Epidemic characteristics, high-risk townships and space-time clusters of human brucellosis in Shanxi Province of China, 2005-2014. BMC Infect Dis. 2016;16:760.
Li YJ, Li XL, Liang S, et al. Epidemiological features and risk factors associated with the spatial and temporal distribution of human brucellosis in China. BMC Infect Dis. 2013;13:547.
Jia P, Joyner A. Human brucellosis occurrences in Inner Mongolia, China: a spatio-temporal distribution and ecological niche modeling approach. BMC Infect Dis. 2015;15:36.
Yang L, Bi ZW, Kou ZQ, et al. Time-series analysis on human brucellosis during 2004-2013 in Shandong Province, China. Zoonoses Public Health. 2015;62(3):228–35.
Yang GH, Stroup DF, Thacker SB. National public health surveillance in China: implications for public health in China and the United States. Biomed Environ Sci. 1997;10(1):1–13.
Zhang WY, Guo WD, Sun SH, et al. Human brucellosis, Inner Mongolia, China. Emerg Infect Dis. 2010;16(12):2001–3.
Ministry of Health of the People's Republic of China. Brucellosis diagnostic criteria and treatment principles. Beijing: China People's Health Publishing House; 2008. p. 71–3.
Zhou L, Fan M, Hou Q, et al. Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China. Math Biosci Eng. 2018;15(2):543–67.
Lou P, Wang L, Zhang X, et al. Modelling seasonal brucellosis epidemics in Bayingolin Mongol Autonomous Prefecture of Xinjiang, China, 2010-2014. Biomed Res Int 2016;2016:5103718.
Chen JD, Ke CW, Deng X, et al. Brucellosis in Guangdong Province, People’s Republic of China, 2005-2010. Emerg Infect Dis. 2013;19(5):817–8.
Nematollahi S, Ayubi E, Karami M, et al. Epidemiological characteristics of human brucellosis in Hamadan Province during 2009-2015: results from the National Notifiable Diseases Surveillance System. Int J Infect Dis. 2017;61:56–61.
Liu Q, Cao L, Zhu XQ. Major emerging and re-emerging zoonoses in China: a matter of global health and socioeconomic development for 1.3 billion. Int J Infect Dis. 2014;25:65–72.
Chen S, Zhang H, Liu X, et al. Increasing threat of brucellosis to low-risk persons in urban settings, China. Emerg Infect Dis. 2014;20(1):126–30.
Wang Y, Zhang W, Ke Y, et al. Human brucellosis, a heterogeneously distributed, delayed, and misdiagnosed disease in China. Clin Infect Dis. 2013;56(5):750–1.
Ducrotoy MJ, Bertu WJ, Ocholi RA, et al. Brucellosis as an emerging threat in developing economies: lessons from Nigeria. PLoS Negl Trop Dis. 2014;8(7):e3008.
Pappas G, Papadimitriou P, Akritidis N, et al. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.
Galinska EM, Zagorski J. Brucellosis in humans—etiology, diagnostics, clinical forms. Ann Agric Environ Med. 2013;20(2):233–8.
Blackburn JK, Kracalik IT, Fair JM. Applying science: opportunities to inform disease management policy with cooperative research within a one health framework. Front Public Health. 2015;3:276.
Godfroid J, DeBolle X, Roop RM, et al. The quest for a true one health perspective of brucellosis. Rev Sci Tech. 2014;33(2):521–38.
Craighead L, Meyer A, Chengat B, et al. Brucellosis in West and Central Africa: a review of the current situation in a changing landscape of dairy cattle systems. Acta Trop. 2018;179:96–108.
This study was supported by the grant from the National Natural Science Foundation of China (Grant No. 71573275) and Science Foundation of Liaoning Provincial Department of Education (LFWK201719 and LK201644).
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
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
Guan, P., Wu, W. & Huang, D. Trends of reported human brucellosis cases in mainland China from 2007 to 2017: an exponential smoothing time series analysis. Environ Health Prev Med 23, 23 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12199-018-0712-5
- Human brucellosis
- Time series analysis
- Exponential smoothing method