Article Text

Download PDFPDF

Association between meteorological factors and bacillary dysentery incidence in Chaoyang city, China: an ecological study
  1. Yang Zhao,
  2. Yaxin Zhu,
  3. Zhiwei Zhu,
  4. Bo Qu
  1. Faculty of Health Statistics, School of Public Health, China Medical University, Shenyang, China
  1. Correspondence to Dr Bo Qu; qubo6666{at}163.com

Abstract

Objectives To quantify the relationship between meteorological factors and bacillary dysentery incidence.

Design Ecological study.

Setting We collected bacillary dysentery incidences and meteorological data of Chaoyang city from the year 1981 to 2010. The climate in this city was a typical northern temperate continental monsoon. All meteorological factors in this study were divided into 4 latent factors: temperature, humidity, sunshine and airflow. Structural equation modelling was used to analyse the relationship between meteorological factors and the incidence of bacillary dysentery.

Material Incidences of bacillary dysentery were obtained from the Center for Disease Control and Prevention of Chaoyang city, and meteorological data were collected from the Bureau of Meteorology in Chaoyang city.

Primary outcome measures The indexes including χ2, root mean square error of approximation (RMSEA), comparative fit index (CFI), standardised root mean square residual (SRMR) and goodness-of-fit index (GFI) were used to evaluate the goodness-of-fit of the theoretical model to the data. The factor loads were used to explore quantitative relationship between bacillary dysentery incidences and meteorological factors.

Results The goodness-of-fit results of the model showing that RMSEA=0.08, GFI=0.84, CFI=0.88, SRMR=0.06 and the χ2 value is 231.95 (p=0.0) with 15 degrees of freedom. Temperature and humidity factors had positive correlations with incidence of bacillary dysentery, with the factor load of 0.59 and 0.78, respectively. Sunshine had a negative correlation with bacillary dysentery incidence, with a factor load of −0.15.

Conclusions Humidity and temperature should be given greater consideration in bacillary dysentery prevention measures for northern temperate continental monsoon climates, such as that of Chaoyang.

  • Bacillary dysentery
  • Meteorological factors
  • Structural equation modeling

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • This study included multiple meteorological factors which can produce a better forecast of high-risk environments for dysentery epidemics than using a single factor alone.

  • The effects of meteorological variables on the incidence of bacillary dysentery were not comprehensively analysed as an all-inclusive model in previous studies, and we used structural equation modelling to clearly quantify the integrated effect of various meteorological factors on bacillary dysentery incidences.

  • We only used data from 1981 to 2010 in Chaoyang city, and further study should expand the time span and study area to reduce errors and improve accuracy in order to make better the applicability of these results feasible.

  • Using cross-sectional data in ecological study limited us to confirm the causal relationship and control the confounding factors.

Introduction

Bacillary dysentery, a type of dysentery, is a severe form of shigellosis, with the characteristics of blood in stool. Patients often develop bloody diarrhoea, fever and stomach cramps with an incubation of 1 or 2 days.1 The infection is spread from person to person via oral-feces, food or drinking water.2 Bacillary dysentery is one of the most common causes of diarrhoea, and the estimated annual incidence of Shigella episodes worldwide is about 164.7 million, of which 163.2 million cases have occurred in developing nations, resulting in 1.1 million deaths.3 Although episodes of bacillary dysentery in China have decreased from 1.3 million in 1991 to 0.24 million in 2011, there is still a considerable risk of bacillary dysentery, especially among the youngest and the oldest population groups and those in low economic regions.4 Although there had been a dramatic drop in overall infectious disease outbreaks between 1975 and 1995, it has since then reverted and maintained a gradual upward trend.5 It is reported that bacillary dysentery, one of the top four infectious diseases from 2005 to 2010, presents nearly 500 000 cases each year in China.6 In addition to high incidence rates, bacillary dysentery treatment is severely threatened by the episodic emergence of new antibiotic resistance strands.7 Thus, prevention heavily relies on public health measures, such as water and food sanitation, since there is a lack of safe and efficacious vaccines.

In addition to the heavy burden of bacillary dysentery on public health, a greater worry is the change in the incidence and distribution of infectious diseases that is concomitant with meteorological effects.8 According to the Intergovernmental Panel on Climate Change, there is estimated to be a 1.1–6.4°C increase in average surface temperature by 2100, an increase 2–9 times higher than the average global warming change during the last century.9 In addition to this temperature increase, it is projected that there will be increase in the frequency and severity of extreme weather conditions such as floods and droughts. Studies show that increase in disease outbreaks and variations can be attributed to this climate change.10–14 For example, in 2000, the increased global burden of diseases, such as diarrhoea and respiratory diseases, due to climate change attributed to over 150 000 deaths worldwide.15

Bacillary dysentery is a foodborne and waterborne infectious disease. The incidences of bacillary dysentery in Chaoyang were not stable in the past few years. From 1981 to 2000, the general trend of incidences was decreased, and the incidence showed a slight increase around 2000, then it reduced and fluctuated, to around 1.5 per million.16 The incidence of bacillary dysentery is related to many factors, including economical status, hygiene, food and water.17 Recently, studies have found associations between climate and bacillary dysentery in light of global warming.18–20 Meteorological factors such as temperature and humidity can affect bacteria reproduction and survival.18 Although a number of studies have addressed these issues, they have focused on a bivariate relationship between observed meteorological variables and incidence of bacillary dysentery rather than the integrated effect of latent meteorological variables. More importantly, few studies have quantified the impact of these meteorological variables. It is important to perform comprehensive quantitative analyses on the effects of meteorological variables on infectious diseases such as bacillary dysentery in order to illustrate trends in the correlations among the different meteorological factors such as temperature, humidity and sunshine. Chaoyang is located in the northwest of Liaoning province and has a northern temperate continental monsoon climate. Affected by the warm moisture from the southeast sea, and the dry cold air from the northern plateau, it has a semiarid subhumid climate, which easily becomes dry. The increasing temperature and decreasing precipitation makes the warm desiccation tendency quite obvious in Chaoyang city.21 Clear assessment of influencing factors on infectious diseases such as bacillary dysentery are necessary to shed light on adaptation strategies, policies and measures to lessen the projected adverse impacts.22

Regression and correlation analysis have been used to study the transmission of bacillary dysentery.23–25 Previous research found that the incidence of bacillary dysentery showed significant correlation with precipitation, relative humidity, pressure and temperature.23 ,24 Researchers also believed that meteorological indicators could be used to predict incidences of bacillary dysentery.19 ,25 However, effects of those meteorological endogenous variables on the incidence of bacillary dysentery were not comprehensively analysed as an all-inclusive model, and the impact sizes were not measured in these studies. The limitations of these methods are that some meteorological factors like temperature, sunshine and humidity cannot be clearly quantified. Structural equation modelling (SEM) is a statistical technique for testing and estimating causal relationships using a combination of statistical data and qualitative causal assumptions. SEM has gained much popularity due to its flexibility in experimental design and rigorous approach to model testing. SEM has the advantages of handling multiple dependent variables at one time, analysing indices that cannot be directly measured, etc.26 SEM has recently been applied to studies of infectious diseases.27 ,28

In this study we used SEM to analyse meteorological factors like temperature, sunshine, airflow and humidity and their impacts on bacillary dysentery incidence and to understand the integrated effect of meteorological factors on the incidence of bacillary dysentery so as to provide valuable suggestions to public health departments when they need to take decisive steps in health awareness and bacillary dysentery prevention.

Materials and methods

Study area

Chaoyang is located in the northwest part of Liaoning province with a population of 3 342 624 in 2000 and 3 399 665 in 2010 (figure 1). The climate is a typical northern temperate continental monsoon climate with four distinctive seasons: spring, March to May; summer, June to August; autumn, September to November; and winter, December to February. Its annual average temperature is about 9.5°C, sunshine time is 2850–2950 hours per year, and the average annual precipitation is limited to 450–580 mm because of its landlocked location. It extends between longitude 118′50″–121′17″ and latitude 40′25″–42′22″, stretching 165 km horizontally and 216 km vertically. The dry cold air from the Mongolian Plateau often descends on the area, forming a semidry and semihumid climate with little precipitation. Drought often occurs in the spring and autumn.

Figure 1

Location of Chaoyang city. Chaoyang is located in the northwest part of Liaoning province, and the climate is a typical northern temperate continental monsoon climate.

Data collection

The time frame chosen for this study was from 1981 to 2010. We obtained data on monthly incidences of bacillary dysentery through the Center for Disease Control and Prevention of Chaoyang city. Monthly meteorological data were collected from the Bureau of Meteorology in Chaoyang city. First, a logarithmic transformation of bacillary dysentery incidence was performed to normalise the distribution. Second, Pearson's correlation was conducted to explore the relevancy between meteorological variables and the monthly incidence of bacillary dysentery with a lag of 0–3 months. SEM was used to assess the independent effects of meteorological variables on the incidence of bacillary dysentery.

We built the hypothesis model which included four latent meteorological factors: temperature, humidity, sunshine and airflow. The temperature factor included average air temperature, maximum air temperature, minimum air temperature, average ground temperature, maximum ground temperature and minimum ground temperature. The humidity factor included monthly average evaporation amount, absolute humidity, relative humidity, maximum frozen soil depth, air pressure, number of non-precipitation days, maximum precipitation, precipitation, maximum depth of snow cover and maximum number of snow cover days. The sunshine factor included monthly average intensity of sunshine, average percentage of sunshine and average duration of sunshine, and the airflow factor included monthly average wind speed and wind direction.

Given that independency of observed variables should be a vital precondition for the application of SEM, the collinearity diagnostic was carried out. The results showed that some observed variables were not independent of each other. In order to adjust for this, the principal component factor analysis was conducted. Two factors from temperature observed variables were extracted, which were defined as ground temperature and air temperature, respectively. And the exploratory factor analysis identified four latent factors accounting for 72.89% of the total observed variables.

Construction of the SEM

SEM includes both the measurement model and the structural model.29–31 The whole model is based on relationships among observed variables and latent factors. In this study, the latent factors such as temperature and sunshine are listed in circles, while the observed variables are in rectangles. All of the observed variables were fit to the linear analysis. There were six main steps in the establishment of SEM: developing a theoretical model, constructing a path diagram of causal relationships, converting the path diagram into a set of structural equations and specifying a measurement model, fitting the proposed model, evaluating goodness-of-fit criteria, and interpreting and modifying the model.

Statistical analysis

SPSS V.13.0 software was used to predetermine meteorological factors in the study. Factor analysis was carried out to test for structural validity. LISREL V.8.5 was used to construct SEM. LISREL provided different goodness-of-fit measures. Goodness-of-fit in the present study was usually measured on the basis of χ2 test, goodness-of-fit index (GFI), comparative fit index (CFI), standardised root mean square residual (SRMR) and root mean square error of approximation (RMSEA). The χ2 test had traditionally been used to test the hypothesis of the relationships suggested in the hypothetical model. The rest of the values range between 0.0 and 1.0, with GFI, CFI values closer to 1.0 and RMSEA, SRMR values closer to 0.0 indicating good fit.

Results

The data of incidences and meteorological factors were collected from 1981 to 2010, which means that for the incidences and each factor there are 360 samples. Annual and monthly incidences of bacillary dysentery from 1981 to 2010 are shown in figure 2. Pearson's correlation analysis showed that temperature, humidity and sunshine factors were related to the monthly bacillary dysentery incidence at the highest level with a lag of 1 month. This result is consistent with other studies on meteorological variables and intestinal infectious diseases that have found an average lagged effect of 1 month.20 The time lag would capture the virus incubation period within the human body and from the onset of the bacillary dysentery to notification to the healthcare system.27

Figure 2

Annual and monthly incidence of bacillary dysentery (Chaoyang city, China, 1981–2010).

The initial diagram of causal relationships is presented in figure 3. The apostrophe represents the meteorological factors that are not listed here. The circles represent the latent variables, while the rectangles show the observed variables. A certain variable affects others with the direction of the arrow.

Figure 3

A hypothetical structural equation modelling of meteorological factors and bacillary dysentery. The apostrophe represents the meteorological factors that are not listed here. The circles represent the latent variables, while the rectangles show the observed variables. A certain variable affects others with the direction of the arrow.

The SEM fitting results showed that RMSEA=0.08, GFI=0.84, CFI=0.88 and SRMR=0.06. The χ2 value was 231.95 (p=0.0) with 15 degrees of freedom. The model produced an acceptable fit to the sample data. The results of SEM confirmed that all relevant meteorological indicators in this study were divided into three latent factors, temperature, humidity and sunshine. Factors like temperature and humidity had positive correlations with incidence of bacillary dysentery, with the factor load of 0.59 and 0.78, respectively. Sunshine had a negative correlation with bacillary dysentery incidence, with a factor load of −0.15 (figure 4). The effects of meteorological factors shown in figure 4 were confirmed by our study on bacillary dysentery. All the results here were completely standardised solutions.

Figure 4

Structural equation modelling of meteorological factors and bacillary dysentery. Temperature and humidity factors had positive correlations with incidence of bacillary dysentery, with the factor loads of 0.59 and 0.78, respectively. The sunshine factor had a negative correlation with the incidence of bacillary dysentery, with a factor load of −0.15.

For the temperature factor, air temperature was the most influential factor, with item loading between temperature of 0.79. This was followed by ground temperature, with respective item load of 0.71. This suggested that air temperature may be the most effective to predict the impact of temperature on bacillary dysentery. For the humidity factor, monthly average evaporation and precipitation influenced the incidence of bacillary dysentery, with respective item loads of 0.65 and 0.61, which suggested an emphasis on evaporation. Sunshine factor was mainly composed of only one item: intensity of sunshine, the item load of which was 0.93.

Figure 4 shows that meteorological indicators correlated with the incidence of bacillary dysentery. The effect of meteorological variables on bacillary dysentery incidence was calculated by multiplying the load between the observed variable and latent factor and the load between the latent factor and disease incidence. For example, monthly evaporation has a positive effect on the incidence of bacillary dysentery, with load of 0.51, which was calculated by multiplying 0.65 and 0.78. The results of calculation showed that the observed variables with greater effects on bacillary dysentery were monthly average evaporation, precipitation and air temperature, with item effects of 0.51, 0.48 and 0.47, respectively. The factor loads between meteorological variables and bacillary dysentery incidence are shown in table 1.

Table 1

The factor loads between meteorological variables and bacillary dysentery incidence in Chaoyang city from 1981 to 2010

Discussion

Global warming has far-reaching effects on public health and infectious diseases outbreaks.32–37 Recent studies of global warming have forecasted an increase in heat wave patterns, droughts and inclement weather cycles in many regions of the world.33 ,34 Incremental heat may give rise to increased likelihood of food pathogens and waterborne diseases.35 ,36 The results of Gilbreath37 showed that the risk of diarrhoea caused by global warming was ∼10% higher in some regions in Africa and Southeast Asia compared with when no climate change occurred. Therefore, public health is particularly vulnerable to climate change, and health policymakers need a variety of approaches in controlling infectious diseases especially in drought-affected regions.

In the present study, we used SEM to evaluate the influence of meteorological factors on the incidence of bacillary dysentery in Chaoyang city, a drought-prone region, China. There have been many studies on the impact of meteorological variables on bacillary dysentery. However, these studies have overlooked the interaction between different meteorological factors and incidences of bacillary dysentery. The results of SEM in this study are based on an all-inclusive model and take into consideration the comprehensive effects of meteorological factors and meteorological indicators on the incidence of bacillary dysentery. The goodness-of-fit statistics mean that SEM was suitable for understanding such relational data in multivariate systems. The final SEM indicated that three meteorological factors: temperature, humidity and sunshine affected the incidence of bacillary dysentery, but factors like humidity and temperature had greater effect and should be paid closer attention in the study area. The observed meteorological variables like monthly average evaporation, precipitation and air temperature had greater impact on the incidence of bacillary dysentery.

In this study, we found a positive correlation between temperature and incidence of bacillary dysentery. This finding is consistent with other studies on bacillary dysentery, as well as studies on other foodborne diseases such as salmonellosis and typhoid.18–20 ,38 ,39 Studies in northern cities of China have shown that a 1°C increase in maximum or minimum temperature is associated with an up to 12% increase in cases of bacillary dysentery.19 ,20 In this study, we examined two specific temperature parameters as well as their quantitative effects on bacillary dysentery. Of these two parameters, air temperature was the most influential factor with item effect of 0.47. Temperature could affect the reproduction of bacteria especially water reservoirs. It could cause food shortage for humans as well as an increased risk of food contamination.38 In addition, increase in heat and temperature can strongly influence the immune system and individual behaviour.40 ,41

Our results also showed a positive relationship between humidity and bacillary dysentery. This was consistent with some other studies.18–20 ,39 ,42 Monthly evaporation and precipitation have major effect on bacillary dysentery, with factor loads of 0.51 and 0.48. A study in two coastal cities in China showed that the monthly incidence of dysentery was positively correlated with rainfall and humidity.18 A study in Taiwan showed that daily precipitation levels were significantly correlated with all eight notifiable infectious diseases, including waterborne and vectorborne infectious diseases.39 Thus, the local meteorological conditions should be seriously considered for public health interventions on reducing future risks of bacillary dysentery. And the 1-month lagged effect also should be taken into account.

Studies in Europe, Australia, the USA and Peru have examined the association between diarrhoeal diseases and climate variation.43 ,44 The results showed that weather conditions, such as temperature, rainfall and humidity, could result in an increased incidence of foodborne and waterborne diseases. However, the effects of these meteorological variables are not consistent in some published studies, which could be due to local meteorological conditions, population and ecological characteristics in different regions. Drought, excess rainfall and flooding can also contribute to epidemics of waterborne infectious diseases due to poor sanitation resulting from runoff from overwhelmed sewage lines or the contamination of water by livestock.45 Owing to the arid climate in Chaoyang, drought is relatively common. Large amount of evaporation may affect access to safe drinking water, increase contaminant concentration in the water, or increase ways of usage of a particular body of water.15 ,44 The results of our study suggest that due attention should be paid to the humidity factor and its effects on bacillary dysentery. Further studies should focus on the effect of humidity indicators on bacillary dysentery transmission in drought-prone areas.

In addition to temperature and humidity, we also found a negative correlation between sunshine and bacillary dysentery. The intensity of sunshine had a negative effect on the incidence of bacillary dysentery, with an item effect of 0.14, and plays an exclusive role in the impact of sunshine on bacillary dysentery. This implies that intensity of sunshine should be taken into account in the prediction of bacillary dysentery. Though there are few studies focusing on sunshine and incidences of infectious diseases, our results indicated that sunshine could have possible effects on bacillary dysentery in the study area. Further understanding of sunshine on bacillary dysentery outbreaks may be needed.

This study provides data ranges to inform public health departments prior to bacillary dysentery outbreaks so that effective measures could be taken to mitigate a bacillary dysentery epidemic. Since this study includes multiple meteorological indicators, it can produce a better forecast of high-risk environments for dysentery epidemics than using a single factor alone. In addition, this understanding of meteorological factors' association with bacillary dysentery can contribute to the current health policy assessment and the development of effective bacillary dysentery control strategies.

Our study also had some limitations. First, the study was conducted in only one Chinese city, which may partly limit the representation of this study sample. Further studies in more areas of China will determine whether these findings can be generalised. Second, using cross-sectional data in ecological study limited us to confirm the causal relationship and control the confounding factors. Additionally, GFIs can only reveal the reality that the hypothesised SEM is compatible with collected data, and the possibility of more plausible SEMs cannot be ruled out.

Conclusion

Temperature, sunshine and humidity are factors that have correlations with the incidence of bacillary dysentery. More specifically, monthly average evaporation, precipitation and air temperature have greater effect on the incidence of bacillary dysentery. They should be given greater consideration in bacillary dysentery prevention measures for northern temperate continental monsoon climates, such as that of Chaoyang.

Acknowledgments

The authors would like to thank the Center for Disease Control and Prevention of Chaoyang city for access to the survey data. We gratefully acknowledge Peng Guan and Desheng Huang for the theoretical discussions as well as practical advice on the analysis of the data.

References

Footnotes

  • Contributors BQ designed and carried out the studies and performed the statistical analysis. YZ participated in the data analysis and drafted the manuscript. YxZ participated in the data collection. ZZ participated in coordination and helped polish the article.

  • Funding Supported by the National Natural Science Foundation of China, No. 81273186 and No. 30700690.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.