Key PointsQuestion
Is the general vegetation level of a residential area, referred to as residential greenness, associated with cardiovascular disease among adults, and does the presence of cardiometabolic disorders mediate or modify the association between residential greenness and cardiovascular disease?
Findings
In this cross-sectional study of 24 845 adults in China, residential areas with higher greenness levels were associated with a lower likelihood of cardiovascular disease. The presence of cardiometabolic disorders partially mediated the association between residential greenness and cardiovascular disease.
Meaning
The study’s findings may be helpful for health care professionals and policy makers in the development of strategies, such as planning for green spaces in residential areas, to mitigate the burden of cardiovascular disease.
Importance
Living in areas with more vegetation (referred to as residential greenness) may be associated with cardiovascular disease (CVD), but little data are available from low- and middle-income countries. In addition, it remains unclear whether the presence of cardiometabolic disorders modifies or mediates the association between residential greenness and CVD.
Objective
To evaluate the associations between residential greenness, cardiometabolic disorders, and CVD prevalence among adults in China.
Design, Setting, and Participants
This analysis was performed as part of the 33 Communities Chinese Health Study, a large population-based cross-sectional study that was conducted in 33 communities (ranging from 0.25-0.64 km2) in 3 cities within the Liaoning province of northeastern China between April 1 and December 31, 2009. Participants included adults aged 18 to 74 years who had resided in the study area for 5 years or more. Greenness levels surrounding each participant’s residential community were assessed using the normalized difference vegetation index and the soil-adjusted vegetation index from 2010. Lifetime CVD status (including myocardial infarction, heart failure, coronary heart disease, cerebral thrombosis, cerebral hemorrhage, cerebral embolism, and subarachnoid hemorrhage) was defined as a self-report of a physician diagnosis of CVD at the time of the survey. Cardiometabolic disorders, including hypertension, diabetes, dyslipidemia, and overweight or obese status, were measured and defined clinically. Generalized linear mixed models were used to evaluate the association between residential greenness levels and CVD prevalence. A 3-way decomposition method was used to explore whether the presence of cardiometabolic disorders mediated or modified the association between residential greenness and CVD. Data were analyzed from October 10 to May 30, 2020.
Main Outcomes and Measures
Lifetime CVD status, the presence of cardiometabolic disorders, and residential greenness level.
Results
Among 24 845 participants, the mean (SD) age was 45.6 (13.3) years, and 12 661 participants (51.0%) were men. A total of 1006 participants (4.1%) reported having a diagnosis of CVD. An interquartile range (1-IQR) increase in the normalized difference vegetation index within 500 m of a community was associated with a 27% lower likelihood (odds ratio [OR], 0.73; 95% CI, 0.65-0.83; P < .001) of CVD prevalence, and an IQR increase in the soil-adjusted vegetation index within 500 m of a community was associated with a 26% lower likelihood (OR, 0.74; 95% CI, 0.66-0.84; P < .001) of CVD prevalence. The presence of cardiometabolic disorders was found to mediate the association between residential greenness and CVD, with mediation effects of 4.5% for hypertension, 4.1% for type 2 diabetes, 3.1% for overweight or obese status, 12.7% for hypercholesterolemia, 8.7% for hypertriglyceridemia, and 11.1% for high low-density lipoprotein cholesterol levels.
Conclusions and Relevance
In this cross-sectional study, higher residential greenness levels were associated with lower CVD prevalence, and this association may be partially mediated by the presence of cardiometabolic disorders. Further studies, preferably longitudinal, are warranted to confirm these findings.
Cardiovascular disease (CVD) remains the primary cause of death and disability worldwide, especially in low- and middle-income countries.1 Between 1990 and 2016, the estimated number of individuals with CVD in China increased from 40.6 million to 93.8 million.2 Interventions and control strategies are therefore needed to mitigate this increase in CVD prevalence.
Accumulating data indicate that a higher amount of neighborhood vegetation (referred to as residential greenness) is associated with a range of beneficial health outcomes.3,4 A large number of studies have examined the cardiovascular associations of greenness,5-37 but more than 50% of those focused on CVD mortality. In addition, most of the previous studies were conducted in middle- and high-income countries, and only a few studies were performed in low-income countries.23,27,30,31 It is unclear if the results of studies from high-income countries are generalizable to lower-income settings because a number of cultural, sociodemographic, and environmental factors differ between high-income nations vs lower-income nations.
Cardiometabolic disorders, including hypertension, type 2 diabetes, overweight or obese status, and dyslipidemia, are documented preclinical factors associated with CVD risk.38 In environmental epidemiologic studies, these cardiometabolic disorders have often been suggested as intervening variables in the association between environmental exposures and CVD. A systematic review and meta-analysis3 as well as previous findings from the 33 Communities Chinese Health Study39-41 have indicated that higher residential greenness levels were associated with a lower risk of cardiometabolic disorders. Two studies within the past 5 years have reported that the presence of cardiometabolic disorders partially mediated the association between greenness and CVD prevalence.28,36 However, it is also plausible that the presence of cardiometabolic disorders may modify the health outcomes associated with residential greenness because individuals with these disorders may change their lifestyles (eg, engage in physical activity or visit parks more frequently) as part of a treatment plan. Therefore, the associations between residential greenness, cardiometabolic disorders, and CVD remain unclear. To address these research gaps, we aimed to (1) examine the association between residential greenness levels and CVD prevalence in a general population of Chinese adults and (2) explore whether the presence of cardiometabolic disorders mediated or modified the association between residential greenness levels and CVD using an advanced 3-way decomposition method.
Participants and Procedures
The 33 Communities Chinese Health Study, a large population-based cross-sectional investigation, was conducted in the Liaoning province of northeastern China between April 1 and December 31, 2009. The province is a large industrial center that is highly urbanized. The prevalence rates of CVD and cardiometabolic disorders have been reported to be high in this region.42-44 This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline for cross-sectional studies. The Human Studies Committee of Sun Yat-sen University approved the study protocols. All participants provided written informed consent.
We used a 4-stage stratified cluster sampling strategy to randomly recruit study participants39-41 (Figure). First, 3 cities (Shenyang, Anshan, and Jinzhou) were randomly selected among 14 total cities in the Liaoning province. Second, 3 residential communities were randomly selected from each of the 11 total districts (5 districts in Shenyang and 3 districts each in Anshan and Jinzhou), generating 33 communities that ranged in size from 0.25 km2 to 0.64 km2. Third, 700 to 1000 households were randomly selected from each community. Fourth, 1 adult was randomly selected from each household. A total of 28 830 individuals were initially invited to participate in the study. Of those, 3985 individuals (13.8%) were excluded for the following reasons: (1) 3634 individuals did not complete the study questionnaire; (2) 233 individuals had resided in the study community for less than 5 years; (3) 79 individuals were pregnant; (4) 27 individuals were younger than 18 years or older than 74 years; and (5) 12 individuals had a serious preexisting disease, such as terminal cancer. After exclusions, 24 845 participants (86.2%) were included in the analysis, and 15 477 participants (53.7%) provided venous blood samples after fasting. Sociodemographic characteristics were similar between the participants who provided blood samples and the participants who did not (n = 9368) (eTable 2 in the Supplement).
Cardiovascular Disease and Cardiometabolic Disorders
We assessed CVD status using a study questionnaire. A CVD case was defined as an affirmative response to the question: “Has a physician ever diagnosed you with myocardial infarction, heart failure, coronary heart disease, cerebral thrombosis, cerebral hemorrhage, cerebral embolism, or subarachnoid hemorrhage?”45
Using a standard mercuric-column sphygmomanometer, trained and certified nurses measured systolic blood pressure and diastolic blood pressure according to procedures recommended by the American Heart Association.46 We defined hypertension as mean systolic blood pressure higher than 140 mm Hg, mean diastolic blood pressure higher than 90 mm Hg, and/or reported receipt of antihypertensive medication in the 2 weeks before the interview.47
Using protocols recommended by the World Health Organization,48 we measured height and weight and calculated body mass index (BMI; calculated as weight in kilograms divided by height in meters squared) for each participant. Based on Asian criteria,48 obese status was defined as a BMI of 27.5 or higher, overweight status as a BMI of 23.0 to 27.5, and normal weight (including underweight) status as a BMI of less than 23.0.
After an overnight fast, peripheral venous blood samples were obtained. The levels of total cholesterol, triglycerides, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol were determined using a Hitachi Autoanalyzer, type 7170A (Hitachi). We defined hypercholesterolemia as a total cholesterol level of 240 mg/dL or higher, hypertriglyceridemia as a triglyceride level of 200 mg/dL or higher, low HDL cholesterol as an HDL cholesterol level of 40 mg/dL or higher, and high LDL cholesterol as an LDL cholesterol level of 160 mg/dL or higher.49
We performed a standard 75-g oral glucose tolerance test and obtained blood samples at 0 hours and 2 hours after glucose intake. We measured fasting and 2-hour glucose levels using an enzymatic colorimetric method. Type 2 diabetes was defined based on American Diabetes Association guidelines as a fasting blood glucose level of 126 mg/dL or higher (to convert to mmol/L, multiply by 0.0555), a 2-hour glucose level of 200 mg/dL or higher, or the reported receipt of antidiabetic medication.50
Residential Greenness Levels
We applied 2 satellite-based vegetation indices, the normalized difference vegetation index (NDVI)51 and the soil-adjusted vegetation index (SAVI),52 to characterize residential greenness levels. Both indices are derived based on the difference of surface reflectance and absorbance in 2 vegetation-informative light bands, visible red and near-infrared. The NDVI is calculated by subtracting visible red light from near-infrared light and dividing the difference by the sum of near-infrared light and visible red light. For SAVI, a correction factor was added to minimize factors associated with the soil background. Therefore, the SAVI was calculated by subtracting visible red light from near-infrared light and dividing the difference by the sum of near-infrared light, visible red light, and the correction factor; the dividend was then multiplied by the sum of the correction factor and 1.0. The resulting NDVI and SAVI values were unitless and ranged from –1.0 to 1.0, with higher values indicating greener areas, negative values indicating bodies of water, and values close to 0 indicating barren areas.
The 2 indices were calculated using Landsat 5 Thematic Mapper (National Aeronautics and Space Administration) satellite images at a spatial resolution of 30 m by 30 m, which were obtained during the greenest month (ie, August, which has the highest level of growing vegetation) and the year closest to the collection of health data (ie, 2010). Residential greenness was defined as the mean NDVI and SAVI levels in circular buffers of 500 m and 1000 m around each community’s centroid. Greenness level was defined using the NDVI, with low greenness level defined as less than 0.40 m of vegetation per 500 m, and high greenness level defined as 0.40 m or more of vegetation per 500 m. All calculations were performed using ArcGIS, version 10.4 (Esri).
Confounding and Mediating Factors
We developed a directed acyclic graph to select confounding variables and potential mediating factors using the DAGitty package, version 2.16.3 for R software (R Foundation for Statistical Computing) (eFigure 1 in the Supplement). The following confounding variables were retained: age, sex, ethnicity, annual mean household income, educational level, physical activity level, district-level gross domestic product,39 and particulate matter with an aerodynamic diameter of 2.5 μm or less.53
Cardiometabolic disorders were selected as potential mediating factors (eFigure 1 in the Supplement). We developed an additional directed acyclic graph to identify potential confounding variables of exposure-mediator associations (eFigure 2 in the Supplement) and mediator-outcome associations (eFigure 3 in the Supplement). We incorporated the common confounding variables (ie, age, sex, ethnicity, household income, educational level, physical activity level, district-level gross domestic product, and particulate matter with an aerodynamic diameter of 2.5 μm or less) into the mediation analysis.
We used generalized linear mixed models54 to assess the association between residential greenness levels and CVD prevalence, with communities considered a random effect. The effect estimates were expressed as odds ratios (ORs) and 95% CIs per interquartile range (IQR) increase in NDVI and SAVI levels. We estimated unadjusted models and main models that were adjusted for the confounding variables selected using the directed acyclic graph (including age, sex, ethnicity, household income, educational level, physical activity level, district-level gross domestic product, and particulate matter with an aerodynamic diameter of 2.5 μm or less) (eFigure 1 in the Supplement). Informed by findings from previous studies,55,56 we focused on a 500-m exposure buffer for the main analysis.
We also performed several sensitivity analyses. First, we estimated NDVI and SAVI levels in a larger exposure buffer of 1000 m to assess the impact of exposure misclassification. Second, we categorized NDVI and SAVI levels per 500 m into quartiles (for NDVI levels, quartile 1 comprised 0.18-0.23, quartile 2 comprised 0.24-0.29, quartile 3 comprised 0.30-0.40, and quartile 4 comprised 0.41-0.80; for SAVI levels, quartile 1 comprised 0.10-0.13, quartile 2 comprised 0.14-0.16, quartile 3 comprised 0.17-0.24, and quartile 4 comprised 0.25-0.48), and we used a natural spline-smoothing function to test for nonlinear associations. Third, we built models using sequential adjustments, and we adjusted the main models for alcohol consumption, cigarette smoking, physical activity level, controlled diet with low calories and low fat, consumption of sugar-sweetened soft drinks, and family history of CVD to assess their potential impact. We performed stratified analyses by age (<50 years vs ≥50 years), sex, annual household income (<10 000 yuan vs ≥10 000 yuan [$1 was equivalent to 6.84 yuan in 2009]), and educational level (<9 years vs ≥9 years) to examine whether the associations were consistent among different subpopulations.
We used a 3-way decomposition method to apportion the association between residential greenness and each cardiometabolic disorder with which the greenness level may have interacted with CVD prevalence. The 3 nonoverlapping components included (1) variables associated with greenness that were directly controlled for potential mediating factors (direct effect), (2) variables associated with additive interaction but not mediation (interactive effect), and (3) variables associated with mediation but not interaction (pure indirect effect) (eFigure 4 in the Supplement).57 Standard errors were calculated using the delta method.58 We also estimated the mediation effect of combined cardiometabolic disorders using the lavaan package, version 3.4.3 for R software (R Foundation for Statistical Computing). All other statistical analyses were performed using SAS software, version 9.4 (SAS Institute Inc). All tests were 2-tailed and paired, with a significance threshold of P < .05.
Among 24 845 adult participants, the mean (SD) age was 45.6 (13.3) years, and 12 661 participants (51.0%) were men (Table 1). Most of the participants (94.5%) were of Han ethnicity and had a high educational level (78.0% of participants had >9 years of education) and high household income (76.8% of participants had a mean household income of >10 000 yuan per year). A total of 1006 participants (4.1%) reported having a diagnosis of CVD. Participants with CVD were more likely to be male (743 participants [73.9%]), 50 years or older (749 participants [75.5%]), of Han ethnicity (980 participants [97.4%]), have an educational level of 9 years or less (836 participants [93.1%]), and have an annual household income of more than 10 000 yuan (640 participants [63.6%]) (eTable 3 in the Supplement). Participants living in communities with higher greenness levels were more likely to be female (3106 participants [59.4%]), younger than 50 years (3662 participants [70.0%]), have an educational level of 9 years or less (3932 participants [75.2%]), have an annual household income of more than 10 000 yuan (4193 participants [80.2%]), exercise less than 180 minutes per week (3547 participants [67.8%]), and live in communities with lower concentrations of particulate matter with an aerodynamic diameter of 2.5 μm or less (71 participants [30.5%]). The prevalence of cardiometabolic disorders ranged from 1333 participants (8.6%) with high LDL cholesterol to 15 459 participants (62.2%) with overweight or obese status.
Greenness levels varied across the communities. For instance, NDVI levels per 500 m ranged from 0.18 to 0.80, with a median of 0.29 (IQR, 0.23-0.40), and SAVI levels per 500 m ranged from 0.10 to 0.48, with a median of 0.16 (IQR, 0.13-0.24) (eTable 4 in the Supplement). In addition, NDVI levels were consistent with SAVI levels. For example, the lowest Spearman correlation coefficient was 0.88 for the correlation between NDVI levels per 500 m and SAVI levels per 1000 m, and the highest Spearman correlation coefficient was 0.98 for the correlation between NDVI levels per 500 m and SAVI levels per 500 m (eTable 4 in the Supplement).
In the adjusted models, an IQR increase in both NDVI and SAVI levels per 500 m was associated with a lower likelihood of CVD prevalence, with ORs of 0.73 (95% CI, 0.65-0.83; P < .001) for NDVI levels and 0.74 (95% CI, 0.66-0.84; P < .001) for SAVI levels (Table 2). Similar results were observed for NDVI levels per 1000 m (OR, 0.79; 95% CI, 0.71-0.89; P < .001) and SAVI levels per 1000 m (OR, 0.78; 95% CI, 0.69-0.87; P < .001). The effect estimates were similar to those of the main models when we adjusted for alcohol use (for NDVI: OR, 0.74; 95% CI, 0.65-0.83; P < .001; for SAVI: OR, 0.74; 95% CI, 0.66-0.84; P < .001), cigarette smoking (for NDVI: OR, 0.73; 95% CI, 0.64-0.86; P < .001; for SAVI: OR, 0.73; 95% CI, 0.65-0.83; P < .001), low-calorie and low-fat diet (for NDVI: OR, 0.73; 95% CI, 0.65-0.82; P < .001; for SAVI: OR, 0.74; 95% CI, 0.65-0.83; P < .001), consumption of sugar-sweetened soft drinks (for NDVI: OR, 0.73; 95% CI, 0.64-0.82; P < .001; for SAVI: OR, 0.73; 95% CI, 0.65-0.83; P < .001), and family history of CVD (for NDVI: OR, 0.73; 95% CI, 0.65-0.82; P < .001; for SAVI: OR, 0.74; 95% CI, 0.65-0.83; P < .001) and when we incorporated sequential adjustments (eg, in model 1, which was adjusted for age, sex, ethnicity, and community size, the OR for NDVI was 0.71; 95% CI, 0.63-0.80; P < .001;for SAVI, 0.72; 95% CI, 0.64-0.81; P < .001) (eTable 5 in the Supplement).
We found similar linear dose-response patterns when using NDVI and SAVI levels per 500 m categorized into quartiles (eg, for quartile 2, the adjusted OR for NDVI was 1.03; 95% CI, 0.86-1.23; for quartile 3, the adjusted OR for NDVI was 0.66; 95% CI, 0.55-0.80; for quartile 4, the adjusted OR for NDVI was 0.63; 95% CI, 0.50-0.79; with P < .001 for trend) (eTable 6 in the Supplement) and when using a restrictive cubic spline analysis (for NDVI, P for nonlinear association = 0.74; for SAVI, P for nonlinear association = 0.62) (eFigure 5 in the Supplement). We found no significant interaction by age, sex, household income, or educational level (eTable 7 in the Supplement).
The mediation analyses examining the association between NDVI level per 500 m and CVD prevalence indicated mediation effects of 4.5% (95% CI, 0.6%-8.3%; P = .02) for hypertension, 3.1% (95% CI, 1.1%-5.1%; P = .002) for overweight or obese status, 4.1% (95% CI, 0.7%-7.4%; P = .02) for type 2 diabetes, 12.7% (95% CI, 4.6%-20.8%; P = .002) for hypercholesterolemia, 8.7% (95% CI, 2.9%-14.4%; P = .003) for hypertriglyceridemia, and 11.1% (95% CI, 3.8%-18.4%; P = .003) for high LDL cholesterol. No significant mediation effect was found for low HDL cholesterol (1.9%; 95% CI, −0.3% to 4.0%; P = .09) (Table 3). Combined, these cardiometabolic disorders jointly mediated 21.2% (95% CI, 6.4%-35.9%; P = .004) of the association between NDVI level per 500 m and CVD prevalence (eTable 8 in the Supplement). Similar results were obtained when we used SAVI level per 500 m as the exposure measure, with a significant mediation effect for all cardiometabolic disorders (eg, for hypercholesterolemia, 13.0%; 95% CI, 5.0%-21.0%; P = .001) with the exception of low HDL cholesterol (2.0%; 95% CI, −0.2% to 4.2%; P = .08). No significant interactive effect was found between residential greenness and any cardiometabolic disorder (Table 3).
In this study, we found that participants living in greener communities had a lower likelihood of receiving a diagnosis of CVD. The association between residential greenness and CVD was robust in a series of sensitivity analyses. The presence of cardiometabolic disorders (with the exception of low HDL cholesterol) partially mediated the association between residential greenness levels and CVD prevalence.
Many previous studies have documented an association between greenness levels and CVD (eTable 1 in the Supplement),5-37 which is consistent with our current findings. For example, a cross-sectional study of 11 404 Australian adults reported that higher neighborhood greenness levels (measured as the mean NDVI level per 1600 m) and variability (measured as the SD of the NDVI level per 1600 m) were associated with a lower prevalence of heart disease or stroke.13 Another cross-sectional study of Chinese adults reported a substantial reduction in the likelihood of coronary heart disease and stroke among populations residing in areas with higher NDVI levels.30 Similar results were reported by cross-sectional studies from the US,36 Israel,28 New Zealand,14 and Brazil.23
Although our study was cross-sectional and used NDVI levels to estimate exposure, our findings may support those from 2 prospective cohort studies15,32 that used different green space measures and 1 natural experiment study.18 The first cohort study followed up 5112 Lithuanian adults for 4.4 years. The investigators found that a greater distance from green spaces was associated with an increased risk of CVD.15 The second cohort study found that a larger tree canopy percentage was associated with a lower risk of incident CVD among 46 786 Australian adults.32 In the natural experiment study, Donovan et al21 analyzed longitudinal data from 156 146 women and examined the association between the loss of trees owing to an invasive forest pest and incident CVD. They reported that women living in areas with greater tree loss had an increased risk of CVD.21
We observed that all of the explored cardiometabolic disorders (with the exception of low HDL cholesterol) substantially mediated the association between residential greenness and CVD prevalence, indicating that the disorders might be involved in the pathways by which greenness is associated with CVD. Our findings were consistent with those of a US study,36 which found that the associations between NDVI levels and CVD were attenuated when regression models were adjusted for the presence of hypertension, type 2 diabetes, and dyslipidemia. Partially consistent with our findings, another study reported that the association of myocardial infarction with NDVI levels was partially mediated by the presence of dyslipidemia and type 2 diabetes but not hypertension.28 Combined with our findings, the data extended our understanding of potential biological mechanisms underlying the association between residential greenness and CVD prevalence. Because cardiometabolic disorders are reported to be precursors to clinical CVD, the results of our mediation analysis suggest that future research and health care decisions that focus on the association of greenness with these earlier disorders may be useful for reducing future CVD events.
An earlier study suggested several pathways, such as reduction in air pollution and heat, encouragement of physical activity, increases in social cohesion, and decreases in mental fatigue and stress,59 to explain the association between greenness and health. However, accumulating data do not provide support for an association between substantial air pollution reductions or increases in physical activity and a greater number of urban greenness areas.60,61 Because data on social cohesion and mental health were not available for our study, we did not investigate those potential factors. Future studies are needed to explore potential mechanisms, preferably through the inclusion of more mediators and the performance of serial mediation tests.
To our knowledge, this is one of a few epidemiologic studies to report an association between residential greenness and CVD in a low- or middle-income country. Our results may have public health importance, especially considering the fact that many low-income countries are experiencing rapid urbanization and subsequent decreases in greenness. A major strength of our study is its large population-based sample, which was randomly selected from communities across 3 different cities, and its high response rate, which minimized selection bias and provided sufficient statistical power to detect modest effects. In addition, we used 2 greenness exposure measures and 2 different exposure buffers, and we adjusted for a parsimonious yet comprehensive set of covariates to control confounding without overadjustment. Furthermore, we performed mediation analyses using a 3-way decomposition method, which incorporates interactive effects and typically provides greater validity and higher statistical power than traditional approaches.57
This study has several limitations. First, the cross-sectional design did not allow us to assess temporality between greenness exposure and CVD, although the likelihood of reverse-association is low. In addition, the cross-sectional design may have overestimated the mediation effects.62 Second, we assessed greenness exposure based on communities and not individuals, which might have produced measurement error (ie, Berkson error63). Although this error did not bias our effect estimates, it could have reduced statistical power. In addition, we have not specifically collected data on limited mobility. Third, although we objectively measured greenness exposure using NDVI and SAVI levels, these 2 measures reflect general vegetation levels and cannot differentiate the structure, type, and quality of greenness. Thus, we were unable to examine aspects of greenness that are most relevant to CVD prevalence.
Fourth, data on antidyslipidemia medication receipt were not available and were not used to define dyslipidemia, which might have reduced the prevalence of dyslipidemia and biased the results of our mediation analysis. Fifth, CVD diagnosis and several potential confounding variables were self-reported; therefore, false-positive findings and recall bias are possible. In addition, we assessed socioeconomic status using a limited number of variables, and other environmental factors, such as noise, walkability, and air temperature, were not considered in our model; thus, our results might have been affected by residual confounding.
Higher community greenness levels were associated with a lower likelihood of CVD prevalence in an industrial and highly urbanized setting within China. The association might be partly mediated by the presence of cardiometabolic disorders. Although the findings require confirmation in longitudinal studies, they contribute to the limited data available to policy makers who are interested in designing population-level interventions to mitigate the prevalence of CVD in China and other low- and middle-income countries.
Accepted for Publication: June 16, 2020.
Published: September 21, 2020. doi:10.1001/jamanetworkopen.2020.17507
Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2020 Yang B-Y et al. JAMA Network Open.
Corresponding Authors: Guang-Hui Dong, MD, PhD, Guangdong Provincial Engineering Technology Research Center of Environmental and Health Risk Assessment, Department of Occupational and Environmental Health, Sun Yat-sen University School of Public Health, 74 Zhongshan Second Rd, Yuexiu District, Guangzhou 510080, China (donggh5@mail.sysu.edu.cn); Yunjiang Yu, PhD, State Environmental Protection Key Laboratory of Environmental Pollution Health Risk Assessment, South China Institute of Environmental Sciences, Ministry of Environmental Protection, Guangzhou 510655, China (yuyunjiang@scies.org).
Author Contributions: Dr Dong had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Drs Yang and Hu contributed equally to this work.
Concept and design: Yang, Hu, Gao, Guo, Y. Yu, Dong.
Acquisition, analysis, or interpretation of data: Yang, Jalaludin, Knibbs, Markevych, Heinrich, Bloom, Morawska, Lin, Jalava, Roponen, Gao, Chen, Zhou, H.-Y. Yu, Liu, Zeng, Zeeshan, Guo, Dong.
Drafting of the manuscript: Yang, Gao, Y. Yu, Dong.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Yang, Gao, H.-Y. Yu, Liu, Guo, Dong.
Obtained funding: Yang, Hu, Zeng, Guo, Dong.
Administrative, technical, or material support: Yang, Hu, Knibbs, Bloom, Lin, Chen, Zeeshan, Y. Yu, Dong.
Supervision: Heinrich, Guo, Y. Yu, Dong.
Conflict of Interest Disclosures: Dr Bloom reported receiving personal fees and nonfinancial support from Sun Yat-sen University during the conduct of the study. No other disclosures were reported.
Funding/Support: This study was funded by grants 81673128, 81703179, 81803196, 81872582, and 91543208 from the National Natural Science Foundation of China; grants 2016A030313342 and 2017A050501062 from the Guangdong Province Natural Science Foundation; grant 2016YFC0207000 from the National Key Research and Development Program of China; and grants 201803010054 and 201807010032 from the Science and Technology Program of Guangzhou.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Additional Contributions: The authors acknowledge the cooperation of the participants in this study, who have been very generous with their time and assistance.
8.Richardson
E, Pearce
J, Mitchell
R, Day
P, Kingham
S. The association between green space and cause-specific mortality in urban New Zealand: an ecological analysis of green space utility.
BMC Public Health. 2010;10:240. doi:
10.1186/1471-2458-10-240
PubMedGoogle ScholarCrossref 9.Dominguez-Berjon
MF, Gandarillas
A, Segura del Pozo
J,
et al. Census tract socioeconomic and physical environment and cardiovascular mortality in the region of Madrid (Spain).
J Epidemiol Community Health. 2010;64(12):1086-1093. doi:
10.1136/jech.2008.085621
PubMedGoogle ScholarCrossref 23.Massa
KHC, Pabayo
R, Lebrao
ML, Chiavegatto Filho
ADP. Environmental factors and cardiovascular diseases: the association of income inequality and green spaces in elderly residents of Sao Paulo, Brazil.
BMJ Open. 2016;6(9):e011850. doi:
10.1136/bmjopen-2016-011850
PubMedGoogle Scholar 24.Ngom
R, Gosselin
P, Blais
C, Rochette
L. Type and proximity of green spaces are important for preventing cardiovascular morbidity and diabetes—a cross-sectional study for Quebec, Canada.
Int J Environ Res Public Health. 2016;13(4):423. doi:
10.3390/ijerph13040423
PubMedGoogle ScholarCrossref 25.Picavet
HSJ, Milder
I, Kruize
H, de Vries
S, Hermans
T, Wendel-Vos
W. Greener living environment healthier people?: exploring green space, physical activity and health in the Doetinchem cohort study.
Prev Med. 2016;89:7-14. doi:
10.1016/j.ypmed.2016.04.021
PubMedGoogle ScholarCrossref 26.Vienneau
D, de Hoogh
K, Faeh
D, Kaufmann
M, Wunderli
JM, Roosli
M; SNC Study Group. More than clean air and tranquillity: residential green is independently associated with decreasing mortality.
Environ Int. 2017;108:176-184. doi:
10.1016/j.envint.2017.08.012
PubMedGoogle ScholarCrossref 30.Jia
X, Yu
Y, Xia
W,
et al. Cardiovascular diseases in middle aged and older adults in China: the joint effects and mediation of different types of physical exercise and neighborhood greenness and walkability.
Environ Res. 2018;167:175-183. doi:
10.1016/j.envres.2018.07.003
PubMedGoogle ScholarCrossref 31.da Silveira
IH, Junger
WL. Green spaces and mortality due to cardiovascular diseases in the city of Rio de Janeiro.
Rev Saude Publica. 2018;52:49.
PubMedGoogle ScholarCrossref 32.Astell-Burt
T, Feng
X. Urban green space, tree canopy and prevention of cardiometabolic diseases: a multilevel longitudinal study of 46 786 Australians.
Int J Epidemiol. 2020;49(3):926-933. doi:
10.1093/ije/dyz239
PubMedGoogle ScholarCrossref 33.Orioli
R, Antonucci
C, Scortichini
M,
et al. Exposure to residential greenness as a predictor of cause-specific mortality and stroke incidence in the Rome longitudinal study.
Environ Health Perspect. 2019;127(2):27002. doi:
10.1289/EHP2854
PubMedGoogle ScholarCrossref 34.Kim
S, Kim
H, Lee
J-T. Interactions between ambient air particles and greenness on cause-specific mortality in seven Korean metropolitan cities, 2008-2016.
Int J Environ Res Public Health. 2019;16(10):1866. doi:
10.3390/ijerph16101866
PubMedGoogle ScholarCrossref 38.Forouzanfar
MH, Alexander
L, Anderson
HR,
et al; GBD 2013 Risk Factors Collaborators. Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks in 188 countries, 1990-2013: a systematic analysis for the Global Burden of Disease study 2013.
Lancet. 2015;386(10010):2287-2323. doi:
10.1016/S0140-6736(15)00128-2
PubMedGoogle ScholarCrossref 44.Zhang
L, Wang
Z, Wang
X,
et al; China Hypertension Survey Investigators. Prevalence of abdominal obesity in China: results from a cross-sectional study of nearly half a million participants.
Obesity (Silver Spring). 2019;27(11):1898-1905. doi:
10.1002/oby.22620
PubMedGoogle ScholarCrossref 45.Qin
X-D, Qian
Z, Vaughn
MG,
et al. Gender-specific differences of interaction between obesity and air pollution on stroke and cardiovascular diseases in Chinese adults from a high pollution range area: a large population based cross sectional study.
Sci Total Environ. 2015;529:243-248. doi:
10.1016/j.scitotenv.2015.05.041
PubMedGoogle ScholarCrossref 46.Pickering
TG, Hall
JE, Appel
LJ,
et al. Recommendations for blood pressure measurement in humans and experimental animals: part 1: blood pressure measurement in humans: a statement for professionals from the Subcommittee of Professional and Public Education of the American Heart Association Council on High Blood Pressure Research.
Circulation. 2005;111(5):697-716. doi:
10.1161/01.CIR.0000154900.76284.F6
PubMedGoogle ScholarCrossref 47.Chobanian
AV, Bakris
GL, Black
HR,
et al; Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure; National Heart, Lung, and Blood Institute; National High Blood Pressure Education Program Coordinating Committee. Seventh report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure.
Hypertension. 2003;42(6):1206-1252. doi:
10.1161/01.HYP.0000107251.49515.c2
PubMedGoogle ScholarCrossref 49.Joint Committee for Developing Chinese Guidelines on Prevention and Treatment of Dyslipidemia in Adults. [Chinese guidelines on prevention and treatment of dyslipidemia in adults].
Zhonghua Xin Xue Guan Bing Za Zhi. 2007;35(5):390-419.
PubMedGoogle Scholar 55.Dadvand
P, Villanueva
CM, Font-Ribera
L,
et al. Risks and benefits of green spaces for children: a cross-sectional study of associations with sedentary behavior, obesity, asthma, and allergy.
Environ Health Perspect. 2014;122(12):1329-1335. doi:
10.1289/ehp.1308038
PubMedGoogle ScholarCrossref 63.Keogh
RH, Shaw
PA, Gustafson
P,
et al. STRATOS guidance document on measurement error and misclassification of variables in observational epidemiology: part 1—basic theory and simple methods of adjustment.
Stat Med. 2020;39(16):2197-2231. doi:
10.1002/sim.8532
PubMedGoogle ScholarCrossref