Article Text
Abstract
Objectives Increased mortality from asthma, chronic bronchitis and emphysema has previously been reported among workers in the silicon carbide (SiC) industry. The objective of the present study was to evaluate the influence of specific exposure factors on mortality from obstructive lung diseases (OLD), using a newly revised job-exposure matrix.
Materials and methods 1687 long-term workers employed in 1913–2003 in the Norwegian SiC industry were characterised with respect to cumulative exposure to quartz, cristobalite, SiC particles and SiC fibres. Standardised mortality ratios (SMRs) for underlying causes of death, 1951–2007, were calculated stratified by category of cumulative exposure, and Poisson regression analyses of OLD were performed using cumulative exposure variables.
Results An increased total mortality (SMR 1.1, 95% CI 1.0 to 1.2) and increased mortality from cancer, non-malignant respiratory diseases and external factors, were observed. The SMR of OLD was increased at the highest level of cumulative exposure to all investigated exposure factors. In the internal analyses, a twofold increased risk of OLD was observed with increasing levels of cumulative exposure to SiC particles. In a multivariate model, SiC particles showed the most stable increased risk estimate when controlled for other exposure factors, among workers with less than 15 years of employment. Among workers with more than 15 years of employment, crystalline silica, primarily cristobalite, seemed to be the most important exposure factor.
Conclusion Exposure to SiC and crystalline silica may contribute to OLD development among SiC industry workers in different time windows, and possibly through different mechanisms.
- Epidemiology
- respiratory
- mortality studies
- inorganic dusts
- inorganic fibres
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.
Statistics from Altmetric.com
What this paper adds
Workers in the silicon carbide industry have increased mortality from non-malignant respiratory diseases.
The specific causal factors of this increase have not previously been known.
The present study shows that several exposure factors, including silicon carbide dust and crystalline silica, may contribute to the increased mortality.
Dust exposure control in the silicon carbide industry is essential.
Introduction
The question of adverse health effects among workers in the silicon carbide (SiC) industry was already being addressed in the 1920s. Winslow et al described an excess of tuberculosis among the workers, and pointed to dust in the working atmosphere as a possible cause.1 In 1949, Bruusgaard showed that workers in the SiC industry had a high prevalence of silicosis, and that many of the workers who had developed silicosis seemed not to have been exposed to crystalline silica but only to SiC dust.2 Two mortality studies have been performed in this industry. In 1994, Infante-Rivard et al found a significantly increased risk of non-malignant respiratory diseases (standardised mortality ratio (SMR) 2.0, 95% CI 1.2 to 3.2) among workers in the Canadian SiC industry.3 The risk increased with increasing exposure to total dust. In 2002, Romundstad et al published a mortality study in the Norwegian SiC industry, with similar findings—an increased risk of total mortality (SMR 1.1, 95% CI 1.1 to 1.2) and mortality from non-malignant respiratory diseases (SMR 1.4, 95% CI 1.1 to 1.7).4 Specifically, the combined risk of asthma, chronic bronchitis and emphysema was elevated, with an SMR of 2.2 (95% CI 1.6 to 3.0), and an increased risk with increasing level of exposure to total dust was also observed.4 However, the exposure data in these two studies were limited.
In the present follow-up of the Norwegian cohort, a revised and extended, retrospective job-exposure matrix (JEM) has been developed, based on more than 8500 samples of total dust and dust fractions (Føreland et al, 2011, submitted). This paper presents an update of the previous Norwegian mortality study,4 with an additional 11 years of follow-up and improved exposure estimates from the revised JEM. As the risk of non-malignant pulmonary diseases was increased in both previous mortality studies,3 4 we will address this issue specifically, and study the association of mortality from obstructive lung disease (OLD) with exposure to fibres and crystalline components in the dust. The effect of long-term cumulative exposure is also addressed.
Materials and methods
SiC production
SiC is produced by the Acheson method.5 A mixture of finely ground quartz sand and petroleum coke, often with additives such as sawdust and aluminium oxide, is placed in open furnaces. The mixture is electrically heated by a graphite core to over 2500°C close to the core and about 1000°C on the periphery. In the region close to the core, silicon dioxide is in the gas phase and reacts with coke to form SiC and carbon monoxide. The thermal process lasts about 48 h, and at the end the zone closest to the core consists of a thick layer of loosely knitted SiC crystals, whereas the periphery consists of unreacted material. The zone between those two layers contains partially reacted material, and here fibrous SiC is found.6 The unreacted and partially reacted material is separated from the SiC fraction, and is reused as raw material. Through several crushing and washing processes the crude SiC is refined into different grain size fractions, according to desired end use. The production of SiC is energy intensive, and Norway, with its abundant supply of hydro electric power, got its first SiC plant in 1913. In the 1960s, two additional plants were established.
Cohort design and follow-up
The study population was based on a previously established cohort of men employed in 1942 and later, with ≥6 months of total employment in the Norwegian SiC industry between 1913 and 2003, and alive after 1 January 1953, giving 2631 individuals in total.7 As data on causes of death were available from 1 January 1951, up to two extra person-years of follow-up were added for all workers employed before 1953 and alive between 1951 and 1953. In addition, six workers who died between 1951 and 1953 were added to the cohort. The information available for this extended cohort of 2637 men included name, date of birth, employment periods with job codes, and smoking status (ever/never-smoker and unknown). For workers alive in 1960 or born later, the Norwegian unique 11-digit personal identification number was also available.
Data on causes of death were available from 1 January 1951. For the period 1960–2007, the cohort was linked, via the 11-digit identification number, with the Norwegian National Population Register and the Norwegian Causes of Death Register to obtain dates of death and emigration, and causes of death. To obtain causes of death before 1960, linkage with the individual's name and date of birth was performed. The Regional Committee for Medical Research Ethics and the Norwegian Data Inspectorate endorsed the study.
Eight workers registered as dead in the Norwegian National Population Register had no death diagnoses in the Norwegian Causes of Death Register, and were thus excluded from the cohort. A recent cancer study in the Norwegian SiC industry found different risk profiles between short- and long-term workers, probably due to a combination of occupational and lifestyle factors.7 It was therefore decided to base the present follow-up on workers with a total employment of ≥3 years, giving 1687 men in total.
The diagnoses from the Causes of Death Register were coded in International Classification of Diseases (ICD)-6 through ICD-10, depending on the date of diagnosis. Due to changes in the coding system, and changes in diagnostic trends concerning OLD during the years,8 we chose to pool diagnoses related to asthma, chronic bronchitis, emphysema and chronic obstructive pulmonary disease (COPD) into OLD.9 In the Norwegian Causes of Death Register, both underlying and contributing causes of death are collected from the death certificate. The underlying cause of death is a disease or injury that initiated the train of events leading to death. A contributory cause of death is a condition that contributes to the fatal outcome, but which is not related to the disease or condition directly causing death.
The follow-up of mortality among the long-term workers started after 3 years of total employment duration, or from 1 January 1951 if the 3 years of employment duration was reached before that date, and lasted until 31 December 2007 or the time of death or emigration. The long-term workers contributed 42 689 person-years at risk to the study.
The JEM
A historical JEM was developed for all job groups in each of the three plants, for all years since the start of production. The development of the JEM was based partly on previous dust measurements, performed both for surveillance by the three plants and in connection with studies performed by the Norwegian National Institute of Occupational Health. Altogether more than 3300 total dust measurements from 1967 to 2005 were available, but few measurements (about 300) of specific agents. Hence, the modelling of historical exposure levels was based on the total dust measurements. To predict exposure to specific dust components, a comprehensive comparative exposure assessment study was performed in connection with the present project during 2002–2003, and included 600–700 samples of total dust collected in parallel with fibres and respirable dust. The respirable dust samples were analysed for quartz, cristobalite and non-fibrous SiC content. This study added almost 5000 measurements to the exposure data.10 Mixed effect models were constructed to establish the relationships between total dust and specific components, and these models were used to compute historical job type-specific exposure estimates of the dust components. Information from the plants about technological changes and changes in working hours was used to estimate relative changes in total dust exposure during the earliest years of production, when no measurement data were available. Based on the assumption (with a few exceptions) that the composition of the dust did not change over time, group mean exposures to the specific components could be estimated from total dust, in about 80 job groups, distributed over the three Norwegian SiC plants, for each year of production in the period 1913–2003. The exposure estimation process is described in detail elsewhere (Føreland et al, 2011, submitted). High correlation between exposure to the different agents in the previous JEM made it difficult to separate potential effects posed by different types of exposures.4 In the revised JEM the correlation coefficients were moderate to low except between total dust and respirable dust, providing the potential for better differentiation of any harmful effects from exposure factors.
Through linkage via job code and year between the cohort database and the JEM, we were able to estimate each person's time weighted average exposure to each investigated dust component during each year of employment in the SiC industry. The individual annual exposure was summed, giving estimates of cumulative exposure to total and respirable dust, respirable quartz, cristobalite and SiC particles, and SiC fibres, for each person-year of follow-up.
Departments and job groups
The main departments of SiC plants are as follows. (i) The furnace department, where workers are exposed to raw materials (coke and quartz), agents formed during the furnace process (cristobalite (a high temperature crystal phase of silica) and SiC fibres) and the end product (non-fibrous SiC). In addition, other exposure factors are likely to be present: organic compounds such as polycyclic aromatic hydrocarbons and graphite, carbon monoxide and sulphur dioxide gases, and in earlier periods, asbestos. (ii) The processing department, where the SiC is refined into size fractions. Here, non-fibrous SiC is the main exposure factor, but agents from the partly reacted material contaminate the crude SiC and are also present. (iii) The maintenance department, which mainly consists of electricians and mechanics, working in both the furnace and processing departments. (iv) Other workers, including low exposed personnel, for example laboratory workers and stockhouse workers. (v) Office workers. As the two latter groups are quite small, we have combined them in the analyses.
For the analyses stratified by department, we chose to allocate department affiliation to those who according to our knowledge had worked in one department only. Workers with employment in more than one department, or with unknown affiliation, were allocated to a ‘mixed’ group. The years of employment in this mixed group were distributed almost equally in furnace (2810 years) and processing (2675 years) departments, in addition to 2181 years in maintenance and 1305 years in ‘other, low exposed’ jobs.
Figure 1 shows the mean total dust levels, per 10-year period, experienced by the long-term workers in the present cohort, by department.
The cumulative exposures to total and respirable dust, respirable SiC, quartz and cristobalite, and fibres among workers in the different departments are shown in table 1. Arithmetic means are estimated from geometric means as described by Seixas et al.11
Data analysis
SMR were calculated as the ratio between observed and expected numbers of deaths. Expected numbers were based on national mortality rates, calculated from national incidences of underlying causes of death and population numbers (available for 1951–2006), summed in 5-year age and period groups. For 2007 we used the same national rates as for 2006, as an approximation. Cumulative exposure to all six exposure factors was categorised in three groups (low, medium and high) ensuring equal numbers of person-years of follow-up in each category. We performed SMR analyses of the main diagnostic groups overall, and of OLD stratified by department and by category of cumulative exposure. We calculated 95% CIs assuming a Poisson distribution of the observed numbers.
The correlations between exposure variables were examined using Pearson's correlation coefficients between the log-transformed cumulative exposures.
In order to study possible exposure–response associations, we performed Poisson regression analyses of OLD related to categories of exposure, adjusted for age (dichotomised at 75 years, due to the high age of those dying with OLD), smoking (never, ever and unknown) and period of diagnosis (before/after 1990). In the internal analyses we, as an alternative to OLD as underlying cause of death (OLD-u), used OLD as underlying or contributing cause of death (OLD-uc), to increase the number of cases and thereby the power of the analyses.12 Tests of trend were done by fitting a Poisson regression model for the linear effect of the exposure variables, assuming a constant increase in the log rate ratio per exposure category level.13 Confounding between the four specific exposure factors (respirable quartz, cristobalite and SiC, and fibres) was studied by modelling Poisson regression analyses with two exposure variables at a time. Continuous cumulative exposure to each exposure factor was log-transformed, due to skewed distribution of the exposure. In order to test the null hypothesis that there was no effect of the second exposure variable after controlling for the first, likelihood ratio testing was performed comparing the model containing one exposure factor with a model containing two exposure factors. These analyses were done among ever-smokers only, stratified by total employment duration (</≥15 years) in order to study the effect of long-term cumulative exposure, and adjusted for age and period of diagnosis. In these analyses, follow-up of the group with ≥15 years of employment started after 15 years of total employment.
All statistical analyses were performed using Stata 11.
Results
Cohort characteristics and the distribution of deaths due to OLD for the 1687 long-term workers are presented in table 2.
We observed both an increased total mortality risk (SMR 1.1, 95% CI 1.0 to 1.2; 788 cases) and increased risks of cancer (SMR 1.2, 95% CI 1.0 to 1.4; 201 cases), respiratory diseases (SMR 1.6, 95% CI 1.3 to 2.0; 91 cases) and external causes (SMR 1.5, 95% CI 1.1 to 2.0; 44 cases). In addition, excess risks were found for the respiratory subgroups OLD (SMR 2.0, 95% CI 1.5 to 2.7; 45 cases), pneumonia (SMR 1.4, 95% CI 1.0 to 1.9; 38 cases) and pneumoconiosis (SMR 15, 95% CI 7.0 to 31; 7 cases). For circulatory diseases the SMR was 1.0 (95% CI 0.9 to 1.2; 347 cases), for digestive diseases 1.1 (95% CI 0.7 to 1.7; 17 cases) and for other diagnoses 0.9 (95% CI 0.7 to 1.1; 88 cases). (Not shown in table.)
The increased risk of mortality from OLD was linked to employment in the two production departments: furnace (SMR 2.7, 95% CI 1.6 to 4.6; 13 cases) and processing (SMR 2.2, 95% CI 1.2 to 4.3; 9 cases). There was no increase in mortality from OLD among workers with affiliation to the maintenance (SMR 1.0, 95% CI 0.4 to 2.7; 4 cases) or ‘other, low-exposed’ (SMR 0.4, 95% CI 0.1 to 2.6; 1 case) departments. Workers with employment in more than one department (SMR 2.5, 95% CI 1.6 to 4.0; 18 cases) had a risk almost as high as furnace workers. (Not shown in table.)
Table 3 presents the analyses of OLD related to categories of cumulative exposure to total dust, respirable dust, respirable SiC, quartz and cristobalite, and fibres. The SMR was significantly increased at the two higher levels of exposure to total dust, cristobalite, SiC dust and respirable dust, at all levels of exposure to fibres, and at the lowest and highest level of exposure to quartz. Internal analyses in models of ever-smokers using OLD as underlying cause of death showed significantly increased risks at the two higher levels of exposure to SiC (incidence rate ratio (IRR) 3.1 and 3.4, respectively), with a significant trend (ptrend=0.03). Using OLD as underlying or contributing cause of death, the IRR estimates for SiC were somewhat reduced, still with a significant trend (ptrend=0.02). This last analysis also showed a significant exposure–response trend for total dust (ptrend=0.03) and cristobalite (ptrend=0.05). In analyses stratified by total employment duration in the SiC industry, there was a significant exposure–response trend for SiC dust among workers with less than 15 years of employment. In contrast, among workers with more than 15 years of employment, there were significant exposure–response trends for total dust and cristobalite, and a borderline significant trend for quartz.
Table 4 shows the results from analyses including two log-transformed cumulative exposure variables in models of ever-smokers only. Correlations between cumulative exposure to the investigated specific exposure variables were at a moderate level (≤0.75), except the correlations between quartz and fibre exposure (CorrPearson 0.76–0.83). In the subcohort of ever-smokers, SiC was the exposure factor with the highest risk estimate when studied as an isolated log-transformed exposure factor. This risk estimate remained at the same level when fibres, quartz or cristobalite were included in the model. Similarly, among workers with less than 15 years of employment, the risk estimate was highest for SiC, and was not changed when other exposure factors were included in the model. SiC was, however, highly correlated to total dust (CorrPearson 0.88–0.93), and by including total dust in the model, the effect of SiC was reduced, both among all ever-smokers (IRR 1.1) and among workers with less than 15 years of employment (IRR 1.5; not shown in the table). Among workers with more than 15 years of employment, risk estimates were highest and most stable for the two crystalline silica agents, cristobalite and quartz. When including both cristobalite and quartz in the model, the effect of quartz was reduced, and cristobalite seemed to be the most important factor. The inclusion of total dust did not change these estimates.
These analyses were performed using OLD as underlying or contributing cause of death (OLD-uc). Analyses using OLD as underlying cause of death (OLD-u) gave the same results overall, but with less precision (results not shown).
Discussion
In a 57-year follow-up of mortality among 1687 long-term workers in the Norwegian SiC industry, we have found increased risks of mortality from all causes, cancer, external factors and non-malignant respiratory diseases, including OLD. SMR analyses of OLD stratified by category of cumulative dust exposure showed significantly increased SMR with higher levels of exposure to all investigated exposure factors. Internal analyses indicated that SiC was the exposure factor with the highest risk estimate among workers with less than 15 years of employment, and cristobalite seemed to be the most important factor among workers with more than 15 years of employment.
In the previous mortality study in the Norwegian SiC industry, an excess mortality from asthma, emphysema and chronic bronchitis combined was found, with an exposure–response association to total dust exposure. High correlations between exposure to the various agents specified in the JEM precluded differentiation between the potential effects of these agents.4 The present study confirmed the main results from the previous study, and the revised and extended JEM made it possible to better distinguish between the effects of different exposure factors.
The improved JEM used in this study provided estimates of exposure to specific exposure factors in more than 80 job groups during almost 100 years of industrial activity. Cumulative exposures to total dust, respirable dust and components in the respirable fraction (quartz, cristobalite and SiC dust, and fibres) were estimated and used in analyses of exposure related mortality from OLD. However, other exposures in the SiC industry that could be important in this connection, for example, sulphur dioxide and carbonaceous materials such as coke, graphite and polycyclic aromatic hydrocarbons, were not measured in sufficient detail to be included in the JEM. These exposure factors are mainly present in the furnace department. In the processing department, SiC is the most abundant exposure factor and represents 60–80% of the respirable dust. In the furnace department, SiC represents 20–40% of the respirable dust.10
A priori we applied categories of cumulative exposure by equal number of person-years in each group in our analyses. Post hoc, we have tested categories based on equal numbers of OLD cases in each group. The results are quite similar. Another alternative could be to choose cut-points from the cumulative exposure data, based on presumed low, medium and high risk levels. This was not attempted because we had no such presuppositions, and arbitrarily chosen cut-points would give a greater risk of collapse in the analyses because of uneven distribution of the cases.
Development of the JEM used in the present study was based on dust measurements performed since 1967. The exposure estimates for the period before 1967 were based on information about changes in processes and working hours. Hence, the exposure estimates before 1967 were less reliable than the later estimates. However, the exposure estimates were developed independently of the outcomes from the epidemiological study, giving a non-differential misclassification of exposure, which in most cases will give a bias towards the null.14
The cohort used in this study has a cross-sectional start in 1942, with start of follow-up in 1951. Thus, the oldest part of the cohort consists of survivors from the period 1913–1951. As exposure intensities were higher in earlier times, one would expect that work related death was also higher in earlier times. In the Norwegian SiC industry, many silicosis cases were reported in the 1940s, and there is reason to believe that quite a few workers died from work related lung diseases, primarily silicosis, with and without tuberculosis, in the 10–20 years before 1950.1 2 This would result in a healthy survivor population at start of follow-up in 1951, probably giving an underestimation of work related mortality in the present study. The fact that we observe only two cases of OLD-u (three cases of OLD-uc) in the period 1951–1969 (table 2) may also point in the direction of a healthy survivor cohort. On the other hand, we do not have any reason to believe that the left-censoring caused an overestimation of the risk. We have also performed analyses with an inception subcohort starting in 1951, and the overall trends are the same as presented in table 3. Further subdivision of the cohort, at 15 years of employment, results in small groups, with fragile and unreliable results (results not shown).
One interesting finding in this study was the contrast between exposure–response associations among workers with less than and more than 15 years of employment in the SiC industry. Occupational exposure to minerals and fibres can elicit symptom-producing responses in the airways, often clinically interpreted as OLD,15 through different mechanisms, depending on the characteristics of the dust, for example, chemical composition, particle size and surface properties.16 17 Other risk factors for OLD are familiar disposition, smoking and earlier infections.18 Subjects with a higher risk of OLD from a combination of causes19 20 may have an increased risk of early airway symptoms in a dusty work environment. This could, eventually, lead to selection out of exposed work.21 The toxicity of SiC dust has previously been studied in in vitro and in vivo experiments, and it has been classified as ‘inert’.22 23 Our finding of an association between SiC and OLD mortality among workers who left the industry before 15 years of employment, and also the close correlation between SiC and total dust may, however, be consistent with SiC being part of a nuisance dust load which may induce OLD symptoms among disposed individuals. The limitation of this association to workers with less than 15 years of employment may be suggestive of a healthy survivor effect.12
Crystalline silica, on the other hand, is present in lower concentrations, and represents only a minor fraction of the total dust exposure in the SiC industry. Crystalline silica will, however, accumulate in the lung,24 and may induce lung injury, with symptoms associated with OLD that do not appear until many years after first exposure.25 This could explain our findings of exposure–response associations with cristobalite among the workers with more than 15 years of employment. Quartz seemed to be a less important exposure factor than cristobalite. Previous studies have indicated that cristobalite is a more potent fibrogen than quartz.26 27 The results of our study suggest that cristobalite is also more potent than quartz when OLD is the outcome. When we summed the exposure to cristobalite and quartz, generating a log-transformed cumulative crystalline silica exposure variable, we found a stronger exposure–response association (IRR 3.1, 95% CI 1.5 to 6.3) than with cristobalite alone, indicating an additive effect of quartz and cristobalite.
As smoking is the most important cause of OLD, adjustment for smoking is necessary in studies of other aetiological factors. Among the long-term workers in the present cohort, 69% were ever-smokers, and only seven of the cases (11%) were among never-smokers or people with unknown smoking status. The small numbers of cases in the latter two groups caused a collapse in the multivariate analyses with adjustment for smoking, and we therefore chose to perform these analyses in the subcohort of ever-smokers. Our findings may still be biased by a possible correlation between smoking and degree of dust exposure, or even an interaction between these factors, but we do not have the information needed to address this question.
Mortality from pneumoconiosis, which is exclusively a work related disease, was increased to a high level (SMR 15, 95% CI 7.0 to 31). It would have been of interest to study exposure–response associations with pneumoconiosis using the JEM, but the numbers were too small to obtain meaningful results (seven cases). All those who died from pneumoconiosis had a diagnosis of silicosis, had worked in the SiC industry for more than 20 years before 1975, and had first employment before 1950. They had been exposed in different departments, including the furnace, processing and maintenance sections.
The mortality from pneumonia was significantly increased with 38 deaths and an SMR of 1.4. A review of the contributing causes of death gave no indication of work relatedness—only one person had COPD as contributing cause of death, and many had dementia and prostate cancer, associated with high age (mean age at death: 84 years). We do not believe that this increased incidence of pneumonia was related to occupational exposure in the SiC industry.
Conclusions
We have in the present study performed analyses on associations between OLD and cumulative exposure to certain dust species among workers in the SiC industry. Mortality from OLD was increased compared to the general population, and significantly increased risks were seen at the highest level of cumulative exposure to all investigated exposure factors. Exposure–response relationships were found with SiC and crystalline silica. However, other unmeasured components in the furnace hall, such as carbonaceous dust and sulphur dioxide, may also have contributed to the observed excess of OLD. Different exposure factors may contribute to OLD development among SiC workers in different time-windows, possibly through different mechanisms.
Acknowledgments
The authors thank Pål Romundstad and Berit Bakke for commenting on the manuscript draft.
References
Footnotes
Funding The project has been financed with the aid of EXTRA funds from the Norwegian Foundation for Health and Rehabilitation, and with support from the Ministry of Labour and Social Inclusion.
Competing interests None.
Ethics approval This study was conducted with the approval of the Regional Committee for Medical Research Ethics South-East Norway, Oslo University, Norway.
Provenance and peer review Not commissioned; externally peer reviewed.