Research | Open | Published:
Correlations between meteorological parameters and prostate cancer
International Journal of Health Geographicsvolume 9, Article number: 19 (2010)
There exists a north-south pattern to the distribution of prostate cancer in the U.S., with the north having higher rates than the south. The current hypothesis for the spatial pattern of this disease is low vitamin D levels in individuals living at northerly latitudes; however, this explanation only partially explains the spatial distribution in the incidence of this cancer. Using a U.S. county-level ecological study design, we provide evidence that other meteorological parameters further explain the variation in prostate cancer across the U.S.
In general, the colder the temperature and the drier the climate in a county, the higher the incidence of prostate cancer, even after controlling for shortwave radiation, age, race, snowfall, premature mortality from heart disease, unemployment rate, and pesticide use. Further, in counties with high average annual snowfall (>75 cm/yr) the amount of land used to grow crops (a proxy for pesticide use) was positively correlated with the incidence of prostate cancer.
The trends found in this USA study suggest prostate cancer may be partially correlated with meteorological factors. The patterns observed were consistent with what we would expect given the effects of climate on the deposition, absorption, and degradation of persistent organic pollutants including pesticides. Some of these pollutants are known endocrine disruptors and have been associated with prostate cancer.
Approximately one in six men will develop prostate cancer in their life-time . However, the risk of prostate cancer is not the same across the United States; northern counties tend to have a higher incidence of the disease than southern counties [2–5] (Figure 1). This north-south pattern in prostate cancer has also been reported in other areas of the world . The current hypothesis for this distribution is that lower exposure to ultraviolet (UV) radiation in the northern states, especially during the winter months, results in lower vitamin D synthesis [2, 3, 5–8]. This vitamin regulates transcription in cells with vitamin D receptors and therefore insufficient levels may increase the risk of prostate cancer [6, 9]. A recent U.S. study on prostate cancer found approximately 5.5% of the variation in this disease could be explained by the UV index ; however, this study did not control for potential confounders.
Several other spatially-distributed factors may contribute to the north-south disease pattern. For example, meteorological parameters such as air temperature, snowfall, and rainfall all vary spatially and it is well documented that these parameters affect the deposition, absorption, degradation of persistent organic pollutants (POPs) [10–13]. Cold trapping and snow scavenging are believed to be the reason why some POPs are found at higher concentration with increasing northern latitudes [10, 13, 14], and may therefore play a significant role in the level of pollutants to which individuals in different geographical areas are exposed. The purpose of this study was to determine whether there was a correlation between meteorological parameters and county-level incidence rates of prostate cancer in the U.S., controlling for exposure to local pesticide use, air pollution, and other known risk factors for prostate cancer that may vary spatially.
The average annual incidence rate of prostate cancer in this study ranged between 39.5 and 311.1 cases per 100,000. Comparison of the mean incidence rate of prostate cancer for counties in the upper and lower quartiles for different meteorological parameters suggested there was a difference in cancer rates between these groups for individual climate parameters not controlling for other variables (Table 1).
Various biologically relevant models were developed using ordinary least squares (OLS) regressions (Table 2). These were developed by building on previous published models with one parameter, UV radiation. The last model we constructed, which included all significant variables available to us (Table 3) suggested radiation and temperature were best modeled using a quadratic term, and several parameters beside UV radiation were correlated with the incidence of prostate cancer (Table 3). In all our models that included meteorological parameters, UV radiation, rainfall, and temperature were always negatively correlated with prostate cancer (Table 2; Figure 2 and 3). Note that HDD is positively associated with prostate cancer, which reflects a negative correlation between temperature and the disease. The higher the HDD value the colder the county. Our index for pesticide use (acres of land used to grow crops) was positively correlated with prostate cancer, but only in counties where there was snow (Figure 4). The potential confounders in our model included premature death from heart disease and unemployment rate. These were both negatively correlated with prostate cancer in all our models (Table 2). Variables that were not significant in our models included EPA permitted air emissions for various pollutants, number of individuals residing in each county, and all interaction terms evaluated between meteorological parameters and pollution indices except acres used to grow crops crossed with snow.
Of the models developed using OLS analyses, the best fit model in a geographically weighted regression (GWR) analysis, based on the Akaike's Information criteria (AIC), included meteorological parameters (shortwave radiation, radiation2, HDD, HDD2, snowfall, and rainfall), confounders (premature mortality associated with heart disease and unemployment rate), the pollution index for pesticide use (acres of land used to grow crops), and the interaction term acres of land used to grow crops crossed with annual snowfall (last model in Table 2). This model explained approximately 43% of the variation in the county level incidence rate of prostate cancer (R2 = 0.43) (Table 2). In comparison, the GWR model with only shortwave radiation explained approximately 31% of the variation in prostate cancer (R2 = 0.31) (Table 2). The modeling assumptions for the final model with the lowest AIC were satisfied, the residuals were approximately normal and there were fewer than 0.58% (16/2571) of the counties with standardized residuals greater than 3 standard deviations above or below the mean (Figure 5). Further, the counties with these extreme residuals were scattered throughout the U.S. and did not cluster in a particular area.
Our analyses suggest meteorological conditions, including daily shortwave radiation, heating degree days (HDD), which is defined as the annual sum of degrees Celsius required to attain 18.3 °C when the air temperature is less than 18.3°C, and average annual snowfall and rainfall, were significantly correlated with the average annual county-level incidence rates of prostate cancer (Tables 2 and 3). This study confirmed the negative correlation between shortwave radiation and the incidence of prostate cancer (Tables 2 and 3). This was consistent with previous analyses and with the hypothesis that lower exposure to UV radiation results in lower Vitamin D synthesis [2, 3, 5–8]. UV radiation may also reduce the risk of cancer by increasing the photodegradation of some chemicals, including pesticides [15, 16]. We improved the previously described UV model for prostate cancer  by including a quadratic term, which suggests there may be an upper threshold effect to the benefit of UV radiation (Figure 2). However, even with this parameter in our model, other meteorological variables appear to be significantly correlated with this cancer.
Temperature was negatively correlated with prostate cancer, after controlling for shortwave radiation, local pesticide use, rainfall, snowfall, premature mortality from heart disease, and unemployment rate. At 3000 degree days, the median value for HDD in our dataset, or higher our final OLS regression model suggests a positive relationship between HDD and the incidence of prostate cancer: the higher the HDD the colder the county (Table 3 and Figure 3). A quadratic term best described the relationship between HDD and prostate cancer (Table 3). Based on our OLS model the correlation between temperature and prostate cancer was biologically negligible for counties with less than 3000 HDD, but after the HDD reached this threshold it was positively correlated with prostate cancer (Figure 3). Interestingly, the model with only HDD and the confounders (premature mortality from heart disease and unemployment) had a lower AIC than the model with only shortwave radiation and the same confounders (Table 2).
We hypothesize that temperature may be associated with the incidence of prostate cancer by modulating exposure to POPs, some of which have been linked to the disease. Temperature affects POPs in a number of ways. For example, cold temperature increases the solid phase portioning of POPs . Organic chemicals, especially semi-volatile organic contaminants such as polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), and organophosphate and organochlorines pesticides, favor a solid phase rather than a gaseous phase at cold temperatures, which causes them to precipitate to the earth's surface [10, 13, 14]. Cold trapping of chemicals partially explains the presence of PCBs and other pollutants in pristine areas at high altitude and latitude [10, 12, 17]. Some semi-volatile compounds are known to be endocrine disruptors (i.e. PCBs, Alpha HCH, gamma HCH, PeCB HCB, and alpha endosulphans) [18, 19], and their increased deposition at colder temperatures may predispose these places to endocrine responsive diseases (i.e. prostate cancer). Similar volatilization occurs with some persistent organic pesticides . Several pesticides have been identified as endocrine disruptors  and have been associated with prostate cancer [22–26].
Temperature also affects the degradation of POPs in the soil and the atmosphere [27, 28]. Experiments have demonstrated that the biodegradation of certain organic compounds by microorganisms is temperature-dependant and slower at colder temperatures . Chemical reactions in general are slower at colder temperatures. Lower degradation of POPs at northerly latitudes suggests that environmental bioaccumulation of some pollutants may be greater in the northern part of the U.S. than in the south where temperature-dependant biodegradation processes are more productive.
Humidity also plays an important role in absorption and degradation of POPs. In general, the higher the humidity the greater the absorption and degradation of non-polar semi-volatile compounds, such as PCBs, PAHs, and the less volatilization of these compounds [27–31]. We observed a strong negative correlation between the incidence of prostate cancer and rainfall (Tables 1 and 2), which may reflect the increase in absorption and degradation of organic pollutants in moist soils and the decrease in volatilization of these compounds in humid environments [21, 28].
In all our models the amount of land used to grow crops was significantly correlated with prostate cancer (Table 2). This is consistent with several studies that have found some types of pesticides associated with this cancer [22–25]. Interestingly this relationship was stronger in the counties with a high average annual snowfall (> 40 cm/year). Areas with low average annual snowfall did not have a significant relationship between land used to grow crops and prostate cancer (Figure 4).
There may be several possible explanations for the interaction between acres of land used to grow crops and average annual snowfall. Snowflakes have a high surface area-to-volume ratio and, as they fall through the atmosphere, they scavenge and collect small particulate matter, including PCBs and PAHs (suspected endocrine disruptors) . Snow trapping of atmospheric pollutants may compound the effect of pesticides by increasing the deposition of POPs. Also associated with snow are cold temperatures, which reduced biodegradation of chemicals [27, 28].
It is also likely that areas with different climates grow different crops that require different pesticides. Refining the crop variable to include different types of crops may help clarify this relationship. Obtaining specific measures on type and quantity of pesticides used in each county would also clarify this association. Currently this information is not collected for all counties in the U.S.
We initially controlled for the effect of the local permitted air emissions and the number of individuals living in each county; however, these variables were never significantly correlated with the incidence of prostate cancer (p values were always greater than 0.185 for all OLS models). It is possible that our measure of air pollution, which was an aggregation of over 350 permitted chemicals, was too crude. Refining the air emissions variable to test individual chemicals or groups of compounds that have similar mechanisms of action and chemical properties may identify specific air contaminants that are problematic for prostate cancer.
Another limitation of this study was the fact that we could not include many known risk factors for prostate cancer because the data were not available at a county level for all of the U.S. Of particular concern to us was the omission of ethnicity and obesity because both of these variables are known to cluster spatially and are associated with prostate cancer [32, 33] thus they had the potential to distort the correlations between the meteorological parameters and the incidence of prostate cancer. Despite not being able to control for these potential confounders directly other parameters such as race, premature mortality from heart disease, and unemployment rates were controlled. Controlling for these other confounders may have inadvertently controlled for the effects of ethnicity and obesity.
Although premature mortality from heart disease and unemployment rates were only included in this study to control for confounding bias, their correlation with prostate cancer is noteworthy and suggests these variables should be included in future models, as both were negatively correlated with the incidence of prostate cancer. We found the more premature heart disease in a county the less prostate cancer there was and, likewise, the more unemployment (lower socioeconomic status) the less prostate cancer in a given county.
Besides controlling for socioeconomic status by including unemployment rate in our model we also tried to minimize the effect of this variable on our outcome parameter (prostate cancer) by using incidence rates instead of mortality rates. We believe mortality rates are influenced by the treatment an individual receives, and this is influenced by their socioeconomic status. The diagnosis of prostate cancer is initially dependant on screening, which is also associated with the individual's socioeconomic status, but presumably an individual with advanced stages of prostate cancer will be diagnosed regardless of their socioeconomic status.
Because all variables were measured and analyzed at the county level in this study, there was potential for bias if individuals in the counties were not actually exposed to the factors included in the model. Given the long incubation period of prostate cancer it is possible that some individuals with this disease migrated between counties during this period, which would not be reflected in our measure of exposure to meteorological parameters. This may have been problematic in our study because older individuals, who are at higher risk of prostate cancer , are more likely to emigrate in one direction: north to south. If anything, this misclassification of individuals would have biased our findings towards the null. It should be noted that to properly control for this type of bias would require conducting studies using data that are collected at the individual level.
Despite the limitations of this study, for example, the fact that some variables were only crude measures of the parameters of interest, that we may have left out some confounding variables, and that all variables were aggregated at the county level, it provides preliminary data suggesting there are correlations between meteorological parameters and prostate cancer. Regardless of the biological parameters included in our models, temperature, shortwave radiation and rainfall were always significant (Table 2). Although it is not possible to determine why meteorological conditions are correlated with prostate cancer using an ecological study, the trends detected in this study are consistent with the literature on environmental chemistry, which suggests meteorological parameters may predispose northern climates to higher levels of pollutants. The transportation and deposition of global sources of POPs to areas with colder temperatures, the reduced efficiency of degradation of compounds in cold dry climates, and the increased volatilization of POPs at low humidity may expose northern regions to higher levels of pollutants. This study, therefore, provides an additional hypothesis for the north-south distribution of prostate cancer, which builds on the existing supposition that individuals at northern latitudes are deficient in Vitamin D because of the low exposure to UV radiation during the winter months. Our study suggests that other meteorological conditions may also significantly affect the incidence of prostate cancer in a county. The findings from this study warrant further investigation using a study design that can more definitively measure the associations between meteorological parameters, and their effects on pollution and prostate cancer.
Methods and Results
We extracted Caucasian average age-adjusted annual incidence rates (cases per 100,000 population per year) of prostate cancer between 2000 and 2004 from the National Cancer Institute (NCI) , for each county in the United States. Analyses were only performed on data for Caucasians (of Hispanic and non-Hispanic origin combined) to control for the effect of race. All data were adjusted to the 2000 U.S. standard population. For six states, including Illinois, Maryland, Minnesota, Mississippi, Tennessee, and Virginia, we obtained data from individual State Cancer Registry websites, as their data were not available through the NCI. For the state of Illinois where data were only available for all races combined, only data from counties where more than 95% of the population was Caucasian were included. We assumed the rates were representative of Caucasians in these cases. We excluded counties with average annual prostate cancer counts of less than 5 from the analysis because stable accurate age-adjusted rates were not available for these counties. The time block used to calculate the average annual age-adjusted incidence rate varied slightly by states (i.e. 2000-2004 or 2001-2005, and in one case 1999-2003); however, it was always an average for a 5 year period. Data indicated a north-south spatial distribution (Figure 1).
Given the increasing frequency of studies reporting associations between different types of pesticides and prostate cancer [22, 22, 24, 25], we felt it was necessary to control for this variable in our models. We used acres of land used to grow crops as a proxy for pesticide use; these data were available through the U.S. Census Bureau. We also acquired population demographics from the Census Bureau  for the 3109 counties in the continental U.S. County level data included total population in 2000, household income for Caucasians in 1999, and annual average unemployment rate between 2000 and 2004. The annual average age-adjusted mortality rate for male Caucasians between 1 and 65 years of age in 2000 and 2004 was acquired through the Centers for Disease Control and Prevention .
Environmental information on average shortwave radiation, average temperature, mean heating degree days (HDD), mean cooling degree days (CDD), average number of frost days, and mean precipitation between 1980 to1997 was obtained from Daymet U.S. Data Center . The spatial reference for these data was defined in ArcGIS (v. 9.3.1) using a projection file provided by the Utah State University Spatial Data Group . Data were then re-projected to allow for the calculation of means by county using zonal statistics.
Average monthly snowfall data from U.S. weather stations for 2000, 2001, 2002, and 2003 were obtained from the National Climatic Data Center . Stations with missing data between the months of October and May were excluded from the analysis. The average snowfall for each year was calculated for all remaining weather stations and subsequently an average was calculated for the 4 years of data. Once this summary statistic was available the weather stations were georeferenced, and using zonal statistics the average snowfall was calculated for the period between 2000 and 2003 for each county with weather station information. We then interpolated the data and estimated snowfall values for the 442 counties with missing data. These counties were randomly distributed in the Eastern and Southeastern United States. Average annual rainfall was calculated by subtracting one tenth of the average annual snowfall (converted from inches/10 to cm) from the mean 18 year average annual precipitation.
Permitted air emissions data for 2002 was obtained from the U.S. Environmental Protection Agency . Emissions were reported for over 350 chemicals. We aggregated the chemicals and determined the sum of permitted emissions for each county using zonal statistics in ArcGIS (v 9.3.1).
Prior to creating models for prostate cancer we screened the variables for correlation because several of the variables of interest measured similar parameters. For example, average temperature, HDD, CDD and mean frost days were well correlated (Pearson r was always greater than 0.89). Given we were most interested in the effect of cold temperature on prostate cancer we chose to include HDD, which is defined as the annual sum of degrees Celsius required to attain 18.3 °C when the air temperature is less than 18.3°C.
Snowfall was positively correlated with HDD (Pearson r = 0.730) and negatively correlated with rainfall (Pearson r = -0.435), and rainfall was negatively correlated with HDD (Pearson r =-0.572). Despite the correlation between these variables we retained all of them for initial testing in our models because they measured different biological parameters.
Once we selected the potential variables to be included in our OLS regression analyses we created several biologically relevant candidate models for explaining the incidence of prostate cancer. These models included various levels of complexity. The first model was similar to what has been published by Schwartz and Hanchette  and included only shortwave radiation. This was used as a comparison for other models. We subsequently added potential confounders such as premature mortality from heart disease and county unemployment rate, as well as higher order terms for shortwave radiation and HDD to account for curvature (Table 2). We also tested a model that included all meteorological parameters, potential confounders, and pollution indices (air emissions, acres of land used for crops, and population). The most extensive model tested included all meteorological variables, pollution indices, confounders, and biologically relevant interaction terms between meteorological parameters and pollution indices. Only variables with p values less than 0.05 were considered significant and maintained in any models.
Subsequently, all candidate models were fitted using GWR analyses . These analyses use information from surrounding counties to build a model where the relationship between the dependent variable and independent variables varies spatially. An adaptive kernel type function using 10% of the U.S. counties as our distance criterion was used for all GWR analyzes. This large distance criteria was required because of the co-linearity between the numerous variables in our model, however the influence of the variables on the outcome was weighted by distance . These models were conducted in Spatial Analysis in Macroecology (v 3.0) . The best fit model in our GWR analyses was determined using the AIC . The residuals from this model were standardized (subtracted from the mean and divided by the number of observations) and mapped in ArcGIS (v 9.3.1).
To clarify the relationship between prostate cancer and the statistically significant interaction term as well as the quadratic terms in the best fit GWR model we used the OLS regression equation from our best fit model and graphed the relationships. These plots were generated by introducing the median value for all parameters except those of interest and determining the incidence of prostate cancer associated with the upper and lower quartile range of values for the parameter(s) of interest (Figures 2, 3, and 4). The figures were generated in Excel (2007 Microsoft® Office Excel® 2007).
Surveillance, Epidemiology and End Results (SEER)Program. [accessed Nov 23, 2009],http://seer.cancer.gov/statfacts/html/prost.html
Rhee van der HJ, de Vries E, Coebergh JW: Does sunlight prevent cancer? A systematic review. Eur J Cancer. 2006, 42: 2222-2232. 10.1016/j.ejca.2006.02.024.
Schwartz GG, Hanchette CL: UV, latitude, and spatial trends in prostate cancer mortality: All sunlight is not the same (United States). Cancer Causes Control. 2006, 17: 1091-1101. 10.1007/s10552-006-0050-6.
Mandal R, St-Hilaire S, Kie JG, Derryberry D: Spatial and temporal correlation between breast and prostate cancers. Int J Health Geographics. 2009, 8: 53-10.1186/1476-072X-8-53.
Boscoe FP, Schymura MJ: Solar ultraviolet-B exposure and cancer incidence and mortality in the United States 1993-2002. BCM Cancer. 2006, 6: 264-10.1186/1471-2407-6-264.
Schwartz GG: Vitamin D and the epidemiology of prostate cancer. Sem Dialysis. 2005, 18: 276-89. 10.1111/j.1525-139X.2005.18403.x.
Schwartz GC, Hulka BS: Is vitamin D a risk factor for prostate cancer? (hypothesis). Anticancer Res. 1990, 10: 1307-1312.
Corder EH, Guess HA, Hulka BS, Friedman GD, Sadler M, Vollmer RT, Lobaugh B, Drezner MK, Vogelman JH, Orentreich N: Vitamin D and prostate cancer: A prediagnostic study with stored sera. Cancer epidemiology, Biomarkers & prevention. 1993, 2: 467-472.
DeLuca HF: Evolution of our understanding of vitamin D. Nutrition reviews. 2008, 66: S73-S87. 10.1111/j.1753-4887.2008.00105.x.
Wania F, Mackay D: Global fractionation and cold condensation of low volatility organochlorine compounds in polar regions. Ambio. 1993, 22: 10-18.
MacDonald RW, Bidleman TF, Diamond ML, Gregor DJ, Semkin RG, Strachan WMJ, Li YF, Wania F, Alaee M, Alexeeva LB, Backus SM, Bailey R, Bewers JM, Gobeil C, Halsall CJ, Harner T, Hoff JT, Jantunen LMM, Lockhart WL, Mackay D, Muir DCG, Pudykiewicz J, Reimer KJ, Smith JN, Stern GA, Schroeder WH, Wagemann R, Yunker MB: Contaminants in the Canadian Arctic: 5 years of progress in understanding sources, occurrence, and pathways. Sci Total Environ. 2000, 254: 93-234. 10.1016/S0048-9697(00)00434-4.
Wania F, Westgate JN: On the mechanisms of mountain cold-trapping of organic chemicals. Environ Sci Technol. 2008, 42: 9092-9098. 10.1021/es8013198.
Franz TP, Eisenreich SJ: Snow scavenging of polychlorinated biphenyls and polycyclic aromatic hydrocarbons in Minnesota. Environ Sci Technol. 1998, 32: 1771-1778. 10.1021/es970601z.
Simonich S, Hites RA: Global distribution of persistent organochlorine compounds. Science. 1995, 269: 1851-1854. 10.1126/science.7569923.
Burrows HD, Canle ML, Santaballa JA, Steenken S: Reaction pathways and mechanisms of photodegradation of pesticides. Journal of photochemistry and photobiology B: Biology. 2002, 67: 71-108. 10.1016/S1011-1344(02)00277-4.
Leech DM, Snyder MT, Wetzel RG: Natural organic matter and sunlight accelerate the degradation of 17 β-estradiol in water. Science of the Total Environment. 2009, 407: 2087-2092. 10.1016/j.scitotenv.2008.11.018.
Wania F: Assessing the potential of persistent organic chemicals for long-range transport and accumulation in polar regions. Environ Sci Technol. 2003, 37: 1344-1351. 10.1021/es026019e.
Juberg DR: An evaluation of endocrine modulators: Implications for human health. Ecotoxicology and Environmental Safety. 2000, 45: 93-105. 10.1006/eesa.1999.1851.
Propper CR: The study of endocrine-disrupting compounds: past approaches and new directions. Integr Comp Biol. 2005, 45: 194-200. 10.1093/icb/45.1.194.
Berg Van den F, Kubiak R, Benjey WG, Majewski MS, Yates S, Reeves GL, Smelt JH, Linden van der AMA: Emissions of pesticides into the air. Water Air Soil Poll. 1999, 115: 195-218. 10.1023/A:1005234329622.
McKinlay R, Plant JA, Bell JNB, Voulvoulis N: Endocrine distrupting pesticides: Implications for risk assessment. Environmental International. 2008, 34: 168-183. 10.1016/j.envint.2007.07.013.
Morrison H, Saitz D, Semenciw R, Hulka B, Mao Y, Morison D, Wigle D: Farming and prostate cancer mortality. Am J Epidemiol. 1993, 137: 270-280.
Gulden van der JW, Kolk JJ, Verbeek AL: Work environment and prostate cancer risk. Prostate. 1995, 27: 250-257. 10.1002/pros.2990270504.
Settimi L, Masina A, Andrion A, Alexson O: Prostate cancer and exposure to pesticides in agricultural settings. Int J Cancer. 2003, 104: 458-461. 10.1002/ijc.10955.
Lynch SM, Mahajan R, Beane Freeman LE, Hoppin JA, Alavanja MCR: Cancer incidence among pesticide applicators exposed to butylate in the agricultural health study (AHS). Environmental Research. 2009, 109: 860-868. 10.1016/j.envres.2009.06.006.
Prins GS: Endocrine distruptors and prostate cancer risk. Endocrine-Related Cancer. 2008, 15: 649-656. 10.1677/ERC-08-0043.
Noyles PD, McElwee MK, Miller HD, Clark BW, VanTiem LA, Walcott KC, Erwin KN, Levin ED: The toxicology of climate change: Environmental contaminants in a warming world. Environmental International. 2008, 35: 971-986.
Sanscartier D, Zeeb B, Koch I, Reimer K: Bioremediation of diesel-contaminated soil by heat and humidified biopile system in cold climates. Cold Regions Science and Technology. 2009, 55: 167-173. 10.1016/j.coldregions.2008.07.004.
Goss K-U, Eisenreich SJ: Sorption of volatile organic compounds to particles from a combustion source at different temperatures and relative humidities. Atmospheric Environment. 1997, 31: 2827-2834. 10.1016/S1352-2310(97)00048-4.
Storey JME, Lou W, Isabelle LM, Pankow JF: Gas/solid partitioning of semivolatile organic compounds to model atmospheric solid surfaces as a function of relative humidity. Environ Sci Technol. 1995, 29: 2420-2428. 10.1021/es00009a039.
Xuan R, Blassengale AA, Wang Q: Degradation of estrogenic hormones in a silt loam soil. J Agric Food Chem. 2008, 56: 9152-9158. 10.1021/jf8016942.
Lopez-Otin C, Diamandis EP: Breast and prostate cancer: an analysis of common epidemiological, genetic, and biochemical features. Endocr Rev. 1998, 19 (14): 365-396. 10.1210/er.19.4.365.
Williams H, Powell IJ: Epidemiology, pathology, and genetics of prostate cancer among African Americans compared with other ethnicities. Methods Mol Biol. 2009, 472: 439-53. full_text.
Crawford ED: Epidemiology of prostate cancer. Urology. 2003, 62 (supplemental 6A): 3-11. 10.1016/j.urology.2003.10.013.
National Cancer Institute State Cancer Profiles. [accessed September 2008],http://statecancerprofiles.cancer.gov/incidencerates/index.php
TIGER/Line® Shapefiles, U.S. Census Bureau, Geography Division, last revised 11/5/2009. accessed August 2009.http://www.census.gov/geo/www/tiger/tgrshp2009/tgrshp2009.html
Center for Disease Control and Prevention. [accessed August 2009],http://wonder.cdc.gov/cmf-icd10.html
DAYMET U.S. Data Center. [accessed August 2009],http://www.daymet.org/climateSummary.jsp
USU Spatial Data Group. [accessed August 2009],http://groups.google.com/group/ususmac/browse_thread/thread/7c34c5000de6d4e0
National Climatic Data Center. [accessed August 2009],http://cdo.ncdc.noaa.gov/pls/plclimprod/poemain.accessrouter
U.S. Environmental Protection Agency 2002 National-Scale Air Toxics Assessment. [accessed August 2009],http://www.epa.gov/ttn/atw/nata2002/tables.html#table1
Brunsdon C, Fotheringham S, Charlton M: Geographically Weighted Regression -- Modelling Spatial Non-stationarity. Journal of the Royal Statistical Society, Series D: The Statistician. 1998, 47: 431-443. 10.1111/1467-9884.00145.
Fernadndo TRLVB, Diniz-Filho JAF, Bini LM: Towards an integrated computational tool for spatial analysis in macroecology and biogeography. Global Ecology and Biogeography. 2006, 15: 321-327. 10.1111/j.1466-822X.2006.00237.x.
Burnham KP, Anderson DR: Model Selection and Multimodel Inference. 2002, New York: Springer-Verlag
This project was funded by an Idaho State University faculty research grant (#LFR021). We would like to thank W. Chalmers and C. Evilia for their editorial comments.
The authors declare that they have no competing interests.
SS provided the idea for the project, assisted with the data interpretation, and helped write the manuscript. RM extracted the cancer and health data. SM and AC extracted and geo-referenced all the environmental data. DD conducted all the statistical analyses and helped with the interpretation of the data. All authors participated in the review and final approval of the manuscript.