Investigation and Health Risk Assessment of Volatile Organic Compounds (VOCs) Emitted from Industrial Estate, Rayong City, Thailand in 2024
Susira Bootdee* and Suganya PhantuAbstract Volatile organic compounds (VOCs) in ambient air remain a public health concern due to their potential short and long-term health effects. This study aimed to measure VOC concentrations and assess related health risks in the industrial estate of Rayong city, Thailand. Daily VOC levels were continuously monitored at three Air Quality Monitoring stations in Maptaphut City Municipality from January to December 2024. Mean concentrations of vinyl chloride, benzene, and ethylbenzene varied significantly across the monitoring sites (P < 0.05), while total VOC levels did not differ between the dry and wet seasons (P > 0.05). Meteorological factors, particularly air temperature and relative humidity, appeared to influence VOC levels. Notably, ambient concentrations of 1,3-butadiene (1.13 ± 1.06 to 1.49 ± 1.91 µg/m³) exceeded the annual standard of 0.33 µg/m³ set by the Pollution Control Department. Principal component analysis (PCA) identified several emission sources, including synthetic rubber and thermoplastic resin production, wastewater treatment, vehicle exhaust, industrial solvent use, and emissions from PVC-related processes. Carcinogenic equivalency (CEQ) values were highest near petrochemical facilities producing plastic resins or polymers and areas linked to wastewater treatment, reflecting the presence of carcinogenic VOCs. Despite elevated levels of 1,3-butadiene, both carcinogenic and non-carcinogenic risk assessments suggested no immediate adverse health effects for the surrounding population. However, the findings highlight the need for continued monitoring, particularly due to the long-term cancer risk associated with 1,3-butadiene.
Keywords: Volatile organic compounds (VOCs), Health risk assessment, Industrial area, Principal component analysis (PCA), Carcinogenic equivalency (CEQ)
Funding: The authors are grateful for the research funding provided by King Mongkut’s University of Technology North Bangkok, contract no. KMUTNB-68-BASIC-48.
Citation: Bootdee, S. and Phantu, S. 2026. Investigation and health risk assessment of volatile organic compounds (VOCs) emitted from industrial estate, Rayong city, Thailand in 2024. Natural and Life Sciences Communications. 25(2): e2026032.
Graphical Abstract:

INTRODUCTION
Volatile organic compounds (VOCs) are emitted as gases from certain solids or liquids. VOCs form a large, diverse group of organic chemicals that are gaseous at room temperature. The VOCs are a group of organic compounds occurring in the air. Both NOx and VOCs are significant precursors in photochemical reactions and the generation of secondary air pollutants, which are associated with ozone (O₃) levels in ambient air (Berenjian et al., 2012; Zheng et al., 2023). VOCs are released from various sources, including both those generated by humans and those of natural origin. VOCs generally originate from anthropogenic sources, including automobile exhaust emissions resulting from both complete and incomplete burning of gasoline, gasoline evaporation from engines and tanks, and solvents utilized in industrial and indoor commercial and home applications (Saleem et al., 2022; Wang et al., 2023; Zheng et al., 2023; Lv et al., 2024). VOCs are usually emitted from natural sources such as livestock farms, vegetation, and wildfires (Sharma et al., 2016; Kim et al., 2024). Suwattiga and Limpaseni (2005) reported the source of VOC emissions, including 36% of VOCs from gasoline vehicles, 6% from diesel vehicles, 12% from fuel oil boilers, 12% from gasoline vapor, 5% from paint and thinner vapor, 10% from burning biomass, 3% from cooking food, 8% from municipal waste, and 8% from sources that weren't explained. Therefore, VOCs contaminated in ambient air are emitted from various sources such as vehicles (Suwattiga and Limpaseni, 2005), gas stations (Allahabady et al., 2022), agricultural wastes (Odekanle et al., 2022), industrial estates (Chen et al., 2024), and liquefied petroleum gas (Mihajlović et al., 2021).
VOCs are important atmospheric species that contribute to alterations in atmospheric chemistry. The emissions of VOCs from biomass burning and the VOC chemistry in fire plumes are relatively uncertain (Uttamang et al., 2023). Saeaw and Thepanondh (2015) reported that concentrations of total VOC emissions from the Maptaphut estate in Rayong Province were 24.4-70.9 µg/m³, while the source apportionment of total VOC emissions from the industrial area of Rayong Province in Thailand indicated that 42 to 57% originated from mobile sources, 15 to 45% from industrial processes, and 3 to 10% from household chemical usage. Additionally, the source of VOCs in the community area in a petrochemical industrial area of Maptaphut estate in Rayong province was benzene, total xylene, and ethylbenzene emitted from gasoline combustion, identified as vehicle emissions. Vinyl chloride and 1,2-dichloroethane are emitted from the polyvinyl chloride (PVC) industry's raw materials, while toluene, styrene, and 1,3-butadiene are released from industrial origins (Pinthong et al., 2022). Zhang et al. (2017) revealed that the total ambient VOC emissions from the Jinan petroleum refinery in China amounted to 50.7 µg/m³. The significant VOCs identified include ethane (24.58 µg/m³), ethene (3.94 µg/m³), ethylbenzene (0.82 µg/m³), 1,3-butadiene (0.82 µg/m³), benzene (0.24 µg/m³), and toluene (0.08 µg/m³), among others. Kim et al. (2024) showed that the VOCs emitted from the multi-industrial city of Ulsan, South Korea, included 18.3% toluene (6.37 µg/m³), 13.4% m, p-xylene (4.69 µg/m³), 7.7% ethylbenzene (2.70 µg/m³), and 4.8% benzene (1.69 µg/m³), indicating significant effects of industrial emissions. Moreover, the concentrations of 1,2-dichloroethane, benzene, and 1,3-butadiene emissions from the industrial estate of Rayong city were 4.38 ± 1.99 µg/m³, 1.99 µg/m³, and 2.48 ± 0.61 µg/m³, respectively. Then, the levels of 1,2-dichloroethane, benzene, and 1,3-butadiene exceeded the annual VOC standards by 2.71, 2.58, and 38.00 times (Nakyai et al., 2025).
VOCs include a variety of chemicals, some of which may have short- and long-term adverse health effects. Concentrations of many VOCs are consistently higher indoors (up to ten times higher) than outdoors. VOCs such as benzene, vinyl chloride, 1,2-dichloroethane, trichloroethylene, dichloromethane, 1,2-dichloropropane, tetrachloroethylene, chloroform, and 1,3-butadiene are listed by the Pollution Control Department (PCD) of Thailand as hazardous pollutants that require monitoring (PCD, 2025). Exposure to VOCs can lead to both acute and chronic toxicity in humans. Short-term effects of VOC exposure commonly include irritation of the eyes, nose, throat, and skin, along with symptoms such as headaches, nausea, dizziness, fatigue, and shortness of breath (US-EPA, 2007). Increased concentrations of VOCs can adversely affect health. Long-term exposure to low levels at or exceeding legal limits may be linked to respiratory system problems, including lung tissue damage, respiratory tract irritation, and the exacerbation of pulmonary conditions such as asthma, as well as DNA strand breaks, chromosomal aberrations, mutations, and disruptions in blood cell and vessel function, leading to hematological disorders and injuries (Zhou et al., 2020; Saeedi et al., 2024).
Based on the information submitted, the researcher chose to measure VOC concentrations in ambient air and conduct a health risk assessment within the industrial estate during the rainy, cold, and hot seasons.
MATERIALS AND METHODS
Sampling sites and VOC concentrations
Researchers studied and selected sampling sites in an industrial estate in Rayong province, Thailand. Rayong province's industrial estate lies on the eastern coast of the Gulf of Thailand. The city plays a crucial role in Thailand's industry. The Maptaphut industrial estate covers 8,012.39 acres in Rayong province and is the main industrial estate. Petrochemical and chemical manufacturing facilities, steel and metal industries, machinery and equipment production, coal-powered electricity plants, petroleum refineries, residential districts, and transportation infrastructure dominate the region.
The VOC concentrations provided secondary data for this investigation from the Air Quality Monitoring Station (AQM) of Maptaphut City Municipality. Three AQMs in Maptaphut City Municipality in Rayong city continuously recorded VOC concentrations with the represented IA 1-3 symbol. The IA1 is located in the community area (12º 43' 50.2674"N and 101º 9' 20.2314"E), approximately 2.53 km from the petrochemical industry, which mostly produces plastic resins or polymers molded into products used for food packaging, automotive components, medical devices, electrical appliances, and infrastructure, including pressure-resistant pipes and telecommunication cables, among others. And then, the IA2 is in the public health center and the nearby government office (12º 42' 31.2114" N and 101º9'57.168" E). It is close to the roadside. Furthermore, this sampling site took place 1.38 km away from the oil refineries. The IA3 has been included in the school and community (12º 41' 8.3688" N and 101º 7' 1.218" E). The site is situated approximately 600 m from the liquefied natural gas (LNG) industry and near a high-traffic crossroads during peak times.
VOC concentrations were continuously monitored and recorded at three AQM stations located within Maptaphut City Municipality. The continuous sampling utilized networked online gas analyzers, specifically equipped with Gas Chromatography-Flame Ionization Detectors (GC-FIDs, Baseline 9100, Ametek, USA). The daily real-time measurement focused on the continuous collection of data for vinyl chloride, benzene, toluene, ethylbenzene, p-xylene, and 1,3-butadiene, spanning an entire year from January to December in 2024. The daily meteorological parameters, including wind speed (WS), relative humidity (RH), total precipitation (Rain), air temperature (Temp), and air pressure (PS), were obtained from AQM of Maptaphut City Municipality. The industrial estate remained from January to December 2024, covering both the dry season (January to April and November to December) and the wet season (May to October). The daily sample data started at 00:00 and ended at 23:00 the following day (24 hours). In that time, about 910 data points were recorded from three AMQ stations. Figure 1 depicts the sampling sites.

Figure 1. Sampling sites of Air Quality Monitoring Station (AQM) in Map Ta Phut Town Municipality.
Inhalation of health risk assessment
Non-carcinogenic risk assessment
The hazard quotient (HQ) is a guideline used to evaluate the non-carcinogenic risk associated with human exposure to air contaminants. It was estimated that both children and adults could inhale VOCs through the inhalation exposure pathway (Abd Hamid et al., 2019; Kawichai et al., 2022). Equation (1) calculates the HQ as the ratio of toxicological effects based on the exposure concentration via inhalation (ECInh) and the reference concentration (RfC) for each VOC, as shown in Table 1.

When the HQ >1.0 signifies an opportunity of non-carcinogenic effects on health, whereas HQ<1.0 means no significant risk or negligible hazards. The ECInh for each VOC exposure was determined utilizing Equation (2).

Where EC is the EC of pollutants (µg/m³), C is individual VOC concentrations (µg/m³), ED is the exposure duration (days), EF is the exposure frequency (days/year), ET is the exposure time (hours/day), and AT is the average time (hours). The values of these parameters were shown in Table 2, which is referred to by Morakinyo et al. (2017) and Prasertsin and Nathapindhu (2020).
Table 1. The reference concentration (RfC), inhalation unit risk (IUR), and cancer slope factors (SF) for health risk assessment from the inhalation route for VOCs.
|
VOCs |
RfC (mg/m3) |
IUR (per µg/m3) |
SF (mg/kg-day) |
|
Vinyl chloride |
5.11×10-2 |
4.40×10-6 |
7.20×10-1 |
|
1,3-Butadiene |
2.00×10-3 |
3.00×10-5 |
6.00×10-1 |
|
Benzene |
3.00×10-2 |
7.80×10-6 |
1.50 |
|
Toluene |
5.00 |
|
|
|
Ethylbenzene |
1.00 |
2.50×10-6 |
1.10×10-2 |
|
p-Xylene |
1.00×10-1 |
|
|
Carcinogenic risk assessment
The cancer risk (CRInh) associated with each carcinogenic VOC was determined by assessing the long-term risk of cancer development due to ongoing exposure to these substances, applying an exact method (Abd Hamid et al., 2019; Kawichai et al., 2022). To determine the cancer risk, we investigated the risks linked to each carcinogenic VOC, including vinyl chloride, benzene, ethylbenzene, and 1,3-butadiene, as illustrated by Equation (3).

Table 2. Parameters for health risk assessment for the inhalation route for VOCs.
|
Parameters |
Age category |
||||||
|
Children |
|
Adult |
|||||
|
1-5 years |
6-11 years |
12-19 years |
|
(≥20 years) |
|||
|
Exposure frequency (EF) (days/year) |
350 |
350 |
350 |
|
350 |
||
|
Exposure time (ET) (hour/day) |
24 |
24 |
24 |
|
24 |
||
|
Exposure duration (ED) (year) |
2 |
3 |
3.5 |
|
30 |
||
|
Averaging time (AT) = 70 yearsÍ365 days/year × 24 hour/day (hour) |
613,200 |
613,200 |
613,200 |
|
613,200 |
||
Carcinogenic equivalency (CEQ) is a risk assessment method applied to determine the potential carcinogenicity of volatile organic compounds (VOCs) in relation to reference material. This parameter defines the combined carcinogenic potential of a mixture compared to a particular amount of reference VOC carcinogen, generally benzene (Chen et al., 2024). The CEQ for a VOC has been calculated by multiplying the VOC concentration by the ratio of the cancer slope factors (SF) of the individual VOC to that of benzene, as illustrated in Equation (4).

Data analysis
The VOC concentrations and meteorological parameters were statistically analyzed using SPSS Statistics (Version 28). The VOC concentrations at three sampling sites within the industrial estate in Rayong City were tested for normality using the Kolmogorov–Smirnov test (P < 0.05), which indicated a non-normal distribution. The Mann-Whitney U and Kruskal-Wallis tests were utilized to statistically assess individual VOC concentrations across three sampling sites and seasons at the industrial estate in Rayong City. Moreover, the association between meteorological parameters and VOC concentrations was assessed using Pearson's correlation analysis.
Finally, we statistically analyzed the source identification of VOC concentrations in ambient air at the industrial estate in Rayong City using principal component analysis (PCA), which widely tests the sources of air pollution, agriculture, and the environment (Kawichai et al., 2020; Teravecharoenchai et al., 2021; Chen et al., 2024). In the PCA for source identification, daily VOC concentrations were standardized to achieve simulations of results from data sets of only the concentrations of detection data that were evaluated. The PCA is performed by decomposing the variance matrix of the dataset into a collection of eigenvalues and eigenvectors. The eigenvalues indicate the amount of variance along each principal component (PC), whereas the eigenvectors indicate the positions of the PCs. The varimax method is the most effective and commonly employed orthogonal rotation method for conducting rotation. Components possessing eigenvalues over 1 were preserved. Furthermore, loadings above 0.5 are considered significant for source attribution (Maskey et al., 2018; Mabel and Olayemi, 2020).
RESULTS
VOCs concentrations
Table 3 displays the concentration of volatile organic compounds (VOCs) released from the industrial estate in Rayong city, Thailand. The average concentrations of VOCs for all sampling sites in descending order were ethylbenzene (2.17 ± 4.48 to 11.0 ± 11.7 µg/m³), toluene (1.78 ± 1.66 to 4.06 ± 4.42 µg/m³), 1,3-butadiene (1.13 ± 1.06 to 1.49 ± 1.91 µg/m³), p-xylene (0.85 ± 0.36 to 1.33 ± 1.38 µg/m³), benzene (0.60 ± 0.63 to 1.24 ± 0.43 µg/m³), and vinyl chloride (0.03 ± 0.03 to 0.87 ± 1.64 µg/m³). The highest values of benzene (1.24 ± 0.43 µg/m³) and toluene (4.06 ± 4.42 µg/m³) were found at the IA1 site, while the highest mean annual values of ethylbenzene (11.0 ± 11.7 µg/m³) and p-xylene (1.33 ± 1.38 µg/m³) were emitted from the IA2 site. IA1 is located within a community area and close to a petrochemical facility that produces plastic resins and polymers. The ambient VOCs at this site are likely influenced by both traffic emissions and nearby industrial activity. IA2 is situated near several oil refineries, and this monitoring station is also positioned close to a roadside, which may contribute additional traffic-related VOCs. The IA3 site reported that the vinyl chloride (0.87 ± 1.64 µg/m³) and 1,3-butadiene (1.49 ± 1.91 µg/m³) concentrations were the highest values. The IA3 is situated in the school and community area near a high-traffic crossroads during peak times and is close to the liquefied natural gas (LNG) industry. Previous study of ambient VOC concentration at the Map Ta Phut industrial estate in Rayong province. They found that the major VOCs in the community area were benzene, total xylene, and ethylbenzene. In contrast, samples from the industrial estates showed higher levels of vinyl chloride and 1,2-dichloroethane, which are raw materials used in the PVC industry, as well as toluene, styrene, and 1,3-butadiene released from industrial processes (Pinthong et al., 2022). Furthermore, the total volatile organic compounds (VOCs) during the dry seasons at the IA1, IA2, and IA3 sites were 8.88, 22.33, and 12.05 µg/m³, whereas they were 11.45, 9.20, and 8.45 µg/m³ during the wet season, respectively (Figure 2).
The VOC concentrations at three sampling sites within the industrial estate in Rayong City displayed a non-normal distribution. The Mann-Whitney U and Kruskal-Wallis tests were applied to statistically evaluate individual VOC concentrations across three sampling sites within the industrial estate in Rayong City (Table 3). The results indicated that the mean concentrations of vinyl chloride, benzene, and ethylbenzene differed significantly among the sampling sites (P < 0.05). The vinyl chloride values emitted from the IA1, IA2, and IA3 sites were significantly different (P < 0.05), while the concentrations of benzene and ethylbenzene in the IA2 and IA3 sites were not significantly different (P > 0.05). Moreover, 1,3-butadiene, toluene, and p-xylene levels were not significantly different among sampling sites (P > 0.05). However, the total VOC concentrations during the dry season were not significantly different from the wet season in all sites (P > 0.05).
Table 3. The VOC concentrations emitted from an industrial estate in Rayong city in 2024.
|
VOCs |
Site |
Concentrations of VOCs (µg/m3 ± SD) |
||||||||||||
|
Jan. (n=31) |
Feb. (n=29) |
Mar. (n=31) |
Apr. (n=30) |
May (n=31) |
Jun. (n=30) |
Jul. (n=31) |
Aug. (n=31) |
Sep. (n=30) |
Oct. (n=31) |
Nov. (n=30) |
Dec. (n=31) |
Annual (Mean ± SD) |
||
|
Vinyl Chloride |
IA1 |
0.07 ± 0.08 |
0.07 ± 0.06 |
0.06 ± 0.02 |
0.22 ± 0.42 |
0.09 ± 0.16 |
0.08 ± 0.08 |
0.10 ± 0.18 |
0.16 ± 0.45 |
0.07 ± 0.05 |
0.33 ± 0.59 |
0.38 ± 0.46 |
0.25 ± 0.36 |
0.16 ± 0.11a |
|
IA2 |
0.03 ± 0.08 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.07 ± 0.13 |
0.09 ± 0.24 |
0.04 ± 0.11 |
0.04 ± 0.09 |
0.01 ± 0.00 |
0.03 ± 0.03b |
|
|
IA3 |
0.09 ± 0.00 |
0.09 ± 0.00 |
5.92 ± 14.40 |
0.46 ± 0.47 |
0.09 ± 0.00 |
0.09 ± 0.01 |
0.80 ± 2.46 |
1.31 ± 0.75 |
0.46 ± 0.31 |
0.50 ± 0.24 |
0.32 ± 0.07 |
0.30 ± 0.15 |
0.87 ± 1.64c |
|
|
1,3-Butadiene |
IA1 |
0.33 ± 0.01 |
0.43 ± 0.22 |
0.58 ± 0.31 |
0.52 ± 0.27 |
0.37 ± 0.10 |
0.37 ± 0.11 |
0.41 ± 0.23 |
0.44 ± 0.24 |
0.39 ± 0.21 |
11.5 ± 21.0 |
1.53 ± 1.84 |
0.60 ± 1.22 |
1.46 ± 3.19 a |
|
IA2 |
1.00 ± 0.82 |
1.31 ± 0.39 |
0.70 ± 0.53 |
0.26 ± 0.11 |
0.24 ± 0.13 |
1.54 ± 1.64 |
3.84 ± 4.54 |
0.83 ± 1.78 |
2.36 ± 2.88 |
0.34 ± 0.51 |
0.32 ± 0.60 |
0.76 ± 2.73 |
1.13 ± 1.06 a |
|
|
IA3 |
0.21 ± 0.02 |
0.22 ± 0.18 |
4.72 ± 3.20 |
4.06 ± 3.55 |
0.40 ± 0.46 |
0.22 ± 0.11 |
0.21 ± 0.04 |
0.22 ± 0.03 |
0.52 ± 0.75 |
5.03 ± 3.17 |
0.96 ± 0.74 |
1.14 ± 1.99 |
1.49 ± 1.91 a |
|
|
Benzene |
IA1 |
1.00 ± 0.34 |
1.15 ± 0.44 |
1.30 ± 1.21 |
1.10 ± 0.82 |
0.89 ± 0.56 |
2.25 ± 1.51 |
1.29 ± 0.83 |
1.67 ± 1.27 |
1.66 ± 1.18 |
0.76 ± 0.12 |
0.94 ± 0.69 |
0.85 ± 0.19 |
1.24 ± 0.43a |
|
IA2 |
ND |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00 |
0.19 ± 0.00b |
|
|
IA3 |
0.54 ± 0.40 |
0.09 ± 0.11 |
0.17 ± 0.25 |
0.06 ± 0.00 |
0.61 ± 0.63 |
0.12 ± 0.28 |
0.82 ± 3.79 |
0.17 ± 0.17 |
0.23 ± 0.40 |
0.78 ± 0.45 |
2.01 ± 1.86 |
1.63 ± 0.86 |
0.60 ± 0.63b |
|
|
Toluene |
IA1 |
12.60 ± 15.30 |
2.51 ± 3.92 |
1.43 ± 2.15 |
0.83 ± 0.14 |
3.51 ± 4.48 |
0.87 ± 0.15 |
1.06 ± 1.08 |
1.40 ± 1.40 |
1.61 ± 1.02 |
9.00 ± 9.62 |
1.99 ± 2.54 |
11.9 ± 18.1 |
4.06 ± 4.42a |
|
IA2 |
6.09 ± 3.62 |
3.20 ± 3.61 |
0.76 ± 0.74 |
0.45 ± 0.92 |
1.74 ± 2.57 |
0.33 ± 0.16 |
0.34 ± 0.11 |
0.68 ± 0.67 |
1.36 ± 2.17 |
2.86 ± 2.67 |
1.93 ± 2.04 |
1.60 ± 1.46 |
1.78 ± 1.66 a |
|
|
IA3 |
1.82 ± 1.63 |
ND |
ND |
0.58 ± 0.19 |
3.58 ± 2.49 |
2.04 ± 2.83 |
3.11 ± 2.35 |
2.25 ± 3.00 |
1.63 ± 1.65 |
5.61 ± 3.44 |
5.38 ± 9.18 |
2.90 ± 1.76 |
2.89 ± 1.61 a |
|
|
Ethylbenzene |
IA1 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
0.01 ± 0.00 |
9.74 ± 26.7 |
13.4 ± 30.7 |
0.15 ± 0.00 |
0.20 ± 0.11 |
0.18 ± 0.09 |
0.23 ± 0.12 |
1.47 ± 6.48 |
0.67 ± 2.84 |
2.17 ± 4.48 a |
|
IA2 |
25.40 ± 3.80 |
26.6 ± 1.75 |
23.8 ± 1.24 |
21.5 ± 5.78 |
23.0 ± 0.03 |
7.09 ± 10.1 |
0.77 ± 0.14 |
0.88 ± 0.42 |
0.88 ± 0.24 |
0.82 ± 0.28 |
0.76 ± 0.19 |
0.82 ± 0.33 |
11.0 ± 11.7 b |
|
|
IA3 |
7.70 ± 0.59 |
6.77 ± 1.28 |
2.69 ± 0.35 |
1.98 ± 0.04 |
1.96 ± 0.04 |
1.98 ± 0.10 |
1.95 ± 0.00 |
1.97 ± 0.05 |
1.96 ± 0.02 |
2.03 ± 0.08 |
2.10 ± 0.17 |
2.02 ± 0.14 |
2.93 ± 2.03 b |
|
|
p-Xylene |
IA1 |
0.71 ± 0.00 |
0.81 ± 0.51 |
0.71 ± 0.00 |
0.71 ± 0.01 |
0.71 ± 0.00 |
0.71 ± 0.00 |
0.71 ± 0.00 |
0.82 ± 0.24 |
0.84 ± 0.53 |
0.71 ± 0.0 |
1.98 ± 6.48 |
0.73 ± 0.05 |
0.85 ± 0.36a |
|
IA2 |
2.90 ± 2.81 |
4.93 ± 2.11 |
2.37 ± 4.35 |
0.51 ± 0.01 |
0.52 ± 0.05 |
0.56 ± 0.10 |
0.82 ± 0.90 |
0.70 ± 0.27 |
0.68 ± 0.19 |
0.68 ± 0.33 |
0.60 ± 0.16 |
0.71 ± 0.44 |
1.33 ± 1.38a |
|
|
IA3 |
0.83 ± 0.30 |
0.61 ± 0.20 |
1.02 ± 0.38 |
0.41 ± 0.22 |
3.29 ± 3.40 |
ND |
ND |
0.64 ± 0.69 |
0.52 ± 0.44 |
0.83 ± 0.65 |
1.75 ± 1.24 |
1.25 ± 1.49 |
1.12 ± 0.86a |
|
Note: ND: no data. Mann-Whitney U test for individual VOCs in the industrial estate. a, b, c is a significant difference in each site at the industrial estate in Rayong City for using the Kruskal-Wallis test (P-value > 0.05).

Figure 2. Total VOC concentrations during wet and dry seasons emitted from an industrial estate in Rayong city in 2024.
Table 4 illustrates the correlation between VOC concentrations and meteorological parameters in the Rayong Industrial Estate in 2024. Research results revealed that toluene showed a significant negative correlation with relative humidity (r = -0.166 to -0.266, P < 0.01) and air temperature (r = -0.263 to -0.315, P < 0.01), while benzene and ethylbenzene were positively correlated with air temperature and relative humidity in some sites. Moreover, a moderately significant correlation was seen between toluene and pressure (r = -0.488, P < 0.01) at the IA2 site. The correlation between benzene and atmospheric pressure in the IA2 site was strong and significantly positive (r = 0.713), while benzene and air temperature (r = 0.220) and wind speed (r = 0.253) were positively correlated (P < 0.01).
Table 4. Pearson’s correlation coefficients (r) between VOC concentrations and meteorological parameters in a Rayong Industrial Estate (2024).
|
Parameters |
Vinyl Chloride |
1,3-Butadiene |
Benzene |
Toluene |
Ethylbenzene |
P-Xylene |
|
|
IA1 (n=363) |
|||||||
|
WS (m/s) |
0.237** |
0.000 |
0.040 |
-0.021 |
-0.004 |
0.044 |
|
|
RH (%) |
0.016 |
-0.062 |
0.156** |
-0.266** |
0.163** |
-0.049 |
|
|
Rain (mm) |
-0.073 |
0.147** |
-0.087 |
-0.079 |
0.028 |
-0.018 |
|
|
Temp (ºC) |
0.016 |
-0.062 |
0.155** |
-0.263** |
0.162** |
-0.050 |
|
|
PS (hPa) |
0.043 |
0.072 |
-0.156** |
0.185** |
-0.013 |
0.067 |
|
|
IA2 (n=351) |
|||||||
|
WS (m/s) |
-0.092 |
0.156** |
0.253** |
-0.246** |
-0.082 |
-0.069 |
|
|
RH (%) |
0.045 |
0.003 |
0.175** |
-0.166** |
0.119* |
-0.076 |
|
|
Rain (mm) |
0.050 |
0.014 |
0.096 |
-0.084 |
-0.103 |
0.040 |
|
|
Temp (ºC) |
-0.017 |
-0.038 |
0.220** |
-0.315** |
0.207** |
-0.200** |
|
|
PS (hPa) |
-0.040 |
0.073 |
0.713** |
-0.488** |
-0.208** |
-0.263** |
|
|
IA3 (n=196) |
|||||||
|
WS (m/s) |
0.093 |
-0.056 |
-0.076 |
-0.109 |
-0.161 |
0.111 |
|
|
RH (%) |
0.030 |
0.112* |
0.117* |
-0.181** |
-0.128** |
-0.055 |
|
|
Rain (mm) |
-0.017 |
0.082 |
0.061 |
-0.028 |
-0.097 |
-0.009 |
|
|
Temp (ºC) |
0.020 |
0.089 |
-0.052 |
-0.047 |
-0.138** |
0.116 |
|
|
PS (hPa) |
0.024 |
0.047 |
-0.124* |
0.160* |
0.116 |
-0.160 |
|
Note: **Correlation is significant at the 0.01 level (2-tailed); * Correlation is significant at the 0.01 level (2-tailed)
Source apportionment of VOCs by principal component analysis (PCA)
The results of the PCA classifying the sources of VOCs at an industrial estate in Rayong city are shown in Table 5. A threshold of 0.50 for factor loadings was determined to assess statistical significance. Principal component analysis (PCA) accounted for 58 to 77% of the cumulative variance, with eigenvalues from the varimax-rotated factor analysis exceeding 1. At the IA1 site, four components were identified and explained, with cumulative variance ranging from 23% to 77% depending on the site and eigenvalues ranging from 1.028 to 1.376. The component includes vinyl chloride and p-xylene for PC1, 1,3-butadiene for PC2, ethylbenzene for PC3, and toluene for PC4. These components were slurry, open equipment, vessel, and raw material in PVC industrial emissions, synthetic rubber and thermoplastic resin production, wastewater treatment, and vehicle exhaust (Chen et al., 2016; Hosaini et al., 2017; Sukiit et al., 2021; Malakan et al., 2022; Pinthong et al., 2022; Kim et al., 2024). At the IA2 site, the cumulative variance ranged from 37% to 71%, and the eigenvalues ranged from 1.027 to 2.217. PC1 was comprised significantly of toluene, ethylbenzene, and p-xylene, which indicates industrial solvent usage. The compounds PC2 and PC3 are vinyl chloride and 1,3-butadiene, which are classified as belonging to the oil refinery, slurry, open equipment, vessel, and raw material in PVC industrial emissions, as well as the synthetic rubber and thermoplastic resins industry and wastewater treatment (Chen et al., 2016; Hosaini et al., 2017; Sukiit et al., 2021; Kim et al., 2024). Moreover, the IA3 site indicated that three factors, including 1,3-butadiene (PC1), benzene (PC2), and vinyl chloride (PC3), for the cumulative variance ranged from 24 to 58%, and the eigenvalues ranged from 1.003 to 1.449. The components were classified as sources from the synthetic rubber and thermoplastic resins industry and wastewater treatment, vehicle exhaust, and slurry, open equipment, vessel, and raw material in PVC industrial emissions, respectively (Chen et al., 2016; Hosaini et al., 2017; Malakan et al., 2022; Pinthong et al., 2022; Kim et al., 2024).
Health risk assessment of VOCs
Non-carcinogenic risk assessment for VOCs
Non-carcinogenic risks are health problems that occur in an organism due to VOC exposure, excluding cancer, which is calculated by the Hazard Quotient (HQ) in the study of inhalation routes to VOCs for children and adult groups, as shown in Table 6. In the Rayong City industrial estate, inhalation exposure to all VOCs posed no significant non-carcinogenic risk, as indicated by HQs that remained below 1.0. The highest of HQ values for 1,3-butadiene for children and adult groups, which showed no effects on non-carcinogenic risk.
Table 5. Rotated principal component loadings of VOC concentrations emitted from an industrial estate in Rayong city in 2024.
|
VOCs |
IA1 (n=363) |
|
IA2 (n=351) |
|
IA3 (n=196) |
||||||||||||
|
PC1 |
PC2 |
PC3 |
PC4 |
|
PC1 |
PC2 |
PC3 |
|
PC1 |
PC2 |
PC3 |
|
|||||
|
Vinyl Chloride |
0.747 |
0.258 |
-0.090 |
0.060 |
|
-0.015 |
0.957 |
-0.005 |
|
-0.027 |
-0.072 |
0.917 |
|
||||
|
1,3-Butadiene |
0.409 |
0.715 |
-0.078 |
-0.259 |
|
-0.019 |
0.000 |
0.982 |
|
0.648 |
-0.188 |
-0.220 |
|
||||
|
Benzene |
-0.061 |
-0.558 |
-0.504 |
0.215 |
|
-0.824 |
-0.035 |
-0.001 |
|
0.118 |
0.890 |
-0.064 |
|
||||
|
Toluene |
-0.093 |
0.238 |
0.496 |
0.763 |
|
0.678 |
0.195 |
-0.090 |
|
-0.305 |
0.453 |
0.031 |
|
||||
|
Ethylbenzene |
0.165 |
-0.324 |
0.666 |
-0.534 |
|
0.717 |
-0.274 |
-0.181 |
|
-0.776 |
0.020 |
-0.096 |
|
||||
|
p-Xylene |
0.666 |
-0.313 |
0.337 |
0.209 |
|
0.739 |
-0.103 |
0.181 |
|
0.474 |
0.111 |
0.358 |
|
||||
|
Eigenvalues |
1.376 |
1.149 |
1.071 |
1.028 |
|
2.217 |
1.034 |
1.027 |
|
1.449 |
1.032 |
1.003 |
|
||||
|
% of Variance |
22.940 |
19.144 |
17.847 |
17.130 |
|
36.950 |
17.229 |
17.119 |
|
24.148 |
17.192 |
16.725 |
|
||||
|
Cumulative % |
22.940 |
42.084 |
59.930 |
77.060 |
|
36.950 |
54.179 |
71.298 |
|
24.148 |
41.340 |
58.064 |
|
||||
|
Estimate sources |
Slurry, Open equipment, Vessel, and raw material in the PVC industry |
Synthetic rubber and thermoplastic resin production and Waste water treatment |
Vehicles exhaust
|
Vehicles exhaust
|
|
Mixed vehicles exhaust and Industrial solvent usage |
Slurry, Open equipment, Vessel, and raw material in the PVC industry |
Synthetic rubber and thermoplastic resin production and Waste water treatment |
|
Synthetic rubber and thermoplastic resin production and Waste water treatment |
Vehicles exhaust |
Slurry, Open equipment, Vessel, and raw material in the PVC industry |
|
||||
Table 6. Hazard quotient (HQ) associated with VOCs exposure affects both children and adults in an industrial estate located in Rayong city.
|
Site |
Age group |
Mean HQ value |
|||||
|
Vinyl chloride |
1,3-Butadiene |
Benzene |
Toluene |
Ethylbenzene |
P-Xylene |
||
|
IA1 |
1-5 years |
8.47×10-5 |
2.00×10-2 |
1.13×10-3 |
2.22×10-5 |
5.96×10-5 |
2.32×10-4 |
|
6-11 years |
1.27×10-4 |
3.00×10-2 |
1.70×10-3 |
3.33×10-5 |
8.94×10-5 |
3.47×10-4 |
|
|
12-19 years |
1.48×10-4 |
3.50×10-2 |
1.98×10-3 |
3.89×10-5 |
1.04×10-4 |
4.05×10-4 |
|
|
Adult (>20years) |
1.27×10-3 |
3.00×10-1 |
1.70×10-2 |
3.33×10-4 |
8.94×10-4 |
3.47×10-3 |
|
|
IA2 |
1-5 years |
1.50×10-5 |
1.54×10-2 |
1.74×10-4 |
9.74×10-6 |
3.02×10-4 |
3.64×10-4 |
|
6-11 years |
2.24×10-5 |
2.31×10-2 |
2.60×10-4 |
1.46×10-5 |
4.53×10-4 |
5.47×10-4 |
|
|
12-19 years |
2.62×10-5 |
2.70×10-2 |
3.04×10-4 |
1.71×10-5 |
5.28×10-4 |
6.38×10-4 |
|
|
Adult (>20years) |
2.24×10-4 |
2.31×10-1 |
2.60×10-3 |
1.46×10-4 |
4.53×10-3 |
5.47×10-3 |
|
|
IA3 |
1-5 years |
4.67×10-4 |
2.04×10-2 |
5.50×10-4 |
1.58×10-5 |
8.02×10-5 |
3.06×10-4 |
|
6-11 years |
7.00×10-4 |
3.07×10-2 |
8.25×10-4 |
2.37×10-5 |
1.20×10-4 |
4.59×10-4 |
|
|
12-19 years |
8.17×10-4 |
3.58×10-2 |
9.62×10-4 |
2.77×10-5 |
1.40×10-4 |
5.35×10-4 |
|
|
Adult (>20years) |
7.00×10-3 |
3.07×10-1 |
8.25×10-3 |
2.37×10-4 |
1.20×10-3 |
4.59×10-3 |
|
Carcinogenic risk assessment for VOCs
The cancer risk (CRInh) for inhalation exposure routes was estimated to assess the potential for cancer of VOCs on humans, including vinyl chloride, 1,3-butadiene, benzene, and ethylbenzene. Table 7 presents the cancer risk (CR) of VOCs for children and adults exposed within the industrial estate in Rayong city. The results of CRInh for inhalation exposure routes for children and adult groups calculated that the risk was below 10⁻⁶, indicating that the no risk from all VOCs related to cancer. The CRInh values of 1,3-butadiene for both age groups in the industrial estate of Rayong City were between 10⁻⁶ and 10⁻⁴, indicating that the cancer risk from VOCs was seen as acceptable over a lifetime.
Table 7. Cancer risk (CRInh) associated with VOCs exposure affects both children and adults in an industrial estate located in Rayong city.
|
Site |
Age group |
Mean CRInh value |
|||
|
Vinyl chloride |
1,3-Butadiene |
Benzene |
Ethylbenzene |
||
|
IA1 |
1-5 years |
1.90×10-8 |
1.20×10-6 |
2.65×10-7 |
1.49×10-7 |
|
6-11 years |
2.86×10-8 |
1.80×10-6 |
3.97×10-7 |
2.23×10-7 |
|
|
12-19 years |
3.33×10-8 |
2.10×10-6 |
4.63×10-7 |
2.61×10-7 |
|
|
Adult (>20years) |
2.86×10-7 |
1.80×10-5 |
3.97×10-6 |
2.23×10-6 |
|
|
IA2 |
1-5 years |
3.36×10-9 |
9.25×10-7 |
4.06×10-8 |
7.55×10-7 |
|
6-11 years |
5.04×10-9 |
1.39×10-6 |
6.09×10-8 |
1.13×10-6 |
|
|
12-19 years |
5.88×10-9 |
1.62×10-6 |
7.11×10-8 |
1.32×10-6 |
|
|
Adult (>20years) |
5.04×10-8 |
1.39×10-5 |
6.09×10-7 |
1.13×10-5 |
|
|
IA3 |
1-5 years |
1.05×10-7 |
1.23×10-6 |
1.29×10-7 |
2.00×10-7 |
|
6-11 years |
1.57×10-7 |
1.84×10-6 |
1.93×10-7 |
3.01×10-7 |
|
|
12-19 years |
1.84×10-7 |
2.15×10-6 |
2.25×10-7 |
3.51×10-7 |
|
|
Adult (>20years) |
1.84×10-7 |
2.15×10-6 |
2.25×10-7 |
3.51×10-7 |
|
The carcinogenic equivalency (CEQ) is a method of risk evaluation applied to evaluate the possible carcinogenicity of volatile organic compounds (VOCs) such as vinyl chloride, 1,3-butadiene, benzene, and ethylbenzene. Figure 2 illustrates that the CEQ values in the industrial estate, Rayong city, were 1.91 ± 1.21 µg/m³, 0.72 ± 0.41 µg/m³, and 1.64 ± 1.38 µg/m³ for IA1, IA2, and IA3, respectively. The highest CEQ value was found at the IA1 site, while the lowest values of CEQ were found in the IA2 site, as shown in Figure 3.

Figure 3. Carcinogenic equivalency (CEQ) of volatile organic compounds (VOCs) in an industrial estate located in Rayong city.
DISCUSSION
We performed an analysis of the average VOC concentrations emitted from Rayong City, Thailand's major petrochemical industrial estate. It was found that the concentrations of VOCs in all sampling sites in descending order were ethylbenzene (2.17 ± 4.48 to 11.0 ± 11.7 µg/m³)> toluene (1.78 ± 1.66 to 4.06 ± 4.42 µg/m³) >1,3-butadiene (1.13 ± 1.06 to 1.49 ± 1.91 µg/m³) > p-xylene (0.85 ± 0.36 to 1.33 ± 1.38 µg/m³) > benzene (0.60 ± 0.63 to 1.24 ± 0.43 µg/m³) > vinyl chloride (0.03 ± 0.03 to 0.87 ± 1.64 µg/m³). Moreover, the results revealed mean levels of vinyl chloride, benzene, and ethylbenzene were significantly different among their mean levels in each sampling site (P < 0.05). The total VOC concentrations during the dry season were not significantly different from the wet season in all sites (P > 0.05). The investigation of VOC emissions from the petrochemical industrial complex in Taiwan found that the amounts released each year were 117.7 tons of benzene, 44.6 tons of vinyl chloride, 42.4 tons of ethylbenzene, and 22.3 tons of 1,3-butadiene (Chen et al., 2016). Moreover, Pinthong et al. (2022) revealed that the VOCs emitted from the Maptaphut industrial area included benzene, total xylene, and ethylbenzene as main indicators of gasoline burning in automobiles, whereas industrial sources emitted toluene, styrene, vinyl chloride, and 1,3-butadiene. Similar to this study, Kim et al. (2024) showed that the volatile organic compounds (VOCs) emitted from the multi-industrial city of Ulsan, South Korea, included 18.3% toluene (6.37 µg/m³), 13.4% m, p-xylene (4.69 µg/m³), 7.7% ethylbenzene (2.70 µg/m³), and 4.8% benzene (1.69 µg/m³), indicating significant effects of industrial emissions. Malakan et al. (2022) reported that VOC emission inventories from industrial sources and mobile sources in the Maptaphut industrial estate included benzene, 1,3-butadiene, and vinyl chloride at 5.11, 3.19, and 37.78 tons/year, respectively, while mobile sources emitted benzene and 1,3-butadiene at 73.89 and 9.96 tons/year, respectively. In 2024, the ambient concentrations of 1,3-butadiene (1.13 ± 1.06 to 1.49 ± 1.91 µg/m³) at three sites in the industrial estate of Rayong city exceeded the annual standard (0.33 µg/m³) set by the Pollution Control Department (PCD). In a similar study by Nakyai et al. (2025), concentrations of 1,3-butadiene emitted from industrial estates in Rayong Province during the rainy and cold seasons exceeded the annual standards of PCD.
Moreover, the result of the association between VOC concentrations and meteorological parameters at the Rayong industrial estate in 2024 was obtained using Pearson’s correlation. Toluene had a substantial negative connection with relative humidity (r = -0.166 to -0.266) and air temperature (r = -0.263 to -0.315) across all sites. Volatile organic compound concentrations decreased with increasing relative humidity and air temperature. The increase in temperature enhances vertical mixing and the overall dispersion of the pollutants. Furthermore, increased relative humidity enhances the wet scavenging or washout effect, facilitating the removal of VOCs from the atmosphere. Moreover, wind speed and total rainfall were increased to reduce VOCs in ambient air (Nakyai et al., 2025).
Benzene and ethylbenzene were positively correlated with air temperature and relative humidity in some sites. The correlation between benzene and atmospheric pressure in the IA2 site was strong and significantly positive (r = 0.713), while benzene and air temperature (r = 0.220) and wind speed (r = 0.253) were positively correlated. The IA2 site is located close to the high traffic density, various government buildings, and the oil refineries. The main results are that BTEX compounds are hydrophobic compounds whose adsorption effects and solubility decrease in water vapor. Consequently, elevated temperature and high relative humidity increased benzene and ethylbenzene levels in the ambient air. Meteorological circumstances influence the various components due to their high Henry's law constant, low water solubility, and low polarity (Xu and Zhang, 2011; Zhou et al., 2017).
The PCA evaluating the sources of VOCs in the industrial estate in Rayong City indicated that the IA1 site comprised vinyl chloride and p-xylene for PC1, signifying the presence of petrochemical companies, including slurry, open equipment, vessels, and raw material in PVC industrial emissions. The IA1 site is located close to the community, approximately 2.53 km from the petrochemical factory, which mainly manufactures plastic resins or polymers that are molded into products used for food packaging, automotive parts, medical instruments, electrical devices, and infrastructure, such as pressure-resistant pipes and telecommunication cables, among others. Moreover, Malakan et al. (2022) revealed that industrial sources of VOC emissions, including stack, wastewater treatment, slurry, open equipment, and vessels, contributed 37.78 tons/year from the Maptaphut petroleum and petrochemical industrial estates. Sass et al. (2005) and Pinthong et al. (2022) confirmed that industrial sources released toluene, styrene, vinyl chloride, and 1,3-butadiene from using raw material in the PVC industry. PC2, PC3, and PC4 indicated 1,3-butadiene, ethylbenzene, and toluene, respectively. The sources were classified as synthetic rubber and thermoplastic resins, and vehicular emissions. Chen et al. (2016) reported that the source of 1,3-butadiene is emissions from plastic and rubber, nitrile rubber, and acrylonitrile copolymer resins. Sukjit et al. (2021) reported that 1,3-butadiene in ambient air at Maptaphut estates emitted from wastewater treatment plants was about 92% from the synthetic rubber industry. Moreover, ethylbenzene and toluene are released from industrial activity, vehicular exhaust, and gasoline evaporation (Hosaini et al., 2017; Kim et al., 2024).
Additionally, PC1 for the IA2 site was predominantly composed of toluene, ethylbenzene, and p-xylene, implying the possibility of mixed vehicle exhaust and industrial solvent usage (Malakan et al., 2022; Pinthong et al., 2022). The IA2 site is located close to the public health center and a nearby government building. It is adjacent to the roadway. Moreover, this sampling site was located 1.38 km from the oil refineries. Malakan et al. (2022) reported that mobile sources, including motorcycles, passenger cars, vans and pickups, trucks, and buses, emitted 73.39 tons/year of benzene and 9.96 tons/year of 1,3-butadiene from the Maptaphut petroleum and petrochemical industrial estates. Similar research from Pinthong et al. (2022) revealed that the high VOC concentrations emitted from Maptaphut petrochemical industrial areas in Thailand were found to be benzene, total xylene, and ethylbenzene emitted from gasoline combustion in mobile motors, while industrial sources released toluene, 1,3-butadiene, styrene, 1,2-dichloroethane, and vinyl chloride. Kim et al. (2024) investigated the source of VOCs released from the multi-industrial city of Ulsan in South Korea using the positive matrix factorization model (PMF). They found that toluene, ethylbenzene, benzene, m, p, o-xylene, ethyl acetate, and vinyl acetate were related to petrochemical and non-ferrous industrial activity and industrial solvent usage. The compounds PC2 and PC3 were vinyl chloride and 1,3-butadiene, which are classified as belonging to the oil refinery and slurry, open equipment, and vessels, as well as the synthetic rubber and thermoplastic resins industry and wastewater treatment (Chen et al., 2016; Sukjit et al., 2021). Moreover, the IA3 site indicated that three factors, including 1,3-butadiene (PC1), benzene (PC2), and vinyl chloride (PC3), were classified as sources from the synthetic rubber and thermoplastic resins industry and wastewater treatment, vehicle exhaust, and the slurry, open equipment, vessels, and raw material in PVC industrial emissions, respectively. The IA3 site is in the school and community area. The location is approximately 600 meters from the liquefied natural gas (LNG) industry and close to a heavily trafficked intersection during peak times. The polymerization of the polyvinyl chloride (PVC) industry is emitted from vinyl chloride (Chen et al., 2016; Pinthong et al., 2022), while 1,3-butadiene is related to plastic and rubber industrial processes (Chen et al., 2016; Sukjit et al., 2021) and oil refineries (Pinthong et al., 2022). Hosaini et al., (2017) reported that principal component analysis (PCA) identified toluene and benzene as the main sources of ambient VOCs in Kuala Lumpur, Malaysia, which were associated with gasoline evaporation and vehicle exhaust emissions.
Non-carcinogenic risk assessments in the study of inhalation routes to VOCs for children and adult groups in the industrial estate of Rayong City are calculated by the Hazard Quotient (HQ). It was found that the inhalation exposure to all VOCs posed no significant non-carcinogenic risk, as indicated by HQs that remained below 1.0, which were no effects of non-carcinogenic risk. Moreover, the cancer risk (CRInh) of VOCs attached to children and adults exposed to the industrial estate, Rayong City. The results of CRInh for inhalation exposure routes for children and adult groups calculated that the risk was below 10⁻⁶, indicating that there is no risk from all VOCs related to cancer. However, the CEQ values in the industrial estate, Rayong city, were 1.91 ± 1.21 µg/m³, 0.72 ± 0.41 µg/m³, and 1.64 ± 1.38 µg/m³ for IA1, IA2, and IA3, respectively. The highest CEQ value was obtained at the IA1 site, which contained high levels of carcinogenic VOCs, including vinyl chloride, 1,3-butadiene, and benzene, associated with the petrochemical industries primarily producing plastic resins or polymers. Although both non-carcinogenic and carcinogenic risks associated with VOCs have shown no adverse impact on human health, previous studies found that exposure to benzene and 1,3-butadiene significantly decreases DNA repair capacity in children, with both pollutants causally linked to DNA damage and the disruption of this essential cellular function. Consequently, these exposure occasions could increase the risk of cancer (Au et al., 2010; Ruchirawat et al., 2010). Exposure to vinyl chloride may be associated with Raynaud's syndrome, acral sclerodermatous transformations, acroosteolysis (AOL), skin cancer, and cutaneous signs of liver disease (Goodman et al., 2023). Moreover, the occupational exposure to 1,3-butadiene is linked to leukemia and cardiovascular disease (Lin et al., 2020; Chen and Zhang, 2022). Furthermore, current evidence indicates that nonoccupational exposure to 1,3-butadiene may be linked to specific reproductive effects and, notably, to various children developing tumors, autism, and asthma (von Ehrenstein et al., 2014; Hall et al., 2019; Kuang et al., 2021). This study focuses on 1,3-butadiene because of its long-term cancer risk.
CONCLUSION
This study performed an analysis of the average VOC concentrations emitted from the petrochemical industrial estate of Rayong city in Thailand. The results revealed mean levels of vinyl chloride, benzene, and ethylbenzene were significantly different among their mean levels in each sampling site. The total VOC concentrations during the dry season were not significantly different from the wet season in all sites. Moreover, the ambient concentrations of 1,3-butadiene (1.13 ± 1.06 to 1.49 ± 1.91 µg/m³) in the industrial estate of Rayong city exceeded the annual standard (0.33 µg/m³) set by the Pollution Control Department (PCD). Although the level of vinyl chloride (0.60 ± 0.63 to 1.24 ± 0.43 µg/m³) and benzene (0.03 ± 0.03 to 0.87 ± 1.64 µg/m³) remained below the particular annual standards, their continued distribution in the atmosphere requires consideration regarding potential human health effects, specifically increased cancer risk. The main effect of meteorological conditions on VOCs might be air temperature and relative humidity. The PCA evaluated the sources of VOCs in the industrial estate in Rayong City, which were classified as sources from the synthetic rubber and thermoplastic resins industry, wastewater treatment, vehicle exhaust, industrial solvent usage, and the petrochemical industry, including slurry, open equipment, vessels, and raw material in PVC industrial emissions. Moreover, the highest CEQ value, which contained high levels of carcinogenic VOCs, might be associated with the petrochemical industries primarily producing plastic resins or polymers and wastewater treatment. However, both non-carcinogenic and carcinogenic risks associated with VOCs indicated no adverse impact on human health. The conclusion clarified that health risks assessed account for exposure duration, frequency, and toxicological potency, which moderate the potential health impact. Elevated concentrations require the risk of dose accumulation from exposure conditions to cause significant adverse health effects. This study focuses on 1,3-butadiene because of its long-term cancer risk.
The study has several limitations. Geographic coverage was restricted to only three monitoring sites, and the assumptions used in the risks modelling may have introduced data gaps. Future studies should expand monitoring to additional areas, incorporate personal exposure assessments, refine seasonal source apportionment, and investigate the potential synergistic health effects of combined VOC exposures. Moreover, optimizing air quality in industrial areas requires a synergistic approach encompassing continuous monitoring, implementation of targeted source controls (with emphasis on high-risk compounds, e.g., 1,3-butadiene), and strategic health risk communication directed at informing relevant policymakers and stakeholders.
ACKNOWLEDGEMENTS
The authors are grateful to the Faculty of Science, Energy, and Environment, King Mongkut’s University of Technology North Bangkok. King Mongkut's University of Technology North Bangkok funded this research under Contract no. KMUTNB-68-BASIC-48. The authors would like to thank Maptaphut City municipality and Thailand’s Pollution Control Department (PCD) for providing VOCs data.
AUTHOR CONTRIBUTIONS
Susira Bootdee: Conceptualization (Lead), Methodology (Lead), Data Curation (Lead), Data Analysis (Equal), Investigation (Lead), Resources (Lead), Writing – Original Draft (Lead), Writing – Review & Editing (Lead), Visualization (Lead), Supervision (Lead), Project Administration (Lead); Suganya Phantu: Data Curation (Supporting), Visualization (Supporting), Statistical & Data Analysis (Equal). All authors have read and approved the final manuscript.
CONFLICT OF INTEREST
The authors declare that they hold no competing interests.
REFERENCES
Abd Hamid, H.H., Latif, M.T., Mohd Nadzir, M.S., Uning, R., Khan, M.F., and Kannan, N. 2019. Air quality. Atmosphere & Health. 12: 341-351. https://doi.org/10.1007/s11869-019-00664-1
Allahabady, A., Yousefi, Z., Mohammadpour Tahamtan, R.A., and Sharif, Z.P. 2022. Measurement of BTEX (benzene, toluene, ethylbenzene, and xylene) concentration at gas stations. Environmental Health Engineering and Management Journal. 9(1): 23-31. https://doi.org/10.34172/EHEM.2022.04
Au, W.W., Giri, A.K., and Ruchirawat, M. 2010. Challenge assay: A functional biomarker for exposure-induced DNA repair deficiency and for risk of cancer. International Journal of Hygiene and Environmental Health. 213: 32-39. https://doi.org/10.1016/j.ijheh.2009.09.002
Berenjian, A., Chan, N., and Malmiri, H.J. 2012. Volatile organic compounds removal methods: A review. American Journal of Biochemistry and Biotechnology. 8(4): 220-229. https://doi.org/10.3844/ajbbsp.2012.220.229
Chen, M.-J., Lin, C.-H., Lai, C.-H., Cheng, L.-H., Yang, Y.-H., Huang, L.-J., Yeh, S.-H., and Hsu, H.-T. 2016. Excess lifetime cancer risk assessment of volatile organic compounds emitted from a petrochemical industrial complex. Aerosol and Air Quality Research. 16: 1954-1966. https://doi.org/10.4209/aaqr.2015.05.0372
Chen, W.-H., Chang C.-F., Lai, C.-H., Peng, Y.-P., Su, Y.-J., and Chen, G.-F. 2024. Multivariate analysis of carcinogenic equivalence (CEQ) to characterize carcinogenic VOC emissions in a typical petrochemical industrial park in Taiwan. Environment International. 186: 108548. https://doi.org/10.1016/j.envint.2024.108548
Chen, W.-Q. and Zhang, X.-Y. 2022. 1,3-Butadiene: A ubiquitous environmental mutagen and its associations with diseases. BMC Genes and Environment. 44: 3. https://doi.org/10.1186/s41021-021-00233-y
Goodman, R.S., Mittal, L., and Parker, E.R. 2023. Public health risks, dermatological manifestations, and environmental justice associated with vinyl chloride exposure: Narrative review. JMIR Dermatology. 6: e48998. https://doi.org/10.2196/48998
Hosaini, P.N., Khan, M.F., Hamizah Mustaffa, N.I., Amil, N., Mohamad, N., Jaafar, A., Mohd Nadzir, M.S., and Latif, M.T. 2017. Concentration and source apportionment of volatile organic compounds (VOCs) in the ambient air of Kuala Lumpur, Malaysia. Natural Hazards. 85: 437-452. https://doi.org/10.1007/s11069-016-2575-7
Hall, C., Heck, J.E., Ritz, B., Cockburn, M., Escobedo, L.A., and von Ehrenstein, O.S. 2019. Prenatal exposure to air toxics and malignant germ cell tumors in young children. Journal of Occupational and Environmental Medicine. 61(6): 529-534. https://doi.org/10.1097/JOM.0000000000001609
Kawichai, S., Prapamontol, T., Chantara, S., Kanyanee, T., Wiriya, W., and Zhang, Y.-L. 2020. Seasonal variation and sources estimation of PM2.5-bound PAHs from the ambient air of Chiang Mai city: An all-year-round study in 2017. Chiang Mai Journal of Science. 47(5): 958-972.
Kawichai, S., Bootdee, S., Sillapapiromsuk, S., and Janta, R. 2022. Epidemiological study on health risk assessment of exposure to PM2.5-bound toxic metals in the industrial metropolitan of Rayong, Thailand. Sustainability. 14: 15368. https://doi.org/10.3390/su142215368
Kim, S.-J., Lee, S.-J., Kim, H., Hong, Y., and Choi, S.-D. 2024. Contribution of individual sources of volatile organic compounds to their cancer and non-cancer risks in the multi-industrial city of Ulsan, South Korea. Air Quality, Atmosphere & Health. 17: 1937-1949. https://doi.org/10.1007/s11869-024-01554-x
Kuang, H., Li, Z., Lv, X., Wu, P., Tan, J., Wu, Q., Li, Y., Jiang, W., Pang, Q., Wang, Y., and Fan, R. 2021. Exposure to volatile organic compounds may be associated with oxidative DNA damage-mediated childhood asthma. Ecotoxicology and Environmental Safety. 210: 111864. https://doi.org/10.1016/j.ecoenv.2020.111864
Lin, C.-Y., Lee, H.-L., Jung, W.-T., Sung, F.-C., and Su, T.C. 2020. The association between urinary levels of 1,3-butadiene metabolites, cardiovascular risk factors, microparticles, and oxidative stress products in adolescents and young adults. Journal of Hazardous Materials. 396: 12275. https://doi.org/10.1016/j.jhazmat.2020.122745
Lv, Z., Liu, X., Bai, H., Nie, L., and Li, G. 2024. Process-specific volatile organic compounds emission characteristics, environmental impact, and health risk assessments of the petrochemical industry in the Beijing-Tianjin-Hebei region. Environmental Science and Pollution Research. 31: 3938-3950. https://doi.org/10.1007/s11356-023-31351-5
Mabel, O.A. and Olayemi, O.S. 2020. A comparison of principal component analysis, maximum likelihood, and the principal axis in factor analysis. American Journal of Mathematics and Statistics. 10(2): 44-54.
Malakan, W., Thepanondh, S., and Kondo, A. 2022. Modeling of inhalation health risk of volatile organic compounds in the vicinity of Maptaphut petroleum and petrochemical industrial estate, Thailand. Sustainability. 14: 12073. https://doi.org/10.3390/su141912073
Maskey, R., Fei, J., and Nguyen, H.-O. 2018. Use of exploratory factor analysis in maritime research. The Asian Journal of Shipping and Logistics. 34(2): 91-111. https://doi.org/10.1016/j.ajsl.2018.06.006
Mihajlović, V., Grba, N., Sudi, J., Eichert, D., Krajinović, S., Gavrilov, M.B., and Marković, S.B. 2021. Assessment of occupational exposure to BTEX in a petrochemical plant via urinary biomarkers. Sustainability. 13: 7178. https://doi.org/10.3390/su13137178
Morakinyo, O.M., Adebowale, A.S., Mokgobu, M.I., and Mukhola, M.S. 2017. Health risk of inhalation exposure to sub-10 µm particulate matter and gaseous pollutants in an urban-industrial area in South Africa: An ecological study. Occupational and Environmental Medicine. 7(3): 013941. https://doi.org/10.1136/bmjopen-2016-013941
Nakyai, T., Polyong, C.P., Kongsombatsuk, M., Boonnuk, T., and Thetkathuek, A. 2025. Seasonal impact and meteorological factors affecting the distribution of volatile organic compounds concentrations and health risk assessments inside and outside industrial estates: A case study of Rayong province, Thailand. Case Studies in Chemical and Environmental Engineering. 11: 101121. https://doi.org/10.1016/j.cscee.2025.101121
Odekanle, E.L., Olubiyo, F.W., Sonibare, O.O., Ogunlaja, O.M., Aremu, C.O., Ojediran, J.O., and Fakinle, B.S. 2022. Health impact assessment of volatile organic compounds (VOCs) emission from the combustion of agricultural wastes. Cogent Engineering. 9: 213049. https://doi.org/10.1080/23311916.2022.2143049
Prasertsin, S. and Nathapindhu, G. 2020. The temporal pattern of ambient PM2.5 and health risk assessment in Thailand. Indian Journal of Public Health Research & Development. 11(03): 1096-1100.
Pinthong, N., Thepanondh, S., and Kondo, A. 2022. Source identification of VOCs and their environmental health risk in a petrochemical industrial area. Aerosol and Air Quality Research. 22(2): 210064. https://doi.org/10.4209/aaqr.210064
Pollution Control Department (PCD). 2025. Standard of annual VOCs in ambient air. [Cited 2025 Jul 8]. Available from: https://www.pcd.go.th/laws/2810/.
Ruchirawat, M., Navasumrit, P., and Settachan, D. 2010. Exposure to benzene in various susceptible populations: Co-exposures to 1,3-butadiene and PAHs and implications for carcinogenic risk. Chemico-Biological Interactions. 184: 67-76. https://doi.org/10.1016/j.cbi.2009.12.026
Saeaw, N. and Thepanondh, S. 2015. Source apportionment analysis of airborne VOCs using positive matrix factorization in industrial and urban areas in Thailand. Atmospheric Pollution Research. 6: 644-650. https://doi.org/10.5094/APR.2015.073
Saeedi, M., Malekmohammadi, B., and Tajalli, S. 2024. Interaction of benzene, toluene, ethylbenzene, and xylene with the human body: Insights into characteristics, sources, and health risks. Journal of Hazardous Materials Advances. 16: 100459. https://doi.org/10.1016/j.hazadv.2024.100459
Saleem, I., Ali, Z., and Hussain, A. 2022. Assessment of indoor air pollutants exposure in hair salons of Lahore. Chiang Mai University Journal of Natural Science. 21(4): e2022065. https://doi.org/10.12982/CMUJNS.2022.065
Sass, J.B., Castleman, B., and Wallinga, D. 2005. Vinyl chloride: A case study of data suppression and misrepresentation. Environmental Health Perspectives. 113(7): 809-812. https://doi.org/10.1289/ehp.7716
Sharma, S., Giri, B., and Patel, K.S. 2016. Ambient volatile organic compounds in the atmosphere of industrial central India. Journal of Atmospheric Chemistry. 73: 381-395. https://doi.org/10.1007/s10874-016-9329-5
Sukjit, N., Thepanondh, S., Kwonpongsagoon, S., Jinsart, W., Punya, L., and Hanma, P. 2021. Source concentration of 1,3 butadiene in the vicinity of a petrochemical industrial area. Applied Environmental Research. 43(4): 29-38. https://doi.org/10.35762/AER.2021.43.4.3
Suwattiga, P. and Limpaseni, W., 2005. Seasonal source apportionment of volatile organic compounds in Bangkok ambient air. ScienceAsia. 31, 395-401. https://doi.org/10.2306/scienceasia1513-1874.2005.31.395
Teravecharoenchai, J., Chalermchaiwat, P., and Thuwapanichayanan, R. 2021. Principal component analysis application on nutritional, bioactive compound, and antioxidant activities of pigmented dough grain. Chiang Mai University Journal of Natural Science. 20(2): e2021036. https://doi.org/10.12982/CMUJNS.2021.036
United States Environmental Protection Agency (US-EPA). 2007. An introduction to indoor air quality (IAQ), volatile organic compounds (VOCs). Washington DC: USEPA, 2007. [Cited 2025 Jul 8]. Available from: https://www.epa.gov/indoor-air-quality-iaq/volatile-organic-compounds-impact-indoor-air-quality.
Uttamang, P., Janta, R., Bran, S.H., Macatangay, R., Surapipith, V., Tala, W., and Chantara, S. 2023. Effects of biogenic volatile organic compounds and anthropogenic NOx emissions on O3 and PM2.5 formation over the northern region of Thailand. Frontiers in Environmental Science. 11: 1146437. https://doi.org/10.3389/fenvs.2023.1146437
von Ehrenstein, O.S., Aralis, H., Cockburn, M., and Ritz, B. 2014. In utero exposure to toxic air pollutants and risk of childhood autism. Epidemiology. 25(6): 851-858. https://doi.org/10.1097/EDE.0000000000000150
Wang, Z., Zhang, P., Pan, L., Qian, Y., Li, Z., Li, X., Guo, C., Zhu, X., Xie, Y., and Wei, Y. 2023. Ambient volatile organic compound characterization, source apportionment, and risk assessment in three megacities of China in 2019. Toxics. 11: 651. https://doi.org/10.3390/toxics11080651
Xu, J. and Zhang, J.S. 2011. An experimental study of the relative humidity effect on VOCs' effective diffusion coefficient and partition coefficient in a porous medium. Building and Environment. 46(9): 1785-1796. https://doi.org/10.1016/j.buildenv.2011.02.007
Zhang, G., Wang, N., Jiang, X., and Zhao, Y. 2017. Characterization of ambient volatile organic compounds (VOCs) in the area adjacent to a petroleum refinery in Jinan, China. Aerosol and Air Quality Research. 17: 944-950. https://doi.org/10.4209/aaqr.2016.07.0303
Zheng, H, Chen, D., Qian, D., Meng, K., Hu, C., Li, S., Feng, X., Xu, H., Peng, Q., Yu, H., et al. 2023. Source analysis and control strategy of VOCs and PM2.5 in the high ozone season in the industrial zone of Shaoxing city, China. Water, Air, & Soil Pollution. 234: 338. https://doi.org/10.1007/s11270-023-06324-6
Zhou, Y., Wang, K., Wang, B., Pu, Y., and Zhang, J. 2020. Occupational benzene exposure and the risk of genetic damage: A systematic review and meta-analysis. BMC Public Health. 20: 1113. https://doi.org/10.1186/s12889-020-09215-1
Zhou, C., Zhan, Y., Chen, S., Xia, M., Ronda, C., Sun, M., Chen, H., and Shen, X. 2017. Combined effects of temperature and humidity on indoor VOCs pollution: Intercity comparison. Building and Environment. 121: 26-34. https://doi.org/10.1016/j.buildenv.2017.04.013
OPEN access freely available online
Natural and Life Sciences Communications
Chiang Mai University, Thailand. https://cmuj.cmu.ac.th
Susira Bootdee1, * and Suganya Phantu2
1 Chemical Industrial Process and Environment program, Faculty of Science, Energy and Environment, King Mongkut’s University of Technology North Bangkok (Rayong campus), Rayong 21120, Thailand.
2 Data Science and Computation for Business and Industry program, Faculty of Science, Energy and Environment, King Mongkut’s University of Technology North Bangkok (Rayong campus), Rayong 21120, Thailand.
Corresponding author: Susira Bootdee, E-mail: susira.b@sciee.kmutnb.ac.th
ORCID iD:
Susira Bootdee: https://orcid.org/0000-0003-4486-9544
Suganya Phantu: https://orcid.org/0000-0003-1446-5610
Total Article Views
Editor: Sirasit Srinuanpan,
Chiang Mai University, Thailand
Article history:
Received: September 25, 2025;
Revised: December 9, 2025;
Accepted: December 16, 2025;
Online First: January 8, 2026