Abstract
Background:
Research has shown that higher altitude is associated with lower risk of lung cancer and improved survival among patients. The current study assessed the influence of county-level atmospheric pressure (a measure reflecting both altitude and temperature) on age-adjusted lung cancer mortality rates in the contiguous United States, with 2 forms of spatial regression.
Methods:
Ordinary least squares regression and geographically weighted regression models were used to evaluate the impact of climate and other selected variables on lung cancer mortality, based on 2974 counties.
Results:
Atmospheric pressure was significantly positively associated with lung cancer mortality, after controlling for sunlight, precipitation, PM2.5 (µg/m3), current smoker, and other selected variables. Positive county-level β coefficient estimates (P < .05) for atmospheric pressure were observed throughout the United States, higher in the eastern half of the country.
Conclusion:
The spatial regression models showed that atmospheric pressure is positively associated with age-adjusted lung cancer mortality rates, after controlling for other selected variables.
Introduction
In 2017, lung cancer accounted for an estimated 222 500 of new cases (14% of all cancer cases) and 155 870 of deaths (25% of all cancer deaths) in the United States. 1 Tobacco smoking is the leading cause of lung cancer, with men and women who smoke 25 times more likely to develop the disease than nonsmokers. 1 Smoking also explains about 80% of deaths from lung cancer, with the percentage even higher for small-cell lung cancer. 2 Other risk factors include exposure to second-hand smoke, asbestos, certain metals, radon gas, air pollution, diesel exhaust, and some organic chemicals. 1 In addition, research has found that lung cancer risk decreases with higher altitude. 3 -5 A decreased risk for other forms of cancer has also been associated with higher altitude, which equates to lower atmospheric oxygen concentration, barometric pressure, and natural background radiation. 6 -10 The association between higher altitude and lower lung cancer risk is stronger than that with other cancers (eg, breast, prostate, or bowels), indicating that the inhalation process plays an important role. 3 One study found that, on average, lung cancer decreases by 7.23 per 100 000 cases for every 1000-m rise in altitude, after adjusting for selected environmental correlates of elevation. 3
Inspired molecular oxygen (O2) is essential for life. About 21% of the atmosphere consists of O2. Oxygen levels are the same at both high and low altitudes. However, at high altitudes air pressure is much lower, allowing air particles to spread farther apart and making oxygen less accessible. Oxygen tends to easily react with other molecules in a cell, forming a type of unstable molecule called reactive oxygen species (ROS). Under normal conditions, cells control ROS levels. However, under oxidative stress conditions, ROS buildup in cells can contribute to carcinogenesis due to the damage of lipids, DNA, RNA, and cellular protein, as well as inhibit apoptosis and activate proto-oncogenes. 11 -16
The purpose of the current study was to assess the influence of atmospheric pressure, along with other selected variables (ie, sunlight, precipitation, sex, ethnicity, race, poverty, population density, smoking, obesity, diabetes, and ambient air pollution), on lung cancer mortality using 2 forms of spatial regression. Because ambient air temperature can influence known risk factors (eg, radon gas and ground level ozone) for lung cancer, our analysis used atmospheric pressure (a composite indicator of altitude and temperature). This is the first study to assess whether atmospheric pressure based on population-weighted, county-level altitude and temperature data is associated with lung cancer mortality.
Methods
Cancer Mortality Data
County-level lung cancer mortality rates (age-adjusted using the 2000 US standard population) were obtained from the CDC WONDER database, 1999 to 2010. 17 These rates were based on death certificate data from the National Center for Health Statistics (NCHS). 18 Annual detailed mortality files from the NCHS are coded according to the 10th revision of the International Statistical Classification of Diseases and Related Health Problems (ICD-10). Counties reporting less than 20 deaths from individual cancers were excluded because of confidentiality and large standard errors. This left 2974 counties for assessment.
Health and Socioeconomic Factors
County-level socioeconomic variables for 1999 to 2010 were included: percentage male, caucasian, and Hispanic from the CDC WONDER database 17 ; percentage below the poverty line from the US Census Bureau 19 ; percentage of current smoker from Dwyer-Lindgren et al 20 ; percentage obese and percentage diabetic (2004-2005) from the Behavior Risk Factor Surveillance System 21 ; and population density, calculated using data from the US Census Bureau and CDC WONDER database. 17
Climate Factors
Climate factors were atmospheric pressure, average daily sunlight, precipitation, and PM2.5. Weighted county-level altitude was calculated to account for population center residence. This was achieved by using a 30-m elevation data set available through ArcGIS Online from the US Geological Survey. 22 Census tract- and county-level shape files were obtained from the US Census Bureau. 23 Altitude values were then aggregated to the census tract-level, which approximates neighborhoods. Using census tract population estimates from 2000 and census tract altitude values, we calculated weighted county-level altitude values. Considering the nonlinear relationship between altitude and atmospheric pressure, and the impact of temperature on atmospheric pressure, we derived county atmospheric pressure using the hypsometric atmospheric pressure formula, weighted county elevation, and average daily maximum temperature. 24 Average daily maximum temperature was used because there is greater county-to-county variation in maximum than minimum temperature. Average daily sunlight, temperature, and precipitation (1999-2010), and PM2.5 (2003-2006) were obtained from the CDC WONDER database. 17
Model Selection
Both ordinary least squares (OLS) regression and geographically weighted regression (GWR) models were computed to evaluate the impact of climate and other selected variables. The GWR is a form of spatial regression that uses a moving bandwidth to calculate regional regression models and is used to determine whether the regression coefficients are consistent across space. These models were estimated using the software packages GeoDa (v 1.12) 25 and GWR v.4.09. 26 In ArcGIS Pro (v 2.0.1), Moran I identified significant positive spatial autocorrelation in the OLS regression model residuals. There was low multicollinearity for atmospheric pressure coefficients (variation inflation factors ≤2.59).
The presence of spatial autocorrelation was identified (P < .0001), meaning there was a clustering of model residuals. Hence, due to spatial influences, the observations were not considered independent. We used 2 forms of spatial regression to account for this spatial dependence. First, we added a spatial lag to the OLS regression. Spatial weights were calculated using a second-order queen contiguity matrix, and the spatial lag was added using GeoDa. Second, GWR was used to evaluate variation in the association between lung cancer mortality rates, climate, and other selected variables across regions. The use of the spatial lag and GWR produced a better-fitting, model-based lower Akaike information criteria (AIC, a measure of model fit) and larger adjusted R 2 values for the spatial regression models compared to the OLS regression model. The GWR bandwidth (the number of points included in the moving bandwidth) was calculated using a technique called Golden Selection Search, which chooses the optimal bandwidth size for each model according to AIC values. 25
Results
Atmospheric pressure values are presented by county for the contiguous United States (Figure 1). The lowest values are in the Rocky Mountain and plateau region of the United States. The highest values are along the Atlantic and Pacific coastal regions and the Gulf of Mexico. On the state level, mean atmospheric pressure is 960.1 hPa (standard deviation = 54.3), ranging from 809.6 to 1010.4. Eight (16%) states had scores below 900 hPa (Colorado, Utah, Wyoming, New Mexico, Nevada, Idaho, Montana, and Arizona). Seven states (Florida, Delaware, Louisiana, Rhode Island, New Jersey, Mississippi, and Connecticut) and the District of Columbia had scores above 1000 hPa. Based on our calculations, at low altitude and at constant temperature, a decrease in atmospheric pressure of 1 hPa is approximately equal to an 8.33-meter increase in elevation. Because of the nonlinear relationship between atmospheric pressure and altitude, there is a greater pressure change with a fixed altitude change at higher altitudes.

Atmospheric pressure and sunlight in the contiguous United States.
Estimates of the association between county-level lung cancer mortality rates and atmospheric pressure and other selected variables are shown in Table 1. Estimated coefficients in the table were simultaneously computed. A spatial lag variable was included in order to account for correlated error terms. The model statistics (adjusted R 2 and AIC) indicate that the variables included in the model do a good job explaining variation in lung cancer mortality rates. County-level, age-adjusted lung cancer mortality rates were significantly positively associated with current smoking (%), atmospheric pressure, and sunlight. Lung cancer mortality rates were significantly negatively associated with being Hispanic (%) or obese (%). In other words, smoking, higher atmospheric pressure, and more sunlight increased the risk of lung cancer, whereas being Hispanic or obese was protective against lung cancer.
Coefficient Estimates Using Ordinary Least Squares Regression.a
Abbreviations: AIC, Akaike information criteria; CI, confidence interval.
aCoefficient estimates were simultaneously computed for data from 2974 counties.
The GWR model provided a better fit to the data compared with the original OLS regression model, according to GWR analysis of variance that compares sum of squares (P < .0001). Additionally, the GWR model had a higher adjusted R 2 and lower AIC (Table 2) than the OLS model, after the addition of the spatial lag. The direction of association for each of the coefficient estimates is the same as seen in the previous table, except for precipitation. However, this variable is not statistically significant. Although diabetes (%) was marginally insignificant in the OLS model, the association with lung cancer mortality was greater in the better-fitting GWR model.
Geographically Weighted Regression Model Summary and Coefficient Estimates.a
Abbreviations: AIC, Akaike information criteria; GWR, geographically weighted regression; SD, standard deviation.
aCoefficient estimates were simultaneously computed for data from 2974 counties.
Figure 2 shows the GWR county-level β coefficient estimates (P < .05) of the association between lung cancer mortality rates and atmospheric pressure, adjusted for the other variables shown in Table 2. All significant β coefficients are positive, with the highest levels in the eastern half of the United States, particularly in Wisconsin, Iowa, Arkansas, Mississippi, and Louisiana.

Significant geographically weighted regression atmospheric pressure coefficients.
Discussion
The purpose of the current study was to assess the relationship between atmospheric pressure and age-adjusted lung cancer mortality rates, using 2 forms of spatial regression. Other selected variables were included in the study. Overall, atmospheric pressure had a significant positive relationship with lung cancer mortality, after adjusting for other selected variables. The positive association observed between atmospheric pressure and lung cancer is consistent with previous studies. 3 -5 County-level significant positive relationships between lung cancer mortality rates and atmospheric pressure were observed throughout much of the country and were stronger in the eastern half of the contiguous United States.
From the spatial lag model, a 100-hPa decrease (or an increase in altitude by approximately 833 meters) leads to a decrease of 3.58 lung cancer deaths per 100 000, after adjusting for sunlight, precipitation, sex, race, ethnicity, poverty, population density, smoking, obesity, diabetes, and PM2.5 (µg/m3). A previous study observed an average decrease in lung cancer incidence of 7.23 cases per 100 000 for every 1000-meter rise in altitude, after adjusting for selected variables associated with elevation. 3 The GWR results showed that the impact of atmospheric pressure on lung cancer mortality may be even greater, especially regionally. In the better-fitting GWR model, we found that the average coefficient for GWR is 0.1498 (median 0.1411), meaning a 100-hPa decrease in atmospheric pressure leads to a decrease of 14.98 lung cancer deaths per 100 000.
One of the primary benefits of GWR is that it allows us to evaluate whether there is any variation in the β coefficients of association across space. The fact that there were no significant negative coefficients for atmospheric pressure indicates consistency in the direction of the association and confirms the results of the spatial lag OLS model. However, the positive county-level association between atmospheric pressure and lung cancer mortality varies considerably by region. The spatial lag in the OLS model was statistically significant. This suggests that other spatially correlated variables (such as diet), which are related to both atmospheric pressure and lung cancer mortality, are unaccounted for in the model. In addition, unaccounted for variables can influence the regional associations between atmospheric pressure and lung cancer mortality.
Atmospheric pressure coefficients in those areas with the highest β coefficients can exceed 0.4. This means that in those areas an increase of 3 hPa (or about 25 meters in attitude) would lead to a decrease in 1.2 cancer deaths per 100 000. This is of particular importance because day-to-day variations in atmospheric pressure may exacerbate the physiological effects of lung cancer and, subsequently, lead to an increase in lung cancer mortality. Additional research is needed to examine the impact of seasonal and daily changes in atmospheric pressure on lung cancer mortality.
Higher altitude is correlated with greater Vitamin D synthesis from sunlight. Higher altitude is related to less atmospheric filters of ultraviolet radiation. Ultraviolet radiation levels increase from 10% to 12% for each 1000 m\ in altitude. 27 Greater natural radiation from cosmic rays at higher altitudes increases vitamin D levels produced by the skin. 28 Vitamin D synthesis from sunlight is also greater at lower latitudes because the zenith angle is smaller. Vitamin D is important because it modulates cell growth and stabilizes chromosomal structure, thus guarding against chromosomal aberrations and the formation of cancer. 29,30 Nevertheless, sunlight had a small positive effect on lung cancer mortality. The positive effect was unexpected, based on previous studies, 31 -33 possibly occurring because of an unaccounted confounding variable. It is important to note that observed positive correlation is between sunlight and lung cancer mortality and not vitamin D concentration. Sunlight may have an effect, such as with temperature, that is independent of sun-synthesized vitamin D. One study involving cardiovascular disease showed that sunlight has an impact on the body independent of serum vitamin D concentration. 34 Therefore, we may not fully understand the impact of sunlight on the body independent of vitamin D concentrations. Additional research is needed to fully examine the impact of sunlight on lung cancer mortality.
The hormetic model states that low doses of radiation reduce cancer risk, whereas high doses above a certain threshold increase cancer risk. 35 -37 Ecologic studies have shown an association between natural background radiation and lower cancer mortality, which supports the possibility of radiation hormesis. 6 -8,38 Although the study results are consistent with the possibility of radiation hormesis, the ecologic data limit drawing definitive conclusions.
The relatively large, positive association between smoking and lung cancer mortality is expected. Smoking has been an established cause of lung cancer since the 1964 Surgeon General’s Report. 39 Thus, it is important to note that with this variable in the model, atmospheric pressure continues to have a positive association with lung cancer mortality.
Obesity was negatively associated with lung cancer mortality, which is consistent with 2 previous meta-analyses that found that overweight and obesity were protective against lung cancer. 40,41 Studies have also shown that lung cancer-induced cachexia may lead to decreases in body weight. 42,43 However, another meta-analysis found that abdominal obesity was positively associated with the development of lung cancer. 44 Yet another study found no association between body mass index and lung cancer risk, regardless of smoking status. 45
Diabetes was independently associated with an increased risk of lung cancer mortality. This result is consistent with previous research. 46 Elevation and temperature both affect blood glucose, 47 so we included this variable in the model.
Hispanics had a lower rate of lung cancer mortality. The lower mortality is despite Hispanic patients tending to be diagnosed at later stage of disease and having higher levels of obesity and diabetes than non-Hispanic whites. 48 Although we adjust for obesity and diabetes in our models, we did not take into account the relatively low rates of lung cancer incidence among Hispanics. 49 The comparatively low incidence of lung cancer among Hispanics likely explains this result.
The present study has 2 primary limitations. First, it is an ecological study and, therefore, does not account for individual-level exposure to atmospheric pressure, sunlight, precipitation, and so on. Second, data are aggregated over a 12-year period, so the results do not capture daily and yearly changes in atmospheric pressure, sunlight, precipitation, pollution levels, and human migration.
Conclusion
This study shows that atmospheric pressure is positively associated with age-adjusted lung cancer mortality rates, after adjusting for other selected variables. This association appears to be greater for mortality than has been observed in previous studies involving incidence. Conclusions cannot be made from our results about the role of sunlight and sunlight-synthesized vitamin D, as they relate to lung cancer mortality rates. Air pollution was not significant in the models. This may be because of the large variation in air pollution within a county and because air pollution is associated with other variables in the models (ie, sunlight, temperature, and diabetes). Because of the influence of atmospheric pressure on lung cancer, individuals at high risk of lung cancer or patients with lung cancer may benefit from living at higher altitude.
Footnotes
Authors’ Note
All authors listed have made substantial, direct, and intellectual contribution to the work and approved it for publication.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
