Skip to main content

Outdoor air pollution, green space, and cancer incidence in Saxony: a semi-individual cohort study



There are a few epidemiological studies that (1) link increased ambient air pollution (AP) with an increase in lung cancer incidence rates and (2) investigate whether residing in green spaces could be protective against cancer. However, it is completely unclear whether other forms of cancer are also affected by AP and if residential green spaces could lower cancer incidence rates in general. Therefore, the objective was to estimate whether AP and green space are associated with several cancer types.


The analysis was based on routine health care data from around 1.9 million people from Saxony who were free of cancer in 2008 and 2009. Incident cancer cases (2010–2014) of mouth and throat, skin (non-melanoma skin cancer - NMSC), prostate, breast, and colorectum were defined as: (1) one inpatient diagnosis, or (2) two outpatient diagnoses in two different quarters within one year and a specific treatment or death within two quarters after the diagnosis. Exposures, derived from freely available 3rd party data, included particulate matter with aerodynamic diameter of less than 10 μm (PM10) and nitrogen dioxide (N02) as well as green space (Normalized Difference Vegetation Index - NDVI). Associations between air pollutants, green space, and cancer incidence were assessed by multilevel Poisson models. Age, sex, physician contacts, short- and long-term unemployment, population density, and having an alcohol-related disorder were considered as potential confounders.


Three thousand one hundred seven people developed mouth and throat cancer, 33,178 NMSC, 9611 prostate cancer, 9577 breast cancer, and 11,975 colorectal cancer during the follow-up period (2010–2014). An increase in PM10 of 10 μg/m3 was associated with a 53% increase in relative risk (RR) of mouth and throat cancer and a 52% increase in RR of NMSC. Prostate and breast cancer were modestly associated with PM10 with an increase in RR of 23 and 19%, respectively. The associations with N02 were in the same direction as PM10 but the effect estimates were much lower (7–24%). A 10% increase in NDVI was most protective of mouth and throat cancer (− 11% RR) and of NMSC (− 16% RR). Colorectal cancer was not affected by any of the exposures.


In addition to the studies carried out so far, this study was able to provide evidence that higher ambient AP levels increase the risk of mouth and throat cancer as well as of NMSC and that a higher residential green space level might have a protective effect for NMSC in areas with low to moderate UV intensity. Nevertheless, we cannot rule out residual confounding by socioeconomic or smoking status.

Peer Review reports


Outdoor air pollution (AP) is ubiquitous with its exposure having effects on a large proportion of the world population [1]. There is strong evidence from experimental and epidemiological studies that AP such as particulate matter (PM), nitrogen dioxide (NO2), and ozone (O3) are major risk factors for cardiovascular and cardiopulmonary diseases, and potentially cancer [2,3,4,5,6,7,8,9]. While the trigger function of AP on cardiac and pulmonary events is considered as causal, the role in cancer onset is only suggestive. In the case of PM and cardiovascular mortality, there is a solid association that fulfils both a temporal and a close exposure-response relationship. There is coherence of results between several scientific disciplines, including experimental studies offering plausible biological mechanisms reviewed in [3]. Both long-term and short-term effects are involved. An increase in PM2.5 long-term exposure per 10 μg/m3 increased the risk for cardiovascular mortality by 11% in a recent meta-analysis [5]. Even a daily rise in mean PM2.5 level per 10 μg/m3 increased the risk for cardiovascular mortality by approximately 0.4 to 1.0% [10]. In the case of cancer, it is much more difficult to prove causality. One reason is the long latency between exposure to carcinogens and the development of cancer. Another reason is that doses are typically low in the environment, and therefore direct causal inferences are hardly realized. Nevertheless, there is suggestive epidemiological evidence that outdoor AP increases incidence rates of some cancer types with strongest evidence for lung cancer [11,12,13,14,15,16,17,18,19]. Molecular epidemiological studies [20,21,22] have shown that the biological mechanisms causing cancer from outdoor AP involve genotoxic effects of the chemical compounds that accumulate over time, including PAH-DNA adducts, chromosome aberrations, sister chromatid exchanges, ras oncogene overexpression, and radically induced (oxidative) DNA damage. PM is listed as Group I carcinogen by the International Agency for Research on Cancer (IARC). Sources of PM, other than household, are numerous including traffic, agricultural and industrial emissions as main anthropogenic origins [23, 24]. We use PM10 (particles with aerodynamic diameter of less than 10 μm) instead of PM2.5, because we have more variation in PM10 data. PM2.5 is far more evenly distributed throughout Saxony. We use NO2as marker for traffic-related AP to investigate the effect of urban agglomeration in Saxony on cancer incidence [25, 26]. We also want to narrow the research gap on the impact of green spaces on health, especially cancer. While numerous studies have analyzed associations between AP and lung cancer, few studies have addressed associations between green space and cancer incidence [27,28,29,30]. Green space is known to be inversely related to AP due to lack of emission sources in green places [31]. Additionally, people residing in greener places might be more motivated to conduct physical activity [31], which is suggested to decrease cancer risk [32] through reductions of adipose tissue volume and endocrine activity (e.g. sex hormones) [33, 34]. On the other hand, some studies have also reported adverse associations with green space that may be context-specific, like a negative effect on the human skin through increased sunlight exposure, while spending time in green space [31]. Since a handful of studies have reported a link between lower all-cause mortality among people residing in green spaces [28, 35,36,37], it is a valid assumption that residential green space might decrease risk of some cancer types.

The aim of the study was to analyze the association between residential exposure to PM10 and NO2 as well as green space on different cancer types. We restricted the analysis to those cancer types, which are less strongly associated with smoking as lung cancer is, but which might be affected by AP exposure either by inhalation or dermal route of exposure to PM10 particles and NO2. Further, we needed a sufficient number of cases to ensure that we have enough power for our analyses. Therefore, we had to use cancers that are frequent in the population. We identified mouth and throat cancer, non-melanoma skin cancer (NMSC) and colorectal cancer in both sexes, as well as prostate cancer in men and breast cancer in women, which fulfilled these requirements. At the beginning of the study, we also wanted to examine the effects of pesticides on various forms of cancer in addition to AP and green spaces. According to the literature [15, 16], glandular tissues like the prostate or the breast may be susceptible to some agricultural used pesticides (hormone disruptive agents). Unfortunately, the data situation was not sufficient and therefore this investigation was not continued. Nevertheless, prostate and breast cancer were still used as outcomes.


AOK PLUS study population and case definitions

In Germany, approximately 90% of the population is covered by statutory health insurances. We used routine health care data from AOK PLUS, a large statutory health insurance in Saxony (area ~ 18,000 km2, population ~ 4 Mio), which covers almost half of the local general population. The data include information from inpatient and outpatient care with respect to diagnosis, procedures and prescriptions as well as socio-demographic information of the insured population such as age, sex and residential district (first four digits of the 5-digit postal code of the residential address). Age distribution and sex-ratio of the AOK PLUS beneficiaries in Saxony are comparable to the Germany-wide population [38]. We used data for the years 2007–2014 from the outpatient as well as from the inpatient sector. All beneficiaries were allocated to 186 four digit postal code districts based on their residential address in 2007. We excluded all cancer cases from the year 2007 for the determination of prevalent cases, because we had only outpatient data for this year. Further, beneficiaries that were diagnosed (outpatient or inpatient) in the years 2008 and 2009 with one of the analyzed cancer types were excluded, to estimate incidence rates for the years 2010 until 2014 (see Additional file 1: Figure S3). Age groups in increments of 10 years were built based on the age of every beneficiary in the year 2012. Therefore, the newborns of the years 2013 and 2014 had to be excluded as well. Further, the ages from 0 to 49 years were collapsed to one group. See Table 1 for an exact partitioning of ages.

Table 1 Baseline characteristics of the study population (AOK PLUS data)

In accordance with the Good epidemiological practice for secondary data analyses [39], cancer diseases were defined as one diagnosis of the corresponding ICD-10 code (International Classification of Diseases) for the inpatient data, while two diagnoses in different quarters within one year were necessary for all outpatient data plus a prescription of a specific treatment (e.g. radiotherapy, cytostatic medication) or death within two quarters after the second diagnosis (see Additional file 2: Table S1, Additional file 3: Table S2, Additional file 4: Table S3). We used the codes C00-C14 for mouth and throat cancer, C44 or L57.0 for NMSC, C61 for prostate cancer, C50 for breast cancer, and C18-C21 for colorectal cancer.

Ethics approval and consent to participate

The present analysis is based on secondary data from the health insurance company AOK PLUS, which were collected for the purpose of billing medical services. The study is supported by the AOK PLUS with which a data use and transfer agreement exists. The data was available for us only in anonymous form, so that no conclusions could be drawn about the individuals. Personal data of participants were anonymized through AOK PLUS before data sharing. Personal identifiers were masked or deleted (clear name to pseudonym; no social security number provided). Quasi-identifiers were generalized (only year of birth used; dropping of last digit of the zip code). Unfortunately, the data are not publicly accessible. There was no influence whatsoever on the policyholders and no intervention was carried out. It is therefore a purely observational study. According to paragraph 75 SGB X (Zehntes Buch Sozialgesetzbuch - German federal law) it is not reasonable to get consent for data sharing and analysis from around 2 million people, like in our investigation. Therefore, we submitted and get granted an application to the Saxon State Ministry for Social Affairs and Consumer Protection for obtaining consent to data transmission and analysis on behalf of the insured.

This study was conducted in accordance with the Helsinki Declaration [40] and follows the principles of Good Epidemiological Practice and Good Practice in Secondary Data Analysis [41]. The study was also registered in the database “Versorgungsforschung Deutschland” under the number VfD_ECo_epi_16_003770.

Air pollution and green space exposure assessment

Annual NO2 and PM10 concentrations for the year 2007 were derived from freely available maps in resolution of 100 m developed for Western Europe [42]. These maps were created by land use regression models based on more than 1500 EuroAirnet monitoring sites. Predictor variables for land use regression models included land use characteristics, population density, road length, altitude, distance to sea, and satellite-derived NO2 and PM10 data [43].

Green spaces were defined by the Normalized Difference Vegetation Index (NDVI), which was derived from freely available MODerate-resolution Imaging Spectroradiometer (MODIS) satellite images at the resolution of 250 m [44]. Briefly, NDVI is a commonly used indicator of vegetation level, ranging from − 1 (water) to + 1 (absolutely vegetated area). The algorithm for NDVI is based on two vegetation-informative bands: near-infrared (841 nm to 876 nm) and visible red (from 620 nm to 670 nm). For current analysis, we averaged 115 16-day composite NDVI images for the years 2005 to 2009 [45, 46].

These PM10, NO2, and NDVI estimates for 386 five digit postal codes were weighted by the population number from the German census 2011 and then averaged to 186 four digit postal code districts, because the address information were only available for the four digit postal codes due to privacy regulations. Freely available postal code vector data were obtained from

Geographic data management and calculations were conducted using the ArcGIS 10.1 Geographical Information System (GIS) (ESRI, Redlands, CA, USA) software program.

Confounders/effect modifiers

As control variables for the individual data analysis, age, sex, alcohol-related disorder, absolute number of physician contacts in the four digit postal code districts (AOK data) and the proportion of short- (up to one year) and long-term unemployment (more than 1 year; statistical office Saxony), ranging from 0 to 1, were considered. Additionally, for each four digit postal code district, population density per km2, all-cause mortality and proportion of persons with an alcohol-related disorder were considered in regression analysis on aggregated data.

Alcohol-related disorder was defined as one inpatient F10 diagnosis without F10.0 (i. e. acute intoxication) or three F10 prescriptions within four quarters of a year in the outpatient sector and was used to adjust the models to correct for its influence on cancer incidence [47, 48].

Statistical analysis

All the analyses were conducted with the program R, version 3.3.2 (Vienna, Austria) R Core Team [49], but data preprocessing was done with the software Stata (StataCorp. 2013. Stata Statistical Software: Release 13. College Station, TX: StataCorp LP). The individual associations between residential air pollutants, green space and cancer incidence rates were assessed by multilevel Poisson models with 95% confidence intervals with the R software package lme4. Due to the fact that not all individuals were fully insured or alive in the five year-long study period (2010–2014), we used the exact observation time in days as offset in the models. Because of high correlations (Pearson correlations above |0.75|) between PM10, NO2, and NDVI, associations with each exposure were analyzed individually (see Additional file 5: Figure S1).

For the aggregated data, cancer incidence rates of the 186 four digit postal codes were age-standardized by the European standard population [50] and adjusted for measurement errors due to the limited observation period using conditional autoregressive models (CAR), which use the first order spatial dependencies (shared borders) with the software BayesX for R [51, 52]. These associations were analyzed with linear models.

We used the best subset approach according to the lowest Akaike Information Criterion (AIC) for variable selection, implemented in the R package glmulti [53], to filter out unassociated variables for each model (adjusted models). We computed crude and adjusted models for each outcome and exposure pair (Additional file 6: Table S5).

We conducted a sensitivity analysis for the influence of the change of residence on the effect estimates by exclusion of movers and by comparing these results with the primary analysis where both, movers and non-movers, were included (Additional file 7: Table S4).


Baseline characteristics of the study population

Of more than 1.9 million initially cancer-free persons included into the analysis (Table 1), 3107 people developed mouth and throat cancer, 33,178 NMSC, 9611 prostate cancer, 9577 breast cancer, and 11,975 colorectal cancer during the follow-up period (2010–2014). Although there were slightly more females in the total population (53%), more males were affected by mouth and throat (74 vs. 26%) and colorectal cancer (53 vs. 47%), respectively. The number of incident cancer cases increased with age as expected, but decreased for the oldest age-groups (80–89, 90+). Mean age for developing cancer of interest ranged from 65 to 76 years. In the mouth and throat cancer subgroup, the proportion of people with an alcohol-related disorder reached almost 36%. Altogether, 69,722 people had a prevalent alcohol-related disorder (~ 3.6% from total) and 302,818 people (~ 15.8% from total) changed their place of residence between different four digit postal code districts at least once within the study period.

Air pollution and green space

PM10 and NO2 levels were differentially distributed over the 186 four digit postal code districts in Saxony (Fig. 1). While NO2 showed highest concentrations in agglomeration areas with high population density (cities Dresden, Leipzig, Chemnitz, and Zwickau), PM10 was much smoother distributed, even in rural areas, with highest concentrations in the north-west part of Saxony around Leipzig and in Middle Saxony to Dresden. Lowest concentrations of PM10 were observed in the Erzgebirge region in the south-west of Saxony. On the contrary, green space, as defined by NDVI, was higher in rural areas and the lowest in the metropolitan areas of Leipzig, Dresden, Chemnitz, and Görlitz and their vicinities. PM10 concentrations ranged from 15.47 to 26.30 μg/m3 (mean: 20.89), NO2 concentrations ranged from 9.32 to 31.55 μg/m3 (mean: 20.44), and NDVI varied between 0.38 and 0.64 (mean: 0.51). The admissible annual averages for PM10 and NO2 in accordance with the Federal Immission Control Act (39. BImSchV) and the European Directive on Air Quality (2008/50/EG) are 40 μg/m3 each.

Fig. 1
figure 1

Exposure and Case maps; Mean concentrations of the exposures PM10 and NO2 (μg/m3) and mean NDVI (0 to + 1) are shown. Case maps of specific cancer types in Saxony, population density and mean physician contacts per year over the years 2010 until 2014 (data source: AOK PLUS) are given. Age-standardized cancer incidences (per 100.000 persons) were smoothed using a Bayesian CAR model

Additional file 5: Figure S1 shows Pearson correlation coefficients for associations between AP and NDVI, as well as physician contacts (PC). Most correlations were at least 75%, except for NDVI and PC (~ − 58%) or NO2 and PC (~ 61%).

Association between air pollution, green space, and cancer incidence

Semi-individual data modeling

An increase in PM10 of 10 μg/m3 was associated with a 53% increase in relative risk (RR) of mouth and throat cancer and a 52% increase in RR of NMSC (Table 2). Prostate and breast cancer were modestly associated with PM10 with an increase in RR of 23 and 19%, respectively. The associations with N02 were in the same direction as PM10 but the effect estimates were much weaker (7–24%) and did not reach statistical significance for prostate cancer. An increase in NDVI by 10% revealed associations with mouth and throat cancer with an 11% decrease in RR and with NMSC with a 16% decrease in RR. No associations were found for prostate or breast cancer and NDVI. Colorectal cancer was not affected by any of the exposures (Table 2).

Table 2 Relative risk (RR) estimates from multilevel Poisson regression models with observation time as offset and controlled for age as cubic term; 1.9 Mio. Persons in 186 postal code districts in Saxony were considered. 95% Wald confidence intervals (CI) are given in brackets; NMSC – non-melanoma skin cancer 

Alcohol-related disorder was associated with a more than 9 times increased risk of mouth and throat cancer a 50% increase of colorectal cancer and a 22% increase of breast cancer. Men were also more affected by any of the cancers which occur in both sexes (RR 1.6 to 2.7×; Table 2).

When only the people who did not change the place of residence were considered, similar results were observed except for the effect of PM10 on mouth and throat cancer which increased by 8% (Additional file 7: Table S4).

Aggregated data analysis

Cancer incidence rates generally showed a scattered distribution over Saxony (Fig. 1). The range of differences between the four digit postal code districts was small for colorectal cancer (86 to 94 cases per 100.000 insured persons) and breast cancer (144 to 158 cases per 100.000 insured persons), but high for mouth and throat cancer (22 to 36 cases), prostate cancer (140 to 240 cases), and especially for NMSC (150 to 400 cases). High correlation of physician contacts with PM10 (~ 75%) prevented using this variable of personal demand as a confounder in the adjusted models.

Crude linear regression analyses showed associations between PM10 and mouth and throat cancer (R2 0.203) or NMSC (R2 0.144). Associations with breast- or prostate cancer were weak (Fig. 2). NO2 was associated with NMSC (R2 0.164), but asssociations with breast-, and prostate cancer, or colorectal cancer were weak (Additional file 8: Figure S2 and Additional file 6: Table S5 in supplement). Moderate negative correlations were observed for NDVI and mouth and throat cancer (R2 0.114) or NMSC (R2 0.137), but again, the association with prostate cancer was weak (Additional file 8: Figure S2 and Additional file 6: Table S5).

Fig. 2
figure 2

Scatter plots of the crude linear regression analysis; Assoziations between PM10 in μg/m3 and age-standardized cancer incidence rates per 100.000 persons. For each diagram, coefficients of determination (R2) are given which correspond to a moderate model fit for mouth and throat cancer and for NMSC, but show only a poor model fit for breast- and prostate cancer

Adjusted linear regression models were in the same direction as the crude analysis. Also PM10 and mouth and throat cancer (R2 0.331) or NMSC (R2 0.258) showed associations, but effect estimates were lower compared to the crude analyses (Additional file 6: Table S5). NO2showed also some association with NMSC (R2 0.259) in the adjusted analysis, while associations to prostate-, breast-, and colorectal cancer remained weak. Negative correlations between NDVI and mouth and throat cancer (R2 0.275) or NMSC (R2 0.250) were still present in the adjusted analysis, but again, effect estimates were lower than in the crude models.


To the best of our knowledge, this study is novel and substantially adds to the current knowledge, because first of all the modeling is based on a very large cohort (1.9 million beneficiaries), which is needed to study associations between environmental exposures and cancer incidence. While earlier work has exclusively focused on the association between AP and lung cancer [1, 12,13,14,15,16,17,18,19, 54,55,56], we observed that PM10 is also related to mouth and throat cancer, NMSC, and maybe also to prostate and breast cancer. It is not so surprising that other types of cancer such as mouth and throat cancer or NMSC are also associated with AP, since chemical compounds are inhaled from the tidal air or even passed through the skin into the human body and accumulate there in the respective tissue over time. Some of these compounds could induce different forms of DNA damage [20,21,22], which could lead to tumor development. On the other hand, no association was found with colorectal cancer. Perhaps the carcinogens of the ambient air cannot accumulate in the intestine as they are already stopped by the lung or skin tissue beforehand. Associations between cancer incidence and green space (NDVI) or NO2 were present, but much smaller than for PM10.

We speculate that the high regional variation in NMSC (Fig. 1) could be due to differences in health care utilization in the population, whereas this is unlikely for more severe cancer types. One Australian study reported a higher skin cancer risk for regions with more vegetation [27], but exposure to UV light is certainly significantly greater in Australia compared to our study setting. Our results indicate that green space might be protective against NMSC in Germany. This observed protective effect might be secondary as the distribution of green spaces and PM10 were highly inversely correlated, and, as aforementioned, there are fewer emission sites around green spaces [31]. Thus, in areas where the UV intensity is much weaker than in Australia like in Germany, residential green space could lower the skin cancer risk. However, more studies are needed to clarify whether residing in greener places is beneficial or detrimental for NMSC and if this effect is causal or only indirectly transmitted through a lack of AP at these sites.

In the age group 18 until 64 years 5.0% of the sample population had an alcohol-related disorder. Survey data from 2012 for Germany and this age group found 3.1% with alcohol abuse and additionally 3.4% with alcohol dependency according to DSM IV. Our 5.0% estimate lies between 3.1 and 6.5% [57]. We expect therefore, that not all cases of alcohol abuse or dependency have already been recognized by physicians, but yet we consider it a valuable confounding variable. Further, a recent meta-analysis [58] and other studies [48, 59] found positive associations with high alcohol consumption on colorectal, mouth and throat, as well as breast cancer, and limited evidence for prostate cancer. Our effect estimate of alcohol on mouth and throat cancer is possibly a bit overestimated, as part of the people with an alcohol-related disorder are also smokers and these risk factors interact with each other [59].

We speculate that the high correlation of physician contacts and PM10 (~ 75%) suggests that PM10 is a wide-ranging risk factor for many diseases and therefore reducing it should be a central focus for community health.

Strength and limitations

One of the strengths of the presented analysis is that the dataset contains a huge population. The population covers almost half of the residential population of Saxony and is regarding age and sex distribution very similar to the general population of Saxony, which allows us to investigate cancer incidence for the federal state of Saxony. Further, the study offers complete information about out- and inpatient treatments and is not affected by sampling or nonresponse bias. Selection bias because of changes of the health insurance due to the cancer disease is also unlikely, since no differences in cancer treatment between statutory health insurances are present in Germany. We have invested a lot of work in exact case definitions and can therefore build our analyses on valid cancer diagnoses.

Nevertheless, several important limitations should be acknowledged. Address information up to four digit postal code districts was available, which did not allow a high spatial resolution of individual outcome and exposure assessments. Therefore, we could only use the exposure of the four digit postal code district as a surrogate for the individual exposure (semi-individual). Additional information on some potentially important individual confounders, such as smoking, socioeconomic status, diet, or physical activity was unfortunately not available, which could have led to residual confounding. Further confounding could occur due to genetic predisposition of patients for certain cancer types, or due to other environmental pollutants, like pesticides. But farm workers, with a potential high exposition with pesticides, are not included in our data set, since they have their own statutory health insurance in Germany. Unfortunately, we had no information about occupation, indoor radon exposure of patients and virus infections possible related to cancer which could also confound our analyses.

The observed associations could therefore be overestimated or be biased due to exposure misclassification, for example, but the magnitude should not be too high. In one study on lung cancer, the hazard ratio of PM10 decreased by only 11% if it was additionally controlled for smoking status, smoking intensity, square of smoking intensity, smoking duration, time since quitting smoking, environmental tobacco smoke, occupation, fruit intake, marital status, education level, and employment status [54]. Our models assumed that the residential population did not move and AP was constant over time. Both assumptions, especially the first, are unlikely. Nevertheless, it is reassuring that our findings remained robust, when analyzing only the population that did not move outside their four digit postal code district during the observation period. Still, it is unknown what their movement history was until 2007. From the literature it is known that migrants have better health because people with low socioeconomic status or severe diseases are less mobile [60]. This could bias our results and therefore, we trust our primary analysis more than the results of the sensitivity analysis with excluded migrants.


Beyond the current study situation we found some evidence that higher ambient AP levels increase the risk of mouth and throat cancer and of NMSC, while a higher residential green space level might have a protective effect. In summary, we assume that our effect estimates are not strongly biased by residual confounding, but we cannot exclude that for sure. Further research should try to measure the environmental exposures through the life course and focus on the effects of relocation together with environmental factors and extent the analyses to other cancer types.



Akaike information criterion


Air pollution


Diagnostic and Statistical Manual of Mental Disorders


Normalized Difference Vegetation Index


Non-melanoma skin cancer

NO2 :

Nitrogen dioxide


Physician contacts


Particulate matter


Relative risk


Uniform Value Scale


German modification of the International Classification of Procedures in Medicine


Pharmaceutical registration numbers


  1. Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, Balakrishnan K, Brunekreef B, Dandona L, Dandona R, et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the global burden of diseases study 2015. Lancet (London, England). 2017;389(10082):1907–18.

    Article  Google Scholar 

  2. Brook RD, Newby DE, Rajagopalan S. The global threat of outdoor ambient air pollution to cardiovascular health: time for intervention. JAMA cardiology. 2017;2(4):353–4.

    Article  PubMed  Google Scholar 

  3. Brook RD, Rajagopalan S, Pope CA 3rd, Brook JR, Bhatnagar A, Diez-Roux AV, Holguin F, Hong Y, Luepker RV, Mittleman MA, et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121(21):2331–78.

    Article  PubMed  CAS  Google Scholar 

  4. Cosselman KE, Navas-Acien A, Kaufman JD. Environmental factors in cardiovascular disease. Nat Rev Cardiol. 2015;12(11):627–42.

    Article  PubMed  CAS  Google Scholar 

  5. Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, Kaufman JD. Long-term air pollution exposure and cardio- respiratory mortality: a review. Environmental health : a global access science source. 2013;12(1):43.

    Article  CAS  Google Scholar 

  6. Lee BJ, Kim B, Lee K. Air pollution exposure and cardiovascular disease. Toxicological research. 2014;30(2):71–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Li J, Sun S, Tang R, Qiu H, Huang Q, Mason TG, Tian L. Major air pollutants and risk of COPD exacerbations: a systematic review and meta-analysis. International journal of chronic obstructive pulmonary disease. 2016;11:3079–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Faustini A, Rapp R, Forastiere F. Nitrogen dioxide and mortality: review and meta-analysis of long-term studies. Eur Respir J. 2014;44(3):744–53.

    Article  PubMed  CAS  Google Scholar 

  9. Atkinson RW, Butland BK, Dimitroulopoulou C, Heal MR, Stedman JR, Carslaw N, Jarvis D, Heaviside C, Vardoulakis S, Walton H, et al. Long-term exposure to ambient ozone and mortality: a quantitative systematic review and meta-analysis of evidence from cohort studies. BMJ Open. 2016;6(2):e009493.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Pope CA 3rd, Dockery DW. Health effects of fine particulate air pollution: lines that connect. Journal of the Air & Waste Management Association (1995). 2006;56(6):709–42.

    Article  CAS  Google Scholar 

  11. WHO. An overview of the evidence on environmental and occupational determinants of Cancer. In. 2011:1–7.

  12. Perera FP. Environment and cancer: who are susceptible? Science. 1997;278(5340):1068–73.

    Article  PubMed  CAS  Google Scholar 

  13. Hamra GB, Laden F, Cohen AJ, Raaschou-Nielsen O, Brauer M, Loomis D. Lung Cancer and exposure to nitrogen dioxide and traffic: a systematic review and meta-analysis. Environ Health Perspect. 2015;123(11):1107–12.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Loomis D, Grosse Y, Lauby-Secretan B, El Ghissassi F, Bouvard V, Benbrahim-Tallaa L, Guha N, Baan R, Mattock H, Straif K, et al. The carcinogenicity of outdoor air pollution. The Lancet Oncology. 2013;14(13):1262–3.

    Article  PubMed  CAS  Google Scholar 

  15. Clapp RW, Jacobs MM, Loechler EL. Environmental and occupational causes of cancer: new evidence 2005-2007. Rev Environ Health. 2008;23(1):1–37.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Blackadar CB. Historical review of the causes of cancer. World journal of clinical oncology. 2016;7(1):54–86.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Pope CA 3rd, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, Thurston GD. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. Jama. 2002;287(9):1132–41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Perez L, Rapp R, Kunzli N. The year of the lung: outdoor air pollution and lung health. Swiss Med Wkly. 2010;140:w13129.

    PubMed  Google Scholar 

  19. Fajersztajn L, Veras M, Barrozo LV, Saldiva P. Air pollution: a potentially modifiable risk factor for lung cancer. Nat Rev Cancer. 2013;13(9):674–8.

    Article  PubMed  CAS  Google Scholar 

  20. Farmer PB, Singh R, Kaur B, Sram RJ, Binkova B, Kalina I, Popov TA, Garte S, Taioli E, Gabelova A, et al. Molecular epidemiology studies of carcinogenic environmental pollutants. Effects of polycyclic aromatic hydrocarbons (PAHs) in environmental pollution on exogenous and oxidative DNA damage Mutation research. 2003;544(2–3):397–402.

    PubMed  CAS  Google Scholar 

  21. Perera FP, Hemminki K, Gryzbowska E, Motykiewicz G, Michalska J, Santella RM, Young TL, Dickey C, Brandt-Rauf P, De Vivo I, et al. Molecular and genetic damage in humans from environmental pollution in Poland. Nature. 1992;360(6401):256–8.

    Article  PubMed  CAS  Google Scholar 

  22. Perera FP, Jedrychowski W, Rauh V, Whyatt RM. Molecular epidemiologic research on the effects of environmental pollutants on the fetus. Environ Health Perspect. 1999;107(Suppl 3):451–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Council NR. Global sources of local pollution: an assessment of long-range transport of key air pollutants to and from the United States. Washington, DC: the National Academies Press; 2010.

    Google Scholar 

  24. Weichenthal S, Farrell W, Goldberg M, Joseph L, Hatzopoulou M. Characterizing the impact of traffic and the built environment on near-road ultrafine particle and black carbon concentrations. Environ Res. 2014;132:305–10.

    Article  PubMed  CAS  Google Scholar 

  25. Parent ME, Goldberg MS, Crouse DL, Ross NA, Chen H, Valois MF, Liautaud A. Traffic-related air pollution and prostate cancer risk: a case-control study in Montreal, Canada. Occup Environ Med. 2013;70(7):511–8.

    Article  PubMed  Google Scholar 

  26. Al-Ahmadi K, Al-Zahrani A. NO(2) and cancer incidence in Saudi Arabia. Int J Environ Res Public Health. 2013;10(11):5844–62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Astell-Burt T, Feng X, Kolt GS. Neighbourhood green space and the odds of having skin cancer: multilevel evidence of survey data from 267072 Australians. J Epidemiol Community Health. 2014;68(4):370–4.

    Article  PubMed  Google Scholar 

  28. Gascon M, Triguero-Mas M, Martinez D, Dadvand P, Rojas-Rueda D, Plasencia A, Nieuwenhuijsen MJ. Residential green spaces and mortality: a systematic review. Environ Int. 2016;86:60–7.

    Article  PubMed  Google Scholar 

  29. Richardson EA, Mitchell R. Gender differences in relationships between urban green space and health in the United Kingdom. Soc Sci Med. 2010;71(3):568–75.

    Article  PubMed  Google Scholar 

  30. Mitchell R, Popham F. Effect of exposure to natural environment on health inequalities: an observational population study. Lancet. 372(9650):1655–60.

  31. Markevych I, Schoierer J, Hartig T, Chudnovsky A, Hystad P, Dzhambov AM, de Vries S, Triguero-Mas M, Brauer M, Nieuwenhuijsen MJ, et al. Exploring pathways linking greenspace to health: theoretical and methodological guidance. Environ Res. 2017;158:301–17.

    Article  PubMed  CAS  Google Scholar 

  32. Kerr J, Anderson C, Lippman SM. Physical activity, sedentary behaviour, diet, and cancer: an update and emerging new evidence. The Lancet Oncology. 2017;18(8):e457–71.

    Article  PubMed  Google Scholar 

  33. Friedenreich CM, Shaw E, Neilson HK, Brenner DR. Epidemiology and biology of physical activity and cancer recurrence. Journal of molecular medicine (Berlin, Germany). 2017;

  34. McTiernan A. Mechanisms linking physical activity with cancer. Nat Rev Cancer. 2008;8(3):205–11.

    Article  PubMed  CAS  Google Scholar 

  35. James P, Hart JE, Banay RF, Laden F. Exposure to greenness and mortality in a Nationwide prospective cohort study of women. Environ Health Perspect. 2016;

  36. Toftager M, Ekholm O, Schipperijn J, Stigsdotter U, Bentsen P, Gronbaek M, Randrup TB, Kamper-Jorgensen F. Distance to green space and physical activity: a Danish national representative survey. J Phys Act Health. 2011;8(6):741–9.

    Article  PubMed  Google Scholar 

  37. van den Berg M, Wendel-Vos W, van Poppel M, Kemper H, van Mechelen W, Maas J. Health benefits of green spaces in the living environment: a systematic review of epidemiological studies. Urban For Urban Green. 2015;14(4):806–16.

    Article  Google Scholar 

  38. Trautmann F, Schuler M, Schmitt J. Burden of soft-tissue and bone sarcoma in routine care: estimation of incidence, prevalence and survival for health services research. Cancer Epidemiol. 2015;39(3):440–6.

    Article  PubMed  Google Scholar 

  39. Swart E, Gothe H, Geyer S, Jaunzeme J, Maier B, Grobe TG, Ihle P. German Society for Social M, prevention, German Society for E: [good practice of secondary data analysis (GPS): guidelines and recommendations]. Gesundheitswesen. 2015;77(2):120–6.

    Article  PubMed  CAS  Google Scholar 

  40. World Medical Association Declaration of Helsinki. ethical principles for medical research involving human subjects. Jama. 2013;310(20):2191–4.

    Article  CAS  Google Scholar 

  41. Swart E, Gothe H, Geyer S, Jaunzeme J, Maier B, Grobe TG, Ihle P. Good practice of secondary data analysis (GPS): guidelines and recommendations. Gesundheitswesen. 2015;77(2):120–6.

    Article  PubMed  CAS  Google Scholar 

  42. Small Area Health Statistics Unit [].

  43. Vienneau D, de Hoogh K, Bechle MJ, Beelen R, van Donkelaar A, Martin RV, Millet DB, Hoek G, Marshall JD. Western European land use regression incorporating satellite- and ground-based measurements of NO2 and PM10. Environmental science & technology. 2013;47(23):13555–64.

    Article  CAS  Google Scholar 

  44. Moderate Resolution Imaging Spectroradiometer Land Product Subsets [].

  45. Tucker CJ. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens Environ. 1979;8(2):127–50.

    Article  Google Scholar 

  46. Weier J, Herring D: Measuring vegetation (NDVI & EVI). 2000.

    Google Scholar 

  47. IARC_Working_Group. Personal habits and indoor combustions. Volume 100. E A review of human carcinogens IARC Monographs on the Evaluation of Carcinogenic Risks in Humans. 2012;100(Pt E):373–472.

    Google Scholar 

  48. IARC_Working_Group. Alcohol consumption and ethyl carbamate. IARC Monographs on the Evaluation of Carcinogenic Risks in Humans. 2010;96:3–1383.

    Google Scholar 

  49. R Development Core team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2010.

    Google Scholar 

  50. Eurostat: Revision of the European Standard Population - Report of Eurostat's task force. 2013.

    Google Scholar 

  51. Umlauf N, Adler D, Kneib T, Lang S, Zeileis A. Structured Additive Regression Models. An R Interface to BayesX: Faculty of Economics and Statistics, University of Innsbruck; 2012.

    Google Scholar 

  52. Brezger A, Kneib T, Lang S. BayesX: analyzing Bayesian structural additive regression models. J Stat Softw. 2005;014(i11)

  53. Calcagno V, de Mazancourt C. Glmulti: an R package for easy automated model selection with (generalized) linear models. J Stat Softw. 2010;34(12):1–29.

    Article  Google Scholar 

  54. Raaschou-Nielsen O, Andersen ZJ, Beelen R, Samoli E, Stafoggia M, Weinmayr G, Hoffmann B, Fischer P, Nieuwenhuijsen MJ, Brunekreef B, et al. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European study of cohorts for air pollution effects (ESCAPE). The Lancet Oncology. 2013;14(9):813–22.

    Article  PubMed  Google Scholar 

  55. Raaschou-Nielsen O, Beelen R, Wang M, Hoek G, Andersen ZJ, Hoffmann B, Stafoggia M, Samoli E, Weinmayr G, Dimakopoulou K, et al. Particulate matter air pollution components and risk for lung cancer. Environ Int. 2016;87:66–73.

    Article  PubMed  CAS  Google Scholar 

  56. Vineis P, Husgafvel-Pursiainen K. Air pollution and cancer: biomarker studies in human populations. Carcinogenesis. 2005;26(11):1846–55.

    Article  PubMed  CAS  Google Scholar 

  57. Kraus L, Pabst A, Piontek D, EGd M. Substanzkonsum und substanzbezogene Störungen: Trends in Deutschland 1980 – 2012. Sucht. 2013;59(6):333–45.

    Article  Google Scholar 

  58. Bagnardi V, Rota M, Botteri E, Tramacere I, Islami F, Fedirko V, Scotti L, Jenab M, Turati F, Pasquali E, et al. Alcohol consumption and site-specific cancer risk: a comprehensive dose-response meta-analysis. Br J Cancer. 2015;112(3):580–93.

    Article  PubMed  CAS  Google Scholar 

  59. Purdue MP, Hashibe M, Berthiller J, La Vecchia C, Dal Maso L, Herrero R, Franceschi S, Castellsague X, Wei Q, Sturgis EM, et al. Type of alcoholic beverage and risk of head and neck cancer--a pooled analysis within the INHANCE consortium. Am J Epidemiol. 2009;169(2):132–42.

    Article  PubMed  Google Scholar 

  60. Eronen J, von Bonsdorff M, Rantakokko M, Portegijs E, Viljanen A, Rantanen T. Socioeconomic status and life-space mobility in old age. J Aging Phys Act. 2016;24(4):617–23.

    Article  PubMed  Google Scholar 

Download references


We thank the health insurance AOK PLUS for cooperation in data utilization and Stefan Baßler for proofreading the manuscript. We thank Beate Brecht for conducting data queries. We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the TU Dresden.


The project was carried out without funding.

Availability of data and materials

The data that support the findings of this study are available from AOK PLUS, but restrictions apply to the availability of these data which were used under license for the current study and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of AOK PLUS.

Author information

Authors and Affiliations



TD and FT2 designed and performed the statistical analyses and drafted the manuscript. IM, FT2 and TD conducted geographical analyses. TD, FT2, IM, FT1, JS, and JH revised the manuscript and contributed to the interpretation of the results. JH and JS initiated and monitored the study. All authors approved the submitted manuscript and agree to be responsible for this work.

Corresponding author

Correspondence to Thomas Datzmann.

Ethics declarations

Ethics approval and consent to participate

The consultation of an ethics committee and consent to participate is not required for analyses based solely on secondary data.

Competing interests

Jochen Schmitt reports institutional funding for IITs from Sanofi, Novartis, ALK, and Pfizer. The other authors do not indicate any financial or non-financial competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S3. STROBE diagram; Flowchart showing selection of incident cancer cases following the guidelines of the STROBE initiative (STrengthening the Reporting of OBservational studies in Epidemiology). (PNG 26412 kb)

Additional file 2:

Table S1. ICD-10 Codes; Analyzed cancer entities and corresponding ICD-10-GM codes are shown. (DOCX 12 kb)

Additional file 3:

Table S2. Validation of outpatient cancer cases other than NMSC; Used OPS-, EBM codes and prescribed medications in outpatient cancer care. (DOCX 14 kb)

Additional file 4:

Table S3. Validation of outpatient cancer cases (NMSC); Used OPS- and EBM codes in outpatient cancer care of NMSC. (DOCX 14 kb)

Additional file 5:

Figure S1. Correlation plot; Correlation between AP (PM10; NO2), NDVI, and physician contacts. Numbers are Pearson correlation coefficients. (PNG 1050 kb)

Additional file 6:

Table S5. Crude and adjusted linear regression of aggregated data; Results of the linear regression analysis are shown. Estimators, 95% confidence intervals (CI) and coefficients of determination (R2) are given for crude and adjusted models. (DOCX 19 kb)

Additional file 7:

Table S4. Non-mover multilevel Poisson regression (CIs are Wald confidence intervals); Relative risk estimates for an increase of environmental exposures (10 μg/m3 for PM10 / N02; 10% for NDVI) on cancer incidence in Saxony for patients with no relocation within study period between 2010 until 2014. (DOCX 15 kb)

Additional file 8:

Figure S2. Scatter plots of the crude linear regression analysis; Associations between PM10, NDVI, and NO2 and different cancer types. For each diagram, coefficients of determination (R2) are given. We see positive but weak associations between NO2 and mouth and throat cancer, prostate cancer and breast cancer, but an elevated positive association with NMSC exists (R2 = 0.164). With increasing vegetation level (NDVI) of the neighborhood, cancer incidence rates decrease for breast cancer, prostate cancer, mouth and throat cancer, and NMSC (in increasing order). (PNG 1295 kb)

Supplemental introduction text

All case definitions were based on the respective coding systems for diagnosis (ICD-10-GM), procedures (Uniform Value Scale (EBM) and German modification of the International Classification of Procedures in Medicine (OPS)), as well as prescriptions (Anatomical Therapeutic Chemical code - ATC and pharmaceutical registration numbers - PZN).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Datzmann, T., Markevych, I., Trautmann, F. et al. Outdoor air pollution, green space, and cancer incidence in Saxony: a semi-individual cohort study. BMC Public Health 18, 715 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: