ORIGINAL RESEARCH article

Front. Mar. Sci., 21 July 2025

Sec. Marine Pollution

Volume 12 - 2025 | https://doi.org/10.3389/fmars.2025.1621429

Multi-matrix contamination by cyclic and linear siloxanes in a highly industrialized estuarine environment of Korea: source identification, seasonal variation, and bioaccumulation potential

  • 1. Ministry of Education Key Laboratory for Ecology of Tropical Islands, Key Laboratory of Tropical Animal and Plant Ecology of Hainan Province, College of Life Sciences, Hainan Normal University, Haikou, Hainan, China

  • 2. Department of Environmental Health Sciences, Soonchunhyang University, Asan, Republic of Korea

  • 3. Department of Marine Science and Convergence Engineering, College of Science and Convergence Technology, Hanyang University, Ansan, Republic of Korea

  • 4. Maritime Safety and Environmental Research Division, Korea Research Institute of Ships and Ocean Engineering, Daejeon, Republic of Korea

Abstract

Cyclic and linear siloxanes are widely used in personal care, industrial, and consumer products. To better understand their fate and ecological risks in estuarine environments, this study investigated spatial distribution, seasonal variations, and the bioaccumulation potential of siloxanes in surface water, sediment, and in benthic organisms in an industrialized bay of Korea. Cyclic siloxanes exhibited higher concentrations than linear siloxanes across multiple environmental matrices, indicating their predominance and persistence. Water and sediment had distinct spatial distributions of siloxanes, which were influenced by local sources and hydrodynamic conditions. Industrial activities and power plant effluents were identified as major contributors to siloxane contamination in the coastal environment. Overall, the concentration of siloxanes in surface water was highest in winter. In particular, cyclic siloxanes showed a greater sensitivity to seasonal variation than linear siloxanes, with concentrations fluctuating significantly across sampling periods. The spatial distribution of siloxanes in sediment was strongly associated with organic carbon. Benthic invertebrates exhibited a strong potential for bioaccumulation of D5 and L9, with the highest bioaccumulation factors and biota-sediment accumulation factors among the detected compounds. These findings highlight the need for continued monitoring and management of siloxane contamination in industrialized estuarine environments.

1 Introduction

Siloxanes, including both cyclic (Dn, where n represents the number of silicon atoms) and linear (Ln) forms, are used in a broad range of products, such as cosmetics and personal care products (CPCPs), cooking utensils, pharmaceuticals, and as intermediates in the production of silicone polymers (; ). The main chain of siloxane molecules is composed of repeated silicon oxygen (Si-O) bonds, where silicon atoms are not only connected to oxygen atoms to form a skeleton, but also typically bonded to organic groups such as methyl groups. This structure endows siloxanes with a range of excellent properties, including high thermal stability (derived from strong Si-O bonds), excellent lubricity, extremely low surface tension, good hydrophobicity, flexibility, physiological inertness, and a wide viscosity range (depending on molecular weight, low molecular weight is in liquid form, high molecular weight can be in rubber or resin form) (; ; ; ). Their essential properties, such as high thermal stability, lubricity, and low surface tension, make them commercially useful and thus in strong demand (; ). Global silicone production has been estimated to be ~8.5 million metric tons in 2021 (), with the market expected to grow from USD 16.7 billion in 2021 to USD 23.4 billion by 2026 (). The large consumption of siloxanes has led to widespread contamination across multiple environmental matrices, including air, water, soil, sediment, and biological samples, on a global scale (; ; ; ; ; ; ).

Previous studies have reported ecotoxicological effects at concentrations below 10 μg/L for several cyclic siloxanes, such as octamethylcyclotetrasiloxane (D4) and decamethylcyclopentasiloxane (D5), in several freshwater and marine species (Sousa et al., 1995). More recent studies have confirmed that siloxanes exert toxic effects on the nervous, immune, endocrine, and reproductive systems in both animals and humans (; ). However, earlier ecotoxicological assessments of cyclic siloxanes–including D4, D5, and D6–suggested no significant health risks to aquatic organisms (; ; ; ; ). Despite these findings, subsequent research has reported a significant potential for bioaccumulation and biomagnification of siloxanes in a variety of aquatic ecosystems (; ; ; ; ; ). Additionally, high concentrations of siloxanes, ~1,000 ng/g lipid weight, have been detected in a range of aquatic species, including benthic invertebrates and fish (; , ; ; ).

Regulatory actions on the use of cyclic siloxanes (D4, D5, and D6) in CPCPs have been implemented by the European Commission (), the United Kingdom (, , ), Canada (, , ), the United States (), and Japan (). However, no specific regulations on the use of siloxanes currently exist in Korea. Although the release of polydimethylsiloxane (PDMS) that contains primary siloxanes into the marine environment is legally prohibited in Korea, PDMS is still used as an antifoaming agent in fossil fuel power plants. In 2016, government investigations revealed that several fossil fuel power plants had discharged siloxanes into the marine environment, raising both social and environmental concerns ().

Our previous studies confirmed that industrial activities are major sources of siloxane contamination in coastal environments (; ; ). Due to their extensive use in CPCPs and a number of industrial applications, siloxanes are continuously released into aquatic systems. Their physicochemical properties, such as hydrophobicity, volatility, and environmental persistence, raise concerns about their potential for long-term accumulation and ecological impact in marine ecosystems (). Given these concerns, a comprehensive understanding of the distribution of siloxanes across different environmental matrices is essential for evaluating their environmental fate and ecological risks. In the present study, the occurrence and bioaccumulation potential of cyclic and linear siloxanes were investigated in multiple environmental matrices present in coastal waters influenced by large-scale industrial complexes and fossil fuel power plants in Korea. The seasonal variability in the concentration of siloxanes in surface water was also examined to gain insight into the environmental dynamics and potential influencing factors in the aquatic environment.

2 Materials and methods

2.1 Sampling

Ulsan Bay and its connected river and stream systems are among the most heavily industrialized coastal zones and host the largest commercial harbors in Korea (). A total of 100 surface water and 25 surface sediment samples were collected from 25 locations that were categorized into streams, rivers, and the bay. Additionally, 80 surface water samples were collected from the outfalls of three power plants between 2017 and 2018 (Figure 1). Surface water sampling was conducted four times–in June, September, and December 2017, and January 2018–while sediment samples were collected in June 2017. Water and sediment were collected using stainless steel baskets and Van Veen grab samplers, respectively. In June 2017, three benthic species, mussels (Mytilus coruscus), conches (Strombus gigas), and sea cucumber (Holothuroidea), were collected near the power plant outfalls. Benthic invertebrates were collected and analyzed due to their sedentary nature and close association with sediments, particularly under strong tidal current conditions. Given the high hydrophobicity of PDMS and its strong affinity for sediment particles, these organisms serve as suitable bioindicators for assessing PDMS bioaccumulation. Their presence near the power plant’s drainage outlet provides an effective means of evaluating PDMS exposure levels in areas directly influenced by point source contamination. The soft tissues of each species were pooled and homogenized for analysis. All samples were stored at –20°C until further processing. In addition, PDMS (n = 5), which was used as an anti-foaming agent in the fossil fuel power plants, was collected for siloxane analysis.

Figure 1

2.2 Sample pretreatment and instrumental analysis

The experimental standards and reagents used in this study were consistent with those described in previous studies (; ). Water samples (~2L) were extracted using liquid-liquid extraction following the addition of mass-labeled internal standards (13C12-D4, D5, and D6; 100ng each; Moravek Biochemicals and Radiochemicals, Brea, CA, USA). Sequential extractions were performed using hexane (ultra-trace residue analysis grade; J.T. Baker, Phillipsburg, NJ, USA), a mixture of hexane and dichloromethane (ultra-trace residue analysis grade; J.T. Baker), and a mixture of hexane and ethyl acetate (HPLC grade; Sigma-Aldrich, St. Louis, MO, USA). All extracts were combined and concentrated to less than 1mL under a nitrogen stream. Freeze-dried sediment (~5g) or biological samples (~0.5g) were placed in polypropylene (PP) tubes and spiked with the same internal standards (100ng each). The samples were extracted by shaking on an orbital shaker at 250rpm for 60min. For sediment, the same solvents and extraction method used for the water samples were applied. Biological samples were extracted three times with hexane. After each extraction, the mixtures were centrifuged at 3000rpm for 5min, and the supernatants were transferred to clean PP tubes. Sulfur in sediment extracts was removed using activated copper pre-washed with hydrochloric acid (Sigma-Aldrich). For biological extracts, 10% of each was used to determine lipid content. A Turbovap evaporator (Classic LV; Vimpelgatan, Uppsala, Sweden) was used to concentrate each extract to 1mL of hexane for instrumental analysis.

A total of four cyclic siloxanes (D4–D7) and 15 linear siloxanes (L3–L17) were analyzed using a gas chromatograph coupled with a tandem mass spectrometer (GC-MS/MS; Agilent 7890/7000C, Wilmington, TX, USA) equipped with a DB-5MS capillary column (30m length, 0.25mm inner diameter, 0.25µm film thickness; J&W Scientific, Palo Alto, CA, USA). Analyses were performed in electron impact ionization mode, and siloxanes were identified and quantified using the multiple reaction monitoring method. Instrumental conditions followed those reported by . Briefly, Oven temperature was programmed to increase from 40°C (2 min) to 220°C (20°C/min) and to 280°C (5°C/min). This latter temperature was held for 10 min followed by 5 min at 300°C. The MS was operated in electron impact ionization mode at 70 eV. D4, D5, and D6 were quantified using the internal standard method, while D7 and linear siloxanes (L3–L17) were quantified using an external standard method.

Water temperature and salinity were measured in real time using a CTD instrument (SBE 19 plus V2, Sea-Bird Electronics, Bellevue, WA, USA). Total organic carbon (TOC) content in sediment samples was analyzed using a CHN analyzer (FLASH 2000 Series, Thermo Scientific, Boston, MA, USA).

2.3 Quality control

To minimize background contamination, the use of CPCPs was strictly avoided during all experimental procedures, following protocols established in previous studies (; ). The overall linearity (R2) of the calibration curves for siloxane standards in GC-MS instrument was 0.9991, with a range of 0.9919–0.9998. Limits of quantification (LOQs), calculated based on a signal-to-noise ratio of 10, ranged from 0.07ng/g dry weight (dw) (L6) to 0.97ng/g dw (L17). For water samples, the recoveries of the internal standards 13C12-D4, 13C12-D5, and 13C12-D6 were 62%±28%, 71%±21%, and 76%±22% (mean ± standard deviation), respectively. In sediment samples, the recoveries were 72%±18%, 82%±15%, and 83%±12%, respectively. For biological samples, recoveries of 13C12-D4, 13C12-D5, and 13C12-D6 were 71%±3.4%, 74%±4.2%, and 79%±4.7%, respectively. Matrix-spiked sample test was performed on water, sediment, and biological samples to assess the matrix effect. The overall recoveries of matrix-spike samples were in the ranges of 69%–118% (mean: 88%) for cyclic siloxanes and 65%–107% (mean: 86%) for linear siloxanes. To ensure instrumental stability and accuracy throughout the analysis, test standards were analyzed after every 10 sample injections. Detailed quality control results are provided in Supplementary Table S1.

2.4 Statistical analysis

Concentrations of individual siloxanes below the LOQs were substituted with zero for the calculation of mean and total concentrations. Spearman correlation analysis was conducted to assess the relationships between the concentration of siloxanes and other parameters. One-way ANOVA was performed to identify significant differences in the concentrations of siloxanes across operational groups (e.g., seasonal variability and spatial distribution). All statistical analyses were conducted using SPSS® version 23.0 (Armonk, NY, USA), with a significance level set at p < 0.05.

3 Results and discussion

3.1 Occurrence and concentrations of siloxanes

The concentrations of cyclic and linear siloxanes in water, sediment, and biological samples collected from Ulsan Bay and its adjacent regions of Korea are summarized in Table 1. With the exception of L4, all siloxanes were detected in at least one matrix, indicating widespread contamination in the coastal environment. Specifically, D5–D7 and L7–L15 in water, D4–D7 and L5–L17 in sediment, and D4–D7, L6, and L8–L12 in biological samples were detected in more than 50% of the total samples. In contrast, D4, L3–L6, and L16–L17 in water, L3–L4 in sediment, and L3–L5, L7, and L13–L17 in biological samples were detected less frequently (< 30% of total samples). The frequent detection of specific siloxanes for multiple environmental matrices suggests their environmental persistence, consistent with previous studies (; ; ; ). Species-specific accumulation patterns of siloxanes were observed, as evidenced by variations in the concentration of siloxanes among target species. This suggests that selective bioaccumulation and species-dependent metabolic processes are crucial in the distribution of siloxanes in biological samples (; ).

Table 1

CompoundsSurface water (ng/L)Sediment (ng/g dry weight)Benthos (ng/g dry weight)
aDR (%)MinMaxMedianMeanDR (%)MinMaxMedianMeanDR (%)MinMaxMedianMean
Temp (°C)8.836.621.422.7TOCb (%)0.1835.781.311.39Lipid (%)1.9820.72.865.89
Cyclic siloxanes
D431< LOQc252-d3.6575< LOQ2411.7413.567< LOQ19.110.69.92
D570< LOQ1804.079.4196< LOQ10347.448.11002341684370234
D684< LOQ12722.4546.296< LOQ15744.453.61002.1627.410.22.16
D763< LOQ3450.715.4196< LOQ38.04.947.561001.022.601.811.02
Linear siloxanes
L30< LOQ–––18< LOQ2.94–0.260< LOQ–––
L40< LOQ–––0< LOQ–––33< LOQ27.9–9.30
L519< LOQ0.81–0.0450< LOQ5.970.030.370< LOQ–––
L624< LOQ1.62–0.09100< LOQ3.610.180.4967< LOQ1.330.760.69
L763< LOQ11.90.230.6789< LOQ7.780.401.030< LOQ–––
L856< LOQ38.81.065.3293< LOQ35.54.046.701000.2010.65.645.48
L980< LOQ46.23.975.3586< LOQ5.070.580.921004.7548.714.122.5
L1091< LOQ48.76.718.9996< LOQ15.53.744.7267< LOQ10.14.184.77
L1198< LOQ18.91.682.7696< LOQ13.22.604.011007.1119.619.115.3
L1299< LOQ14.31.572.241000.5610.42.103.101002.5717.113.110.9
L1397< LOQ10.01.031.5096< LOQ7.301.372.0933< LOQ5.78–1.93
L1487< LOQ3.540.320.5096< LOQ2.780.500.7333< LOQ3.03–1.01
L1570< LOQ0.880.060.1296< LOQ0.860.120.190< LOQ–––
L162< LOQ0.29––89< LOQ0.280.050.090< LOQ–––
L172< LOQ0.14––71< LOQ0.26–0.040< LOQ–––
∑Cyclice< LOQ127511.664.70.16532106123247173338478.0
∑Linearf0.02016120.827.60.7110717.024.714.612972.271.9
∑Siloxaneg1.85131538.592.30.876391261483751805399860

Summary of the concentration of siloxanes in surface water, sediment, and benthic invertebrates collected from Ulsan Bay and its adjacent regions of Korea.

aDR, detection rate of individual siloxanes for all samples; bTOC, Total organic carbon; c< LOQ, less than limit of quantification; d-, not available; e∑Cyclic, total concentrations of cyclic siloxanes; f∑Linear, total concentrations of linear siloxanes; g∑Siloxane, total concentrations of all siloxanes.

Siloxane concentrations varied significantly among sampling locations within three orders of magnitude. The concentration of all types of siloxanes (∑Siloxane) ranged from 1.85 to 1315ng/L (mean: 92.3ng/L) in surface water, 0.87 to 639ng/g dw (mean: 147ng/g dw) in sediment, and 375 to 1805ng/g dw (mean: 860ng/g dw) in biological samples. Significantly higher concentrations of siloxanes in biological samples compared to water and sediment indicated a strong potential for bioaccumulation (; ; ; ). Specifically, the concentration of total cyclic siloxanes (∑Cyclic) in the biological samples ranged from 247 to 1733ng/g dw (mean: 788ng/g dw), significantly higher (p<0.05) than those in sediment (0.16–532ng/g dw; mean: 123ng/g dw). Similarly, the total concentration of linear siloxanes (∑Linear) in biological samples ranged from 14.6 to 129ng/g dw (mean: 71.9ng/g dw), also significantly higher (p<0.05) than those in sediment (0.71–107ng/g dw; mean: 24.7ng/g dw). In the water samples, ∑Cyclic concentrations ranged from below the LOQ to 1274ng/L (mean: 64.7ng/L), while ∑Linear concentrations ranged from 0.020 to 161ng/L (mean: 27.6ng/L). Across all matrices, ∑Cyclic concentrations were consistently higher than ∑Linear, likely due to the greater use and discharge of cyclic siloxanes (; ; ).

Among individual siloxanes, the highest concentrations in both water and sediment were observed for D6, with mean values of 46.2ng/L and 53.6ng/g dw, respectively, followed by D5 (mean: 9.41ng/L in water and 48.1ng/g dw in sediment). These results were consistent with previous findings from industrialized coastal regions, where D5 and D6 were the predominant siloxanes detected in sediments (, ; ). In the biological samples, D5 exhibited the highest mean concentration among cyclic siloxanes (763ng/g dw), followed by D6 (13.3ng/g dw). These results suggested a strong tendency for benthic organisms to selectively bioaccumulate D5 over D6, as was described in previous studies (; ). For linear siloxanes, the highest concentration in water was observed for L10 (mean: 8.99ng/L), followed by L9 (5.35ng/L) and L8 (5.32ng/L). In sediment and biological samples, L8 had the highest concentrations (6.70 and 22.5ng/g dw, respectively), followed by L10 (4.72 and 15.3ng/g dw) and L11 (4.01 and 10.9ng/g dw). Notably, L7–L13 exhibited higher concentrations in the industrialized bays of Korea compared to other linear siloxanes, consistent with previous reports (; ; ).

Spearman correlation analysis was performed on siloxanes detected in water, sediment, and biota samples to further investigate the relationships among these environmental media (Supplementary Table S2). Significant correlations were observed between water and sediment samples collected in different seasons (r = 0.743–0.867, p < 0.01), suggesting the potential for siloxane exchange within the water-sediment system. Additionally, the concentration of siloxanes in biota showed significant correlations with those in water and sediment (r = 0.525–0.704, p < 0.01), indicating that organisms may bioaccumulate siloxanes from both sources. When examining individual siloxanes, only certain linear siloxanes (L10–L14) in surface water displayed significant correlations with one another (r = 0.614–0.847, p < 0.01) (Supplementary Table S3), while correlations among the remaining compounds were weak or absent. This suggests the presence of multiple sources for siloxane contamination in water. In contrast, almost all siloxanes in sediment samples exhibited significant correlations (r = 0.330–0.968, p < 0.05), implying a common source or similar environmental behavior in a sedimentary environment. Furthermore, individual siloxanes were significantly correlated with TOC content in sediments (r = 0.376–0.727, p < 0.05), indicating that TOC may play a crucial role in governing siloxane distribution in sedimentary environments.

3.2 Global comparison of siloxane concentration in water, sediment, and biological samples

The concentration of siloxanes in water, sediment, and biological samples observed in this study was compared with those reported in previous research to better understand the current global status of siloxane contamination. Although the specific types and numbers of siloxanes measured varied across the studies, the primary siloxanes detected were largely comparable. For samples of water and sediment, comparisons were made using data from natural water bodies in rivers, lakes, and coastal environments. For biological samples, benthic invertebrates and small, low-trophic-level fish were considered. Overall, the concentrations of siloxanes in water, sediment, and marine benthos in this study fell within the worldwide ranges reported in previous studies (Supplementary Table S4).

The mean concentration of ΣCyclic siloxanes in water samples (64.7 ng/L) was higher than that reported in some freshwater samples from China (18.3–41.3 ng/L; ; ; ; ) and the Geum River in Korea (49.6 ng/L; ). However, it was lower than concentrations reported in other regions, including China (78.6–291 ng/L; ; ; ), several rivers in Korea (88.7–630 ng/L; ; ), Turkey (96.9 ng/L; ), Japan (221 ng/L; ), Spain (203–1531 ng/L; ), and Vietnam (350 ng/L; ). The mean concentration of ΣLinear siloxanes in water (27.6 ng/L) was higher than those previously reported in China (<LOQ–9.93 ng/L; ; ), but was 4–10 times lower than levels observed in China (180 ng/L; ) and Korea (96.7–304 ng/L; ; ). The detection of siloxanes, especially linear types, in surface waters remains limited and warrants greater attention.

In sediments, the mean concentration of ΣCyclic siloxanes (123 ng/g dw) was higher than that reported in the Llobregat River, Spain (3.39 ng/g dw; ), some rivers in Korea (3.39–106 ng/g dw; ; ; ), China (28.5–58.1 ng/g dw; ; ; ), Canada (50.9 ng/g dw; ), and Japan (71.9 ng/g dw; ). However, these concentrations were 2–20 times lower than those reported in coastal sediments from Korea (219–245 ng/g dw; , ), freshwater sediments in China (170–1035 ng/g dw; ; ; ; ; ; ), Japan (905 ng/g dw; ), Spain (2070 ng/g dw; ), and Vietnam (2518 ng/g dw; ). The mean concentration of ΣLinear siloxanes in sediment (24.7 ng/g dw) was higher than those in coastal sediments from Korea (12.5–22.8 ng/g dw; ; ; ) and rivers in China (12.2 ng/g dw; ). However, it was 3–50 times lower than values reported in sediments from Korean waters (70.1–467 ng/g dw; , ; ), Chinese rivers (84.0–1558 ng/g dw; , ; ; ), and Tokyo Bay, Japan (117–1310 ng/g dw; ).

Studies on siloxanes in biological samples are limited and often report concentrations using different units. Nevertheless, the mean ΣCyclic concentration in marine benthos (788 ng/g dw) observed in our study was within the ranges reported in China (9.96 ng/g dw, 11.9–95 ng/g wet wt; ; ; ; ; ), Canada (18.2–236.7 ng/g wet wt; ), Japan (155 ng/g wet wt; ), Korea (139–4309 ng/g lw; ; ), and Norway (719–1042 ng/g lw; , ). The mean ΣLinear concentration in marine benthos (71.9 ng/g dw) was also within the range reported in Canada (1.4 ng/g wet wt; ), China (11.3–21.9 ng/g wet wt; ; ), and Korea (272 ng/g dw; ).

3.3 Spatial distribution of siloxanes

The spatial distribution of cyclic and linear siloxanes in surface water collected from Ulsan Bay and its adjacent regions of Korea is presented in Figure 2. The mean concentrations of ΣSiloxane in surface water across the different coastal zones were in the following order: Taehwa River (120 ng/L) > Gosa Stream (118 ng/L) > Ulsan Bay (106 ng/L) > plant outfall (73.4 ng/L), with no significant differences observed. The contamination status by ΣSiloxane varied slightly depending on the sampling period. In June, the ΣSiloxane concentrations were in the following order: Gosa Stream (mean: 94.2 ng/L) > Taehwa River (76.7 ng/L) > Ulsan Bay (45.5 ng/L) > plant outfall (36.1 ng/L). In September, the order shifted to: Taehwa River (261 ng/L) > plant outfall (220 ng/L) > Ulsan Bay (164 ng/L) > Gosa Stream (145 ng/L). In December, the order was: Gosa Stream (169 ng/L) > Ulsan Bay (143 ng/L) > Taehwa River (83.0 ng/L) > plant outfall (34.9 ng/L). In January, Gosa Stream (78.9 ng/L) > Taehwa River (60.0 ng/L) > Ulsan Bay (35.2 ng/L) > plant outfall (11.3 ng/L). These results suggested that siloxane contamination in surface water was generally similar across sampling sites when seasonal variability was not considered. Notably, the highest ΣSiloxane concentrations were detected in September at Stations B01 (1315 ng/L) and D04 (1176 ng/L), both located near the power plant outfall. This indicated that power plant discharges may play a crucial role in siloxane contamination in water environments, consistent with previous findings (; ; ).

Figure 2

The concentrations of ΣSiloxane in sediment were in the following order: Gosa Stream (mean: 303 ng/g dw) > plant outfall (134 ng/g dw) > Ulsan Bay (108 ng/g dw) > Taehwa River (103 ng/g dw) (Figure 3). A similar trend was observed for sedimentary TOC content: Gosa Stream (mean: 2.26%) > Ulsan Bay (1.38%) > plant outfall (1.15%) > Taehwa River (0.84%). Unlike the spatial distribution in surface water, the sedimentary concentrations of siloxanes were influenced by both TOC content and source inputs, due to the hydrophobic nature of siloxanes. Previous studies also reported significant correlations between the sedimentary concentrations of siloxanes and TOC content (; , ; ). In our study, the ΣCyclic concentrations in the sediment were 4–7 times higher than those of ΣLinear, indicating that cyclic siloxanes are preferentially enriched in the sedimentary environment, in line with previous findings (, ; ). The highest sedimentary concentrations of ΣSiloxane were observed at Stations G2 (639 ng/g dw) and G3 (393 ng/g dw), both located in Gosa Stream. Previous studies identified Gosa Stream as a major contamination route for the flame retardants, plasticizers, and other contaminants (; , ). Our findings suggested that industrial activities contribute more significantly to sedimentary siloxane contamination than power plant discharges.

Figure 3

3.4 Seasonal variability of siloxanes in surface water

The seasonal variability in the concentrations of ΣSiloxane in surface water collected from Ulsan Bay and its adjacent regions of Korea is presented in Figure 4. The ΣCyclic concentrations in surface water collected in September (169 ± 277 ng/L) were significantly or slightly higher than those measured in December (64.0 ± 147 ng/L; p > 0.05), June (10.3 ± 8.18 ng/L; p < 0.05), and January (8.98 ± 17.4 ng/L; p < 0.05). Similarly, the mean surface water temperature in September (32.0 ± 2.78°C) was significantly higher than in June (21.4 ± 2.71°C; p < 0.05), December (19.4 ± 4.31°C; p < 0.05), and January (18.9 ± 4.31°C; p < 0.05). The ΣLinear concentrations in surface water collected in June (38.3 ± 28.7 ng/L) and September (31.6 ± 16.2 ng/L) were significantly higher than those in December (20.0 ± 27.1 ng/L; p < 0.05) and January (19.8 ± 24.4 ng/L; p < 0.05). In September and December, the ΣCyclic concentrations were 6.1 and 2.4 times higher than the corresponding ΣLinear concentrations, respectively. In contrast, in June and January, the ΣCyclic concentrations were lower than those of ΣLinear. These findings suggested that cyclic siloxanes in surface water may be more sensitive to seasonal variability than linear siloxanes. Previous studies reported higher levels of cyclic siloxanes in lake water during colder seasons, which were attributed to their low volatility and increased partitioning from sediment to water at lower temperatures (; ). However, such patterns were not observed in our study. This discrepancy could be due to the fact that seasonal variability affects not only temperature, but also factors such as salinity, chlorophyll content, suspended solids, surface runoff, and the half-life of organic matter (; ). Further research is needed to better understand the seasonal behavior of siloxanes in aquatic environments.

Figure 4

3.5 Compositional profiles and source tracking of siloxanes

The relative contributions of individual siloxanes to the ΣSiloxane concentrations in water, sediment, and benthic invertebrates collected from Ulsan Bay and its adjacent regions are shown in Figure 5. Additionally, the composition profile of siloxanes in defoamer used in power plant effluent is presented (Supplementary Table S5). Overall, the compositional profiles of siloxanes in defoamer, water, and sediment differed notably. In the defoamer, the contributions of ΣCyclic (mean: 55%) and ΣLinear (45%) were similar. In contrast, the relative contributions in surface water varied considerably by sampling period: ΣCyclic accounted for 23% (June), 63% (September and December), and 28% (January), while ΣLinear accounted for 77%, 37%, 37%, and 72%, respectively. In sediment, ΣCyclic contributed 80%, significantly higher than ΣLinear (20%), which aligned with previous findings from industrialized bays in Korea (, ). For benthic invertebrates, ΣCyclic dominated in mussels (96%), conches (96%), and sea cucumbers (66%), suggesting selective bioaccumulation of cyclic siloxanes from sedimentary environments. Similar patterns were reported in mollusks from the Bohai Sea (China) and biota from the Geum River (Korea) (; ).

Figure 5

The predominant siloxanes in surface water also varied by sampling period. In June, L8 (24%) and L10 (23%) were dominant, followed by L9 (18%) and D5 (15%). In September, D6 (46%) was most abundant, followed by L10 (13%) and D5 (12%). December samples showed a higher contribution of D7 (26%), followed by D6 (20%), D5 (15%), L10 (13%), and L9 (12%). In January, L12 (18%) was dominant, with L11 (15%), L13 (12%), L10 (12%), and D4 (11%) also contributing. These seasonal differences suggested that the composition of siloxanes in surface water is significantly influenced by the sampling period. Moreover, compositional profiles of siloxanes differed between the different water bodies (Ulsan Bay, Gosa Stream, Taehwa River, and power plant outfall), likely reflecting varying contamination sources in coastal zones. Previous studies of siloxanes in wastewater treatment plants (WWTPs) in Korea reported distinct composition profiles between domestic and industrial WWTPs, which was linked to differing usage patterns in household and industrial products ().

In sediment samples, D5 and D6 collectively accounted for 72% of the ΣSiloxane concentrations, consistent with earlier findings (; ; ; ). This is likely reflected the widespread use of D5 and D6 in industrial applications, as well as their association with high TOC content in sediments (; ; ). In contrast to surface water and sediment, D5 was predominant in benthic organisms, comprising 93% of ΣSiloxane concentrations in mussels and conches, and 62% in sea cucumbers, indicating a strong potential for selective bioaccumulation of D5. A previous study reported that the weak binding between D5 and the cytochrome leads to higher bioaccumulation in benthic invertebrates ().

3.6 Potential for bioaccumulation of siloxanes

The bioaccumulation factor (BAF) and biota-sediment accumulation factor (BSAF) are commonly used indicators to assess the bioaccumulation potential of chemicals in aquatic organisms (; ; ). In this study, field-based BAF (L/kg) was calculated as the ratio of lipid-normalized concentrations of siloxanes in aquatic organisms (CB, ng/g lipid weight) to the concentration of siloxanes in water (CW, ng/L). According to the Canadian Environmental Protection Act () and the Stockholm Convention on Persistent Organic Pollutants (), BAF values exceeding 5,000 L/kg are indicative of significant bioaccumulation potential. BSAF was calculated as the ratio of the lipid-normalized concentration of siloxanes in organisms (CB, ng/g lw) to the TOC-normalized concentration in sediment (CS, ng/g TOC), with values > 1.7 indicating preferential partitioning into biological lipids (; ). Among the 19 target siloxanes, 11 compounds (D4–D7 and L8–L14) were available for the BAF and BSAF calculation based on the detection rate (Supplementary Table S6).

Mussels exhibited notable bioaccumulation potential for D4–D6 and L8–L14, with field-based BAFs ranging from 1.30 × 104 to 6.12 × 106. They also showed significant BSAFs for D5 (6.45) and L9 (6.13), indicating strong sediment-associated uptake. Conchs showed high BAFs for D4–D6, L9, and L11–L12 (1.27 × 104–1.35 × 106), and also accumulated L9 from sediment (BSAF = 2.06). Sea cucumbers had bioaccumulation potential for D5–D6 and L8–L12 (BAF = 7.28 × 10³–8.50 × 105), and particularly high BSAF for L9 (21.1). Overall, D5 and L9 showed the highest bioaccumulation potential among the measured siloxanes across all three benthic species. Previous studies reported BAFs for D5 ranging from 7.06 × 10³ to 1.33 × 104 (; ), and BSAF v7alues from 1.5 to 2.1 (; ), confirming its notable bioaccumulation potential. However, data on the bioaccumulation behavior of linear siloxanes remain limited, suggesting the need of further studies.

4 Conclusions

This study provides a comprehensive assessment of the occurrence, distribution, seasonal variability, and bioaccumulation of cyclic and linear siloxanes in an industrialized estuarine environment. The predominance of cyclic siloxanes across water, sediment, and biota indicates their persistence and widespread contamination. The study finds distinct spatial distributions of siloxanes, which are driven by local industrial sources and hydrodynamic conditions. Siloxane contamination in water and sediment appears to be jointly influenced by industrial activities and power plant discharges. Seasonal variation significantly affects siloxane levels, particularly for cyclic siloxanes in surface water. The strong bioaccumulation potential of D5 and L9 in benthic invertebrates suggests ecological health risks. These findings highlight the need for regulatory action and continued monitoring of siloxanes in estuarine ecosystems impacted by industrial and energy production activities.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics statement

Wild benthic organisms were collected for this study; however, no ethical approval was required or obtained for these research activities.

Author contributions

WC: Data curation, Formal Analysis, Methodology, Writing – original draft. SL: Data curation, Formal Analysis, Investigation, Software, Writing – review & editing. H-KL: Data curation, Formal Analysis, Investigation, Methodology, Writing – review & editing. ML: Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing. H-BM: Conceptualization, Funding acquisition, Investigation, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Korea Institute of Marine Science and Technology Promotion (KIMST) through projects entitled “Development of management technology for persistent organic pollutants (POPs) in the marine environment (RS-2024-00417889)”, “Development of source identification and apportionment methods for toxic substances in marine environments (KIMST-20220534)”, and “Development of technology for impact assessment and management of HNS discharged from marine industrial facilities (KIMST-20210660)” funded by the Ministry of Oceans and Fisheries (MOF), Korea.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer SM declared a shared affiliation with the authors H-KL and H-BM to the handling editor at the time of review.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2025.1621429/full#supplementary-material

References

Summary

Keywords

siloxane, cyclic, linear, seasonal variation, benthos, bioaccumulation

Citation

Chen W, Lee S, Lee H-K, Lee M and Moon H-B (2025) Multi-matrix contamination by cyclic and linear siloxanes in a highly industrialized estuarine environment of Korea: source identification, seasonal variation, and bioaccumulation potential. Front. Mar. Sci. 12:1621429. doi: 10.3389/fmars.2025.1621429

Received

01 May 2025

Accepted

30 June 2025

Published

21 July 2025

Volume

12 - 2025

Edited by

Zhiyuan Wu, Changsha University of Science and Technology, China

Reviewed by

Jiashen Tian, Arizona State University, United States

Min-Kyu Park, Jeju National University, Republic of Korea

Sori Mok, Hanyang University, Republic of Korea

Yohei Sato, National Institute of Advanced Industrial Science and Technology (AIST), Japan

Updates

Copyright

*Correspondence: Hyo-Bang Moon,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics