Probabilistic human health risk assessment of 1,3


HomeHome / Blog / Probabilistic human health risk assessment of 1,3

Sep 06, 2023

Probabilistic human health risk assessment of 1,3

Scientific Reports volume 12,

Scientific Reports volume 12, Article number: 22103 (2022) Cite this article

713 Accesses

1 Citations

1 Altmetric

Metrics details

Chemicals containing Volatile Organic Compounds (VOCs) are commonly used in the machine carpet production. 1,3-butadiene and styrene are main components of the carpenter's glue used in carpet factories. Exposition to these chemicals can lead to a number of adverse health effects. This is the first study of the human health risk assessment due to inhalational exposure to 1,3-butadiene (BD) and styrene (ST) performed among workers in the carpet factories in Kashan city, Iran. The importance of the study was related with the fact of high popularity of carpet production in the South Asia countries. Inhalation exposure to BD and ST were measured based on the National Institute for Occupational Safety and Health (NIOSH) 1024 and 1501 methods, respectively. The cancerogenic risk (CR) and non-cancerogenic risk described as Hazard Quotient (HQ) values were calculated based on the United States Environmental Protection Agency (USEPA) method. The sensitivity and uncertainty analysis were performed by the Monte Carlo simulation (MCS) technique. The average concentration measured of BD and ST during work shifts of employees were 0.039 mg m−3 (0.017 ppm) and 12.108 mg m−3 (2.84 ppm), respectively. The mean ± SD value of estimated cancerogenic risk in inhalation exposure to BD and ST were equal to 5.13 × 10–3 ± 3.85 × 10–4 and 1.44 × 10–3 ± 2.36 × 10–4, respectively exceeding the acceptable risk level of 10–6 defined by USEPA. The average non-carcinogenic risk (HQ) values of BD and ST were equal to 8.50 × 100 and 5.13 × 100, respectively exceeding the acceptable risk level of 1. As the results of our studies exceeded both cancerogenic and non-carcinogenic risk values it indicates that adverse health effects due to inhalational exposure to BD and ST for workers in the machine carpet industry are very likely. To avoid negative health effects protective measures for employees in the factories should be introduced immediately and furher detailed research are recommended.

The development of industry, in addition to improving and increasing the level of well-being of human life1, has also created numerous problems for people at the global and regional levels2. Today, the use of chemicals in human life is inevitable. The use of chemicals in many aspects of life and economic activities has brought significant benefits and has changed the quality of human life1. On the other hand, these chemicals can be problematic for human health and the environment1,2. Air quality in workplaces has been one of the main concerns following the growth of industries in recent decades3. Poor air quality at workplaces has been linked to the presence of several hazardous chemicals, including the volatile organic compounds (VOCs) in the air4. Some VOCs, such as benzene, 1,3-butadiene, and styrene, are toxic and might evoke adverse health effects5. Irritation of the eyes, nose, throat, and lungs as well as damage to the liver, kidneys, and central nervous system are acute effects of VOCs’ exposure4. Asthma, respiratory symptoms, cardiovascular disease and different types of occupational cancers are chronic effects of VOCs’ exposure5. Extensive human exposure to these compounds in various industries around the world is a major issue of concern for human health6.

1,3-butadiene (BD) is a synthetic, colorless gas with the formula (CH2=CH)2. BD has been categorized as a human carcinogen by the International Agency for Research on Cancer (IARC) and also, the United States Environmental Protection Agency (USEPA) since 2019 consider BD as a high-priority for risk assessment7. Safety Data Sheet8 indicate that BD is an extremely flammable gas, may cause genetic effects and cancer in inhalational exposure, suspected of damaging fertility or the fetus. Exposure controls for humans regarding the standards for BD were set to be as follows8: ACGIH OEL TWA 2 ppm, OSHA PEL TWA 1 ppm, OSHA PEL STEL 5 ppm, NIOSH ILDH 2000 ppm; where ACGIH is American Conference of Governmental Industrial Hygienists, OSHA is Occupational Safety and Health Administration, NIOSH is National Institute for Occupational Safety and Health, OEL means occupational exposure limit, PEL means permissible exposure limit, TWA means time weighted average, STEL means short term exposure limit, ILDH means immediately dangerous to life or health. Based on the results of a cohort study, associations between inhalation exposure to BD and leukemia and bladder cancer, were found7. Diseases of cardiovascular system were described as the chronic effects of the exposure to BD9. Some metabolites of BD, such as monoepoxide, diepoxide, and epoxydiol, are suspected of causing DNA damage. Diepoxybutane (DEB), that is the most important metabolite of BD, causes adverse effects on DNA by generating reactive oxygen species (ROS) and 8-hydroxydeoxyguanosine (8-OHdG). Irritation of the eyes, nasal passage, and respiratory system as well as fatigue, with great effects on blood pressure, heart rate and damage of central nervous system have been reported as short term effects of BD6,9.

Styrene is a chemical with the formula of C8H8. ST is an aromatic hydrocarbon derived from benzene and has a sweet odor. Safety Data Sheet10 indicate that ST is a flammable liquid and vapor, causes skin irritation and serious eye irritation, is harmful if inhaled, may cause respiratory irritation and drowsiness or dizziness, is suspected of damaging the fetus, is suspected of causing cancer. Defined target organs due to exposure to BD are respiratory system, ears, and central nervous system11. Exposure controls for humans regarding the standards for ST were as follows10: ACGIH TLV TWA 10 ppm, ACGIH TLV STEL 20 ppm, OSHA PEL TWA 50 ppm, OSHA PEL STEL 100 ppm, NIOSH ILDH 700 ppm; where ACGIH is American Conference of Governmental Industrial Hygienists, OSHA is Occupational Safety and Health Administration, NIOSH is National Institute for Occupational Safety and Health, OEL means occupational exposure limit, PEL means permissible exposure limit, TWA means time weighted average, STEL means short term exposure limit, ILDH means immediately dangerous to life or health. Occupational exposure to ST affects negatively the human health, including effects on the peripheral and central nervous system (with symptoms of drowsiness, headache, and imbalance), respiratory system, and liver damage12. Based on the evaluation, absorption of ST is immediate via skin contact and through the lung, mainly disseminated in adipose tissue, and extensively metabolized in the body12. The results of Ruder et al. conducted in plastic injection industry have presented an increase in the occurrence of leukemia and lymphoma due to exposure of workers to ST13. Styrene has been designated as Group 2A (probably carcinogenic to humans) by the IARC in 201914 and also, US National Toxicology Program (US NTP) considered ST as a carcinogenic substance15.

In spite of its toxic properties, both chemicals are used mainly as monomers to produce various types of polymers and copolymers such as styrene-butadiene copolymer and as a chemical intermediate in the manufacturing of some industrial chemicals16. The global BD market volume is over 12 million tons per year17. Carpet manufacturing workers are exposed to BD and ST18,19.

Regarding workers safety, US Occupational Safety and Health Administration (OSHA) limits exposure to BD to not more than 1 ppm (2.21 mg m−3) for eight hours or 5 ppm (11 mg m−3) for 15 min20. US OSHA determined the permissible exposure level (PEL) for ST to be 100 ppm and American Conference of Governmental Industrial Hygienists (ACGIH)—the Threshold Limit Value (TLV) to be 20 ppm21. The main exposure pathway for the BD and ST is inhalation due to volatile properties of these chemicals. One of the most reliable methods for measuring exposure is direct measurement of the concentration of the contaminant in a person's respiratory zone9. Risk assessment is described as determining the potential unfavorable health effects of environmental hazards22, as well as a tool to determine the risks of hazards in the workplace by considering the existing control measures and deciding whether to accept or not23.

Quantitative estimate of carcinogenic risk and non-carcinogenic risk from inhalation exposure to VOCs were developed by various agencies, such as USEPA5. The quantitative method proposed by the United States Environmental Protection Agency (USEPA) is an important and common method in the field of chemical risk assessment6. In this method, to determine the level of carcinogenic risk of exposure to chemical compounds, the cancerogenic risk (CR) index is used6. The use of quantitative risk assessment methods is considered by many international organizations, including the World Health Organization (WHO)6 and the USEPA23, as the basis for legislation on chemical compounds24.

Numerous studies have investigated the health risk of exposure to harmful chemical compounds in various industries like petrochemical6 and oil refinery in Iran25, leather, wooden furniture, printing, dyeing, garment manufacturing in China24, and among hospital workers in Malaysia26 (Table 1). Wani and Jaiswal (2012) reported that carpet weaving in Kashmir, India was related with various health hazards, not only various chemical compounds but also dust27Subedi and Banamala in their research revealed that carpet factory workers in Nepal in more than 50% were young women and for 44% of the carpet factory workers their families worked in the same occupation28. Due to the harmful work conditions in India Wani et al. 2015 pointed the necessity to introduce some provisions to protect workers’ health like masks, earplugs, first aid facilities, and proper personal protective equipment29. Literature search revealed that more investigations were carried on the ready carpet products and their impact on customers due to carpeted floors in their offices, schools, houses, etc.30.

In addition, other research indicated that removal of these pollutants from the environment required sophisticated technologies, like biofiltration34,35, due to the fact the air pollutants do not occur solely and the removal method has to be sufficient for variety of chemicals at the same time36. Thus, based on the literature review, we stated that no study was conducted before to assess the health risk of BD and ST in carpet factories.

As all chemicals present in the air are potentially hazardous to the environment the pressure is put to control their emissions37 .This becomes particularly important during the occupational exposure, when chemicals concentrations and exposure time might be much higher. In some countries these permissible or recommended levels of VOCs in industry might not be in force to protect the employees and some countries might not have restrictive occupational law protecting employees at all33. Due to the above, in our studies we performed Human Health Risk Assessment (HHRA) based on the USEPA methodology of the inhalational exposure to evaluate the general level of risk and to recognize if the further actions will be needed regarding the health safety of investigated subpopulation of workers in the machine carpet industry. The novelty of our research was performing the HHRA based on the measured concentrations of contaminants in the indoor air on the workplaces in the carpet industry factories. As the carpet industry is very popular in South Asia countries, like Iran, this was important to perform this preliminary research, especially that best to our knowledge, such research was not performed before. Moreover, the traditional carpet production process is not performed in the factory shops with modern workers protection measures but more often in family houses in India, China, Turkey, Iran, and Pakistan, with women and children participation38. In our research we have used the USEPA method for calculating health risk as it can be used for various exposure scenarios, not only occupational as is the case of occupational safety agencies. In the HHRA developed by the USEPA the conservative risk assessment principle is used. It means that the slightest adverse health effects are in request. Occupational risk assessment assumes high doses and intensive exposure for workers. In the carpet production industry these conditions not always have to be as such. Workers in the carpet production in South Asia often are not typical employees as whole families might produce them in ordinary buildings and without personal protection measures. Such conditions are more like environmental exposure, which in the USEPA method after adjusting exposure scenario is also used for determining the occupational risk.Thus, based on our preliminary research in this study the same methodology can be further used for health risk assessment investigations for other susceptible subpopulations, like (pregnant) women, children, and customers in the future.

Thus, the objective of our studies was the assessment of the cancerogenic and non-carcinogenic risk for workers of machine carpet industry in Iran due to the inhalational exposure to 1,3-butadiene and styrene using USEPA methodology. The detailed objectives were as follows: (1) to determine the concentrations of BD and ST at workstations in machine carpet industry factories; (2) to investigate exposure levels, health risks, and related uncertainty due to BD and ST presence, and (3) to determine the main exposure factors affecting the risk value. The findings from the present study will provide a baseline data and scientific support needed for environmental pollution management in machine carpet production factories.

This study was conducted in Kashan city, Iran in 2022 to evaluate the health risk related with the inhalational exposure to BD and ST present in the air on the workstation in the machine carpet factories. Based on our previous study33, technological process description, and experts knowledge it was indicated that exposure to BD and ST occured in the finishing shops of carpet factories. A total of 75 employees working in finishing shops of carpet factories were included in the questionaire surveys. The research was performed in compliance with ethical principles of the Declaration of Helsinki and they were approved by the Research Ethics Committee of Kashan University of Medical Sciences, Iran (No 180IR.KAUMS.NUHEPM.REC.1401.004) and informed consent was obtained from all subjects and/or their legal guardian(s). Face-to-face interviews were performed with the employees for collecting data on demographics, past diseases, and occupational history.

Machine-woven carpets woven by carpet loom, in their raw form, are not suitable for supply to the market and delivery to the customer. The presence of dead pile threads floating behind some carpets, the weakness of the roots of the pile thread and the possibility of them coming out of the carpet, the unevenness of the carpet surface, the ugliness of the sides and tarpaulins of the woven carpet, the looseness of the warp yarns and the weft of the stalks and the possibility of them coming out of the carpet texture, the filthy and dirty surface of the carpet and the presence of defects caused by different stages of weaving are among the most important defects that can be observed in raw (unfinished) carpet. To eliminate any of the above-described defects, a finishing process consisting of several stages must be provided during the machine carpet production. As chemicals containing VOCs are used at each stage of the machine carpet finishing process, vapors of BD and ST affect workers when the glue containing styrene-butadiene copolymer is poured and discharged into the glue injection boiler, and when due to the heat the vapors from the glue pan of the sizing part, glued carpets, and from the dryer are released to employees’ workstations19.

Air samples of BD and ST were collected from the worker's breathing zone in two factories located in the Kashan city, Iran. The concentrations of BD and ST were determined in winter and during eight hours work shift (8:00 to 16:00) based on the optimized NIOSH 10246 and 150131 methods. Altogether 247 (3 samples for each worker and 22 blanks) indoor air samples were collected during different times in the shift (beginning, middle and the end) and these results were averaged for further calculations in the study. The samples were pumped by activated coconut charcoal (front (400 mg) and rear (200 mg) for BD and front (100 mg) and rear (50 mg) for ST) made by SKC company which was connected to the worker collar (approximately in the breathing zone) using a pump at the recommended flow rate of 0.2 L min−1. Personal sampling pump of model AirChek TOUCH (5–5000 mL min−1, SKC, Inc.) was used. The time of sampling was adjusted from 70 to 120 min based on the pretest which had been carried out to control the breakthrough volume. In the next step, in order to prevent sample loss during transfer to laboratory, both sides of the adsorbent tubes were sealed with plastic caps and placed in a cool box.

The collected samples of BD and ST were moved to extraction vials. Desorption was carried out using 4-ml methylene chloride (99.95%) (Merck Inc., Germany) and CS2 for BD and ST, respectively. The samples of BD and ST were subjected to ultrasonic waves for 30 min in order to complete extraction. One microliter (1 μl) of the extracted sample was measured using GC–MS (7890 gas chromatograph, and 5975 mass spectrometer, Agilent technologies, CA, USA). Helium was applied as a carrier gas, at a flow rate of 1 mL min−1.

In the present study, the allowable occupational exposure limits for BD and ST vapors were calculated to be 1 ppm (2.21 mg m−3) and 100 ppm (425 mg m−3), respectively, based on the values reported by US OSHA. Given that the amount of TLV-TWA provided is assuming 8 h of work per day and 5 days of work per week. In the cases when working hours exceeded 40 h per week, the amount of TLV-TWA were modified using the Brief & Scala correction model39, which model is used for adjusting non usual work schedules and considers longer workdays and shorter recovery period.

In the present study, the quantitative risk assessment method proposed by the USEPA was used. In this method, to estimate the carcinogenic risk of exposure to BD and ST vapors, the cancerogenic risk (CR) index was used. The value of the index for BD and ST compositions during the present study were calculated using Eq. (1)6:

where CDI is chronic daily intake (mg kg−1 day−1) and SF is slope factor (mg kg−1 day−1)−1. As no Inhalational Unit Risk (IUR) values were available for ST, instead of exposure concentration, we used chronic daily intake (CDI) to have the unified calculating methodology as for our chemicals slope factor (SF) values for both BD and ST were available. SF is an acceptable range in which there is a potential for a response per unit of chemical exposure over a lifetime. SF values for each carcinogenic compound are provided from toxicological databases. In our calculations we used SF values equal to 6 × 100 (mg kg−1 day−1)−1 for BD40 and 5.7 × 10–4 (mg kg−1 day−1)−1 for ST41. Chronic daily intake (CDI) is the dose of particular pollutant taken daily averaged over exposure expressed in years. CDI values in the present study were calculated using Eq. (2)22:

where C = concentration of pollutant (mg m−3), IR = inhalation rate (m3 day−1), ED is exposure duration (years), EF is exposure frequency (days year−1), BW is body weight (kg), AT is averaging time (days). The exposure and toxicological parameters used in the study are presented in Table 2.

According to USEPA guidelines31 in our study the acceptable carcinogenic risk was set to be equal to 1 × 10–6 (one additional cancer risk in 1,000,000 population).

According to the USEPA methodology for calculating the non-carcinogenic risk assessment the Hazard Quotient (HQ) is used. Hazard quotient (HQ) is the ratio between actual exposure to pollutant and its reference concentration (RfC). RfC expresses the continuous inhalation exposure concentration that will not provide adverse health effects during a lifetime. The HQ values were calculated based on the Eq. (3):

where EC is the exposure concentration, RfC reference concentration. The target non-carcinogenic risk value was set to be equal to 1, meaning that HQ values ≥ 1 indicate nonacceptable risk levels43.

Exposure concentration is the daily exposure to pollutant in the inhalational route of exposure. The exposure concentration was calculated using the Eq. 422:

where EC is the exposure concentration (mg m−3), ET is exposure time (hours day−1), ED is exposure duration (years), EF is exposure frequency (days year−1), AT is averaging time (ED in years × 365 days/year × 24 h/day in hours). For values of the parameters see Table 2.

Given that human health is followed by some uncertainties, neglecting these uncertainties may result in the loss of important information. Therefore, unrealistic and incorrect decisions can be made regarding the protection of human health. Monte Carlo simulation (MCS) is a theory based on probabilistic and statistical-mathematical approach applied to study uncertainty using random sampling of each parameter. This technique can reduce uncertainty. By determining statistical indicators or identifying its distribution function, the uncertainty degree of the output variable can be described. The general structure for determining uncertainty by the Monte Carlo technique is a combination of simulations. The calculations were conducted with 10,000 iterations, which eventually yielded the results with a degree of confidence in the range of 1 to 99%44. In the present study, Crystal Ball software (version, Oracle, Inc., USA) was used.

In the present study, blank samples were tested in field sampling and laboratory analysis to check the levels of contamination and potential errors in the stages of sampling, transfer and measurement. The results showed that the concentration of compounds in each blank sample was less than 1% of the values measured in the original samples. The concentration values in the blank samples were also subtracted from the values found in the main samples. To determine the recovery rate of the analytes and to determine the accuracy of the measurements, the spike sample test was used. The results of the measurement accuracy showed that the percentages of analytes recovery were 87% for BD and 89% for ST.

All of the subjects had full consent to participate in the current study.

Cancerogenic and non-carcinogenic human health risk assessment was conducted in our studies due to inhalational individual exposure of carpet factories’ employees to BD and ST. Best to our knowledge, the present study was the first attempt to assess the cancerogenic and non-cancerogenic health risk in occupational exposure to BD and ST for employees in the carpet industry using USEPA method.

Findings from this study revealed that mean personal exposure to BD and ST was equal to 0.039 mg m−3 (0.017 ppm) and 12.108 mg m−3 (2.840 ppm), respectively. Concentrations of BD and ST determined in our studies were lower than the Occupational Exposure Limit (OEL) values recommended by the Environmental and Occupational Health Center of Iran (EOHCI). Occupational Exposure Limit Time-Weighted Average [OEL-TWA] was set to be 4.42 mg m−3 (2 ppm) for BD6 and 86 mg m−3 (20 ppm) for ST43. The temperature of the glue injection boiler and dryer and the percent of styrene-butadiene copolymer entering the boiler impacted the concentration of BD and ST. For instance, in some process the proportion of monomers of 70% to 75% of butadiene versus 20% to 25% of ST are typically used. While, in other process the proportion of monomers of 55% to 65% of butadiene versus 45% to 35% of ST are mainly used. These differences could influence the concentration of BD and ST emitted from the source.

Cancerogenic risk assessment (CR) was applied to determine the cancer risk due to BD and ST exposure. Based on the International Agency for Research on Cancer (IARC), BD is designated as a group 1 human carcinogen (carcinogenic to humans) and ST is listed as a group 2A human carcinogen (probably carcinogenic to humans)43. The mean ± SD values of estimated cancer risk in inhalation exposure to BD and ST were 5.13 × 10–3 ± 3.85 × 10–4 (5.13 additional cases of cancer per 1,000 workers exposed) and 1.44 × 10–3 ± 2.36 × 10–4 (1.44 additional cases of cancer per 1,000 workers exposed), respectively (Table 3). According to the results, the 95th percentile of the carcinogenic risk values obtained for BD and ST (1.22 × 10–2 and 5.77 × 10–3, respectively) exceeded recommended by USEPA the acceptable risk level (1 × 10−6), meaning a significant cancerogenic risk for the exposed employees. The probability distributions of the cancerogenic risk for the pollutants (mean, median, 5th percentile, and 95th percentile) are depicted in Fig. 1. The mean and median values of BD and ST were higher than the USEPA acceptable level of 1 × 10–6. The cancerogenic risk values calculated for 5% of the high-risk population for BD and ST were 1.22 × 10–2 and 5.77 × 10–3, respectively. Even for BD the 5% of low-risk population, the cancerogenic risk values was 1.66 × 10–3 (1.66 additional cases of cancer per 1,000 individuals exposed). These findings revealed that occupational exposure through inhalation to BD and ST posed potential cancer risks for workers of finishing shops in carpet factories.

The cancerogenic risk for (A) 1,3-butadiene (BD), (B) Styrene (ST).

Sadeghi-Yarandi et al. conducted a study in health risk assessment of inhalation exposure to 1,3-butadiene in a petrochemical plant. Their results showed that carcinogenic risk due to 1,3-butadiene at all investigated units was higher than the acceptable level6. The health risk of ST exposure in respiratory zone of workers in an electronics industry assessed by Mohammadyan et al.45 revealed the average CR value equal to 1.4 × 10−3, which was higher than the acceptable carcinogenic risk level. Also, carcinogenic risk values in the inhalation exposure of ST exceeding the acceptable level were reported in the other industries: motor and motorcycle, repair services, ship and boat building, basic chemical and plastic products manufacturing15, petrochemical industry43, and urban environment46.

Non cancerogenic risk in the inhalation exposure to BD and ST was calculated and presented as hazard quotient (HQ) values in Table 4. Our research revealed that all calculated non-cancerogenic risk values for BD and ST were higher than the acceptable risk level of 1, except minimum ST concentration. The probability distributions of the non-cancer risk for the pollutants (mean, median, 5th percentile, and 95th percentile) are presented in Fig. 2. The mean and median risk values for BD and ST were higher than the accepatable non-carcinogeinic risk value (HQ = 1). The non-cancerogenic risk values calculated for 5% of the high-risk population for BD and ST were equal to 2.15 × 101 and 2.01 × 101, respectively.

The non-cancerogenic risk for (A) 1,3-butadiene (BD), (B) Styrene (ST).

Mohammadyan et al., in risk assessment of ST exposure in an electronics industry in Iran45, revealed that non-cancerogenic risk values for ST were higher than the acceptable level of 1 in all industrial units45, which is consistent with the results of the present study. The research of Yimrungruang et al.47 on the health risk assessment of gas service station workers to VOCs exposure indicated that the mean ST concentration was equal to 2.4 × 10–3 mg m−3 and calculated non-cancerogenic risk value was < 147, which is different from this study. As a direct relationship exists between the non-cancerogenic risk value for particular chemical and its concentration, it might be the most important factor of the calculation variations in these studies. For instance, the mean inhalation exposure to ST in the current study was 12.108 mg m−3, which is vastly higher than in the Yimrungruang et al. studies which was reported to be 2.4 × 10–3 mg m−347. This has resulted that in the present study the mean non-cancerogenic risk values were so many times higher than in the Yimrungruang et al. study47.

Sadeghi-Yarandi et al., reported high HQ value for BD exposure in a petrochemical plant6. Furthermore, studies on BD and ST industrial exposure in microenvironments like residential homes, in-office, in-vehicle9, motor and motorcycle, repair services, ship and boat building, basic chemical and plastic products manufacturing15, and petrochemical industry43, reported existing risk, thus application of engineering control methods to decrease occupational exposure to BD and ST seems to be necessary.

The sensitivity analysis by the Crystal Ball software was applied to identify factors influencing the most on the calculated risk values.

The parameter influencing the most on the carcinogenic risk values in inhalation exposure pathway was revealed to be the pollutant concentration. Pollutant concentration affected the calculated risk values in 54.9% for BD and in 81.2% for ST (Fig. 3). The second important factor was revealed to be the exposure frequency with the values of 8% for BD and 3.4% for ST. The third important factor was the exposure duration with the values of 5% for BD and 2.8% for ST. The body weight was revealed to have the negative impact on the calculation of the carcinogenic risk values. In the case of non-carcinogenic risk values in the inhalation exposure to BD and ST, the most significant variables were: the pollutant concentration (92.2% for BD and 74.4% for ST), exposure time (6.2% for BD and 16.5% for ST), exposure duration (1.5% for BD and 5.5% for ST), and exposure frequency (0.2% for BD and 3.6% for ST) (Fig. 4).

Sensitivity analysis of the cancer risk (CR) for (A) 1,3-butadiene, (B) styrene.

Sensitivity analysis of the non-cancer risk (HQ) for (A) 1,3-butadiene, (B) styrene.

Based on the obtained findings it should be highlighted that BD and ST concentrations need to be controlled in the workplaces as they are major influencing factors of the potential risk. Also, inhalation exposure to BD and ST should be controlled via engineering control measures as well as the exposure time of the workers should be decreased via management controls. Also, biomonitoring of the metabolites of these pollutants in the biological fluids is suggested for the future research.

As the uncertainty is an inherent component of the quantitative risk assessment in our studies, we considered its three components. Regarding the environmental part, uncertainty was related with performed measurements. First aspect was the number of samples collected from the worker's breathing zone in the carpet production industry for establishing the range of inhalational BD and ST concentrations. Secondly, samples were taken only during winter season, so it is not known at this stage of the research if these concentrations were only seasonal or could be representative for the whole year. Regarding the population part, we polled 75 male employees working in finishing shops and based on these results population data for HHRA calculations were obtained, like body mass and exposure rates: duration, frequency, and time. Regarding the toxicological part, for non-carcinogenic risk calculations in our research we used the inhalational exposure pathways expressed by the exposure concentration (EC) parameter as RfC values were available in toxicological databases for both, SD and ST. In the case of carcinogenic risk, we used chronic daily intake (CDI) parameter, instead of EC, because for inhalational unit risk (IUR) needed in further calculations values were missing for ST. In carcinogenic risk analysis, as slope factor (SF) values were available for both BD and ST (IUR value was available only for BD), we used exposure pathway expressed by calculating CDI to have the consistency between BD and ST in carcinogenic risk calculations. However, we have used the lowest available SF value for ST obtained from scientific literature41 in accordance with the conservative risk assessment principle. SF value for ST used in our studies did not appear in the toxicological databases as the result of the lack of the agreement on confirmed SF and/or IUR value for ST. In accordance with the used conservative risk assessment principle, our risk results might be overestimated as they were prone to human health protection purposes. Expanded and more detailed future research would allow for more reliable risk calculation based on higher number of participants, different sex and age of participants, higher frequency, and longer time of air pollutants measurements in the future studies.

As it was underlined earlier, best to our knowledge, this was the first studies of the health risk assessment during the process of carpet production. Thus, comparison of our results with the other research is rather limited. The results of a study by Yarandi et al.6, the mean cancerogenic risk (CR) of BD among the petrochemical plant workers was 2.71 × 10–3, which exceeded the acceptable risk value (10–6) recommended by the USEPA43. The average carcinogenic risk value of ST for workers of electronics industry in Iran was equal to 1.4 × 10–3 exceeding the acceptable level of USEPA45. The results of Ahmadi-Moshiran et al.43 studies in Iran, indicated that the cancerogenic risk (CR) values for BD and ST in a petrochemical industry were also higher than the acceptable level of 10–6 as well as the mean HQ values of non-carcinogenic risk for BD and ST were 4.04 and 0.19, respectively exceeding the acceptable risk level of 1. In the study by Hahm al.15, the non-carcinogenic risk (HQ) values for ST in motor vehicle, motorcycle maintenance, repair services, ship and boat building, basic chemical and plastic products manufacturing industries were higher than the acceptable risk level of 1 recommended by the USEPA.

Regarding the toxicological data in the case of BD both RfC and SF values were available in the Integrated Risk Information System (IRIS) database40. In the case of ST only RfC value was confirmed in the toxicological databases. In the case of the SF value for ST we used the value available in other research41, having in mind that so far there is no one agreed value. Banton et al. 201948 reported that RfC value for ST for general population based on ototoxicity was estimated to be 6 ppm (about 25 mg m3). However, as this proposed RfC value was not time-adjusted from the worker occupational exposure limit of 20 ppm, it must be considered that other endpoints might occur at lower exposure duration if based on worker data48. Toxicological values extrapolated from non-ototoxic endpoints from animal studies may identify health values lower than those estimated on occupational ototoxicity when dose-duration adjustments, especially that application of appropriate AFs is required48. Also, having in mind that regarding the cancerogenic risk assessment of ST researchers are leaning to the Margin of Exposure (MOE) approach instead of the dose–response model48, in this study we used the existing SF value for ST to have the consistency with the applied risk calculation methodology. We have this point in mind and in the future research it is planned to apply other risk assessment models using Margin of Exposure (MOE), Point of Departure (POD), or Mode of Action (MOA), depending on the trends in the new generation risk assessment methods49. As the HHRA aim is to indicate the lowest exposure at which risk of the adverse health effects might occur according to the conservative risk assessment principle, in our calculations we have used the lowest SF value available from research to have the idea of the risk values, especially because this research was not performed before. Banton et al.48 also reported that in the fiber-reinforced polymer composites industry, HQ values for workers in open-mold processes without any respiratory protection exceeded the acceptable risk level of 1. Non-carcinogenic risk values were reduced < 1 when appropriate respiratory protection measures were used in open-molding processes, indicating that respiratory protection should be required by risk management measure for high workplace exposures48.

Our research revealed the need of the continuation of this study. We investigated only some of VOCs, namely BD and ST. However in our research only men were investigated, while women and even children, who often work in such industries in Asian countries, were not investigated. Regarding the chemical nature of the substances used in the carpet industry also customers seem to be affected with the vapours that can be released from the new carpet products50. The recent research indicate the need for intensive research in the context of the consumer health safety as carpets are presumed an exposure source of many compounds like emerging contaminats, perfluoroalkylated substances (PFAS), volatile organic compounds (VOCs), and semi-volatile organic compounds (SVOCs), especially that mechanisms like abrasion, diffusion, partitioning to airborne particles and settled dust are not well know51. Biomonitoring research on exposed population would be an important part of the further investigations52.

Our research were the first studies on occupational inhalational exposure to 1,3-butadiene (BD) and styrene (ST) in the finishing shops of the carpet production industry in Kashan city, Iran. Using the USEPA Human Health Risk Assessment methodology carcinogenic (CR) and non-carcinogenic (HQ) risk values were calculated. The concentrations of BD and ST in the respiratory zone of the employees were lower than the permissible occupational exposure limits, 4.42 mg m−3 (2 ppm) for BD (6) and 86 mg m−3 (20 ppm), and were equal to 0.039 mg m−3 (0.017 ppm) and 12.108 mg m−3 (2.840 ppm), respectively. Non-cancerogenic risk (HQ) values for BD and ST were higher than the acceptable level of 1 indicating the possibility of adverse health effects to occur. The mean HQ values for the carpet industry workers were equal to 8.67 × 100 for SD and 5.13 × 100 for ST. Cancerogenic risks (CR) values calculated for BD and ST were higher than acceptable level of 10–6. The mean CR values for employees were equal to 5.13 × 10–3 for SD and 1.44 × 10–3 for ST. The sensitivity analysis revealed that BD and ST concentrations measured in the employees’ respiratory zone were the most significant parameter in the total value of both carcinogenic and non-carcinogenic risk. The application of technical and engineering control measures is recommended to reduce employees’ risk including exposure time reduction performed by management controls. More detailed research, including other susceptible populations and biomonitoring analyses of the metabolites of BD and ST in the biological fluids are required and justified for future research based on our preliminary investigations.

The data sets analysed during the current study are available from the corresponding author upon reasonable request.


American Conference of Governmental Industrial Hygienists

Averaging time

American Thoracic Society


Body weight

Concentration of pollutant

Chronic daily intake

Carcinogenic risk


Deoxyribonucleic acid

Exposure concentration

Exposure duration

Exposure frequency

Human Health Risk Assessment

Hazard quotient

International Agency for Research on Cancer

Immediately Dangerous to Life or Health

Inhalation rate

Integrated Risk Information System

Inhalational unit risk

Mode of action

Margin of exposure

Monte Carlo simulation

National Institute for Occupational Safety and Health

Occupational exposure limit

Occupational Safety and Health Administration

95Th percentile

Permissible exposure limit

Perfluoroalkylated substances

Point of departure

Quality control/quality assurance

Reference concentration

Reactive oxygen species

Standard deviation

Safety data sheet

Slope factor


Short term exposure limit

Semi-volatile organic compounds

Threshold limit value

Time weighted average

United States National Toxicology Program

United States Environmental Protection Agency

Volatile organic compounds

World Health Organization

Lee, C.-W., Dai, Y.-T., Chien, C.-H. & Hsu, D.-J. Characteristics and health impacts of volatile organic compounds in photocopy centers. Environ. Res. 100, 139–149 (2006).

Article CAS Google Scholar

Fermo, P. et al. Improving indoor air quality through an air purifier able to reduce aerosol particulate matter (PM) and volatile organic compounds (VOCs): Experimental results. Environ. Res. 197, 111131 (2021).

Article CAS Google Scholar

Sarkhosh, M. et al. Indoor Air Quality associations with sick building syndrome: An application of decision tree technology. Build. Environ. 188, 107446 (2021).

Article Google Scholar

Soni, V., Singh, P., Shree, V. & Goel, V. Air Pollution and Control 119–142 (Springer, 2018).

Book Google Scholar

Mo, Z., Lu, S. & Shao, M. Volatile organic compound (VOC) emissions and health risk assessment in paint and coatings industry in the Yangtze River Delta. China. Environ. Pollut. 269, 115740 (2021).

Article CAS Google Scholar

Sadeghi-Yarandi, M. et al. Cancer and non-cancer health risk assessment of occupational exposure to 1, 3-butadiene in a petrochemical plant in Iran. Toxicol. Ind. Health. 36, 960–970 (2020).

Article CAS Google Scholar

Sathiakumar, N. et al. 1, 3-Butadiene, styrene and selected outcomes among synthetic rubber polymer workers: Updated exposure-response analyses. Chem. Biol. Interact. 347, 109600 (2021).

Article CAS Google Scholar

Linde. Safety Data Sheet: 1,3-butadiene (2022).

Huy, L. N., Lee, S. C. & Zhang, Z. Human cancer risk estimation for 1, 3-butadiene: An assessment of personal exposure and different microenvironments. Sci. Total Environ. 616, 1599–1611 (2018).

Article ADS Google Scholar

Scientific, T. Safety Data Sheet. (2022).

Roth, C. Safety Data Sheet. (2021).

Marczynski, B., Peel, M. & Baur, X. New aspects in genotoxic risk assessment of styrene exposure: A working hypothesis. Med. Hypotheses 54, 619–623 (2000).

Article CAS Google Scholar

Ruder, A. M. & Bertke, S. J. Cancer incidence among boat-building workers exposed to styrene. Am. J. Ind. Med. 60, 651–657 (2017).

Article CAS Google Scholar

Humans, I. W. G. o. t. E. o. C. R. t. Styrene, Styrene-7, 8-Oxide, and Quinoline. (2019).

Hahm, M., Lee, J., Lee, M.-Y. & Byeon, S.-H. Health risk assessment of occupational exposure to styrene depending on the type of industry: Data from the Workplace Environmental Monitoring Program in Korea. Hum. Ecol. Risk Assess. 22, 1312–1322 (2016).

Article CAS Google Scholar

Organization, W. H. & Cancer, I. A. f. R. o. 1, 3-Butadiene, Ethylene Oxide and Vinyl Halides (Vinyl Fluoride, Vinyl Chloride and Vinyl Bromide) (IARC, 2008).

Google Scholar

Fedotov, A. et al. The production of 1, 3-butadiene from bio-1-butanol over Re-W/α-Al2O3 porous ceramic converter. Catal. Commun. 128, 105714 (2019).

Article CAS Google Scholar

Organization, W. H. & Cancer, I. A. f. R. o. 1, 3-Butadiene, Ethylene Oxide and Vinyl Halides (Vinyl Fluoride, Vinyl Chloride and Vinyl Bromide) (IARC Press, International Agency for Research on Cancer, 2008).

Google Scholar

Mulholland, J. A. et al. Air Emissions from Carpet Manufacturing Processes (EPA, 2002).

Google Scholar

Yoshida, T., Tainaka, H., Matsunaga, I. & Goto, S. The airborne 1, 3-butadiene concentrations in rubber and plastic processing plants. J. Occup. Health. 44, 56–63 (2002).

CAS Google Scholar

Meyer, K. C., Sharma, B., Kaufmann, B., Kupper, A. & Hodgson, M. Lung disease associated with occupational styrene exposure. Am. J. Ind. Med. 61, 773–779 (2018).

Article Google Scholar

Tong, R., Ma, X., Zhang, Y., Shao, G. & Shi, M. Source analysis and health risk-assessment of ambient volatile organic compounds in automobile manufacturing processes. Hum. Ecol. Risk assess. 26, 359–383 (2018).

Article Google Scholar

Wang, H. et al. Characterization and assessment of volatile organic compounds (VOCs) emissions from typical industries. Chin. Sci. Bull. 58, 724–730 (2013).

Article CAS Google Scholar

Xu, Q. et al. Quantitative differences between common occupational health risk assessment models. J. Occup. Health. 62, e12164 (2020).

Article CAS Google Scholar

Hoseini, L. K., Yengejeh, R. J., Rouzbehani, M. M. & Sabzalipour, S. Health risk assessment of volatile organic compounds (VOCs) in a refinery in the southwest of Iran using SQRA method. Public. Health. Front. 10, 452–459 (2022).

Google Scholar

Zain, S. M. S. M., Azmi, W. N. F. W., Veloo, Y. & Shaharudin, R. Formaldehyde exposure, health symptoms and risk assessment among hospital workers in malaysia. J. Environ. Prot. 10, 861–879 (2019).

Article CAS Google Scholar

Ahmad Wani, K. & Jaiswal, Y. K. Health risk factors in different seasons of carpet industry in Kashmir, India. Int. J. Occup. Saf. Ergon. 18, 571–577 (2012).

Article Google Scholar

Subedi, R. K. Occupational safety and health among carpet factory workers in a district of Nepal. Open. Access. Libr. 2, 1 (2015).

Google Scholar

Wani, K. A. & Khan, R. Evaluation of occupational exposure of carpet weavers in northern province of Madhya Pradesh (India) during different seasons. Indian J. Occup. Environ. Med. 19, 110 (2015).

Article Google Scholar

Becher, R. et al. Do carpets impair indoor air quality and cause adverse health outcomes: A review. Int. J. Environ. Res. Public Health. 15, 184 (2018).

Article Google Scholar

Khoshakhlagh, A. H. & Morais, S. Volatile organic compounds in carpet manufacturing plants: Exposure levels and probabilistic risk assessment using Monte-Carlo simulations. Hum. Ecol. Risk Assess. 28, 972–982 (2022).

Article CAS Google Scholar

Khoshakhlagh, A. H., Chuang, K.-J. & Kumar, P. Health risk assessment of exposure to ambient formaldehyde in carpet manufacturing industries. Environ. Sci. Pollut. Res. (2022).

Article Google Scholar

Khoshakhlagh, A. H., Saberi, H. R., Gruszecka-Kosowska, A. & Kumar, V. Respiratory functions and health risk assessment in inhalational exposure to vinyl acetate in the process of carpet manufacturing using Monte Carlo simulations. Environ. Sci. Pollut. Res. (2022).

Article Google Scholar

Natarajan, R. Performance evaluation of biodegradation of isoprene–acetone mixture in integrated biofilter. J. Environ. Eng. 146, 04020030 (2020).

Article CAS Google Scholar

Natarajan, R., Al-Sinani, J., Viswanathan, S. & Manivasagan, R. Biodegradation of ethyl benzene and xylene contaminated air in an up flow mixed culture biofilter. Int. Biodeterior. Biodegrad. 119, 309–315 (2017).

Article CAS Google Scholar

Rajamohan, N. Performance evaluation of hybrid reactor for bio filtration of n-hexane–styrene mixture. Desalin. Water Treat. 1, 6 (2017).

Google Scholar

Moradi, H., Talaiekhozani, A., Kamyab, H., Chelliapan, S. & Nadda, A. K. Development of equations to predict the concentration of air pollutants indicators in Yazd City, Iran. J. Inorg. Organomet. Polym. Mater. (2022).

Article Google Scholar

Wani, K. A. & Jaiswal, Y. Occupational health risk factors in carpet industry: A Review. Asian J. Exp. Biol. Sci. 2, 135–139 (2011).

Google Scholar

Traumann, A., Tint, P., Järvik, O. & Oja, V. Management of health hazards during shale oil handling. Agron. Res. 11, 479–486 (2013).

Google Scholar

Agency, U. E. P. (2022).

Moshiran, V. A. et al. Quantitative and semiquantitative health risk assessment of occupational exposure to styrene in a petrochemical industry. Saf. Health. Work. 12, 396–402 (2021).

Article Google Scholar

Means, B. Risk-Assessment Guidance for Superfund. Volume. 1 Human Health Evaluation Manual. Part A. Interim Report (Final) (Environmental Protection Agency, Office of Solid State, 1989).

Google Scholar

Ahmadi-Moshiran, V. et al. Carcinogenic and health risk assessment of respiratory exposure to Acrylonitrile, 1, 3-Butadiene and Styrene (ABS) in a Petrochemical Industry Using the United States Environmental Protection Agency (EPA) Method. Int. J. Occup. Saf. Ergon. 28, 1–9 (2022).

Article Google Scholar

Badeenezhad, A. et al. Factors affecting the nitrate concentration and its health risk assessment in drinking groundwater by application of Monte Carlo simulation and geographic information system. Hum. Ecol. Risk Assess. 27, 1458–1471 (2021).

Article CAS Google Scholar

Mohammadyan, M. et al. Health risk assessment of occupational exposure to styrene in Neyshabur electronic industries. Environ. Sci. Pollut. Res. 26, 11920–11927 (2019).

Article CAS Google Scholar

Du, Z., Mo, J. & Zhang, Y. Risk assessment of population inhalation exposure to volatile organic compounds and carbonyls in urban China. Environ. Int. 73, 33–45 (2014).

Article CAS Google Scholar

Yimrungruang, D., Cheevaporn, V., Boonphakdee, T., Watchalayann, P. & Helander, H. F. Characterization and health risk assessment of volatile organic compounds in gas service station workers. EnvironmentAsia. 2, 21–29 (2008).

Google Scholar

Banton, M. et al. Evaluation of potential health effects associated with occupational and environmental exposure to styrene: An update. J. Toxicol. Environ. Health. B. Crit. Rev. 22, 1–130 (2019).

Article CAS Google Scholar

Gruszecka-Kosowska, A., Ampatzoglou, A. & Aguilera, M. Integration of omics approaches enhances the impact of scientific research in environmental applications. Int. J. Environ. Res. Public Health 19, 8758 (2022).

Article Google Scholar

Onyshko, J. & Hewlett, R. Toxics in Carpets in the European Union. Published by Anthesis (2018).

Haines, S. R. et al. Ten questions concerning the implications of carpet on indoor chemistry and microbiology. Build. Environ. 170, 106589 (2020).

Article Google Scholar

Pena, A. et al. Human biomonitoring of selected hazardous compounds in Portugal: Part I—lessons learned on polycyclic aromatic hydrocarbons, metals, metalloids, and pesticides. Molecules 27, 242 (2021).

Article Google Scholar

Download references

The authors gratefully acknowledge the participation of the employees and managers of the carpet manufacturing plants in this study. We gratefully acknowledge Dr. Jianfeng li for his valuable comments to improve the paper.

The study is supported by Kashan University of medical sciences.

Department of Occupational Health Engineering, School of Health, Kashan University of Medical Sciences, Kashan, Iran

Amir Hossein khoshakhlagh

Department of Environmental Protection, Faculty of Geology, Geophysics, and Environmental Protection, AGH University of Science and Technology, Al. Mickiewicza 30, 30-059, Krakow, Poland

Agnieszka Gruszecka-Kosowska

Department of Chemistry and Chemical Technology, National University of Lesotho, P. O. Roma 180, Lesotho, South Africa

Abiodun Olagoke Adeniji

Institute of Occupational Medicine, Edinburgh, EH14 4AP, UK

Lang Tran

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

A.H.K.: conceptualization, formal analysis, sample collection, project administration, methodology. A.G.-K.: revised, re-writing- reviewing & editing. A.O.A.: re-writing—reviewing & editing. Lang Tran: re-writing—reviewing & editing.

Correspondence to Amir Hossein khoshakhlagh.

The authors declare no competing interests.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

khoshakhlagh, A.H., Gruszecka-Kosowska, A., Adeniji, A.O. et al. Probabilistic human health risk assessment of 1,3-butadiene and styrene exposure using Monte Carlo simulation technique in the carpet production industry. Sci Rep 12, 22103 (2022).

Download citation

Received: 13 September 2022

Accepted: 15 December 2022

Published: 21 December 2022


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.