Exposure assessment of elemental carbon, ultrafine particles, and crystalline silica at highway toll booths
Article information
Abstract
Highway toll booth workers have been reported to be at an increased risk of occupational lung cancer. Moreover, insufficient studies have been performed on exposure assessment of workers at highway toll booths. Elemental carbon (EC), black carbon (BC), and respirable crystalline silica (RCS) concentrations were measured at highway toll booths in Gyeongsangbuk-do (Republic of Korea). The particle number (PN) concentration and size distribution of ultrafine particle were measured using three SMPS devices. The average concentration of the EC inside the booth was 3.3 μg/m3, with the maximum being 5.8 μg/m3. The concentrations of EC were highest for booths that operated exclusively for trucks. The average PN concentration inside the booth was 3.54 × 104 cm−3, approximately 5-fold higher than the reference indoor background. The average BC concentration in the booths were approximately 8 μg/m3, and the instantaneous peak concentration was 271 μg/m3. The RCS was below the detection limit in all samples. This study revealed that toll workers were most frequently exposed to diesel engine particle less than 100 nm, with an extremely high respiratory deposition rate. Therefore, a respiratory protection program is necessary to safeguard these workers against vehicle-related pollutants.
1. Introduction
As of 2017, the Korea Expressway Corporation, which manages the toll road network in South Korea, had 30 operating routes with a total length of 4,113 km (41 km in 10 lanes, 434 km in 8 lanes, 503 km in 6 lanes, and 3,135 km in 4 lanes) and 604 km in 8 private sector routes.
According to the 2015 data, 5,495 workers were employed at 369 toll stations (335 on the Korea Expressway Corporation and 34 on private roads) in Korea [1].; further, according to 2017 data [2], Roughly 8,000 workers operate in 420 toll stations, most of whom are women. Demographics indicate that most workers are female. Henschke and Miettinen [3] show that women are at a greater risk of lung cancer, rendering the studies on the exposure of workers to various pollutants and carcinogens particularly important. The nationwide count of toll gates includes highways operated by private companies such as the Gyeongchun Expressway, and the general highway operated by the Korea Highway Corporation, and the Namsan 1 and 3 tunnels as well as the Woomyeonsan tunnel in Seoul.
Workers at toll booths are exposed to the smoke generated from numerous vehicles that transit through these tolls throughout the day; further, they are also exposed to the atmospheric particlates matter (PM).
In 2012, the International Agency for Research on Cancer (IARC) designated diesel engine emission as a Group 1 carcinogen, indicating that there is sufficient evidence of its carcinogenicity in humans. PM was also classified as a Group 1 carcinogen in 2013 [4].
The substances emitted from vehicle engines are mixtures of particulate and gaseous substances from incompletely burned fuels and lubricants. Gasoline engines produce higher emissions of poly-cyclic aromatic hydrocarbons (PAHs), carbon monoxide, and carbon dioxide than diesel engines, but exhibit significantly lower emissions of nitrogen oxides and PM [5]. The main components of diesel engine particle (DEP) are elemental carbon (EC) and PAHs. In diesel engine exhaust, a maximum of three benzene rings are generated, alongside low molecular weight PAHs. In gasoline engine exhausts, benzo[a] pyrene and dibenz[a,h]anthracene are generated [6, 7].
Diesel engine exhaust (DEE) is harmful to the human body because of the toxic complex components of many gaseous and particulate matters including PAHs. In addition, the deposition of particulate matter in the respiratory system depends on the particle size. This is explained by the deep penetration and high rate of respiratory deposition (cited in archive toxicity). The particle size of DEP varies depending on engine load, which is determined by vehicle acceleration, the ratio of engine running and stopping, and catalyst use, amongst other factors; however, because the geometric mean diameter of DEP is very small (less than 100 nm), DEP gives rise to concerns about health effects by nanoparticles [5, 7, 8].
Incomplete combustion of diesel fuel in engines is known to generate reactive and resonance-stabilized radicals, along with hydrocarbons such as acetylene and PAHs. Subsequently, these radicals undergo chain reactions to form larger and more stable radicals. These radicals produce a nuclei with a diameter of 1 to 4 nm and the hydrocarbons and radicals in the vicinity are quickly attached to the nuclei, forming a soot nanoparticle with a diameter of 10 to 50 nm. In addition, it is known that hydrocarbons such as PAH, nitrogen oxides, sulfur oxides, and heavy metals adhere to the surface of these soot nanoparticles, and these soot nanoparticles are also aggregated to become larger particles [9]. On the other hand, the lung deposition rate of particles depends on the aerosol concentration in the ambient environment and the mode diameter of the particles. Additionally, the lung deposition particle fraction increases with decreasing diameter of the main mode of a particle size distribution [10]. Therefore, when evaluating exposure to DEP, it is possible to obtain a more in-depth exposure evaluation by analyzing the EC component and the size distribution as well as the number concentration of UFP [10–12]. Elemental carbon are thermally stable carbonaceous particles and normally quantified by an evolved gas analysis such as an OCEC Analyzer. Black carbon is similar to elemental carbon in that both are carbonaceous particles produced from the combustion process in flames and engines. Black carbon is usually measured using a filter-based light absorption technique or photothermal spectroscopy [5].
Ultrafine particles (UPFs) generally refer to particles smaller than 100 nm [13] and account for most of the particle number (PN) concentration measured by a Scanning Mobility Particle Sizer (SMPS) device. Road vehicles are the primary source of UFP emissions in urban pollution hotspots [14].
Toll booth workers are repeatedly exposed to vehicle emission as automobiles slow down and accelerate in front of the toll booth. In addition, booth windows remain open for the majority of the working hours; therefore, workers experience continuous exposure to vehicle emissions and PM in the atmosphere.
Recently, there have been increasing cases of occupational lung cancer due to DEE exposure, and there are concerns about the relationship between PM and diseases amongst outdoor workers, such as municipal household waste workers [15], shipyard transporter signal workers [16], and construction laborers [17]. Previous studies on the exposure of toll workers to hazardous materials have been primarily conducted on the PM2.5 and PM10 [18]. Studies on PN concentration and size distribution have been rarely reported; in particular, few studies have investigated the concentration of respirable crystalline silica (RCS) as a carcinogen for lung cancer [19].
In addition, according to Jancsek-Turóczi et al. [20], the mineral phase composition of resuspended road dust at three kerbside locations in Veszprém, Hungary contained 10 to 20% crystalline silica (as quartz). In another study by Gunawardana et al. [21] the road dust consist primarily of soil derived minerals (60%) with quartz averaging 40–50% and remainder being clay forming minerals of albite, microcline, chlorite and muscovite originating from surrounding soils. Therefore, the toll gate workers may be at risk of exposure to crystalline silica contained in resuspended dust on the road.
Therefore, the objectives of this study are to evaluate the exposure level of RCS and DEP on workers in toll booths through the level of surrogates, EC, black carbon (BC), and PN concentration, and to measure the particle size distribution of airborne UFP.
2. Material and Methods
2.1. Study Sites
This study conducted environmental measurements for DEP and RCS that the toll workers were exposed to at five booths of two tollgates, denoted as Gate A and Gate B, on the Gyeongbu Line highway on February 2, 2015. Interviews with workers and employers were conducted to obtain information on the traffic volume, lane characteristics, and air curtain usage. Workers were allowed to work in three shifts, eight hours a day, which were scheduled as follows: first shift from 06:30 to 14:30, middle shift from 14:30 to 22:00, and last shift from 22:00 to 06:30.
Measurements of the air quality at Gate A were performed both inside and outside for a single exit booth, inside the other three exit booths, and at the entrance gate (boothless gate). Measurements were also made inside and outside the office building environment at Gate A. At Gate B, measurements were performed at the inside of the four exit booths and at the entrance gate, similar to Gate A.
The toll gate entrance is divided into two vehicle types – freights and general vehicles. General vehicles include general passenger cars, vans with less than 16 passengers, freight cars with less than 2.5 tons, and light-weight vehicles.
The measurement position inside the booth was close to the breathing position of the workers, and the measurement positions outside the booth and outside the office building were approximately 1.0 m above the ground (Fig. 1). The PN concentration as a reference was also measured at inside of our research institute located in Incheon.
2.2. Measurement and Analysis Methods
We used 37 mm quartz filters preloaded in three-piece clear plastic cassettes (Cat.No. 225–401, SKC Inc., USA) with respirable dust cyclone (GK2.69, BGI Inc., Waltham, MA, USA) using a sampling pump (Aircheck XR5000, SKC Inc., Eighty Four, Pa., USA). The GK2.69 cyclone has a 50% cut point of 4 μm at a flow rate of 4.2 L/min. The EC concentrations were obtained from the samples using a Lab OCEC Analyzer (Sunset Laboratory Inc., USA) according to the National Institute for Occupational Safety and Health (NIOSH) Method 5040 [22]. The analysis accuracy was confirmed through the analysis of external standard samples (Sunset Laboratory Inc., USA). The limit of detection of the analyzer was 0.15 μg/cm2, and all sample measurements for this study exceeded the detection limit. LOD was calculated directly in our laboratory using the NIOSH guidelines for air sampling and analytical method development and evaluation [23].
The RCS concentration was analyzed according to the NIOSH 7500 guidelines [24]. The repirable dust samples were collected on the polyvinyl chloride filters (diameter 37 mm, pore 5.0 μm) using GK2.69 cyclone in combination with AirChek XR5000 pump. The samples were analyzed by X-ray diffraction (D8 Advance, Bruker Corporation, Ettlingen, Germany).
The analysis took place at the Institute of Occupation and Environment, Korea Workers’ Compensation and Welfare Service, participating in the American Industrial Hygiene Association Proficiency Analytical Testing program.
Realtime, continuous measurement BC mass concentration was measured using the four MicroAeth (AE51, TSI Inc., US) at a total of 4 locations inside and outside booth 1 and outside of the office building at Gate A and inside of booth 1 at Gate B. Data was recorded from approximately 7:30 to 14:00 at Gate A, and from approximately 8:55 to 14:30 at Gate B. The measurement interval for the MicroAeth was 30 s, and the pump flow rate was 0.1 L/min. All BC instruments were replaced with new filters and then measured.
For the measurement of the realtime number concentration and particle size distribution of the PM at both study sites, an SMPS (Model 3910, TSI Inc., US) was used. Particles ranging in size from 0.01 to 0.42 μm were divided into 13 channels. Three SMPSs were used to perform measurements at three locations, which were inside and outside the toll booths at Gate A and outside office building A. The SMPS is calibrated annually by the manufacturer (US TSI Inc.) to maintain the accuracy of the results, and the particle density of all three devices was set to 1.2 g cm−3.
2.3. Meteorological Data
Meteorological factors, including temperature (°C), humidity (%), wind speed (km/h), and wind direction (°), all of which can affect toll booths, were supplied by the Automatic Weather System of the Korea Meteorological Administration. We documented temperature, humidity, wind speed, and the wind direction each minute from 7:30 to 14:00 on the measurement day. The arithmetic means of temperature, humidity, and wind speed were used for statistical analysis to confirm their effects.
2.4. Data Analysis
The minimum, median, and maximum values of PN concentration (SMPS) and BC concentration were calculated using descriptive statistics. A distribution of the SMPS and BC measurement results on a log probability plot resulted in linearity. Thus, the geometric mean (GM) concentration was calculated after the logarithmic conversion. Pearson’s correlation analysis was performed to determine the correlation between BC and PN concentration and a correlation analysis between meteorological factors and BC concentration was also conducted.
3. Results
3.1. Elemental Carbon and Black Carbon Concentrations
The concentration of EC measured inside 8 toll booths was 3.3 μg/m3 on average, with a minimum concentration of 1.85 μg/m3 and maximum concentration of 5.84 μg/m3.
One toll booth was open all day at Gate A, and the concentrations of EC and BC, indicators of DEP, were the highest in Booth 1, open exclusively for trucks. The average EC concentration inside Booth 1 was 3.4 μg/m3, which was 70% of the concentration outside of the Booth (4.8 μg/m3). The GM concentration of BC inside Booth 1 was 7.3 μg/m3 (2.2 – 60.7 μg/m3), which was 88% of the GM concentration of BC outside Booth 1 (Table 1 shows the value as 8.3 μg/m3 [0.2 – 242.3 μg/m3]). The average EC concentration at Booth 2, open all day along with Booth 1, was 2.8 μg/m3, which was the second highest concentration measured after Booth 1. It was expected that the EC concentration at Booth 2 was lower because it was a gate for general vehicles, unlike Booth 1 which was open exclusively to trucks. It is likely that the emissions generated by general vehicles is relatively less compared to that by trucks.
On the other hand, the average EC concentration at Booth 3, open for only one hour in the morning for 347 general vehicles, was 1.6 μg/m3 which was 47% lower than Booth 1. At Booth 4, without any vehicle traffic, the average EC concentration inside the booth was 2.1 μg/m3 which was 62% lower than Booth 1. The EC concentration at the entrance gate without a booth was relatively low at 2.2 μg/m3, despite being open only for trucks and having a large amount of traffic.
Gate B showed a similar trend to Gate A. The average EC concentration was the highest at Booth 2, which was open exclusively to trucks and experienced heavy traffic, at 5.8 μg/m3 and lowest at Booth 4, which had no traffic, at 1.9 μg/m3. The average and maximum BC concentration in Booth 1 was 8.8 ug/m3 and 271 ug/m3, respectively.
3.2. Particle Number Concentration and Particle Size Distribution
Table 2 and Fig. 2 show PN concentration levels for different sampling sites. The average PN concentration inside the booth was 3.54 × 104 cm−3, and the maximum value was 2.48 × 105 cm−3, which was approximately 20 times higher than the minimum value of 1.21 × 104 cm−3. The average PN concentration outside the booth was 6.05 × 104 cm−3, which was 1.7 times higher than the average concentration inside the booth. The minimum value was 1.10 × 104 cm−3, which was 10% lower than inside the booth, and the maximum value was 2.95 × 106 cm−3, which was approximately 12 times higher than inside. The average PN concentration measured outside the office building was 3.13 × 104 cm−3, which was 88% of the average PN concentration inside the booth, but the highest concentration was 1.39 × 105 cm−3, which was much lower, at 56% of the maximum concentration inside the booth. The indoor background PN concentration of our Research Institute (in another region) was approximately 25% of the concentration in the air outside the tollgate office building.
The particle size distribution outside the booth varied in shape depending on the PN concentration. As for the size distribution outside the booth, when the PN concentration was equivalent to the background atmospheric concentration in the absence of vehicles, the mode diameter was 115 nm. Bimodal size distribution was observed when the PN concentration is at its median value. The particle size distribution changed to a bimodal distribution with peaks at particle diameters of 27 and 120 nm. When the PN concentration was at its highest, unimodal size distribution with a mode diameter of 20 nm was observed (Fig. 3).
The particle size distribution inside the booth exhibited similar behavior to the size distribution outside the booth. As for the size distribution inside the booth, when the PN concentration was at a maximum a unimodal size distribution was formed, with a mode diameter of 27 nm. In contrast, when the PN concentration was at the median value, the size distribution changed to a bimodal distribution with a mode diameter of 37 and 87 nm. When the PN concentration was at a minimum, the size distribution became unimodal again with the peak at 87 nm. Interestingly, smaller particles were observed inside the booth, meaning that diesel particles transported inside the booth remained stagnant while the diesel particles outside the booth were scavenged along with air stream due to blowing wind.
The particle size distribution measured outside the office building at low PN concentration was similar to the particle size distribution at atmospheric background level with a mode diameter of 115 nm. When the PN concentration was at its median value, a bimodal distribution with mode diameters of 15 nm and 36 nm was observed. At the highest PN concentration, a multimodal distribution was observed because both diesel particles and atmospheric background particles were present, with peaks at 20 and 115 nm. These results imply that the vehicle emission generated on the highway affected the air quality of office building 30 m away from vehicle traffic.
3.3. Real-Time Particle Number Concentration and Black Carbon
Fig. 4 shows the variation in PN concentration and black carbon levels inside the Booth 1 at Gate A over the sampling period. The PN concentration traced BC concentration well. There was statistically significant correlation between the PN concentration and BC concentration confirmed through Pearson’s correlation analysis. The correlation coefficient was 0.54 (p < 0.01).
3.4. Crystalline Silica
As a result of scanning the collected respirable dust by XRD, trace amounts of quartz and silicate minerals were found in some samples, but they were below the detection limit so could not be used for quantitative analysis.
3.5. Meteorological Influences
During the measurement time, the mean value (standard deviation, SD) of the temperature, humidity, wind speed, and wind direction were −0.76°C (4.46), 55.32% (19.17), 0.416 m/s (0.425) and 136.8° (131.4), respectively.
BC concentration did not show a significant correlation with meteorological factors. The correlation coefficients between BC measured per minute and temperature, humidity, wind speed, and wind direction were 0.047, −0.043, 0.006, and −0.036 respectively (p > 0.05).
4. Discussion
In June 2012, the IARC classified DEE as Group 1 carcinogenic agents, resulting in occupational lung cancer being linked to exposure to DEE in Korea. The number of cases where workers have applied to the government for industrial accident compensation have increased since 2012.
A study by Kim et al. [25] found that 9.7% of the general Korean population are exposed to carcinogens. The population attributable fraction (PAF) of total cancer was 1.1% and lung cancer was 7.0%. In another study by Leem et al. [26], the combined proportion of definite and probable occupational lung cancer among all lung cancers investigated over the four-year (2006–2009) study period was 11.7%. The main carcinogens include asbestos, RCS, radon, PAHs, DEE, chromium, and nickel. In addition, according to Choi et al. [27], the total number of Korean workers exposed to DEE was estimated to be 417,034 in 2013 (2.2% of the total labor force).
However, data on the concentration of DEE and characteristics of vehicle related particles in various workplaces is lacking.
Previous studies conducted at highway tolls showed that the mass median diameter (MMD) of suspended PM in the air was 0.78 μm on average, measured using a high-volume cascade impactor (Anderson Sierra-352, US) at highway tolls in 1999 [28]. However, the MMD measured in downtown Seoul (Bulgwang-dong) was 3.2–3.5 μm, and 4.6 μm during the yellow dust event. This means that highway toll stations are mainly polluted by vehicle emissions with small particle sizes. Additionally, evaluation of the concentration of PAHs in the same study showed that dibenzo[a,h]anthracene, which is known to have high toxicity among PAHs substances, had an average concentration of l.64 μg/m3, and the representative carcinogen, benzo[a]pyrene, had an average concentration of 5.67 μg/m3. This concentration was higher than in some urban areas such as Seoul [29].
Nam and Phi [18] evaluated PM10 concentration at toll gate booths. In 2020, Nam et al. [30] investigated the chemical composition of the total dust using XRF, and measured concentrations of PM10, PM2.5, and PM1 using direct-reading devices. There is a possibility of exposure to the carcinogen crystalline silica because its presence was analyzed by XRF in a total dust sample scattered around the toll station with high SiO2 content. Additionally, a study in India showed that Si was 60% in the urban PM2.5 elemental component [31].
This study uses XRD to analyze the crystal structure of crystalline silica unlike XRF, which analyzes chemical components. However, in our work, RCS was not detected in all samples (Below detection limit).
In contrast, according to Kim et al. [32] who investigated the concentration of RCS in the air in some regions of Korea, the average concentration per year was 0.8 μg/m3, and the average concentration in spring was 1.3 μg/m3. When there was a yellow dust event, the average concentration increased to 8.8 μg/m3. In addition, the concentration of RCS at the time of Yellow Dust Warning with a fine dust concentration of 800 μg/m3 is estimated to be 33 μg/m3, exceeding 25 μg/m3 of the American Conference of Governmental Industrial Hygienists Threshold Limit Value criterion. However, in this study, as part of the measurement of the working environment, it was measured according to the NIOSH 7500 method, so the amount of collected air volume was approximately 1700 L and the lower detection limit (LOD) was 3.4 μg/m3. Therefore, the concentration of crystalline silica in the air on the day of the toll booths survey could not have exceeded the LOD because levels were less than 3.4 μg/m3. However, a respiratory mask is necessary as the workers can be exposed to crystalline silica above the LOD on dusty days or when the wind is strong and the distribution of sediment dust is severe.
In a study by Kim et al. [33], a black carbon assessment was conducted at a toll gate at a university in Seoul. The concentration of BC was found to be 2.44 μg/m3, which was 1.4 times higher than that of the control group at 1.75 μg/m3. It is thought that black carbon concentration is higher in proximity to moving traffic, especially when a truck or an old vehicle using diesel fuel passes by, when the toll gate worker is exposed to an instantaneous high concentration of BC.
A relatively small number of studies have used SMPS to study PN concentration and particle size distribution at highway tolls in Korea. International studies have used SMPS to evaluate factors such as PN concentration and particle size distribution. A study [34] evaluated how these factors varied depending on the distance from major highways in Los Angeles, USA. As the distance from the highway increased, it was observed that the concentration of the smallest particle sizes decreased, and disappeared at 90 m. In addition, the concentration of PN in general decreased with distance. The average concentration of PN at the nearest measurement point, 30 m from the highway, was 1.5×105 cm−3, and the average concentration of BC was 5.4 μg/m3 which was similar to the concentration outside the booth in this study. In our study, the average PN concentration was 1.0×105 cm−3 and the average of BC concentration was 8.3 μg/m3. The PN concentration is considered to be higher in Los Angeles due to the substantially high traffic volume. The concentration of BC was higher in this study because the measurement location was a lane dedicated to freight cars.
In this study, the average BC concentration inside two toll booths was 7 to 9 times higher than the BC general indoor background concentration of 1 μg/m3 [35].
The average BC concentration was 2.2 times higher than the elemental carbon concentration in this study. This is comparable to the results of Lee et al. [15] who found a BC concentration that was 1.99 times higher than the EC concentration. The presence of black carbon may be overestimated in our study because light absorption occurs in organic materials such as brown carbon in addition to elemental carbon components in the measuring device [36].
DEP and RCS flowing into the toll booths may be affected by meteorological conditions because the toll booths are located outdoors. Real-time temperature, humidity, wind speed, and wind direction on the day of the measurement were not related to the BC concentrations for real-time DEP confirmation. The concentration of PM outdoors may increase when the wind speed is slow, and on a day with a high wind speed, the exposure level may be lower due to diffusion [37, 38]. According to Shin et al. [15], although the type of diesel engines are different from this study, it was confirmed that workers exposed to DEP outdoors were exposed to high concentrations on days when wind speeds were lower than 1.3 m/s. In this study, the average wind speed was lower than 1.3 m/s, but the change of BC concentration according to the real-time wind speed change could not be confirmed. Although the toll booth is located outdoors, it has a roof and vehicles continuously move in one direction along the one-lane road next to the booth. This seems to be because DEP was more influenced by the wind direction and speed. On the other hand, on days with strong wind speed, the concentration of DEP can be diluted and lowered as in other studies [37, 38].
This study is the result of only one day of measurements at two toll gates in South Korea. Therefore, this study is limited in its ability to represent the level of diesel exhaust exposure to all toll gate workers globally. The short measurement period in this study also inhibited the detection of crystalline silica. A larger variety of toll gates require additional investigation and measurements should be conducted for longer durations to determine crystalline silica exposure levels.
5. Conclusion
This study is the first to assess the exposure of highway toll workers to DEP and RCS in Korea. The objectives of this study are to evaluate the exposure level of RCS and DEP on workers in toll booths through the level of surrogates, EC, BC, and PN concentration, and to measure the particle size distribution of UFP.
The RCS concentration was below the detection limit. However, caution is required because roadside sedimentary dust is re-dispersed on windy days and may be exposed to RCS during yellow dust periods.
Workers in the toll booths of highways are continuously exposed to vehicle emissions, including diesel exhaust, and intermittently exposed to extremely high concentrations of diesel emissions. Results of the particle size distribution in the toll booths indicated that the most common particles had diameters in the range of 20–70 nm; these particles have a very high respiratory deposition rate.
In addition, the difference in the concentration of DEP exposure to workers inside the booth and outside the booth was not significant. Therefore, improving the work environment of the toll booths for worker’s protection and wearing an N95 respirator is recommended to be incorporated into national policy. Further daily monitoring of vehicle traffic information, meteorological conditions and DEP concentrations that allow for estimates of annual cumulative EC exposure will foster protection of workers’ health. We also recommend establishing an Occupational Exposure Limit and reduction controls against diesel exhaust sources to prevent diseases caused by DEP.
Acknowledgment
We would also like to thank the reviewers and editors of Environmental Engineering Research for their detailed and useful comments.
Notes
Author Contributions
J.S. (Ph.D. Student) and B.K. (Ph.D.) conceptualized. J.S., B.K. and K.L.(M.S.) conducted measurement and all the experiments. J.-S.J. (Ph.D.) and Y.-C.S. (Professor) reviewed the methodology and validated the technical aspects. B.K., J.L. (Professor) and J.S. wrote original draft. B.K., J.L., and J.S. reviewed and revised the draft. All authors have read and agreed to the published version of the manuscript.