
Water pollution prevention and control of the Xiang River has become an issue of great concern to China's central and local governments. To further analyze the effects of central and local governmental policies on water pollution prevention and control for the Xiang River, this study performs a big data analysis of 16 water quality parameters from 42 sections of the mainstream and major tributaries of the Xiang River, Hunan Province, China from 2005 to 2016. This study uses an evidential reasoning-based integrated assessment of water quality and principal component analysis, identifying the spatiotemporal changes in the primary pollutants of the Xiang River and exploring the correlations between potentially relevant factors. The analysis showed that a series of environmental protection policies implemented by Hunan Province since 2008 have had a significant and targeted impact on annual water quality pollutants in the mainstream and tributaries. In addition, regional industrial structures and management policies also have had a significant impact on regional water quality. The results showed that, when examining the changes in water quality and the effects of pollution control policies, a big data analysis of water quality monitoring results can accurately reveal the detailed relationships between management policies and water quality changes in the Xiang River. Compared with policy impact evaluation methods primarily based on econometric models, such a big data analysis has its own advantages and disadvantages, effectively complementing the traditional methods of policy impact evaluations. Policy impact evaluations based on big data analysis can further improve the level of refined management by governments and provide a more specific and targeted reference for improving water pollution management policies for the Xiang River.
Citation: Yangyan Zeng, Yidong Zhou, Wenzhi Cao, Dongbin Hu, Yueping Luo, Haiting Pan. Big data analysis of water quality monitoring results from the Xiang River and an impact analysis of pollution management policies[J]. Mathematical Biosciences and Engineering, 2023, 20(5): 9443-9469. doi: 10.3934/mbe.2023415
[1] | Ya Li, Z. Feng . Dynamics of a plant-herbivore model with toxin-induced functional response. Mathematical Biosciences and Engineering, 2010, 7(1): 149-169. doi: 10.3934/mbe.2010.7.149 |
[2] | Xinru Zhou, Xinmiao Rong, Meng Fan, Josué-Antonio Nescolarde-Selvaa . Stoichiometric modeling of aboveground-belowground interaction of herbaceous plant. Mathematical Biosciences and Engineering, 2019, 16(1): 25-55. doi: 10.3934/mbe.2019002 |
[3] | Alhadi E. Alamir, Gomah E. Nenaah, Mohamed A. Hafiz . Mathematical probit and logistic mortality models of the Khapra beetle fumigated with plant essential oils. Mathematical Biosciences and Engineering, 2015, 12(4): 687-697. doi: 10.3934/mbe.2015.12.687 |
[4] | Jeong-Mi Yoon, Volodymyr Hrynkiv, Lisa Morano, Anh Tuan Nguyen, Sara Wilder, Forrest Mitchell . Mathematical modeling of Glassy-winged sharpshooter population. Mathematical Biosciences and Engineering, 2014, 11(3): 667-677. doi: 10.3934/mbe.2014.11.667 |
[5] | Yang Kuang, Jef Huisman, James J. Elser . Stoichiometric Plant-Herbivore Models and Their Interpretation. Mathematical Biosciences and Engineering, 2004, 1(2): 215-222. doi: 10.3934/mbe.2004.1.215 |
[6] | Guangyu Sui, Meng Fan, Irakli Loladze, Yang Kuang . The dynamics of a stoichiometric plant-herbivore model and its discrete analog. Mathematical Biosciences and Engineering, 2007, 4(1): 29-46. doi: 10.3934/mbe.2007.4.29 |
[7] | Zhilan Feng, Wenzhang Huang, Donald L. DeAngelis . Spatially heterogeneous invasion of toxic plant mediated by herbivory. Mathematical Biosciences and Engineering, 2013, 10(5&6): 1519-1538. doi: 10.3934/mbe.2013.10.1519 |
[8] | Yang Li, Jia Li . Stage-structured discrete-time models for interacting wild and sterile mosquitoes with beverton-holt survivability. Mathematical Biosciences and Engineering, 2019, 16(2): 572-602. doi: 10.3934/mbe.2019028 |
[9] | Sunmi Lee, Chang Yong Han, Minseok Kim, Yun Kang . Optimal control of a discrete-time plant–herbivore/pest model with bistability in fluctuating environments. Mathematical Biosciences and Engineering, 2022, 19(5): 5075-5103. doi: 10.3934/mbe.2022237 |
[10] | Luca Galbusera, Sara Pasquali, Gianni Gilioli . Stability and optimal control for some classes of tritrophic systems. Mathematical Biosciences and Engineering, 2014, 11(2): 257-283. doi: 10.3934/mbe.2014.11.257 |
Water pollution prevention and control of the Xiang River has become an issue of great concern to China's central and local governments. To further analyze the effects of central and local governmental policies on water pollution prevention and control for the Xiang River, this study performs a big data analysis of 16 water quality parameters from 42 sections of the mainstream and major tributaries of the Xiang River, Hunan Province, China from 2005 to 2016. This study uses an evidential reasoning-based integrated assessment of water quality and principal component analysis, identifying the spatiotemporal changes in the primary pollutants of the Xiang River and exploring the correlations between potentially relevant factors. The analysis showed that a series of environmental protection policies implemented by Hunan Province since 2008 have had a significant and targeted impact on annual water quality pollutants in the mainstream and tributaries. In addition, regional industrial structures and management policies also have had a significant impact on regional water quality. The results showed that, when examining the changes in water quality and the effects of pollution control policies, a big data analysis of water quality monitoring results can accurately reveal the detailed relationships between management policies and water quality changes in the Xiang River. Compared with policy impact evaluation methods primarily based on econometric models, such a big data analysis has its own advantages and disadvantages, effectively complementing the traditional methods of policy impact evaluations. Policy impact evaluations based on big data analysis can further improve the level of refined management by governments and provide a more specific and targeted reference for improving water pollution management policies for the Xiang River.
Hand-foot-and-mouth disease (HFMD) causes a substantial disease burden in Asian regions including China, mainly in children below 5 years old and partly in adolescent students [1,2]. Enterovirus 71 (EV71), coxsackievirus A16 (CA16) are the most common enterovirus serotypes causing HFMD in China [3,4]. Although most cases of HFMD are mildly self-limiting, characterized by fever, flat spots or bumps on the hands, feet and mouth, possible transmission routes of HFMD contain casual contact with contaminated discharges, fluid of blisters or stool from infected persons, or contaminated objects, severe cases may develop into heart and lung failure and even death [5,6,7].
In the last decade, many severe HFMD outbreaks have occurred in East Asia and Southeast Asia, which has become a serious public health problem in the affected countries and regions. Despite generation and approval of an EV71 vaccine in China in 2016, recently several HFMD outbreaks were caused by other pathogens, making it important to take additional actions to control HFMD epidemics [8]. The increasing burden of HFMD has become a serious public health problem in China, causing the government to pay greater attention to the risk factors of this disease, and focus on the effective prevention and control strategies.
A distinct seasonal pattern with annual peaks in the summer was reported in temperate regions [9], whereas in tropical and subtropical regions, more than two peaks may occur within a year [10,11,12,13,14]. Many studies have confirmed that meteorological factors have an important influence on the incidence of HFMD despite the inconsistent results from different studies [9,12,15,16,17,18,19]. Relative-humidity and temperature are regarded as the top two meteorological factors affecting the HFMD transmission [15,17,18,19]. Attending kindergarten, child-care center or school and public gathering were also pointed out as risk factors for HFMD [20,21,22]. However, it is not clear whether above factors are related to multiple-peak pattern, and when and how these factors play a critical role in determining the HFMD transmission in one year.
Mathematical epidemiology is the study of the spread of disease in space and time, with the aim to identify the factors that are responsible for or contributing to their occurrence. The basic reproduction number (R0) is often used to assess the transmission potential or transmissibility of an infectious disease, which indicates the average number of expected secondary infections resulting from a single infectious case [23]. In general, if it is greater than one, the disease will continue to spread, while if it is less than one, the disease will eventually disappear [24,25]. Both statistically robust analytical methods (model-free estimation) and dynamic modelling of infectious diseases (model-based estimation) have contributed greatly to estimate R0, which may be used to guide infectious disease prevention, control and eradication [10,26,27].
Now, mathematical models have become one of the important tools for analyzing the spread and control of HFMD [3,28,29,30,31,32,33]. The aim of this study is to investigate the transmission dynamics and examine when and how its potential factors play a decisive role in the annual multiple-peak pattern based on the surveillance data of HFMD in Wenzhou city, Zhejiang province, China. We first estimated daily and annual model-free basic reproduction number Rfree0(day) and Rfree0(year) based on the infected data and the infection characteristics of HFMD, in which Rfree0(day) means the average number of secondary cases that each infected individual would infect if the conditions remained as they were at the time of a day, and Rfree0(year) is the annual average of Rfree0(day). Then we developed a mathematical model incorporating the potential factors, which are screened and identified by the correlation with the obtained daily model-free basic reproduction number, including school opening or not, and meteorological factors. After then, reported HFMD infection data were further used to estimate the undetermined model parameters. Thus, the annual model-based basic reproduction number Rmodel0(year) can be calculated after some reasonable approximation, and sensitivity analysis of the model was performed to estimate the resultant effects of the potential factors. In the end, some feasible control strategies of HFMD were proposed based on the above analyses.
Infected data of HFMD in Wenzhou between 1 January 2010 to 31 December 2015 were collected by Wenzhou Center for Disease Control and Prevention [34], including basic demographic characteristics of HFMD cases, and daily incident cases. According to the obtained records, all the cases were children aged between 0 and 14 years old, and there were 203,488 cases in all.
School summer and winter vacation as well as school opening duration was determined according to the guidelines of the school calendar issued by the Education Bureau [35].
Daily meteorological factors, including minimum temperature, mean temperature, maximum temperature, minimum relative humidity and mean relative humidity, were collected from the China Meteorological Data Sharing Service System [36]. A weather monitoring station is located in Ruian County, Wenzhou City (N27o', E120o39'), from which the data are usually used to represent the meteorological conditions of the whole Wenzhou city.
To quantify the intensity of transmission over time, we first suppose that serial interval (the time between the onset of symptoms in a primary case and the onset of symptoms in secondary cases) is equal to the sum of L and D, where L represents the duration of incubation and D denotes the duration of infection [37]. Then serial interval distribution was obtained by calculating the joint probability density function of L and D. Thus, according to the method in [26], we can estimate the daily model-free basic reproduction number Rfree0(day) by using "epiEstim" Excel software package [26]. According to the annual average, we can obtain the annual model-free basic reproduction number Rfree0(year). See Supplementary A for more details.
According to the primary epidemiological properties of HFMD and the implemented public health interventions (quarantine, isolation and hygiene precaution), we divide the total population (N) into five compartments, named as susceptible (S), exposed but not yet infectious (E), infectious (I), quarantine (Q) and recovered (R). Figure 1 illustrates the flowchart of the five compartments, which describes how individuals can move among the states.
Based on the flowchart in Figure 1, the dynamic model can be described by the following non-autonomous differential equations:
{dSdt=λ−(d+c)S−β(t)SIN+ηR,dEdt=β(t)SIN−(d+σ)E,dIdt=σE−(d+δ1+k+γ1)I,dQdt=kI−(d+δ2+γ2)Q,dRdt=γ1I+γ2Q−(d+η+c)R,dNdt=λ−dN−cS−cR−δ1I−δ2Q. | (2.1) |
The meaning and values of parameters, and initial conditions in model (2.1) are given in Table 4 (See Supplementary B). Note that once individuals are infected with the virus, even though they are growing in age, they can only leave the system after evolving into recovered compartment. Thus, the parameter of transfer rate due to age growth c only appears in the first, fifth and last equations in (2.1).
Here β(t) is the contact rate between susceptible (S) and infectious (I), which is a time-varying, continuous, positive function. Concretely, β(t) can be divided into five sections: (1) the interaction factor between school opening and meteorological effect during spring semester, βs×m(t); (2) the interaction factor between school opening and meteorological effect during fall semester, βf×m(t); (3) the meteorological factor during spring semester, βms(t); (4) the meteorological factor during fall semester, βmf(t); (5) the meteorological factor during winter and summer vacation, βmr(t). Therefore, β(t) can be expressed as
β(t)=βs×m(t)+βf×m(t)+βms(t)+βmf(t)+βmr(t). | (2.2) |
For detailed expressions of (2.2) and more details, see Supplementary B.
Based on the reported HFMD cases and the dynamic model (2.1), an adaptive Metropolis-Hastings algorithm was conducted to carry out the Markov Chain Monte Carlo (MCMC) procedure to estimate the undetermined parameters in contact rate and their standard deviations (see Supplementary C for details). The algorithm runs for 1000000 iterations with a burn-in of 500000 iterations, with the Geweke convergence diagnostic method employed to assess the convergence [38]. The mean values of the estimated parameters and their standard deviations are listed in Table 5. Note that the fluctuations of meteorological factors in different years can be approximately regarded as periodic changes. i.e., the contact rate β(t) can be approximately regarded as annual periodic function. Thus, based on the results in [39,40], the annual model-based basic reproduction number Rmodel0(year) can be calculated by estimated parameters and simulation numerical integration (See Supplementary D for details).
Note that there are 19 undetermined parameters in contact rate, each with different mean value and variance, which means that it is tedious and extremely complex to study their influence on the outcomes. In order to detect the underlying factors more simply, a better idea is to do the sensitivity analysis for different categories of parameters. Thus, according to the expression (2.2) of contact rate, we classified all parameters into five factors, also marked as βs×m, βf×m, βms, βmf and βmr (see Supplementary B for details). Concretely, we use Latin hypercube sampling (LHS) and partial rank correlation coefficients (PRCCs) [41,42] to study the dependence of the above five factors on daily contact rate and annual model-based basic reproduction number and daily contact rate from 2010 to 2015.
Firstly, to identify key underlying factors that influenced the HFMD dynamics, we used LHS and PRCCs to examine the sensitivity of the annual model-based reproduction number to the five factors [42,43]. Normal distribution for all input parameters with the mean values and standard deviations was given in Table 5. By substituting the sampled parameters' values into the formulas (7.2-7.6) in Supplementary B, all five factors' values were obtained. Then the annual global PRCCs were calculated between the five factors and the corresponding output variable, i.e., the annual model-based reproduction number.
Secondly, because the correlation between daily model-free basic reproduction number and daily contact rate has statistical significance, time-varying sensitivity was done between daily contact rate and five factors. Concretely, we still choose all five factors as the input variables, but use the estimated daily contact rate as the output variable. Using similar sampling parameter method as above, time-varying PRCCs were calculated for each factor. Here the sample size Ns is 10000 in both global and time-varying sensitivity analyses. Similar to [41,42], we considered absolute values of PRCC >0.4 as indicating an important or strong correlation between input parameters and output variables, values between 0.2 and 0.4 as moderate or weak correlations, and values between 0 and 0.2 as not significantly different from zero.
For non-normal distribution data, such as all meteorological factors, median with the corresponding first and third quartiles (Q1-Q3) was reported. Otherwise, Mean ± standard deviation (SD) was presented, such as annual model-free and model-based basic reproduction number. Spearman correlation coefficient was used to measure the correlation between meteorological factors, and the correlation between daily model-free basic reproduction number and contact rate. According to whether the daily model-free basic regeneration number is less than one, we divide the meteorological factors into two groups (Rfree0(day)≥1 and Rfree0(day)<1), and Mann-Whitney test was used to analyze the difference between the two groups. Furthermore, independent-samples t test was performed to analyze the difference between annual model-free and model-based basic reproduction numbers. In order to compare with the threshold value (one), one-sample t test was performed for both annual model-free and model-based basic reproduction numbers. A two-tail P<0.05 was considered statistically significant.
Figure 2 showed the estimated daily model-free basic reproduction number (Rfree0(day)) and the reported daily HFMD infectious cases in Wenzhou between 1 January 2010 to 31 December 2015. In every year, we can see that the basic reproduction number of HFMD (blue solid line in Figure 2) began to increase above unit during early March and it usually remained above one for around four months from early March to late June. After that, a rebound of basic reproduction number occurred between September and October in autumn, then after wards fluctuated below one for most of the time during winter. This is consistent with the multiple-peak pattern of HFMD outbreaks, which can be seen from the reported data (orange and green histogram in Figure 2).
Figure 3 revealed that daily model-free basic reproduction number changed with school opening days during spring semester and fall semester, respectively. In the early stage of the spring semester, the basic reproduction number began to rise exponentially and quickly exceeded threshold one in no more than half a month (Figure 3A), and it usually peaked around the 30th day after school beginning and then remained above one until about the 120th day. Similarly, at the beginning of the fall semester, the basic reproduction number also grew rapidly and exceeded one in half a month (Figure 3B), and it usually peaked around the 30th day after the beginning of the fall semester and then fell down to less than one in a month. After two months of school beginning, the basic reproduction number began to increase again and reached a peak slightly above one, and then continued to decline through winter.
Various meteorological factors were found to be inter-related with each other. Table 1 summarized the correlation coefficients between all these factors. It was found that mean temperature, maximum temperature and minimum temperature were highly correlated with each other, with the correlation coefficients larger than 0.90. Moreover, mean relative humidity and minimum relative humidity were also highly interrelated with a correlation coefficient 0.910.
Mean TEM | Maximum TEM | Minimum TEM | Mean RHU | |
Maximum TEM | 0.977 | |||
Minimum TEM | 0.988 | 0.943 | ||
Mean RHU | 0.190 | 0.114 | 0.265 | |
Minimum RHU | 0.206 | 0.087 | 0.295 | 0.910 |
Note: TEM represents temperature; RHU represents relative humidity. |
Table 2 showed the comparison of each meteorological factors between the group of daily model-free reproduction number Rfree0(day)≥1 and Rfree0(day)<1. There were significant differences in all factors between the two groups (P<0.05). Hence, combining the previous correlation analysis results, we included only mean temperature and mean relative humidity in the dynamic model, and excluded other highly correlated variables.
Parameters | Rfree0(day)≥1 | Rfree0(day)<1 | P value |
Mean TEM | 19.3 (14.4, 23.2) | 20.0 (10.1, 27.7) | <0.001 |
Maximum TEM | 23.3 (18.3, 27.2) | 23.8 (14.1, 31.2) | <0.001 |
Minimum TEM | 16.8 (11.6, 20.8) | 17.1 (7.7, 25.0) | <0.001 |
Mean RHU | 78.5 (68.0, 86.0) | 76.0 (68.0, 84.0) | <0.001 |
Minimum RHU | 59.0 (46.0, 70.0) | 57.0 (45.0, 68.0) | 0.005 |
Note: TEM represents temperature; RHU represents relative humidity. |
We employed an adaptive Metropolis-Hastings algorithm to carry out extensive MCMC simulations [38], and to estimate mean values of parameters in the contact rate expression (2.2). Using large sample realization, we fitted the model (2.1) to the data of reported HFMD cases and the fitted parameters are shown in Table 5 in Supplementary C.
Figure 4 illustrates the best-fit numerical simulation of the model with the fitted parameters, which shows that the model was well fitted to the reported cases, and can provide the overall trend of HFMD infected population in Wenzhou city. Furthermore, the annual model-based basic reproduction number Rmodel0(year) from 2010 to 2015 were given in Table 3. Comparing with the threshold one, we find that all annual model-based and model-free basic reproduction numbers were significantly larger than threshold one, which indicates the consistency of two estimation methods. This suggests that our findings may be a possible mirror of the real transmission of HFMD i.e., HFMD is an endemic disease in Wenzhou city, which is consistent with the relevant studies in other regions of China [29,30,31,33]. However, we can also find that the model-based estimates are usually less than model-free estimates (Table 3), which should arouse our attention and further research in the future.
Year | Rfree0(year) | Rmodel0(year) |
2010 | 1.034±0.012♯ | 1.026±0.012∗♯ |
2011 | 1.153±0.012♯ | 1.125±0.012∗♯ |
2012 | 1.133±0.012♯ | 1.119±0.012∗♯ |
2013 | 1.123±0.011♯ | 1.112±0.012∗♯ |
2014 | 1.149±0.013♯ | 1.123±0.012∗♯ |
2015 | 1.068±0.012♯ | 1.041±0.011∗♯ |
∗: P<0.05 represents the comparison with Rfree0(year); ♯: P<0.05 represents the comparison with threshold one |
Taking annual model-based reproduction number as the target output variable, we calculated the PRCCs of all five factors (βs×m, βf×m, βms, βmf and βmr) in contact rate. From Figure 5, we found that all five factors were sensitive to the annual model-based basic reproduction number every year. As expected, the stronger the factor, the greater annual basic reproduction number. In particular, we find that the PRCCs of each factor change little every year, which indicates that each factor has better robustness to the impact of HFMD outbreak. By comparing PRCCs of each factor, we can further conclude that the influence of factors (βs×m, βf×m and βms) are greater than that of factors (βmf and βmr). These may provide potential guidance for HFMD prevention and control in strategy formulation.
Firstly, the correlation between the daily model-free basic reproduction number and the daily model-based contact rate was demonstrated. Figure 6A shows the time series diagram of estimated daily model-free basic reproduction number and model-based contact rate. Furthermore, the results of Spearman correlation analysis indicates that there is a significant correlation between them (ρ=0.446, P<0.05).
Due to the positive correlation between them, time-varying analysis is performed through the time-varying PRCCs between the above-mentioned five factors and daily contact rate β(t) with time (Figure 6B). Because the PRCCs of meteorological factor (βmr) is very large during summer vacation, we can naturally conclude that the cases during this period are mainly caused by the effect of meteorological factors. In addition, during the spring semester, the PRCCs of meteorological factor (βms) is increasing gradually, but the PRCCs of interaction factor between school opening and meteorological effect (βs×m) is decreasing gradually, and those PRCCs almost maintain a high level throughout the spring semester. This may be the reason for the observed persistent peak of HFMD outbreaks in spring semester. In autumn semester, the variation trend of PRCCs of βmf and βf×m is similar to βms and βs×m, but the time for the interaction factor to maintain a higher level is obviously shorter than that for the meteorological factors to maintain a higher level. This may be the reason for the observed multiple small peaks in autumn semester. As a result, these may provide a quantitative characterization of the underlying factors in the annual multi-peak pattern of HFMD outbreak.
Similar to the pattern in other subtropical and tropical regions, HFMD epidemics in Wenzhou occurred at least twice a year. The underlying factors for the annual pattern of HFMD multiple outbreaks have been investigated widely [10,11,12,13,14]. Some studies suggested a strong association between HFMD and weather [11,12]. In Hong Kong, researchers found that HFMD transmission was mainly associated with depletion of susceptible and limited impact were suggested from meteorological factors and school holidays [44]. However, in our study, it is clear that high contact rate was synchronized with the rapid increase of the daily model-free basic reproduction number, especially at the beginning of the school semester (Figure 6A). This pattern may be driven by a new cohort of susceptible children and few infected students entering kindergartens as well as the increase in contact rates among children in the new school semester [44]. Nevertheless, further confirmatory investigation is needed to demonstrate the impact of new students in kindergartens on HFMD transmission.
Though many studies indicated that meteorological factors play an important role in HFMD transmission [9,12,15,16,17,18,19], this study further find that there are different impact effect in different seasons and their interaction with school opening. Meteorological factors in spring semester, the interaction factor between school opening and meteorological effect, made more significant contributions to HFMD transmission (Figure 5 and 6). Meanwhile, both daily model-free basic reproduction number and contact rate in summer holidays were obviously higher than that in winter holidays (Figure 6A). Those results were consistent with the previous studies that temperature and humidity may jointly affect this childhood disease [45], and such interactive impact needs to be considered when evaluating at different school semesters.
Based on the whole analysis, we can identify two feasible strategies for efficiently lessening the basic reproduction number, i.e., for containing HFMD epidemics. The first measure is to increase social distance. The transmission of enterovirus and coxsackievirus is most efficient in crowded settings and, therefore, most countries and regions, including Malaysia, Singapore, Taiwan, Hong Kong, and China, have adopted social distancing measures, such as closures of childcare facilities and schools, and cancellation of public functions involving children [46,47]. In our study, the optimum timing for implementation was at the beginning of every school semester when the peer-group at school started to gather together. At this time, as soon as an infected was detected in the childcare facilities and schools, social distancing measures must be immediately implemented. Such control measure may lower the basic reproduction number and decrease the peak HFMD incidence. Secondly, health education focusing on personal hygiene and good sanitation, including frequent hand washing, proper disposal of soiled nappies, and disinfection of soiled surfaces with chlorinated (bleach) disinfectants, should be addressed in the spring semester. There are two possible reasons why we highlight this measure. On the one hand, physical activity of adolescence is lower during winter and increases during warmer months and therefore, weather conditions may be associated with behavioral patterns that lead to increased contact among children, thereby facilitating the spread of HFMD infection [9,48,49]. On the other hand, weather changes have been credited with improving the duration of viral survival and virulence in the environment [50,51,52], and such changes increase the number of opportunities for hosts to become infected as a result of contact with contaminated surfaces, aspiration of contaminated fomites, or inhalation of aerosols produced by the coughs of infected individuals [53,54]. More evidence on the underlying biological mechanism of HFMD transmission is needed to improve the understanding of the potential effect of meteorological factors.
A few limitations should be considered when interpreting the findings from this study. First, the estimation of model-free basic reproduction number was based on hospitalized HFMD cases and may not reflect the trend of population incidence reliably. This may be the reason why model-free annual basic reproduction number had a minor difference with model-based annual reproduction number. Second, we adopted the dates of summer and winter vacation from the Education Bureau, but some kindergartens started schools earlier. This may partly explain the rebound of model-free basic reproduction number slightly earlier than the start of the school year in September. Finally, despite that we developed a mathematical model to establish causation between HFMD transmission, we may not incorporate all potential driving factors.
The authors are very grateful to the anonymous referees for their valuable comments and suggestions. This research was supported by the National Natural Science Foundation of P. R. China (Grant Nos. 11771448, 61672013, 61772017) and Huaian Key Laboratory for Infectious Diseases Control and Prevention (HAP201704).
The authors declare that they have no competing interests.
Supplementary A: Estimation of the daily and annual model-free basic reproduction number
Similarly to [37], we assumed that serial interval (the time between the onset of symptoms in a primary case and the onset of symptoms in secondary cases) is equal to the sum of L and D, where L represents the duration of incubation and D denotes the duration of infection. As a result, [55] obtained the probability density function of incubation period of HFMD whose best-fitted function was gamma distribution. Concretely, since the estimated median incubation period was 5.4(95% CI 4.4-6.5) days, based on the parameters in [55], we know the probability density function of incubation L was fIncubation(x)=ba11Γ(a1)xa1−1e−b1x if x≥0, and fIncubation(x)=0 if x<0, in which a1=6.9637,b1=1.2217.
In addition, the disability length for symptomatic HFMD or infection was estimated at 5.7(95% CI 5.4-6.0) in [56,57]. Suppose that HFMD infection is also obeys a gamma distribution. Then, based on the parameters in [56,57], the probability density function of infection D was fInfection(x)=ba22Γ(a2)xa2−1e−b2x if x≥0, and fInfection(x)=0 if x<0, in which a2=1.0264,b2=0.1901.
Hence the joint probability density function of the sum of incubation period and infection duration (L+D), described as the probability density function of serial interval [26], was given by
fSI(x)=∫+∞0fIncubation(x−s)fInfection(s)ds=∫x0fIncubation(x−s)fInfection(s)ds=ba11ba22Γ(a1)Γ(a2)∫x0(x−s)a1−1e−b1(x−s)sa2−1e−b2sds. | (4.1) |
Due to the difficulty in solving the above integral, we used numerical integration method to obtain the result (see Figure 7). After that, the daily model-free basic reproduction number (Rfree0(day)) can be calculated based on the infected data and the probability density distribution function of serial interval through the method and Ecxel software package "epiEstim" in [26]. Taking the annual average of Rfree0(day), we can obtain the annual model-free basic reproduction number Rfree0(year).
Supplementary B: Model parameters, initial conditions and contact rate
B1. Model parameters and initial conditions in model (2.1)
According to Statistical Year book and Census data from Wenzhou Statistical Bureau [58], the average daily number of new babies from 2010 to 2015 is λ=334.459, and average daily mortality rate is d=1.904×10−5. Assume that the average daily deaths of people under 14 years old is d⋅N(t), the average daily number of people transferred out due to age growth is c⋅N(t), we can get dNdt=λ−(d+c)N. By solving above equation, we obtained that N(t)=λd+c+(N(0)−λd+c)e−(d+c)t, and N(0) is the initial susceptible. Fortunately, according to the Statistical Year book and Census data from Wenzhou Statistical Bureau, we find that the total population under 14 years in January 1, 2011 and December 31, 2015 is 1305373 and 1355698, respectively. Using these information, we can estimate the values of the parameters with the help of least square method. As a result, c=2.113×10−4 and N(0)=1292553. Because the number of susceptible is almost equal to total population, we have S(0)=N(0).
Overall, the definitions and values of all parameters and initial values of model (2.1) are listed in Table 4.
Parameter | Definition | Source | Mean value |
λ | Recruitment rate of susceptible (per day) | Calculated | 334.459 |
c | Transfer rate due to age growth (day) | Calculated | 2.11312×10−4 |
d | Per capita natural mortality rate (per day) | Calculated | 1.9040×10−5 |
1/σ | Mean incubation period (day) | [55] | 5.4 |
k | Quarantined rate (per day) | [3] | 0.007177 |
1/γ1 | mean recovery time of an infected person (day) | [56] | 5.7 |
1/γ2 | mean quarantined time (day) | [3] | 16 |
δ1 | Disease-induced mortality from infectious (per day) | [3] | 0 |
δ2 | Disease-induced mortality from quarantine | [3] | 0.024887 |
η | Progression rate of the recovered (per day) | [3] | 0.015 |
S(0) | Initial value of susceptible | Calculated | 1292553 |
E(0) | Initial value of exposed | [3] | 30 |
I(0) | Initial value of infectious | Reported | 81 |
Q(0) | Initial value of quarantine | [3] | 46 |
R(0) | Initial value of recovered | [3] | 46 |
B2. Contact rate
The contact rate between susceptible and infected, β(t), is a time-varying, continuous, positive function, which is mainly determined by school opening or not, and meteorological factors. The former can be divided into three phases, spring semester, fall semester and vacation (summer and winter vacation), while the mean temperature and relative humidity were selected to represent meteorological factors. For the sake of narrative convenience, several variables are expressed as bellow: δr(t) indicates whether it is the school vacation or not (1 represents school vacation while 0 denotes school opening); δs(t) indicates whether it is the spring semester or not (1 represents spring semester (including the weekend and other short-term holiday during spring semester) while 0 means not); δf(t) indicates whether it is the fall semester or not (1 represents fall semester (including the weekend and other short-term holiday during fall semester) while 0 means not); xs(t) refers to the cumulative opening days of spring semester in one year; xf(t) refers to the cumulative opening days of fall semester in one year. Since we can see from Figure 3 that the both variation between daily model-free basic reproduction number Rfree0(day) and school opening days in spring and fall semesters are similar to the form of probability density function similar of gamma distribution, the effect of the cumulative school opening days on the daily contact rate is defined as ys(t)=xs(t)b1−1e−b2xs(t) in spring semester and yf(t)=xf(t)b3−1e−b4xf(t) in fall semester.
Since [59,60] points out that there is no time lag effect between temperature and disease outbreak, to normalize the mean temperature per day, T(t) is defined as the average temperature per day divided by 40 because the highest daily mean temperature of Wenzhou city is no more than 40oC. However, [59] indicates that there exists time lag effect between humidity and disease outbreaks. Note that the daily model-free basic reproduction number Rfree0(day) is a measure of the severity of HFMD outbreaks. We first calculate the correlation coefficient ρˉHR(t−τ),Rfree0(day) between Rfree0(day) and the measured mean relative humidity with different time lag days ˉHR(t−τ), τ=0,1,⋅⋅⋅,30 (Figure 8). Then, according to the magnitude and stationarity of the correlation coefficients, we get that the outbreak of disease at time t is mainly driven by the average relative humidity between t−12 and t−23 days. Thus, the average relative humidity at time is defined as HR(t)=23∑τ=12r(τ)ˉHR(t−τ), in which r(τ)=ρˉHR(t−τ),Rfree0(day)/23∑s=12ρˉHR(t−s),Rfree0(day) is corresponding weight coefficient.
In summary, according to the seasonal variation, school opening and their interaction effect, we divide contact rate β(t) approximately into the following five parts:
(1) the interaction factor between school opening and meteorological effect during spring semester,
βs×m(t)=ys(t)(b5+b6T(t)+b7HR(t)). | (4.2) |
(2) the interaction factor between school opening and meteorological effect during fall semester,
βf×m(t)=yf(t)(b8+b9T(t)+b10HR(t)). | (4.3) |
(3) the meteorological factor during spring semester,
βms(t)=δs(t)(b11+b12T(t)+b13HR(t)). | (4.4) |
(4) the meteorological factor during fall semester,
βmf(t)=δf(t)(b14+b15T(t)+b16HR(t)). | (4.5) |
(5) the meteorological factor during winter and summer holiday,
βmr(t)=δr(t)(b17+b18T(t)+b19HR(t)). | (4.6) |
Therefore, β(t) can be expressed as (2.2) in the main text.
Supplementary C: Parameter estimation in contact rate
The variance of measured components, I(t), was given by an inverse gamma distribution with hyper-parameters (0.01,4), where 0.01 is the initial error variance which is updated by the inverse gamma distribution [61], and the small MCMC package provided in this website was used to estimate the parameters. When estimating unknown parameters and initial conditions for model (2.1), the following prior information were given: bi∈[0,+∞),i=1,2,3,4 and bj∈(=∞,+∞),j=5,6,⋅⋅⋅,19; the proposed density was chosen to be a multivariate normal distribution. These ranges were used to ensure good convergence of the MCMC chain.
By fitting model (2.1) to the reported HFMD cases from 2010 to 2015 in Wenzhou city, we estimate mean values of bi,i=1,2,⋅⋅⋅,19, and their standard deviations (SD) in Table 5. Figure 9 is the distribution diagram of each parameter.
Parameter | Mean±SD | Parameter | Mean±SD |
b1 | 0.3257±2.3079×10−3 | b11 | 0.2683±1.5734×10−3 |
b2 | 7.9764×10−9±6.0361×10−9 | b12 | 0.3240±2.6312×10−3 |
b3 | 2.7990±1.0528×10−2 | b13 | −0.3737±1.0570×10−3 |
b4 | 0.3086±1.0272×10−3 | b14 | 0.3377±3.0434×10−4 |
b5 | 2.5272±1.1446×10−2 | b15 | 0.0527±2.7698×10−4 |
b6 | −5.9314±2.0221×10−2 | b16 | −0.3014±3.8947×10−4 |
b7 | 1.7517±1.6270×10−2 | b17 | 0.0788±4.1117×10−4 |
b8 | −0.5895±7.7061×10−3 | b18 | 0.3952±2.1027×10−4 |
b9 | −0.4353±4.8126×10−3 | b19 | 0.0788±4.1117×10−4 |
b10 | 1.2471±1.5149×10−2 |
Supplementary D: Annual model-based basic reproduction number
According to [40], the annual model-based basic reproduction number Rmodel0(year) of model (2.1) can be defined by the spectral radius of the next infection operator. Note that the disease-free steady state of model (2.1) is E0=(λ/(d+c),0,0,0,0,λ/(d+c)), the contact rate β(t) can be approximated as an annual periodic function and the period is marked as ω. Using the notation and definition in [40], as well as the general calculation procedure in [39], for the approximate periodic model (2.1) we have
F(t)=(0β(t)0000000),V(t)=(d+σ00−σd+δ1+k+γ100−kd+δ2+γ2). |
In order to characterize model-based basic reproduction number of model (2.1), we consider the following linear ω-periodic system
dWdt=(−V(t)+F(t)θ)W,t∈R | (4.7) |
with parameter θ∈(0,∞). Let W(t,s,θ),t≥s,s∈R, be the evolution operator of the system (7.7) on R3. For each θ∈(0,∞), the matrix −V(t)+F(t)θ is cooperative. Thus, the linear operator W(t,s,θ) is positive in R3 for each t≥s,s∈R. By Theorem 2.1 in [40], Rmodel0(year) can be defined as the unique solution of ρ(W(ω,0,θ))=1 and Rmodel0(year)=θ in this case, in which ρ(W(ω,0,θ)) is the spectral radius of W(ω,0,θ).
[1] | S. Shao, X. Li, J. Cao, L. Yang, China's economic policy choices for governing smog pollution based on spatial spillover effects, Econ. Res. J., 51 (2016), 73–88. |
[2] |
D. Ghernaout, M. Aichouni, A. Alghamdi, Applying big data in water treatment industry: A new era of advance, Int. J. Adv. Appl. Sci., 5 (2018), 89–97. https://doi.org/10.21833/ijaas.2018.03.013 doi: 10.21833/ijaas.2018.03.013
![]() |
[3] |
J. Wu, S. Guo, J. Li, D. Zeng, Big data meet green challenges: Big data toward green applications, IEEE Syst. J., 10 (2016), 888–900. https://doi.org/10.1109/JSYST.2016.2550530 doi: 10.1109/JSYST.2016.2550530
![]() |
[4] |
Y. Liu, P. Failler, Z. Liu, Impact of environmental regulations on energy efficiency: A case study of China's air pollution prevention and control action plan, Sustainability, 14 (2022). https://doi.org/10.3390/su14063168 doi: 10.3390/su14063168
![]() |
[5] |
T. Wen, X. Niu, M. Gonzales, G. Zheng, Z. Li, S. L. Brantley, Big groundwater data sets reveal possible rare contamination amid otherwise improved water quality for some analytes in a region of Marcellus shale development, Environ. Sci. Technol., 52 (2018), 7149–7159. https://doi.org/10.1021/acs.est.8b01123 doi: 10.1021/acs.est.8b01123
![]() |
[6] | H. Wang, Can voluntary environmental policy approaches be effective in Chinese context, China Popul., Resour. Environ., 20 (2010), 89–94. |
[7] |
Z. Li, F. Zou, B. Mo, Does mandatory CSR disclosure affect enterprise total factor productivity, Econ. Res.-Ekonomska Istraživanja, 35 (2021), 4902–4921. https://doi.org/10.1080/1331677x.2021.2019596 doi: 10.1080/1331677x.2021.2019596
![]() |
[8] | H. Zhang, G. Liu, H. Hao, Review, assessment and recommendations on environmental policy in western China, China Popul., Resour. Environ., 23 (2013), 44–51. |
[9] |
M. H. Hansen, H. Li, R. Svarverud, Ecological civilization: Interpreting the Chinese past, projecting the global future, Global Environ. Change, 53 (2018), 195–203. https://doi.org/10.1016/j.gloenvcha.2018.09.014 doi: 10.1016/j.gloenvcha.2018.09.014
![]() |
[10] |
R. D. Mohr, Technical change, external economies, and the Porter hypothesis, J. Environ. Econ. Manage., 43 (2002), 158–168. https://doi.org/10.1006/jeem.2000.1166 doi: 10.1006/jeem.2000.1166
![]() |
[11] | P. West, P. Senez, Environmental assessment of the NAFTA: The mexican environmental regulation position, report prepared for the Province of British Columbia, Ministry of Economic Development, Small Bus. Trade, 1992 (1992), 69–70. |
[12] |
J. Liu, J. Ren, Y. Zhang, X. Peng, Y. Zhang, Y. Yang, Efficient dependent task offloading for multiple applications in MEC-cloud system, IEEE Trans. Mob. Comput., 22 (2021), 2147–2162. https://doi.org/10.1109/TMC.2021.3119200 doi: 10.1109/TMC.2021.3119200
![]() |
[13] |
M. Greenstone, R. Hanna, Environmental regulations, air and water pollution, and infant mortality in India, Am. Econ. Rev., 104 (2014), 3038–3072. https://doi.org/10.1257/aer.104.10.3038 doi: 10.1257/aer.104.10.3038
![]() |
[14] |
B. Laplante, P. Rilstone, Environmental inspections and emissions of the pulp and paper industry in Quebec, J. Environ. Econ. Manage., 31 (1996), 19–36. https://doi.org/10.1006/jeem.1996.0029 doi: 10.1006/jeem.1996.0029
![]() |
[15] |
D. Zavras, Healthcare access as an important element for the EU's socioeconomic development: Greece's residents' opinions during the COVID-19 pandemic, Natl. Account. Rev., 4 (2022), 362–377. https://doi.org/10.3934/nar.2022020 doi: 10.3934/nar.2022020
![]() |
[16] | H. Wang, Comparison and selection of environmental regulation policy in China: based on Bayesian model averaging approach, China Popul., Resour. Environ., 26 (2016), 132–138. |
[17] |
A. J. Hedley, C. M. Wong, T. Q. Thach, S. Ma, T. H. Lam, H. R. Anderson, Cardiorespiratory and all-cause mortality after restrictions on sulphur content of fuel in Hong Kong: an intervention study, Lancet, 360 (2002), 1646–1652. https://doi.org/10.1016/S0140-6736(02)11612-6 doi: 10.1016/S0140-6736(02)11612-6
![]() |
[18] |
Z. Li, J. Zhu, J. He, The effects of digital financial inclusion on innovation and entrepreneurship: A network perspective, Electron. Res. Arch., 30 (2022), 4697–4715. https://doi.org/10.3934/era.2022238 doi: 10.3934/era.2022238
![]() |
[19] |
S. Tanaka, Environmental regulations on air pollution in China and their impact on infant mortality, J. Health Econ., 42 (2015), 90–103. https://doi.org/10.1016/j.jhealeco.2015.02.004 doi: 10.1016/j.jhealeco.2015.02.004
![]() |
[20] |
Y. Chen, A. Ebenstein, M. Greenstone, H. Li, Evidence on the impact of sustained exposure to air pollution on life expectancy from China's Huai River policy, PNAS, 110 (2013), 12936–12941. https://doi.org/10.1073/pnas.1300018110 doi: 10.1073/pnas.1300018110
![]() |
[21] |
Y. Liu, C. Ma, Z. Huang, Can the digital economy improve green total factor productivity? An empirical study based on Chinese urban data, Math. Biosci. Eng., 20 (2023), 6866–6893. https://doi.org/10.3934/mbe.2023296 doi: 10.3934/mbe.2023296
![]() |
[22] |
H. Wang, R. Zhang, Effects of environmental regulation on CO2 emissions: An empirical analysis of 282 cities in China, Sustainable Prod. Consumption, 29 (2022), 259–272. https://doi.org/10.1016/j.spc.2021.10.016 doi: 10.1016/j.spc.2021.10.016
![]() |
[23] |
L. Yang, K. L. Wang, Regional differences of environmental efficiency of China's energy utilization and environmental regulation cost based on provincial panel data and DEA method, Math. Comput. Model., 58 (2013), 1074–1083. https://doi.org/10.1016/j.mcm.2012.04.004 doi: 10.1016/j.mcm.2012.04.004
![]() |
[24] |
H. Tanaka, C. Tanaka, Sustainable investment strategies and a theoretical approach of multi-stakeholder communities, Green Finance, 4 (2022), 329–346. https://doi.org/10.3934/gf.2022016 doi: 10.3934/gf.2022016
![]() |
[25] |
T. Li, J. Wen, D. Zeng, K. Liu, Has enterprise digital transformation improved the efficiency of enterprise technological innovation? A case study on Chinese listed companies, Math. Biosci. Eng., 19 (2022), 12632–12654. https://doi.org/10.3934/mbe.2022590 doi: 10.3934/mbe.2022590
![]() |
[26] |
S. Ren, X. Li, B. Yuan, D. Li, X. Chen, The effects of three types of environmental regulation on eco-efficiency: A cross-region analysis in China, J. Clean Prod., 173 (2018), 245–255. https://doi.org/10.1016/j.jclepro.2016.08.113 doi: 10.1016/j.jclepro.2016.08.113
![]() |
[27] |
J. Liu, Y. Zhang, J. Ren, Y. Zhang, Auction-based dependent task offloading for IoT users in edge clouds, IEEE Internet Things, 2022 (2022). https://doi.org/10.1109/JIOT.2022.3221431 doi: 10.1109/JIOT.2022.3221431
![]() |
[28] |
Kanupriya, Indian textile sector, competitiveness, gender and the digital circular economy: A critical perspective, Natl. Account. Rev., 4 (2022), 237–250. https://doi.org/10.3934/nar.2022014 doi: 10.3934/nar.2022014
![]() |
[29] | Z. Chu, C. Bian, C. Liu, J. Zhu, Evolutionary game analysis on haze governance in Beijing-Tianjin-Hebei: Based on a simulation tool for proposed environmental regulation policies, China Popul., Resour. Environ., 28 (2018), 63–75. |
[30] | B. Ma, X. Lv, X. Chen, X. Chen, Analysis of atmospheric emission monitoring big data of thermal power plants and study on the policy impact, China Popul., Resour. Environ., 29 (2019), 73–79. |
[31] |
Y. Liu, L. Chen, L. Lv, P. Failler, The impact of population aging on economic growth: A case study on China, AIMS Math., 8 (2023), 10468–10485. https://doi.org/10.3934/math.2023531 doi: 10.3934/math.2023531
![]() |
[32] | Y. Xu, Changes and developing trend of China's marine governance policy: An empirical research based on 161 policy texts from 1982 to 2015, China Popul., Resour. Environ., 28 (2018), 165–176. |
[33] | Y. Wang, M. Li, Study on local government attention of ecological environment governance: Based on the text analysis of government work report in 30 provinces and cities (2006–2015), China Popul., Resour. Environ., 27 (2017), 28–35. |
[34] |
W. Guo, B. Xi, C. Huang, J. Li, Z. Tang, W. Li, et al., Solid waste management in China: Policy and driving factors in 2004–2019, Resour. Conserv. Recycl., 173 (2021), 105727. https://doi.org/10.1016/j.resconrec.2021.105727 doi: 10.1016/j.resconrec.2021.105727
![]() |
[35] |
T. Li, X. Li, G. Liao, Business cycles and energy intensity. Evidence from emerging economies, Borsa Istanbul Rev., 22 (2022), 560–570. https://doi.org/10.1016/j.bir.2021.07.005 doi: 10.1016/j.bir.2021.07.005
![]() |
[36] |
P. Zweifel, Expanding insurability through exploiting linear partial information, Data Sci. Finance Econ., 2 (2022), 1–16. https://doi.org/10.3934/dsfe.2022001 doi: 10.3934/dsfe.2022001
![]() |
[37] |
Y. Liu, Z. Li, M. Xu, The influential factors of financial cycle spillover: Evidence from China, Emerg. Markets Finance Trade, 56 (2019), 1336–1350. https://doi.org/10.1080/1540496x.2019.1658076 doi: 10.1080/1540496x.2019.1658076
![]() |
[38] | H. Yin, Z. Xu, Discussion on China's single-factor water quality assessment method, Water Purif. Technol., 27 (2008), 1–3. |
[39] |
Y. Wang, J. Yang, D. Xu, Environmental impact assessment using the evidential reasoning approach, Eur. J. Oper. Res., 174 (2006), 1885–1913. https://doi.org/10.1016/j.ejor.2004.09.059 doi: 10.1016/j.ejor.2004.09.059
![]() |
[40] |
Y. Zhang, X. Deng, D. Wei, Y. Deng, Assessment of E-Commerce security using AHP and evidential reasoning, Expert Syst. Appl., 39 (2012), 3611–3623. https://doi.org/10.1016/j.eswa.2011.09.051 doi: 10.1016/j.eswa.2011.09.051
![]() |
[41] |
J. Gorin, Assessment as evidential reasoning, Teach. Coll. Rec., 116 (2014), 1–26. https://doi.org/10.1177/016146811411601101 doi: 10.1177/016146811411601101
![]() |
[42] |
Z. Li, Z. Huang, Y. Su, New media environment, environmental regulation and corporate green technology innovation: Evidence from China, Energy Econ., 119 (2023). https://doi.org/10.1016/j.eneco.2023.106545 doi: 10.1016/j.eneco.2023.106545
![]() |
[43] |
X. Si, C. Hu, J. Yang, Q. Zhang, On the dynamic evidential reasoning algorithm for fault prediction, Expert Syst. Appl., 38 (2011), 5061–5080. https://doi.org/10.1016/j.eswa.2010.09.144 doi: 10.1016/j.eswa.2010.09.144
![]() |
[44] |
Y. Liu, P. Failler, Y. Ding, Enterprise financialization and technological innovation: Mechanism and heterogeneity, PLoS One, 17 (2022), e0275461. https://doi.org/10.1371/journal.pone.0275461 doi: 10.1371/journal.pone.0275461
![]() |
[45] |
J. Lein, Applying evidential reasoning methods to agricultural land cover classification, Int. J. Remote Sens., 24 (2003), 4161–4180. https://doi.org/10.1080/0143116031000095916 doi: 10.1080/0143116031000095916
![]() |
[46] |
C. Yu, Z. Li, Z. Yang, A universal calibrated model for the evaluation of surface water and groundwater quality: Model development and a case study in China, J. Environ. Manage., 163 (2015), 20–27. https://doi.org/10.1016/j.jenvman.2015.07.011 doi: 10.1016/j.jenvman.2015.07.011
![]() |
[47] | G. Li, X. Liu, Z. Liu, W. Guo, Water quality assessment of main rivers in Tianjin based on principal component analysis and water quality identification index, J. Ecol. Rural Environ., 27 (2011), 27–31. |
[48] |
J. B. Yang, D. L. Xu, Evidential reasoning rule for evidence combination, Artif. Intell., 205 (2013), 1–29. https://doi.org/10.1016/j.artint.2013.09.003 doi: 10.1016/j.artint.2013.09.003
![]() |
[49] |
J. B. Yang, M. G. Singh, An evidential reasoning approach for multiple-attribute decision making with uncertainty, IEEE Trans. Syst. Man Cybern., 24 (1994), 1–18. https://doi.org/10.1109/21.259681 doi: 10.1109/21.259681
![]() |
[50] | H. Chu, W. Lu, L. Zhang, Application of artificial neural network in environmental water quality assessment, J. Agric. Sci. Technol., 15 (2013), 343–356. |
[51] |
Y. Liu, J. Liu, L. Zhang, Enterprise financialization and R & D innovation: A case study of listed companies in China, Electron. Res. Arch., 31 (2023), 2447–2471. https://doi.org/10.3934/era.2023124 doi: 10.3934/era.2023124
![]() |
[52] | G. Shafer, A Mathematical Theory of Evidence, Princeton University Press, 1976. |
1. | Fan Xia, Feng Deng, Hui Tian, Wei He, Yanni Xiao, Xiaodan Sun, Estimation of the reproduction number and identification of periodicity for HFMD infections in northwest China, 2020, 484, 00225193, 110027, 10.1016/j.jtbi.2019.110027 | |
2. | Zuqin Ding, Yong Li, Yongli Cai, Yueping Dong, Weiming Wang, Optimal Control Strategies of HFMD in Wenzhou, China, 2020, 2020, 1076-2787, 1, 10.1155/2020/5902698 | |
3. | Jinyu Wang, Sheng Li, Nonlinear effect of temperature on hand, foot, and mouth disease in Lanzhou, China, 2020, 99, 0025-7974, e23007, 10.1097/MD.0000000000023007 | |
4. | Pussadee Laor, Tawatchai Apidechkul, Siriyaporn Khunthason, Vivat Keawdounglek, Suntorn Sudsandee, Krailak Fakkaew, Weerayuth Siriratruengsuk, Association of environmental factors and high HFMD occurrence in northern Thailand, 2020, 20, 1471-2458, 10.1186/s12889-020-09905-w | |
5. | Chenxi Dai, Kaifa Wang, Adaptive Weighted Neighbors Method for Sensitivity Analysis, 2022, 14, 1913-2751, 652, 10.1007/s12539-022-00512-4 | |
6. | I. A. Moneim, G. A. Mosa, A realistic model for the periodic dynamics of the hand-foot-and-mouth disease, 2022, 7, 2473-6988, 2585, 10.3934/math.2022145 | |
7. | Zafer BEKİRYAZICI, Saime HASİMOGLU, Comparison of a random model of Hand-Foot-Mouth Disease model with Gaussian and Laplacian parameters, 2021, 2146-538X, 10.17714/gumusfenbil.973990 | |
8. | Chenxi Dai, Dongsheng Zhou, Bo Gao, Kaifa Wang, Eric H. Y. Lau, A new method for the joint estimation of instantaneous reproductive number and serial interval during epidemics, 2023, 19, 1553-7358, e1011021, 10.1371/journal.pcbi.1011021 | |
9. | Aili Wang, Duo Bai, Jingming He, Stacey R. Smith, Optimal control of bi-seasonal hand, foot and mouth disease in mainland China suggests transmission from children and isolating older infected individuals are critical, 2024, 89, 0303-6812, 10.1007/s00285-024-02141-5 | |
10. | Yingying Zhang, Ruohan Wang, Yanan Cai, Atila Bueno, Relaxation Oscillation in SEIR Epidemic Models with the Intrinsic Growth Rate, 2024, 2024, 1076-2787, 10.1155/2024/5373794 |
Mean TEM | Maximum TEM | Minimum TEM | Mean RHU | |
Maximum TEM | 0.977 | |||
Minimum TEM | 0.988 | 0.943 | ||
Mean RHU | 0.190 | 0.114 | 0.265 | |
Minimum RHU | 0.206 | 0.087 | 0.295 | 0.910 |
Note: TEM represents temperature; RHU represents relative humidity. |
Parameters | Rfree0(day)≥1 | Rfree0(day)<1 | P value |
Mean TEM | 19.3 (14.4, 23.2) | 20.0 (10.1, 27.7) | <0.001 |
Maximum TEM | 23.3 (18.3, 27.2) | 23.8 (14.1, 31.2) | <0.001 |
Minimum TEM | 16.8 (11.6, 20.8) | 17.1 (7.7, 25.0) | <0.001 |
Mean RHU | 78.5 (68.0, 86.0) | 76.0 (68.0, 84.0) | <0.001 |
Minimum RHU | 59.0 (46.0, 70.0) | 57.0 (45.0, 68.0) | 0.005 |
Note: TEM represents temperature; RHU represents relative humidity. |
Year | Rfree0(year) | Rmodel0(year) |
2010 | 1.034±0.012♯ | 1.026±0.012∗♯ |
2011 | 1.153±0.012♯ | 1.125±0.012∗♯ |
2012 | 1.133±0.012♯ | 1.119±0.012∗♯ |
2013 | 1.123±0.011♯ | 1.112±0.012∗♯ |
2014 | 1.149±0.013♯ | 1.123±0.012∗♯ |
2015 | 1.068±0.012♯ | 1.041±0.011∗♯ |
∗: P<0.05 represents the comparison with Rfree0(year); ♯: P<0.05 represents the comparison with threshold one |
Parameter | Definition | Source | Mean value |
λ | Recruitment rate of susceptible (per day) | Calculated | 334.459 |
c | Transfer rate due to age growth (day) | Calculated | 2.11312×10−4 |
d | Per capita natural mortality rate (per day) | Calculated | 1.9040×10−5 |
1/σ | Mean incubation period (day) | [55] | 5.4 |
k | Quarantined rate (per day) | [3] | 0.007177 |
1/γ1 | mean recovery time of an infected person (day) | [56] | 5.7 |
1/γ2 | mean quarantined time (day) | [3] | 16 |
δ1 | Disease-induced mortality from infectious (per day) | [3] | 0 |
δ2 | Disease-induced mortality from quarantine | [3] | 0.024887 |
η | Progression rate of the recovered (per day) | [3] | 0.015 |
S(0) | Initial value of susceptible | Calculated | 1292553 |
E(0) | Initial value of exposed | [3] | 30 |
I(0) | Initial value of infectious | Reported | 81 |
Q(0) | Initial value of quarantine | [3] | 46 |
R(0) | Initial value of recovered | [3] | 46 |
Parameter | Mean±SD | Parameter | Mean±SD |
b1 | 0.3257±2.3079×10−3 | b11 | 0.2683±1.5734×10−3 |
b2 | 7.9764×10−9±6.0361×10−9 | b12 | 0.3240±2.6312×10−3 |
b3 | 2.7990±1.0528×10−2 | b13 | −0.3737±1.0570×10−3 |
b4 | 0.3086±1.0272×10−3 | b14 | 0.3377±3.0434×10−4 |
b5 | 2.5272±1.1446×10−2 | b15 | 0.0527±2.7698×10−4 |
b6 | −5.9314±2.0221×10−2 | b16 | −0.3014±3.8947×10−4 |
b7 | 1.7517±1.6270×10−2 | b17 | 0.0788±4.1117×10−4 |
b8 | −0.5895±7.7061×10−3 | b18 | 0.3952±2.1027×10−4 |
b9 | −0.4353±4.8126×10−3 | b19 | 0.0788±4.1117×10−4 |
b10 | 1.2471±1.5149×10−2 |
Mean TEM | Maximum TEM | Minimum TEM | Mean RHU | |
Maximum TEM | 0.977 | |||
Minimum TEM | 0.988 | 0.943 | ||
Mean RHU | 0.190 | 0.114 | 0.265 | |
Minimum RHU | 0.206 | 0.087 | 0.295 | 0.910 |
Note: TEM represents temperature; RHU represents relative humidity. |
Parameters | Rfree0(day)≥1 | Rfree0(day)<1 | P value |
Mean TEM | 19.3 (14.4, 23.2) | 20.0 (10.1, 27.7) | <0.001 |
Maximum TEM | 23.3 (18.3, 27.2) | 23.8 (14.1, 31.2) | <0.001 |
Minimum TEM | 16.8 (11.6, 20.8) | 17.1 (7.7, 25.0) | <0.001 |
Mean RHU | 78.5 (68.0, 86.0) | 76.0 (68.0, 84.0) | <0.001 |
Minimum RHU | 59.0 (46.0, 70.0) | 57.0 (45.0, 68.0) | 0.005 |
Note: TEM represents temperature; RHU represents relative humidity. |
Year | Rfree0(year) | Rmodel0(year) |
2010 | 1.034±0.012♯ | 1.026±0.012∗♯ |
2011 | 1.153±0.012♯ | 1.125±0.012∗♯ |
2012 | 1.133±0.012♯ | 1.119±0.012∗♯ |
2013 | 1.123±0.011♯ | 1.112±0.012∗♯ |
2014 | 1.149±0.013♯ | 1.123±0.012∗♯ |
2015 | 1.068±0.012♯ | 1.041±0.011∗♯ |
∗: P<0.05 represents the comparison with Rfree0(year); ♯: P<0.05 represents the comparison with threshold one |
Parameter | Definition | Source | Mean value |
λ | Recruitment rate of susceptible (per day) | Calculated | 334.459 |
c | Transfer rate due to age growth (day) | Calculated | 2.11312×10−4 |
d | Per capita natural mortality rate (per day) | Calculated | 1.9040×10−5 |
1/σ | Mean incubation period (day) | [55] | 5.4 |
k | Quarantined rate (per day) | [3] | 0.007177 |
1/γ1 | mean recovery time of an infected person (day) | [56] | 5.7 |
1/γ2 | mean quarantined time (day) | [3] | 16 |
δ1 | Disease-induced mortality from infectious (per day) | [3] | 0 |
δ2 | Disease-induced mortality from quarantine | [3] | 0.024887 |
η | Progression rate of the recovered (per day) | [3] | 0.015 |
S(0) | Initial value of susceptible | Calculated | 1292553 |
E(0) | Initial value of exposed | [3] | 30 |
I(0) | Initial value of infectious | Reported | 81 |
Q(0) | Initial value of quarantine | [3] | 46 |
R(0) | Initial value of recovered | [3] | 46 |
Parameter | Mean±SD | Parameter | Mean±SD |
b1 | 0.3257±2.3079×10−3 | b11 | 0.2683±1.5734×10−3 |
b2 | 7.9764×10−9±6.0361×10−9 | b12 | 0.3240±2.6312×10−3 |
b3 | 2.7990±1.0528×10−2 | b13 | −0.3737±1.0570×10−3 |
b4 | 0.3086±1.0272×10−3 | b14 | 0.3377±3.0434×10−4 |
b5 | 2.5272±1.1446×10−2 | b15 | 0.0527±2.7698×10−4 |
b6 | −5.9314±2.0221×10−2 | b16 | −0.3014±3.8947×10−4 |
b7 | 1.7517±1.6270×10−2 | b17 | 0.0788±4.1117×10−4 |
b8 | −0.5895±7.7061×10−3 | b18 | 0.3952±2.1027×10−4 |
b9 | −0.4353±4.8126×10−3 | b19 | 0.0788±4.1117×10−4 |
b10 | 1.2471±1.5149×10−2 |