Secondary Logo

Mortality Update of a Cohort of Canadian Petroleum Workers

Schnatter, A. Robert, DrPH; Wojcik, Nancy C., MS; Jorgensen, Gail

Journal of Occupational and Environmental Medicine: March 2019 - Volume 61 - Issue 3 - p 225–238
doi: 10.1097/JOM.0000000000001523

Objective: This study updates the mortality experience of over 25,000 workers in a large Canadian petroleum company through December 31, 2006.

Methods: Standardized mortality ratios were generated for all-cause and specific cause mortality.

Results: All cause and all cancer mortality were favorable compared with the general Canadian population. Cancers of previous interest were largely consistent with expectation. There is a continuing excess of mesothelioma, which is of similar magnitude as the previous update, although based on larger numbers. This excess is mostly attributable to men who died in their 50s and 60s and who worked in the refining sector.

Conclusion: Most causes of death show mortality rates lower than the Canadian general population. Given the excess of mesothelioma observed, this study supports ongoing vigilance in asbestos exposure control programs, as refineries continue to remove asbestos from their facilities.

ExxonMobil Biomedical Sciences, Inc., Annandale, New Jersey (Drs Schnatter, Wojcik), and Palisades Research, Inc., Bernardsville, New Jersey (Ms Jorgensen).

Address correspondence to: Nancy C. Wojcik, MS, ExxonMobil Biomedical Sciences, Inc., 1545 U.S. Highway 22 East, Annandale, NJ 08801 (

Clinical Significance: The results of this study provide information about potential chronic health patterns observed among a Canadian petroleum worker cohort compared to the general Canadian population. This information can be used by occupational health care professionals for planning and/or implementing facets of their employee health programs.

A. R. Schnatter was employed by ExxonMobil Biomedical Sciences, Inc. at the time of the study. N. C. Wojcik is currently employed by ExxonMobil Biomedical Sciences, Inc. G. Jorgensen was previously employed by ExxonMobil Biomedical Sciences, Inc. and is a paid consultant.

Supplemental digital contents are available for this article. Direct URL citation appears in the printed text and is provided in the HTML and PDF versions of this article on the journal's Web site (

Funding for this study was provided by Imperial Oil Limited.

The authors report no conflicts of interest.

This is an open access article distributed under the terms of the Creative Commons Attribution-Non Commercial-No Derivatives License 4.0 (CCBY-NC-ND), where it is permissible to download and share the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.

The International Agency for Research on Cancer (IARC) has classified exposures in petroleum refining as probably carcinogenic to humans (Group 2A),1 citing more than one significant elevation for both skin cancer and leukemia. Subsequently, excesses of both of these endpoints continue to be reported2–5 in the refining and petrochemical industry, although not consistently. Other endpoints such as non-Hodgkin lymphoma (NHL) are also under investigation in these workers.6,7 The complex exposure scenario in the petroleum industry includes carcinogens such as benzene, polycyclic aromatic hydrocarbons (PAHs), and asbestos. Thus, to assess potential risks from these and other outcomes, continued scrutiny of petroleum industry workers is warranted through occupational cohort mortality studies.

Several long-term health studies have been conducted to assess potential risks in a large Canadian petroleum cohort.8–13 In general, these studies have shown an overall lower mortality for the cohort than that of the Canadian population, but not for every cause of death. The first cohort study8 examined mortality of refinery workers employed between 1964 and 1973. The study found a two-fold risk of cancer mortality of the intestines (including rectum) and other digestive organs compared with nonrefinery workers.

The cohort was expanded to include all operating segments and updated mortality through 1983.9 This update found statistically significant excesses for malignant melanoma [SMR = 1.87, 95% confidence interval (95% CI) = 1.07 to 3.04] and malignant mesothelioma (SMR = 5.65, 95% CI = 2.27 to 11.63). The melanoma excess was concentrated among upstream workers (SMR = 6.00, 95% CI = 2.19 to 13.06). Analyses of specific substances and hydrocarbon streams did not indicate an association with malignant melanoma. Mortality analysis of marketing/transportation workers showed nearly a two-fold but nonsignificant excess of multiple myeloma (SMR = 1.81, 95% CI = 0.73 to 3.73), which was also related to employment duration, latency, and starting employment before 1950. A subsequent study addressing this finding did not reveal an impact of occupational exposures to benzene or total hydrocarbons.10

The mortality experience of this cohort was again updated, extending the follow-up period through 1994.11 The mesothelioma excess previously reported persisted in this cohort update (SMR = 8.68, 95% CI = 5.51 to 13.03). Most of the mesothelioma decedents were long-term employees with jobs that included presumed exposure to asbestos (mechanical and pipefitters). Deaths from multiple myeloma among marketing and distribution workers remained increased (SMR = 2.08, 95% CI = 0.95 to 3.95) within the 11-year update period, but there was no clear pattern by duration of employment or latency.

In 2003, the mortality experience of an inception cohort (those hired in 1964 to 1994) was analyzed.12 Thus, unlike the previous studies, long-term survivors were excluded and the results were more applicable to modern work environments (previous studies included some workers employed before the 1920s). The size of this inception cohort was 25,292, and was established as the basis for all future updates. In addition, cancer morbidity was also assessed for the first time, using the Canadian Cancer Data Base (CCDB) to identify incident cancers diagnosed among employees between 1969 (the CCDB's start date) and 1994.

Analysis of the inception cohort found that gall bladder cancer mortality was increased among males. However, the increase was based on less than five deaths, none of which had common job assignments, and all had worked fewer than 10 years. Observed mesothelioma incidence in this inception cohort was reduced as expected (less than five cases) but was still greater than the number of mesothelioma cases expected from the general population (1.32). In addition to the overall cohort, exposure-specific analyses were conducted and showed associations between hydrogen sulfide (H2S) and accidental deaths and petroleum coke and lung cancer. The petroleum coke/lung cancer association was further followed with additional smoking data as well as more specific exposure data, including asbestos exposure.13 These results showed a stronger correlation between lung cancer and employment in the maintenance trades, rather than coke dust or asbestos, suggesting that multiple exposures experienced by workers in maintenance jobs were more predictive of lung cancer. It was anticipated that future updates of this inception cohort, with additional years of follow-up (viz, the present study), would be increasingly useful in monitoring the disease trends of more recent employees, especially for mesothelioma.

This present study extends monitoring of the mortality experience of employees who started work from 1964 through 2006. Mortality patterns for specific causes of death (viz, lymphohematopoietic (LH) cancers, lung cancer, mesothelioma, malignant melanoma, kidney cancer, gall bladder cancer and angiosarcoma of the liver) are the focus, to determine whether occupational factors influence these causes, and if so, the characterization and extent of the occupational factors’ influence. Twelve additional years of mortality follow-up increase the statistical power to detect elevated health risks, if they are present, or likewise reveal the absence of such risks within the study population. Specific exposures were not captured in this update due to transitions in the exposure tracking software system. Benzene exposure in this population was previously studied in detail using a case–control approach with individual-based exposure assessment.14 In the current study, the total work environment, which includes exposure to several hydrocarbon streams, H2S, benzene, polycyclic aromatic compounds, etc, was the focus. Any excesses found are assessed for occupational influences by duration of employment and latency.

Back to Top | Article Outline


Research Design and Strategy

The design of this study is that of a retrospective cohort. The cohort is followed for mortality and rates of various diseases are compared with the rates expected in a comparable portion of the general Canadian population. Several health outcomes are of a priori interest based on the earlier studies’ findings and/or on the epidemiology literature for petroleum workers. For statistical analyses, these causes of death are evaluated by gender, operating segment of the petroleum industry, and broad job classification. Any suggestive findings are further analyzed by typical occupational influences such as duration of employment and latency (eg, elapsed time between first hire and cancer death).

Back to Top | Article Outline

Data Sources and Definitions

Overall Cohort Source and Definition

The present cohort added 4087 employees, bringing the total cohort size to 29,379 workers. The cohort includes all employees from the previous investigation12 as well as all employees hired between January 1, 1995, and December 31, 2006, who worked for at least a year.

Identification of the overall cohort was based upon the previous study's database and was supplemented by additional computerized work records from the company HR databases. The personnel and payroll systems are integral to the company's function, and are therefore a complete source for identification of company employees and their company-specific work assignments. The analysis was further restricted to employees with a minimum of 1-year employment.

Back to Top | Article Outline

Mortality Status Determination and Definition

The mortality experience of the cohort was updated through December 31, 2006; the last date for complete ascertainment when the study was initiated, using mortality tracing information from the Statistics Canada Mortality Database. Statistics Canada (SC) or provincial registrars obtained death certificates and the underlying cause of death. Underlying cause of death was determined by SC via the International Classification of Disease (ICD) classification scheme according to the revision in effect at the time of death. This assured comparability with the mortality rates used to compute the expected number of deaths for different cause of death categories. Records for employees were also matched to the U.S. National Death Index (NDI) to identify deaths that may have occurred among cohort members who may have transferred or immigrated to the United States.

For angiosarcoma of the liver and mesothelioma of the pleura and peritoneum, previous versions of the ICD scheme did not have a unique code. Thus, we requested that SC review all death certificates for mention of these causes. The specific ICD codes that may relate to both of these diagnoses have been published previously.9 SC identified five mesotheliomas and less than five angiosarcomas via these manual searches.

Records not matched to either the SC mortality database nor the NDI, and not actively working or receiving benefits at the end of the study period (n = 6171) were potentially lost to follow-up (LTF). However, a further “alive tracing” by SC via Tax Summary files from 1984 to 2007 confirmed that only 3564 were LTF as of 2006. Privacy legislation in Canada prevented the identification of these individuals; hence, 3564 of the 6171 workers were selected using a stratified random sampling approach (see statistical analysis section).

Back to Top | Article Outline

Demographic and Work History Sources and Definition

All demographic data (eg, age, gender) for employees were obtained from the company's Human Resource payroll and personnel computerized systems. This information is continuously validated by the company. An operating segment was defined as a part of the industry where different operations are present (eg, “refining” produces products from crude oil, “upstream” consists of the exploration, drilling and production of crude oil or gas, etc). Operating segments were documented well in each employee's work history, and consisted of refinery, agricultural chemicals, petrochemical, marketing, transportation, marine, exploration/drilling/production (upstream), pipeline, coal and minerals, office, and fabricated building products subsidiary. This scheme was developed by three experienced company industrial hygienists. Employees were classified as “ever employed” in an operating segment, such that an employee could be defined in more than one operating segment. In addition, the “Exposure Group” (EG) field, which was linked to the work history (98% linked) from the company's Exposure Tracking and Health Information System (ETHIS),12 was used in analyses. The EG field (hence: job group) describes the jobs among employees who perform similar work tasks and contains the following entries: operations, maintenance, transportation, lab/research, as well as other smaller groups. The job group field was also used to identify employees who were nonexposed or minimally exposed to company operations. For all nonexposed job groups, employees were only classified as such if they never worked in an exposed job group. For exposed job groups, employees were classified as “ever employed”; thus, an employee could be represented in more than one exposed job group, although overlap is minimal.

Back to Top | Article Outline

Data Quality

Data quality assurance was built into designated phases of the project. Company staff performed quality checks on fields abstracted from the payroll and personnel systems. These included out-of-range checks as well as consistency checks between fields. All dates were tested for completeness and logical progressions. Frequency counts were obtained before and after data were received from SC to identify any transcription errors. Within SC, procedures for matching potential deaths and cancer cases to records of cohort members are well developed, and use the Generalized Record Linkage System (GRLS).

Back to Top | Article Outline

Statistical Analysis

Eligible employees were entered into follow-up for mortality 1 year after their date of hire within the cohort entry requirements discussed above. Person-years at risk were accumulated by age, gender, and within quinquennia until the earlier of date of death or the end of study (December 31, 2006). Persons of unknown vital status (LTF) were assumed to be alive until their date last observed (DLO), which was usually date last employed. As the 3564 LTFs were selected from the 6171 original LTF (before additional SC “alive tracing”), we ensured that these workers had the same distribution for DLO as the original 6171 LTF employees by randomly selecting the 3564 LTF within 5-year DLO strata for the 6171 workers. We also performed sensitivity analyses on different LTF strategies to examine the impact of these assumptions.

Expected numbers of cancer cases and deaths were computed using a modified life table approach with the software package OCMAP-Plus, version 3.10 (University of Pittsburgh, Pittsburgh, PA). Expected deaths were calculated by multiplying person-years by mortality rates for the general Canadian population according to age, sex, and calendar time. Mortality rates were obtained from SC databases. Observed deaths were divided by expected deaths to generate standardized mortality ratios (SMRs). Males and females in each of these employment categories were analyzed separately. SMRs and 95% CIs were computed for several causes of death and diagnostic codes. To limit presentation of less meaningful results, SMRs were computed when there are at least five observed outcomes. Although tabulated data for operating segments and job groups represent “ever worked” and are not mutually exclusive, the observed numbers of death for any cause were tabulated only if there were zero deaths or when they reached five or more deaths. This is consistent with disclosure controls for mortality data obtained from the Canadian Vital Statistics Death Database. These disclosure rules were also applied to other tabulated data in this report.

Back to Top | Article Outline


A total of 29,379 workers contributed 702,906 person-years at risk between 1964 and 2006, for an average of about 24 years of follow-up per employee. The distribution of persons by gender, vital status, and other attributes is displayed in Table 1. The mean age of the cohort at the beginning and end of follow-up was 27 and 54 years, respectively. At the end of follow-up, 75% of the cohort was at least age 47, 50% of the cohort was age 54, 25% of the cohort was at least age 60, and 5% of the cohort was at least age 71. Compared with the initial publication of this inception cohort, this population consists of 4087 (16%) additional workers and 309,094 (78%) additional person-years, while average follow-up time increased from 17 to 27 years.



Back to Top | Article Outline


A total of 1669 observed deaths occurred compared with 2418.5 expected deaths, for an all-cause SMR of 0.69 (95% CI = 0.66 to 0.72) (Table 2 ). All-cause SMRs were similar among men, SMR of 0.70 (95% CI = 0.66 to 0.74) and women, SMR of 0.66 (95% CI = 0.59 to 0.74). For the total cohort, nearly all major categories of death were significantly reduced compared with national background rates (Table 2 ). The only significant excess was for mesothelioma, where nine deaths were observed, based on 10th revision coding, compared with 2.8 expected for a significantly raised SMR of 3.27 (95% CI = 1.50 to 6.21).





The 2003 update of this cohort reported a statistical excess of gall bladder cancer based on less than five deaths.12 For this update, no additional gall bladder cancer deaths were found, while there are now a total of 4.9 expected cases (Table 2 ).

Findings for other a priori causes of death were generally unremarkable for the total cohort (Table 2 ). Cancers of the colon (large intestine), lung (bronchus and lung), and kidney showed lower observed versus expected numbers of death. Melanoma and leukemias also showed reduced SMRs.

There was a moderate excess of amyotrophic lateral sclerosis (ALS) based on 15 observed deaths and 8.8 expected deaths for men and women combined, though it was not statistically significant (SMR = 1.70, 95% CI = 0.95 to 2.80) (Table 2 ).

For aortic aneurysm (AA), and the more inclusive category of “Diseases of the arteries, arterioles and capillaries” (DAAC), the SMRs for men versus women were in the opposite direction. Women showed nonsignificant excesses of both diseases (SMR for DAAC = 1.89, 95% CI = 0.87 to 3.59; SMR for AA = 2.48, 95% CI = 0.80 to 5.78). However, men showed deficits for both diseases (SMR for DAAC = 0.44, 95% CI = 0.22 to 0.80; SMR for AA = 0.62, 95% CI = 0.28 to 1.18). There were no other notable differences in male versus female mortality for the entire cohort.

Back to Top | Article Outline

Mortality by Operating Segment

The excess of mesothelioma mortality was largely due to an excess in refinery workers (six observed, SMR = 9.93, 95% CI = 3.64 to 21.60) (Table 3 ). Most a priori cancers (viz. melanoma, leukemia, colon cancer) had observed numbers of death consistent with expectation in refinery workers. The SMR for cancers of unknown sites was 1.60 based on 11 cases but was not statistically significant.







There were no statistically significant increased SMRs in petrochemical workers (Table 3 ). When refinery and petrochemical workers were combined, the findings noted above in refinery workers did not increase in magnitude. Of note, lung cancer showed 40 versus 59 deaths in these two segments, which together showed slight lung cancer excesses within normal statistical variation as previously reported.13

There were no statistically significant excess mortality findings in upstream workers (Table 3 ), including for malignant melanoma, which was elevated for the larger older cohort. A previous observation of an increased incidence of motor vehicle accidents was found to be similar to expected numbers in this current update.

Transportation workers showed a nonsignificant elevation due to ALS, n = 5, SMR = 2.15, 95% CI = 0.70 to 5.03 (Table 3 ).

Back to Top | Article Outline

Mortality by Job Group

Job groups with sufficient numbers included maintenance, operations, and distribution, and are summarized in Table 4 .





As expected, mesothelioma deaths were more evident in maintenance and operations jobs. Workers classified as “operations” had no other notable findings for mortality, with an overall SMR of 0.72 (95% CI = 0.65 to 0.80). On the contrary, workers in the maintenance category showed a slightly higher overall cancer mortality rate than expected [95 observed vs 86.5 expected, SMR = 1.10 (95% CI = 0.89 to 1.34)]. However, the total SMR was only 0.77 (95% CI = 0.67 to 0.88). A nonstatistically significant excess of bladder cancer (5 observed vs 1.7 expected) was observed in maintenance workers.

Distribution jobs showed an overall SMR of 1.00; thus, the expected “healthy worker effect” was not displayed in these workers. This was due to small but notable excesses in the most common causes of death. The SMR for all circulatory diseases was 1.19, while moderate excesses that were not typical of healthy employed populations were present for cerebrovascular disease: 6 versus 4.3, SMR = 1.39 (95% CI = 0.51 to 3.02) and acute myocardial infarction: 20 versus 15.7, SMR = 1.27 (95% CI = 0.78 to 1.96). However, neither SMR was statistically significant. Also, noteworthy was the SMR for motor vehicle accidents in these jobs (SMR = 1.87, 95% CI = 0.86 to 3.55), many of whom are drivers.

Back to Top | Article Outline

Mortality by Exposure Status

Table 5 summarizes SMRs according to a general indicator of exposure to any operations present across jobs and operating segments. Persons classified as nonexposed in Table 5 were never in an exposed job (ie, exclusively nonexposed). A worker who had any portion of their work history in any exposed job was classified as exposed.



There were still significant deficits in all-cause (SMR = 0.74) and all-cancer (SMR = 0.84) in exposed workers, but each SMR was higher than that in nonexposed workers (SMRs = 0.62 and 0.73 for all causes and all cancers, respectively). Mesothelioma was significantly elevated in exposed workers (SMR = 3.53, 95% CI = 1.30 to 7.69), while other nonsignificant elevations were seen for bladder cancer (SMR = 1.25, 95% CI = 0.60 to 2.29), rectal cancer (SMR = 1.12, 95% CI = 0.60 to 1.91), and air and space transport accidents (SMR = 1.61, 95% CI = 0.59 to 3.49). Notably, leukemia, colon cancer, NHL, and lung cancer were all below expectation. There were five AML cases resulting in an SMR of 1.00.

SMRs in nonexposed workers were generally unremarkable, with the exception of 10 ALS cases, which was a significant excess (SMR = 2.48, 95% CI = 1.19 to 4.56). Other nonsignificant elevations were observed for Hodgkin disease based on five cases (SMR = 2.42, 95% CI = 0.79 to 5.64) and follicular lymphoma based on nine cases (SMR = 1.30, 95% CI = 0.60 to 2.47).

Back to Top | Article Outline

Mortality by Hire Date

In order to assess time trends for different diseases, we also examined the cohort by hire date (1964 to 1974) and 1975+ (see Table S-1, Supplemental Digital Content 1, These two periods resulted in a similar number of expected deaths. Overall mortality SMRs decreased from 0.70 to 0.68 for the earlier versus later higher groups, likely to be a reflection of the healthy worker effect wearing off, as most of the decrease was due to circulatory diseases (SMR = 0.70 and 0.62) for earlier and later hire periods, respectively. Malignant neoplasms showed a relatively constant SMR over the earlier and later hire periods (SMR = 0.79 and 0.77, respectively). Mesothelioma risk had a significantly raised SMR of 5.14, based on six deaths in the 1975+ hire group. Some a priori’ causes of death showed a modestly raised (nonsignificant) SMR in the earlier time period only, including follicular NHL (SMR = 1.27 based on 11 cases), AML (SMR = 1.36 based on seven cases), and multiple myeloma/immunoproliferative diseases (SMR = 1.14 based on eight cases). In addition, all central nervous system tumors (the vast majority being brain cancer) showed a higher SMR in the early hire period (SMR = 1.24 based on 25 cases). In contrast, a few causes of death showed a moderately raised, nonsignificant SMR only in the later hire period, including malignant neoplasms of the rectum and rectosigmoid junction (SMR = 1.15 based on 10 cases) and ovary (SMR = 1.62, based on nine cases).

Back to Top | Article Outline


The tenth revision of the ICD, in effect from 2000, designates a single code for mesothelioma; thus, mortality and incidence reporting are less complicated. As excess mesothelioma has been noted previously in this cohort, we also examined death certificates that could contain the diagnosis in previous versions of the ICD, as we did in the previous updates of this cohort.9,11,12 This yielded five additional cases of mesothelioma, four of which were noted in the last update. In addition, the cohort was matched to cancer registry files maintained by SC. This yielded a sixth additional case of mesothelioma. These six cases, along with the nine cases identified from 2000 forward, yielded 15 cases of mesothelioma. Canadian rates assembled in the same manner as case ascertainment techniques yielded an expected number of 5.4 cases for a standardized incidence ratio (SIR) of 2.79 (95% CI = 1.56 to 4.60) (Table 6).15 Subtracting the cases observed in the previously studied observation period yields 11 cases observed versus 4.12 expected for an SIR of 2.67. As previously noted, the mortality-only cases observed from year 1999 forward yielded nine cases and an SMR of 3.27 (95% CI = 1.50 to 6.21).



These 15 cases were further reviewed. As for jobs and duration employed, 10 were in maintenance or operations jobs and worked from 11 to 32 years. As mesothelioma often develops decades after first exposure, the time between death/incidence and first employment was examined. Less than five subjects had a 10 to 14-year interval. Most had either a 15 to 24-year or 25 to 38-year interval. All subjects died or were diagnosed in their 50s and 60s; 10 at ages 50 to 59 years, and five at ages 61 to 63 years. Eleven of the 15 cases were employed in the refinery operating segment. There was only 1.0 expected case in the refinery segment, for a significant SIR of 10.9 (95% CI = 5.45 to 19.5).

We also examined SMRs for mesothelioma by age at hire, as employees hired later may have had exposure in other industries (see Table 6). Of the nine mortality cases, five were hired between ages 25 and 35 years (SMR = 4.89, 95% CI = 1.59 to 11.41). This analysis indicates that mesothelioma cases had some potential for exposures before employment in this company.

Back to Top | Article Outline

Lost to Follow-up

We needed to use a stratified random sampling technique to select 3564 lost to follow-up (LTF) employees. SC had determined that this was the correct number of LTF, but privacy laws prevented their direct identification. Instead, we examined DLO status in the original 6171 LTF employees and selected the 3564 employees to reflect this DLO distribution. We suspected that this strategy would yield the most similar overall SMR to an SMR calculated using the DLO for the 6171 employees. To examine this directly, we computed overall (all causes) and cancer SMRs for two additional scenarios, which we believed represented two extremes: (1) using the original 6171 workers as LTF and removing them from follow-up on their DLO (to produce a minimum number of person-years at risk, hence a lower expected number of deaths and higher SMR), and (2) treating all LTF workers as alive until the end of the study (to produce a maximum number of person-years at risk, hence a larger expected number of deaths and lower SMR). Strategy (1) could be justified, as it represents the usual treatment of LTF, while strategy (2) could be justified by the high sensitivity of the Canadian Mortality Database used by SC to identify deaths.16 Results are summarized in Table 7.



These results show that the all-cause SMR for the baseline scenario (3564 LTF) is within 1.4% of the all-cause SMR, which treats all 6171 employees as LTF on their DLO. On the contrary, treating all LTF employees as alive until the end of the study would have reduced the SMR by 7.7% (from 0.690 to 0.637). Cancer SMRs show similar percent changes.

Finally, as the sensitivity of CMDB has been tested for its ability to identify deaths, one could also assume that this sensitivity (97.6%) applies to this cohort. This would predict that 1710, rather than 1669 deaths occurred. Applying the 1710 deaths to strategy (2) would yield an SMR of 0.653, a 5.4% lower SMR than the baseline approach.

Back to Top | Article Outline


This mortality update for a large Canadian cohort of petroleum workers was based on 16% more workers and 78% additional person-years of observation, allowing a better assessment of risks in this population. However, this is still a young population; thus, the accuracy of risk estimation, especially for rarer causes of death, is still relatively low. The percentage of deaths increased from 2% to 5.7% in this update. We believe it is advantageous to study inception cohorts (in this case, workers hired in 1964 or later) to improve the applicability to current workplace settings and to lessen the possibility of survivor bias influences on the results.17

As in many occupational cohort studies, mortality rates among working employees were compared with general population mortality rates; thus, the well-known healthy worker effect is probable. As this cohort is young, the healthy worker effect is greater, as more workers were close to their “selection date,” which assured a health status sufficient to apply for and be selected for employment. Thus, even moderate, nonsignificant excesses should be viewed more critically for these workers.

While effects of gender and time period were controlled, we had no additional information on lifestyle habits such as drinking and smoking, which can affect mortality from several diseases. The study is also limited by its reliance on surrogates of exposure (eg, duration of employment, operating segment, job type, and general exposure to operations), and as such, excesses restricted to employees only exposed to select compounds may be missed.

Another limitation of this study is the potential LTF for 3564 workers. This is a somewhat larger percentage than usual. We not only used the Canadian Mortality Database, but also used the US NDI to identify deaths to account for people emigrating from Canada to the US who subsequently died. We also benefited from alive tracing performed at SC, which did reduce the number of LTF, but was somewhat compromised by an inability to directly identify which employees were successfully traced. Instead, we needed to select remaining LTF employees utilizing a stratified random sampling approach. We performed sensitivity analyses on this aspect of the study and showed that our estimate of the total SMR was likely to be bounded by an error from -1.4% to +7.7%. We proffer two reasons why the LTF percent is somewhat higher than normal: (1) there could be an increase in movement of workers between this Canadian affiliate and other affiliates from the same company outside of the US, and (2) there could be an increase in the number of short-term workers, as younger persons are frequently more mobile than in years past. Subsequent updates of this cohort will need to employ even more rigorous tracing techniques to limit LTF to more customary levels.

Cancer incidence data were requested for this study but were found to be unreliable, apparently due to the multiple ICD revisions these data cover. We hope to resolve this with Statistics Canada so that future analyses can examine incident cancers.

The mortality findings in the cohort overall continue to show favorable results, with an all-cause SMR of 0.69, slightly higher than the 0.65 measured in the previous update.12 We expect that the all-cause SMR will continue to increase, as the survival advantage of initial employment selection wears off. Decreased mortality was primarily due to respiratory and circulatory diseases (especially acute myocardial infarction), which are commonly noted as most affected by the healthy worker selection effect. There were also very low SMRs for acquired immune deficiency syndrome and infectious and parasitic diseases.

For most a priori causes, results were not significantly elevated. Leukemia, NHL, and multiple myeloma showed unremarkable results, which is in agreement with most of the literature.18 Malignant melanoma, which had been elevated in upstream workers in an earlier update of this cohort,9 showed no excess risk in this update. Kidney, colon, and brain cancer were also not elevated for the full cohort or any subgroup examined.

In the last update, we reported an excess of gall bladder cancer based on less than five cases.12 The present update found no new gall bladder cancer cases, while the expected number increased to 4.9.

The single finding that showed a clear excess in risk was for mesothelioma. With an additional follow-up of this inception cohort, the mesothelioma risk in workers hired post-1964 is clearer (Table 6). Overall, the present cohort found 11 additional cases for an overall SIR of 2.7915 through 2006.

This two to three-fold risk of mesothelioma continues to persist despite continued attention in the company on asbestos exposure control. The risk is higher for maintenance and operations workers, and in the refining sector, where there is more potential for asbestos exposure. Approximately one-third of the mesothelioma cases had latency periods from first employment between 25 and 38 years, which translates to first exposures between 1968 and 1981. Although the risk of mesothelioma was expected to decrease by 2006, we did note that most decedents started work with the company at an age which suggests that they were previously employed elsewhere. If they were employed in prior jobs with asbestos exposure, the latency observed in this cohort could be artificially shortened.

A two-fold nonsignificant elevation of ALS was associated with a previous update of this study.11 This update again reports a nonsignificant elevation of ALS, based on 15 cases (SMR = 1.70), which did not quite reach statistical significance. The etiology of ALS is largely unknown.19,20 It appears that there is no clustering of ALS risk in this study by operating segment. A significant excess of ALS was found in nonexposed workers, which argues against an occupational etiology.

Other findings, some of which did not achieve statistical significance, should not be ignored, as the healthy worker effect is likely operating in this population. A near significant SMR for motor vehicle accidents among distribution workers likely reflects increased time in motor vehicles. The fact that maintenance workers show an all-cancer SMR greater than 1 is also unusual in an employed population. Bladder cancer rates should be monitored in this group, as they are subject to unique work environments (eg, welding, degreasing, etc).

In summary, this study yielded several useful insights on risks for this large Canadian cohort of petroleum workers. The clearest finding was a continuing excess of mesothelioma, which was predicted to decrease in this update. The excess concentrated in refinery workers and maintenance trades, who have a higher potential for asbestos exposure. A nonsignificant excess of ALS was also noted, which achieved significance only in workers not exposed to operations. A nonsignificant excess of motor vehicle traffic accidents in distribution workers may have an occupational link that should be monitored. Continued surveillance of this cohort is expected to yield more insight on these findings.

Back to Top | Article Outline


1. International Agency for Research on Cancer (IARC). Monograph on the evaluation of carcinogenic risks to humans: occupational exposures in petroleum refining, crude oil, and major petroleum fuels. IARC Monogr 1989; 45:39–117.
2. Gun RT, Pratt N, Ryan P, et al. Update of mortality and cancer incidence in the Australian petroleum industry cohort. Occup Environ Med 2006; 63:476–487.
3. Sorahan T. Mortality of UK oil refinery and petroleum distribution workers, 1951-2003. Occup Med 2007; 57:177–185.
4. Wongsrichanalai C, Delzell E, Cole P. Mortality from leukemia and other diseases among workers at a petroleum refinery. J Occup Med 1989; 31:106–111.
5. Huebner WW, Wojcik NC, Rosamilia K, et al. Mortality updates (1970-1997) of two refinery/petrochemical plant cohorts at Baton Rouge, Louisiana, and Baytown. J Occup Environ Med 2004; 46:1229–1245.
6. Steinmaus C, Smith AH, Jones RM, et al. Meta-analysis of benzene exposure and non-Hodgkin lymphoma: biases could mask an important association. Occup Environ Med 2008; 65:371–378.
7. Alexander DD, Wagner ME. Benzene exposure and non-Hodgkin lymphoma: a meta-analysis of epidemiologic studies. J Occup Environ Med 2010; 52:169–189.
8. Hanis NM, Stavraky KM, Fowler JL. Cancer mortality in oil refinery workers. J Occup Med 1979; 21:167–174.
9. Schnatter AR, Theriault G, Katz AM, Thompson FS, Donaleski D, Murray N. A retrospective mortality study within operating segments of a petroleum company. Am J Ind Med 1992; 22:209–229.
10. Schnatter AR, Armstrong TW, Nicolich MJ, et al. Lymphohaematopoietic malignancies and quantitative estimates of exposure to benzene in Canadian petroleum distribution workers. Occup Environ Med 1996; 53:773–781.
11. Lewis RJ, Schnatter AR, Katz AM, et al. Updated mortality among diverse operating segments of a petroleum company. Occup Environ Med 2000; 57:595–604.
12. Lewis RJ, Schnatter AR, Drummond I, et al. Mortality and cancer morbidity in a cohort of Canadian petroleum workers. Occup Environ Med 2003; 60:918–928.
13. Schnatter AR, Nicolich MJ, Lewis RJ, et al. Lung cancer incidence in Canadian petroleum workers. Occup Environ Med 2012; 69:877–882.
14. Schnatter AR, Glass DC, Tang G, Irons RD, Rushton L. Myelodysplastic syndrome and benzene exposure among petroleum workers: an international pooled analysis. J Natl Cancer Inst 2012; 104:1724–1737.
15. Dean AG, Sullivan KM, Soe MN. OpenEpi: Open Source Epidemiologic Statistics for Public Health. Version 2.3. 1; 2011. Available at: Updated April 06, 2013. Accessed June 05, 2017. 2013.
16. Schnatter AR, Acquavella JF, Thompson FS, Donaleski D, Theriault G. An analysis of death ascertainment and follow-up through Statistics Canada's mortality data base system. Can J Public Health 1990; 81:60–65.
17. Saracci R. Survival-related biases survive well. Int J Epidemiol 2007; 36:244–246.
18. Schnatter AR, DeVilbiss EA, Chen M, Lewis RJ, Gallagher EM. Systematic review and meta-analysis of selected cancers in petroleum refinery workers. J Occup Environ Med 2018; 60:329–342.
19. Weisskopf MG, McCullough ML, Morozova N, Calle EE, Thun MJ, Ascherio A. Prospective study of occupational and amyotrophic lateral sclerosis mortality. Am J Epidemiol 2005; 162:1146–1152.
20. Ingre C, Roos PM, Piehl F, Kamel F, Fang F. Risk factors for amyotrophic lateral sclerosis. Clin Epidemiol 2015; 7:181–193.

Canada; Canadian; cause of death; mesothelioma; mortality; occupational cohort; petroleum workers

Supplemental Digital Content

Back to Top | Article Outline
Copyright © 2019 by the American College of Occupational and Environmental Medicine