Spatiotemporal predictions of obesity prevalence in Chinese children and adolescents: based on analyses of obesogenic environmental variability and Bayesian model.

OBJECTIVE: To find variations in Chinese obesogenic environmental priorities from 2000 to 2011, predict spatiotemporal distribution of obesity prevalence aged 7-17 years in 31 provinces, and provide foundations for policy-makers to reduce obesity in children and adolescents.
METHODS: Based on data examination of provincial obesity prevalence aged 7-17 years from three rounds of China Health and Nutrition Surveys (in 9 [2000], 9 [2006], and 12 [2011] provinces) and corresponding years' environments in 31 provinces from China Statistical Yearbooks and other sources, 12 predictors were selected. We used 30 surveyed provinces in three rounds as training samples to fit three analytic models with partial least-square regressions and prioritized predictors by variable importance projection to find variations. And fitted a spatiotemporal prediction model with Bayesian analysis to infer in space-time.
RESULTS: Variations of obesogenic environmental priorities were found at different times. A Bayesian spatiotemporal prediction model with deviance information criterion of 155.60 and statistically significant (P < 0.05) parameter estimates of intercept (-717.0400, 95% confidence intervals [CI]: -1186.0300, -248.0480), year (0.3584, CI: 0.1245, 0.5924), square of food industry level (0.0003, CI: 0.0002, 0.0004), and log (healthcare) (5.3742, CI: 2.5138, 8.2347) was optimized. Totally inferred average obesity prevalence among children and adolescents were 2.23%, 5.11%, 10.77%, 12.20%, 13.99%, and 17.58% in 31 provinces in China in 2000, 2006, 2011, 2015, 2020, and 2030, respectively. Obesity in north and east of China clusters on predicted maps.
CONCLUSIONS: Obesity prevalence in children and adolescents in China is rapidly increasing, growing at 0.3584% annually from 2000 to 2011. From longitudinal observation, prevalence was significantly influenced by food industry ("Amplifier") and healthcare service ("Balancer"). Targeted interventions in north and east of China are pressing. Further researches on the mechanisms underlying the influence of food industry, healthcare service, and so on in children and adolescents are needed.

INTRODUCTION

The prevalence of obesity in children and adolescents has increased substantially over the past three or four decades in developed and developing countries,[1, 2] including China.[3-8] This implies greater reductions in quality of life, greater risk of social isolation, and greater economic disease burden.[3, 9] However, since it is predictable, preventable, controllable, and reversible,[10] the need for preventive action is increasingly being recognized.

Obesity is partly determined by genetic susceptibility, but obesogenic environments are necessary for its phenotypic expression.[11] In industrialized societies, the imbalance between health literacy and rapid food industry, marketing media, and transportation modernization, promotes higher energy density intakes and sedentary behaviors.[12, 13] Given the complexity of obesity in children and adolescents, policy-decision processes, should be more systematic, evidence-based, locally relevant, and stakeholder-informed, for more cohesive and comprehensive actions; but are often non-systematic in practice. In recent years, Analysis Grid for Environments Linked to Obesity (ANGELO) framework was used in Europe by policy-makers at national level to prioritize strategies for effective food and nutrition initiatives.[3, 14] Swinburn et al. suggested this 2 × 4 grid analysis framework to dissect the environments into 2 levels, micro and macro, using 4 dimensions: physical, economic, political, and sociocultural. In this framework, environmental factors influencing food intake and physical activities were characterized as “obesogenic” or “leptogenic” (promoting leanness).[15] Researches on obesogenic environments can help systematically and locally to determine and identify priority areas for population-based prevention of childhood obesity.[3] We need to develop systematic and integrated approaches to facilitate complementary and coordinated actions by various concerned parties, including governments, international organizations, civil society, and corporations.[11, 12]

With rapid industrialization, urbanization, and unbalanced economic development, in the background of complex terrain, distributions of obesity prevalence throughout China are heterogeneous in time-space. Wang et al. found regional differences in the prevalence of obesity among 7–18 year-old Chinese children in 2014.[8] Contrariwise, obesity was shown as closely related to geographic and natural environments,[16] economic activity, dietary patterns, food-production, transportation, social culture, policies and regulations.[17, 18] Often these factors are internally consistent within each of the Chinese provinces. Fighting obesity requires detailed global, national, provincial, and community-level pictures of the problem, especially its severity, size, distribution, and the underlying determinants. Few systematic studies have focused on the relationships between macro-environmental factors and obesity prevalence in Chinese children and adolescents at provincial level. There are perhaps some difficulties, including lack of complete data and suitable methods. Geographic information system (GIS) is increasingly used in thematic mapping and spatiotemporal epidemiology analysis of obesity and other public health problems.[19, 20]

In this study, we collected obesity data for children and adolescents aged 7–17 years from 3 rounds of China Health and Nutrition Surveys (CHNS), and corresponding years’ environmental data at provincial-level from China’s Statistical Yearbooks and other sources. After scanning these data we selected 12 predictors and constructed analytic models with partial least square regression (PLSR) and a predictive model with Bayesian analysis, which have been successfully used in medical fields to address small sample cases.[21-24] Furthermore, obesogenic environmental factors were prioritized for research and intervention, while thematic maps of obesity in 2000, 2006, 2011, 2015, 2020, and 2030 among children and adolescents in 31 Chinese provinces were inferred.

MATERIALS AND METHODS

Study Design, Samples, and Key Variables

We used data from 3 rounds of CHNS; 9 provinces (Liaoning, Heilongjiang, Jiangsu, Shandong, Henan, Hubei, Hunan, Guangxi, and Guizhou) were surveyed in 2000 and again in 2006, then 12 provinces (with Beijing, Shanghai, and Chongqing, in addition to the previous 9 provinces) in 2011. A multistage, stratified, random cluster sampling design was used throughout, to ensure a suitable representation of the population.[25, 26] Specialized and standardized trainings were adopted in advance and data on health, nutrition, and socioeconomic variables at the community, household, and individual levels were collected using questionnaires, dietary surveys, and physical measurements in September (autumn) every round. The Institutional Review Boards of the University of North Carolina at Chapel Hill and the National Institute for Nutrition and Health at Chinese Center for Disease Control and Prevention approved the survey protocols. All subjects completed written informed consent. Further details are in the published profile of this cohort study.[25, 26] We used the data of children and adolescents aged 7–17 with complete height and weight records to assess obesity prevalence at province-level in the corresponding age range; those without height or weight data were excluded from analysis. The final sample included 4 803 participants with 2 249, 1 138, and 1 416 participants in 2000, 2006, and 2011 surveys, respectively. The proportions of male, urban, 7–11 years-old, and total observed obesity prevalence were 53.13%, 23.83%, 41.80%, and 2.45% in 2000 (9 provinces); 52.99%, 26.54%, 49.21%, and 3.51% in 2006 (9 provinces); and 50.71%, 40.40%, 54.38%, and 9.11% in 2011 (12 provinces), respectively. Height was measured without shoes to the nearest 0.1 cm using portable Seca 206 wall-mounted metal tapes (Seca Ltd., Birmingham, UK) while weight was measured without shoes and in light clothing to the nearest 0.1 kg on calibrated Seca 880 beam scales. Body mass index (BMI), the key outcome variable, was defined according to the Body Mass Index Reference Norm for Screening Overweight and Obesity in Chinese Children and Adolescents (aged 7–18 years)[27] developed in 2004 by the Working Group on Obesity in China, organized by the International Life Science Institute Focal Point in China.[27] Relative to sex and age, the definition of obesity, BMI > 95th percentile, has been applied extensively in recent years in China.[27]

The dimensionalities of obesogenic environmental factors were determined by interviewing experts and reviewing literatures. Relevant environmental data, at provincial levels, were found in China Statistical Yearbooks,[28-30] China Health Statistics Yearbooks,[31-33] China Education Yearbooks,[34-36] China Light Industry Yearbooks,[37, 38] China Food industry Yearbooks,[39] the Sixth National Population Census of the People’s Republic of China in 2010,[40] China Electric Power Yearbooks, and other sources, for the corresponding years. We collected data on food consumption, including grains, meat, poultry, vegetables, edible oil, sugar, and dietary behaviors (dining out). We used the provincial annual production of dairy and soft drinks as representative variables for the developmental level of modern food industry. We checked environmental hygiene indicators, such as modified water sources, tap water, and lavatories. We examined economic indicators, including provincial gross domestic product (GDP), GDP per capita, and average incomes. We considered some demographic indicators, including sex ratio, child dependency ratio, urban-rural ratio, and urban population density. Other indicators considered were media and marketing (television coverage), transportation capacity, use of washing machine, healthcare service, mother’s illiteracy rate, sports facilities, built environments (per capita public green space and paved road area, etc.), sedentary life, and school policy scores. From 45 variables of 13 dimensionalities collected and examined, 12 environmental predictors of 11 dimensionalities were selected and analyzed. We did not analyze the data of Hong Kong, Macao, and Taiwan.

Statistical Analysis

In the form of descriptive statistical analysis, tables, and maps, we showed the distributions and variabilities of obesity prevalence and the environmental predictors in time and space at provincial-level among children and adolescents aged 7–17 years.

Compared with principal component analysis, use of PLSR (originally proposed by the econometrician Herman Wold in the 1960s) is more flexible and extensive, especially in medicine.[21, 22, 41, 42] PLSR integrates linear regression, principal component analysis, and canonical correlation analysis to handle multiple co-linear and small sample data, and is even applicable when the sample size is smaller than the number of independent variables with a stable regression. PLSR can use Variable Importance in Projection (VIP) for variable screening and priority ordering conveniently.[43-45]

Firstly, we examined correlations among independent variables (environmental factors) by dimensionality. Secondly, we used independent variables by dimensionality and the dependent variable (observed obesity prevalence among adolescent and children at province-level) as training samples to test for linearity, independence, normality, and equal variance; and to perform the single factor linear regressions. Twelve independent variables (as environmental predictors) of 11 dimensionalities, which eligible for the simple linear regression with statistical significance (P < 0.1), or supported by expert opinions or literatures, were selected. Thirdly, using three rounds of survey as three cross-sectional datasets, we fitted three analytic models with PLSR separately: Y = β0 + β1X1 + β2X2 + ⋯ + βX + ε in which the dependent variable Y was provincial obesity prevalence and the independent variables X1, X2, … X were the 12 environmental predictors we have just selected. We derived three groups of prioritizations of obesity environmental predictors by VIPs.

The analytic PLSR models can explain the observed samples well, but it cannot predict unobserved samples accurately because of overfitting. In addition, a PLSR model cannot give 95% CI of parameter estimates. Bayesian analysis provides robust predictions and inferences that are conditional on the data, without reliance on asymptotic approximation, which treating parameters as unknown random variables and making inferences based on the posterior distributions of the parameters. Moreover, small sample inference proceeds in the same manner as if one had a large sample. It also can estimate any function of parameters directly, without using the “plug-in” method and provide interpretable answers with 95% CI and a model assessment tool of deviance information criterion.

Fourthly, based on analysis of the variations of obesogenic environmental priorities in 3 different rounds, we included all 30 provincial-level observations in 3 survey rounds as a training sample to optimize a space-time predictive model with fewer independent variables by Bayesian analysis. And predicted obesity prevalence for the provinces and for the years not surveyed, as well as a profile of obesity prevalence in 31 provinces in China in 2000, 2006, 2011, 2015, 2020, and 2030. The posterior distribution of the model parameters was estimated using Markov chain Monte Carlo simulations. By defaults in SAS, PROC GENMOD computes three convergence diagnostics: the lag1, lag5, lag10, and lag50 autocorrelations; Geweke diagnostic statistics; and effective sample sizes. We have 95% confidence to say that the intervals contain the overall parameter estimates with statistical significance.

Fifthly, in the Sixth National Population Census of the People’s Republic of China in 2010, we took the children and adolescent proportions as weights, these were the corresponding populations in each of the 31 provinces accounting for the total corresponding population of 31 provinces of China. We calculated the total inferred average prevalence of obesity among the children and adolescents in the 31 provinces.

Sixthly, we made six thematic maps of predicted obesity prevalence among children and adolescents in 31 provinces. We used SAS 9.4 (SAS Institute Inc., Cary, NC, USA) and ArcGIS Desktop 10.3 (Environmental System Research Institute, ESRI) throughout the process.

Finally, we compared the inferred 2015 obesity prevalence with the newest CHNS 2015 round result on the children and adolescents’ obesity prevalence in 15 provinces of China (Shaanxi, Yunnan, and Zhejiang in addition to the 12 provinces surveyed in 2011).

RESULTS

Distributions and Variations for Selected Independent and Dependent Variables in 2000, 2006, and 2011

Table 1 shows the descriptive analysis (minimum, median, maximum, mean, and standard deviation) for selected variables and the number of provinces reported for.

Table 1.

Distribution and variability in the selected independent and dependent variables in years 2000, 2006, and 2011

	2011							2006						2000
Variable	Label	Province	Min	Median	Max	Mean	SD	Province	Min	Median	Max	Mean	SD	Province	Min	Median	Max	Mean	SD
Provincial GDP (100 million yuan)[1]	gdp*##	31	507.46	10368.60	46013.06	14098.13	11401.35	31	251.21	4179.52	22366.54	6380.29	5451.52	31	105.61	1962.98	8464.31	2828.10	2217.83
Children dependency ratio (%)[2]	chilrati*	31	10.41	22.70	38.10	22.75	7.06	31	11.39	29.92	44.65	28.05	8.30	31	18.37	37.47	57.78	35.59	8.29
Modified water source ratio (%)[3]	mdwaterr**	30	78.50	97.15	100.00	94.50	6.13	30	58.00	96.50	100.00	92.80	16.38	31	18.50	93.40	100.00	88.03	16.13
Healthcare (no. doctor/1000 persons)[4]	hlthcare	31	1.04	1.73	5.24	1.99	0.82	31	1.01	1.53	4.28	1.72	0.68	31	1.11	1.65	4.76	1.94	0.80
Illiteracy rate of mothers (%)[5]	m_illi	31	2.75	7.09	38.13	8.71	6.67	31	6.37	15.68	55.76	17.71	10.44	31	9.25	22.27	74.60	23.35	12.69
School policy score (point)[6]	schplcy¶¶	31	0.00	3.00	9.00	2.61	2.43	31	0.00	3.00	6.00	3.10	2.37	31	0.00	2.00	6.00	2.23	2.11
Food industry level (10 thousand tons)[7]	foodindu***##	31	0.46	31.00	345.36	69.61	85.09	31	0.33	16.48	307.53	38.84	59.92	31	0.00	1.24	15.49	2.23	3.26
Transport capability (100 million ton-km)[8]	ft_kilo###¶¶¶	31	38.50	2840.00	18918.20	4210.03	4143.69	31	40.71	1167.38	12592.99	2307.86	3016.64	31	13.80	795.00	5454.70	1147.33	1265.68
Media (television coverage %)[9]	media*	31	91.40	97.70	100.00	97.29	2.13	31	86.00	95.80	100.00	95.34	3.38	31	63.50	90.40	100.00	90.10	7.52
Obesity prevalence (%)[10]	pctn_4[@@@]	12	0.00	11.25	24.66	10.65	6.57	9	0.78	4.03	10.13	4.05	2.88	9	0.51	1.71	8.51	2.59	2.44
Edible oil per capita (g)[11]	oilrpd***	31	8.03	18.11	29.42	17.98	6.25	31	8.94	15.17	29.85	16.69	5.89	31	9.12	16.41	27.67	17.31	5.14
Vegetables per capita (g)[12]	vegerpd#	31	48.11	221.64	411.86	240.63	86.45	31	63.10	247.89	475.56	264.78	102.85	31	70.66	284.68	575.34	280.58	116.51
Washing machines (units/ 100 households)[13]	washmr*	31	8.47	65.75	101.47	62.42	25.74	31	3.89	40.90	97.07	44.18	25.04	31	2.45	20.87	87.60	27.31	22.78

Note.

Gross domestic product by province in units of 100 million yuan;

Ratio of people aged 0–14 to people aged 15–64 by province;

Percentage of improved water source use by province;

Average number of doctors per thousand people by province;

Percentage of illiterate and semi-literate mothers among all those aged 15 and over;

Comprehensive point score for school policies on physical education, health education, and nutrition education;

Dairy products were used to represent the developmental level of the food industry by province;

Cargo transportation by province in units of 100 million freight ton-kilometers;

TV coverage was used to represent the developmental level of the marketing industry by province;

Prevalence of obesity in children/adolescents aged 7–17 in 2011 by province;

Average daily oil consumption in grams per person by province;

Average daily consumption of vegetables in grams per person by province;

Number of washing machines owned per 100 households by province.

The total average obesity prevalence among 2 249 observed participants in 9 provinces was 2.45% in the 2000 round, 3.51% among 1 138 observed participants in 9 provinces in the 2006 round, and 9.11% in 1 416 observed participants in 12 provinces in the 2011 round with statistically significant increasing of obesity prevalence in Cochran-Armitage trend test (Z=−8.8707, p <0.0001).

*, #, and ¶ represent the statistical significance level of the variables in 2011, 2006, and 2000 respectively. 1 indicates P <0.1; 2 P <0.05; and 3 P <0.01 in simple linear regression.

Distribution and Variations in Observed Obesity Prevalence and Environmental Predictors by Provinces in 2000, 2006, and 2011

Figure 1 is a group of maps showing the distribution and variations in the observed obesity prevalence and environmental predictors by provinces.

Figure 1.

Distribution and variations in the observed obesity prevalence and environmental predictors by province in 2000, 2006, and 2011

Note: From A to L, 12 obesogenic environmental predictors are food industry level (in 10 thousand tons), transport capability (in 100 million ton-km), edible oil per capita (g), provincial GDP (in 100 million yuan), modified water source ratio (%), school policy score (point), children dependency ratio (%), media (television coverage in %), washing machine (in units per 100 households), vegetables per capita (g), healthcare service(in doctors/1000 persons), and mother’s illiteracy rate (%).

Analytic Model

We developed parameter estimates, prioritizations of environmental obesity factors (Table 2) and found two environmental factor types. Obesity-promoting factors include edible oil per capita, food industry level, modified water source ratio, provincial GDP, washing machines, media (television coverage), transportation capability, and healthcare service while obesity-impeding factors include child dependency ratio, mother’s illiteracy rate, vegetables per capita, and school policy score. In 2006 and 2011 rounds, the dimensionalities and the direction of the impact of these two factors types on obesity remained consistent. However, in 2000, four predictors (edible oil per capita, washing machines, mother’s illiteracy rate, and school policy score) showed inverse directions of impact on obesity compared with its following two rounds (2006 and 2011).

Table 2.

Parameter Estimates of Analytical Models with PLSR

	2011			2006			2000
Parameter	Original	Standardized	VIP	Original	Standardized	VIP	Original	Standardized	VIP
Intercept	−36.677798250	0.0000000000		−7.469752897	0.0000000000		−2.719622074	0.0000000000
gdp	0.000060460	0.1058703859	1.03560	0.000073516	0.1503879608	1.43884	0.000215251	0.2040779252	1.25047
oilrpd	0.163331260	0.1569353660	1.53510	0.032563828	0.0751809796	0.71929	−0.049534100	−0.1185512105	0.72641
vegerpd	−0.006034430	−0.0692222505	0.67712	−0.003831191	−0.1124128785	1.07551	−0.001228249	−0.0462365202	0.28331
chilrati	−0.079600090	−0.1043472277	1.02070	−0.032019615	−0.0954118052	0.91285	−0.038914199	−0.1032788147	0.63283
foodindu	0.012121670	0.1412824089	1.38199	0.012749799	0.1456202293	1.39322	0.045882704	0.0935850534	0.57343
ft_kilo	0.000103700	0.0853258715	0.83464	0.000328702	0.1747679846	1.67209	0.001207472	0.3444043397	2.11031
media	0.314468060	0.0987461071	0.96591	0.072707207	0.0770338505	0.73702	0.004288565	0.0098473562	0.06034
washmr	0.025566100	0.1012743674	0.99064	0.012011379	0.0896570105	0.85779	−0.001630167	−0.0128682070	0.07885
m_illi	−0.127817320	−0.0845105340	0.82666	−0.014526776	−0.0367419607	0.35153	0.018234197	0.0597247493	0.36596
schplcy	−0.096462360	−0.0415767208	0.40669	−0.022090251	−0.0175449803	0.16786	0.423782344	0.2982890505	1.82774
mdwaterr	0.136225200	0.1158640006	1.13335	0.035757311	0.1024325731	0.98002	0.041515928	0.1747188573	1.07057
hlthcare	0.360857910	0.0659804330	0.64540	0.458255695	0.0566927202	0.54241	0.116509311	0.0200826690	0.12305

Abbreviations: PLSR, partial least square regressions; VIP, Variable Importance in Projection.

In the cross-sectional 3 rounds data, we found variations in Chinese obesogenic environmental priorities at provincial level (Figure 2). In descending order of VIPs, the most important four factors of children obesity respectively were transportation capability, school policy score, provincial GDP, and modified water source ratio in 2000; transportation capability, provincial GDP, food industry level, and vegetables per capita in 2006; and edible oil per capita, food industry level, modified water source ratio, and provincial GDP in 2011.

Figure 2.

Variations of Chinese obesogenic environmental priorities at provincial level in 2000, 2006 and 2011

In the three different survey years, predictors of continuously increasing importance were edible oil per capita, children dependency ratio, food industry level, media, washing machines, mother’s illiteracy rate, healthcare service; predictor of continuously decreasing importance was transportation capability. Predictors lacking steady trend were provincial GDP, vegetables per capita, school policy score, and modified water source ratio.

Optimization of the Bayesian Space-time Predictive Model and the Spatiotemporal Inference

Using all observed 30 provincial-level observations in the 3 rounds as a training sample, we fitted a Bayesian space-time predictive model with a deviance information criterion of 155.60 and with statistically significant parameter estimates (P <0.05) (Table 3). In the posterior distribution test results, no obvious internal correlations among samples and Markov chain Monte Carlo simulations converged rapidly. From longitudinal observation, the prevalence of obesity among children and adolescents increased by 0.3584% annually on average between 2000 and 2011, and was more influenced by the food industry level and healthcare service. The predictive model yielded spatiotemporal inferences on obesity prevalence distributions in the 31 provinces of China in 2000, 2006, 2011, 2015, 2020, and 2030; and inferences in 2015, 2020, and 2030 were based on the obesogenic environmental data of 2011.

Table 3.

Parameter Estimates of the Bayesian Space-time Predictive Model

Parameter	Estimate	Standard Error	Wald 95% Confidence Intervals
Intercept	−717.0400	239.2861	−1186.0300	−248.0480
Year	0.3584	0.1194	0.1245	0.5924
Square of food industry	0.0003	0.0000	0.0002	0.0004
Log(healthcare)	5.3742	1.4594	2.5138	8.2347

Note: The parameters were estimated by maximum likelihood. (Deviance information criterion = 155.60)

Totally Inferred Average Prevalence of Obesity in the 31 Provinces of China

The inferred total average obesity prevalence in all 31 provinces were 2.23%, 5.11%, 10.77%, 12.20%, 13.99%, and 17.58% in 2000, 2006, 2011, 2015, 2020, and 2030 respectively.

Thematic Maps Predicted

Obesity in children and adolescents showed a continuously increasing trend, and obesity in northern and eastern provinces of China was clustered in the predicted maps. In 2000, there were only 3 provinces with the highest obesity levels, from 5% to 10%. In 2006, 3 provinces had a level above 10% and 1 above 15%. In 2011, 9 provinces had a level above 10% and 4 above 15%. In 2015, 14 provinces had a level above 10% and 4 above 15%. In 2020, 20 and 6 provinces would have a level above 10% and 15%, respectively. In 2030, all 31 provinces would have a level above 10% while 15 above 15% (Figure 3).

Figure 3.

Predicted thematic maps of obesity prevalence aged 7–17 in 2000, 2006, 2011, 2015, 2020, and 2030 in 31 provinces of China

Note: We did not analyze the data of Hong Kong, Macao, and Taiwan in this study.

Comparison between Predicted and Newly Observed Results in CHNS 2015 Round

Using the newly observed obesity prevalence among children and adolescents in CHNS 2015 round in 15 provinces, we compared them with their counterparts inferred and found that these two results were almost consistent across regional trends, in provinces on obesity prevalence near or above 10% and in clustering of obesity in the north (Beijing, Liaoning, Heilongjiang, Shaanxi) and east (Shandong, Shanghai, Jiangsu) of China (Supplementary Figure 1).

DISCUSSION

In 3 rounds of surveys, priorities had changed. During the last 10 years, the influence of these 7 factors (edible oil, children dependency ratio, food industry level, media, washing machines, mother’s illiteracy rate, and healthcare service) was continuously increasing. According to the predictive model fitted, we found the following quantitative relations. The prevalence of obesity among children and adolescents increased by 0.3584% annually on average from 2000 to 2011 and was more influenced by food industry level and healthcare service. Food industry level seems to be an “Amplifier” of childhood obesity, with “1” being the critical value, which actually means when the development level of the food industry reaches or exceeds the scale of 10,000 tons of milk products per year in a province the “Amplifier” can accelerate prevalence of obesity more quickly in the form of square of food industry level. Healthcare service level seems to be a “Balancer” of childhood obesity, with “1” being the critical value, which actually means when the number of doctors per thousand persons is less than 1 or more than 1, healthcare service level shows a negative or positive effect on the nutritional status of the population in form of logarithmic relationship. From a qualitative perspective, whether healthcare practitioners focus on “disease treatment” or “nutrition and health integrity” may fundamentally determine the direction, nature and future trends of the role of healthcare in the obesity epidemic. Hence, we can it a “Balancer.” The inferred total average obesity prevalence in 31 provinces of China were 2.23%, 5.11%, 10.77%, 12.20%, 13.99%, and 17.58% in 2000, 2006, 2011, 2015, 2020, and 2030 respectively. Compared with the newly observed CHNS 2015 data in the 15 provinces, predictive results of 2015 were consistent with the observed results, in provinces with obesity rates near or above 10% and in clustering trend of obesity in the north and east of China. We need to understand major contradictions timeously and identify modifiable determinants to control and reverse the epidemic of obesity.

The food industry level combined with media and transportation capacity forms a powerful engine to promote consumption of processed foods and snacks.[15, 46, 47] In this study, food industry level was a positive correlation factor in the 3 rounds with continuously increasing importance. Modern food industry provides abundant inexpensive food, increases the accessibility to food, and plays important roles in reducing protein energy malnutrition. However, with sustained economic growth, increased purchasing power, and decreased physical activity; the excessive oil, salt, and sugar in processed foods, snacks, and beverage contribute to obesity in children and adults.[15, 46] Food marketing and advert directed at children are extensive, widespread, and global; and are focused on fast food, processed food, snacks and beverages. Evidences also imply that television marketing influences children’s food preference, requests, and purchase behaviors. The World Health Organization has provided suggestions in a framework for implementing a set of recommendations on the marketing of foods and non-alcoholic beverages to children (2012), to promote responsible marketing. In our study, media, transportation capability (which determines the access to foods especially during eras or in regions of food shortages), and healthcare service, were positive correlation factors in the 3 rounds with continuously increasing importance.[46]

However, healthcare service not only helps improve residents’ nutritional status but can also help to prevent obesity.[48] Efforts should be directed at setting nutritional guidelines to educate doctors on how to prevent diet-related chronic non-communicable diseases, and not just to treat the resulting diseases only.

We also found other positive correlations, such as provincial-GDP, modified water sources, edible oil, and washing machine use. The studies by Popkin et al. found that when income increases, households from the lower income strata tend to add more oil and sugar to diets that are already high in energy density.[49] Su et al., based on the 1997–2009 Chinese data on children and adolescents found that the percentage of energy intake from fat was positively related to overweight and obesity (odds ratio = 1.01; 95% CI 1.00, 1.02; P = 0.006).[50] Another study found that by 2009 57.0% of children and adolescents aged 7–17 years consumed over 30.0% of their energy from fat, which is the upper limit recommended in Chinese Dietary Reference Intakes.[51] In our study the importance of edible oil in 3 rounds was continuously increasing.

We also found negative correlations between obesity in children and the provincial children dependency ratio, mother’s illiteracy rate, school policy score, and consumption of vegetables. Children dependency ratio reflects the allocation of resources, fertility rate, work participation rate, and change in the amount of housework undertaken.[52, 53] We believe that maintaining a moderately higher fertility rate might lower the obesity prevalence in children. Further research is needed. Mother’s illiteracy rate may affect the availability of food in households and mothers’ practices in feeding children.[54] We consider that the negative effect of mother’s illiteracy rate may be an indirect reflection of the phenomenon that mothers’ high employment rate promotes children’s obesity.[55] The school policy in our study showed a moderately antagonistic effect on child obesity. Children and adolescents are in the key period of knowledge learning, cognitive formation, and eating habits shaping. Education on nutrition, health, and activities will play a positive role, and may also have long-term effects on their offspring.[3, 56, 57]

Wang Shuo et al. in 2014, analyzed obesity data from the Chinese National Survey on Students Constitution and Health among those aged 7–18 years and found regional differences in the prevalence of overweight and obesity. The occurrence of obesity in northern and eastern provinces of China shows obvious clustering. The spatial distribution trend found in our study was consistent with theirs.[8]

In this study, the annual increase in rate of obesity prevalence among Chinese children and adolescents ranged from 0.12%−0.59% from 2000 to 2011. This was close to Wang et al.’s results of the students’ physical examination, in which the annual increase in rate of obesity ranged from 0.10%−0.58% among students aged 7–18 years from 1985 to 2014.[8]

Overall, the association between childhood and adolescent obesity and food industry level and healthcare service were statistically significant. We should create supportive environments for healthier-food choice via continuous systematic research on how the imbalance is formed and how it can be reversed, and by education and the establishment of regulations on food labeling.

Although the data were limited, we adopted PLSR and Bayesian analysis, which was suitable for modeling with small sized sample, and still fitted a stable space-time predictive model. Moreover, as the forecasts for 2015, 2020, and 2030 were based on 2011 environmental data, such predictions may be conservative or underestimated in the context of the overall social developments. In the future, more detailed data should be obtained through routine monitoring to verify the results of this study. Furthermore, “individual-community-province” hierarchical or random-effects models will be adopted in our next research, which can deal with the problem of temporal or spatiotemporal misalignment better and to explore the influencing mechanism and changing rules in obesogenic environments more precisely and clearly. The effects and mechanism of school policy, children dependency ratio, and parents’ educational situation on children obesity also need further explorations.

In summary, there have been obvious changes in obesogenic environmental priorities during the different rounds. Obesity prevalence in children and adolescents aged 7–17 years in China experienced a rapid increase, at about 0.3584% annually between 2000 and 2011. From longitudinal observation, the associations between food industry level (“Amplifier”), healthcare service (“Balancer”) and children and adolescents’ obesity prevalence were statistically significant. Targeted Interventions in northern and eastern provinces of China are necessary and urgent. Further researches on the mechanisms of how levels of food industry, healthcare service, and other environmental factors influence prevalence of obesity are needed. We need to understand the major contradictions timeously and identify modifiable determinants to control and reverse the epidemic of obesity and to achieve enduring health outcomes. Given the complexity of problems, concerted efforts are needed to promote innovation. There should be supply-side reforms in the food industry, healthcare, and educational systems (primary schools, middle schools, medical colleges, and food industry colleges) to support the new generations on how to understand food, nutrition, and activities; and how to produce, select, and consume healthier food.