Abstract
The relationship between heavy metal contamination in sediment and macrofauna in the upper Gulf of Thailand is presented as an indication of the environmental quality of coastal waters. This study aims to monitor five heavy metal elements (As, Cd, Mn, Ni, and Pb) between sediment and the sandworm (Perinereis quatrefagesi (P. quatrefagesi)). The geoaccumulation index (Igeo), enrichment factor (EF), and contamination factor (CF), including the pollution load index (PLI), were used for statistical analyses by one-way analysis of variance (ANOVA), and differences in the data were compared using the least significant difference (LSD) test. The sediment heavy metal concentrations were found to decrease (Al > Mn > Pb > Ni > Cd > As), and the quantity of heavy metal contamination in the sediment was not over the emergency point defined by the Pollution Control Department in Thailand. The Mn at the SH and SP points has Igeo > 0 but <1 and is considered not polluted to moderately polluted. The EF overall is <2 and is deficient in mineral enrichment. The CF of the Samut Sakhon (SH) and Samut Prakan (SP) sites had high Al and Mn values, indicating a low pollution level, but the PLI had an all-site average of 0.0215 ± 0.0252, which is <1. This indicates that the present areas are close to ideal and not polluted. To measure heavy metals in macrofauna, P. quatrefagesi was assessed, and concentrations were found to decrease (Mn > Ni > Cd > Pb > As). The Mn between sediment and the sandworm was directly correlated (r = 0.976; ). The sandworm performs as a bioindicator for the quality of coastal sediment, specifically with Mn; thus, the result present at a baseline level may grow in the future if there is no control measure for defensive measures.
1. Introduction
Since the upper Thai Gulf connects to the central and eastern regions of Thailand, it has become a hub for industry and agriculture. However, the Thai Gulf is a semienclosed topical sea and serves as a boundary to Malaysia, Cambodia, and Vietnam [1]. The upper zone of the Thai Gulf is connected to the industrial zone and major city and has similar characteristics to ก (Kor Kai is the 1st Thai alphabet). Effectively, locals refer to this zone as “Auw Kor Kai” or “Kor Kai Bay,” and its area encompasses seven provinces: Samut Sakhon, Samut Prakan, Samut Songkhram, Chachoengsao, Chonburi, and Rayong. The major rivers of Thailand that empty into the upper Thai Gulf are the Chao Phraya, Mae Klong, Tha Chin, and Bang Pakong. However, the sandy beaches around the country’s tropical regions encourage rapid urban development and tourism. The upper zone of the Thai Gulf also has sandy beaches in the provinces of Chonburi, Rayong, and Phetchaburi and supported over 28 million visitors in 2016 [2, 3]. The sediment comprises the physical dynamics that accommodate and are used by biota, and it serves as a bioindicator of environmental quality [4]. The primary principle of sustainable development is the protection of the environment from human activity, especially in developing countries [5], and one component of protecting the environment is monitoring environmental quality [6]. However, elemental contamination in sediment can also be a secondary source of water quality issues, such as environmental change, and determining the level of heavy metal contamination in sediment is a crucial component of risk evaluation and monitoring in aquatic environments. In the study, we used the geoaccumulation index (Igeo) and enrichment factor (EF) to measure the level of anthropogenic contaminant deposition in sediment and the pollution load index (PLI) to determine the pollutant’s concentration in each individual sediment particle.
Macrofauna is an animal group ranging in size from 0.5 to 5 cm that inhabits soft sediment, sand, and mud, and they are associated with fragmentation and decomposition of organic matter as well as nutrient uptake from organic matter, which promotes growth and biofiltration [7, 8]. Furthermore, macrofauna is among the most suitable organisms for bioindicators and bioaccumulators in environmental monitoring studies as they use nutrients from the beach for maintenance and survival [9, 10]. These nutrients act as specific and appropriate biotic tools that can be used when monitoring environmental changes caused by toxic contamination in an area [11, 12], especially for heavy elements that chronically contaminate the food chain. Nevertheless, this study investigates the macrofauna in the sediment of the upper Thai Gulf to assess whether its ratios of heavy metal absorption reflect those of the environment, and the aim is to evaluate the relationship of heavy metal contamination between the beach sediment and macrofauna of the upper zone of the Thai Gulf. This research selected its environmental monitoring areas using literature reporting environmental contamination in the upper zone of the Gulf of Thailand and focused its survey on five heavy metals (As, Cd, Mn, Ni, and Pb) in the sediment and macrofauna. The macrofauna used in this research was the sandworm Perinereis quatrefagesi (P. quatrefagesi), which is distributed in the sediment around the upper Thai Gulf. P. quatrefagesi is responsible for the decomposition of aquatic animals on the coast, which makes it extremely appropriate for indicating the environmental quality [13–15] of the upper zone of the Thai Gulf. This article’s presence aims to focus on monitoring elemental contamination in the upper zone of the Thai Gulf area, where there are industrial and agricultural zones with high vulnerability to pollution distribution in the environment, as well as the upper zone of the Thai Gulf being very important to the country’s local fishery. The findings of this study suggest that environmental monitoring of sources of pollution will support the continuation of environmental protection policies in the future.
2. Methodology
2.1. The Sample Collection Areas
After 1958, Thailand developed from an agricultural country to an industrial one by opening the public sector and investing in the country. From 1986−1991, the country developed the industrial sector so quickly that Thailand was designated as a newly industrialized country (NIC) of Asia [16]. Presently, the industrial and service sectors represent the majority of the gross domestic product (GDP) of the country [17]. Concurrently, rural areas are being urbanized, which increasingly affects residential areas in the economic zones, and the populations of high-density cities are growing, which is characteristic of urbanization in the country [18]. The sediment may contain some elemental residue, affecting pollutant levels in the environment.
The sediments around the upper zone of the Thai Gulf have collected during the day from 10 : 00 AM to 3 : 00 PM during an ebbing tide from January to March 2022. The collection points and their distribution in the sample site are Phetchaburi (PB) coordinates Lat 13.224112 Long 99.993748; Samut Songkhram (SS) coordinate Lat 13.361663 Long 100.022530; Samut Sakhon (SH) coordinate Lat 13.586446 Long 100.582095; Samut Prakan (SP) coordinate Lat 13.487892 Long 100.814985; and Chachoengsao (CH) coordinate Lat 13.469204 Long 100.979118, as shown in Figure 1. A quadrate was used to create a margin area on the floor and collect sediment at a depth of 0–5 cm, or about 100–200 g, in the red graphic of the quadrate. However, the sample shall have a weight of about 1.3 kilograms per station and be a homogeneous mixture prior to collection in a polypropylene bag.

2.2. Sediment and Macrofauna Collection and Extraction
The sediment was collected from the same sample collection points and mixed homogeneously. Sediment samples weighed approximately 5 kg at each site. The microfauna sample collection was conducted within 0.04 m2 of the quadrat sampler at depths of 0–5 cm using mattock plastic. Twenty samples were collected at each site to ensure the collection of P. quatrefagesi, and samples were soaked and stored in absolute ethanol before analysis in laboratories. However, successful sample collection of microfauna only occurred at three sites in this study: SH, CH, and SP.
The sediment samples were dried in a hot air oven at 105°C for 168 h until dry. They were then ground using a mortar and pestle. Net no 20 sifted soil was selected for extraction and used in ICP-OES analysis. For this, 2 g of each soil sample was mixed with concentrated hydrofluoric acid (HF), concentrated perchloric acid (HClO4), and concentrated nitric acid (HNO3) at a 1 : 1 : 1 ratio, totaling 20 mL. Extractions occurred at 500°C in a SpeedDigester K-425 BUCHI (Switzerland) until dry. Each residue was rinsed with 1% HNO3 and sieved through filter paper. The supernatant was transferred to a 50 mL volumetric flask, and 1% HNO3 was added for continued inductively coupled plasma (ICP) analysis in a PlasmaQuant 9100 series (Germany).
For the sandworm heavy metal analysis, 1 g of wet weight was soaked in 10 mL of concentrated HNO3 in a beaker for 24 h. After initial digestion, concentrated HF, HClO4, and HNO3 were added in a 1 : 1 : 1 ratio, totaling 10 mL. Extractions were performed at 500°C in a SpeedDigester K-425 BUCHI (Switzerland) until dry. Each residue was rinsed with 1% HNO3 and sieved through filter paper. The supernatant was transferred to a 20 mL volumetric flask, and 1% HNO3 was added for ICP analysis in a PlasmaQuant 9100 series (Germany), and the sample was compared with the ICP multielement standard solution of the AccuStandard (USA).
2.3. Geoaccumulation Index (Igeo), Enrichment Factor (EF), and Contamination Factor (CF)
The geoaccumulation index (Igeo) was originally formulated by Muller [19] and is a quantitative measure of pollution in aquatic sediment [20]. It was developed through an understanding of the lithogenic effect. Igeo was derived using the following formula:where Cn is the measured concentration of an element in the sediment and Bn is the background value of the element. The interpreted values of Igeo are <0 = not polluted, 0-1 = not polluted to moderately polluted, 1-2 = moderately polluted, 2-3 = moderately to strongly polluted, 3-4 = strongly polluted, 4-5 = strongly to extremely polluted, and >5 = extremely polluted.
The enrichment factor (EF) was derived by the following formula:where (C/RE) sample is the concentration (C) of an element to a reference element (RE) in the samples, and (C/RE) background is the concentration (C) of an element to a reference element (RE) present in the background. Aluminum (Al) was used as the reference element, as it is a major component of clay, and the interpreted values of EF are <2 = deficiency to mineral enrichment, 2−5 = moderate enrichment, 5−20 = significant enrichment, 20−40 = very high enrichment, and >40=extremely high enrichment.
The concentration factor (CF) is approximated as the ratio of the observed concentration of an element in the sample (Ci) to the background level of the same element (Cb). The concentration factors were calculated as follows:
The interpreted values of CF are <1 = low pollution level, 1−3 = moderate pollution level, 3−6 = considerable pollution level, and >6 = very high pollution level.
The background element values for the Igeo, EF, and CF calculations were formulated from the references for As and Ni [21], Al and Pb [22], Cd [23], and Mn [24].
After calculating the concentration factors, they were used to calculate the pollution load index (PLI), which demonstrates the general contamination level [25]. The formula for PLI is as follows:where CF is the contamination factor of the elements and n is the number of observed elements. The interpreted values of PLI are 0 = perfection, <1 = baseline level, and >1 = polluted.
2.4. Statistical Analysis
Data were analyzed using a one-way analysis of variance (ANOVA), and the differences in the datasets were compared using the least significant difference (LSD) test with a value of <0.05. Correlations of data were assessed using Pearson’s r correlation coefficient (). All analyses were conducted using Statistical Package for the Social Science (SPSS) v.22.
3. Results and Discussion
3.1. The Quantity of Heavy Metals in the Sediment of the Thai Gulf
The heavy metal content in the sediment of the upper Thai Gulf was mainly composed of Al, and the average concentration of Al in the samples was 15,642 ± 2,020 mg/kg dry weight. The heavy metal concentrations of the samples exhibited the following pattern: Al > Mn > Pb > Ni > Cd > As. The sediment concentrations of heavy metals in the upper Thai Gulf are presented in Table 1. The concentration of Al was higher at SH (55,358 ± 802 mg/kg) and SP (12,337 ± 770 mg/kg), but not significantly higher at all the sample collection sites. The average sediment As concentration was 0.0612 ± 0.0225, and the As concentration was significantly higher at the PB site than at the SS collection site (). The Pb concentrations at the SH, CH, and SP collection sites were not significantly different but were significantly higher than those at the PB and SS sites (). The Cd concentration was high at the SP site and significantly different () from the SH, CH, and PB collection sites. The Mn and Ni concentrations were significantly higher at the SP site () than at the SH, CH, PB, and SS collection sites. The concentrations of heavy metals in the sediment of the upper Thai Gulf stand were compared to the emergency points of contamination. In using the scale of land use for agriculture and residential areas by the Pollution Control Department [26], reported metal concentrations are Cd < 67 mg/kg, Mn < 1.7 g/kg, Ni 140.4 mg/kg, Pb < 400 mg/kg, and As < 3.9 mg/kg. Thus, the averages of heavy metals in this series do not exceed the standard for any element. The amount of Al between study sites, however, is not statistically significant, but there are differences in the volume because the nature of the sediment component varies depending on the site. For example, the PB and SS sites have a majority of sand components, whereas the SH, CH, and SP sites almost exclusively have alluvial soil components, which has an impact on the number of elements in the area.
3.2. The Igeo, EF, and CF in Sediment of the Upper Gulf of Thailand
The values of Igeo of the upper Thai Gulf were almost not over 0 (unpolluted) and are presented in Table 2, which aligns with Thongra-ar et al. [27], who presented the geoaccumulation indices (Igeo) of Cd and Ni at values below 0. However, in this study, Mn in the SH and SP collection sites had high Igeo values over 0 but <1 (0.1873 and 0.5936, respectively). This indicates that SH and SP are not polluted or moderately polluted with Mn. The Al at SH was 0.6376 ± 0.0062, which means that the area was not polluted to moderately polluted with AL. The Igeo of Al at SH was significantly different () from those at the CH, SP, PB, and SS sites. The Igeo of As showed no significant differences among any study sites. Mn and Ni Igeo values showed significant differences among all study sites (). Pb Igeo values between the PB and SS sites were not significantly different, but the two sites were significantly different () from the SH, CH, and SP sites. Nevertheless, the Igeo of the sediment of the Thai Gulf indicates that the area is mostly not polluted, and some areas are experiencing light to moderate contamination [28].
The EF values in the upper Thai Gulf area are lower than 2, meaning that the sediment has a deficiency in mineral enrichment (Table 2). The EFs of As, Pb, and Cd between the PB and SS sites were not significantly different, but the EFs of As, Pb, and Cd at the PB and SS sites were significantly different () compared to those at the SH, CH, and SP sites. The EF of Mn at the SS site was significantly higher () than that at all other study sites. The EF of Ni at the SS site was significantly higher () than that at the SH, CH, and SP sites but was not significantly different from that at the PB site. The EF values of the upper Gulf of Thailand indicate low coastal heavy metal contamination. Two sample sites had CF values >1: Al (6.512 ± 0.0943, considerable pollution) and Mn (1.539 ± 0.0166, moderate pollution) at SH, and Mn (3.923 ± 0.0263, considerable pollution) and Al (1.451 ± 0.0906, moderate pollution) at SP.
The three analyses illustrate the sediment quality of the upper Gulf of Thailand. The SH and SP collection sites present Al and Mn as ranging from high to over-recommended levels to no pollution. Therefore, PLI is important to use when conducting holistic environmental monitoring of coastal areas in the Gulf of Thailand. The average PLI of all sites in this study was 0.0215 ± 0.0252, which is <1, and indicates that the areas are ideal and not polluted [29, 30]. The PLI values of the SH and SP collection sites were significantly different () from all other sites, as shown in Table 3.
3.3. P. quatrefagesi in Sediment
Macrofauna is used as bioindicators of coastal environments, as they indicate environmental quality. The sandworm P. quatrefagesi is a bioindicator of acute pollution conditions for heavy metals [31, 32]. The sandworms in this study were collected from three sites: SH, CH, and SP. The features and characteristics of the sandworms from the study sites are presented in Figure 2. As concentrations were significantly higher in worms from the SP collection site () than in those from the SH and CH sites, the Pb and Cd contamination in the sandworms from the SH and SP collection sites were significantly different () compared to the CH site. The Mn concentrations in the worms were significantly different (SP > SH > CH, ), and the Ni concentrations in sandworms from all study sites were not significantly different. The sandworm heavy metal concentrations are presented in Table 4. P. quatrefagesi is the optimal bioindicator for copper, manganese, lead, and cadmium [33, 34]. The sandworm is the primary way heavy metals are introduced into the food chain, as they efficiently absorb heavy metals from the environment [35, 36].

3.4. Correlation between Heavy Metals in Sediment and P. quatrefagesi
The correlation between sample location and heavy metals in sediment and sandworms was evaluated using Pearson’s r correlation coefficient. Positive correlations were detected in the sediment for manganese (r = 0.710; ), lead (r = 0.691; ), and cadmium (r = 0.671; ), and arsenic was positively correlated between sample location and sandworm (r = 0.671; ). Cadmium heavy metal contamination between sediment and the sandworm was significantly positively correlated (r = 0.969; ); similar correlations were also seen with lead and manganese in the sandworm (r = 0.993; ). Lead in the sandworm and manganese in the sediment were significantly positively correlated (r = 0.970; ), and lead in the sandworm and cadmium in the sandworm were also significantly positively correlated (r = 0.954; ). Lead in the sandworm and nickel in the sediment were significantly positively correlated (r = 0.916; ). The correlations of the other components are presented in Table 5.
3.5. Factor Analysis of Heavy Metal Components of Macrofauna and Sediment
The factor analysis of the parameters of the 11 components (heavy metals, sediment, and sandworms) was performed using principal component analysis (PCA). The components found that three PCs had an eigenvalue >1 and explained 88.89% of the cumulative variance in the dataset (Table 6). The component of variance was >10%, as shown in PC1 (62.13%) and PC2 (17.23%) (Table 6, Figure 3). The location was the most important contributor, for which the highest factor load was 0.9394. Location was followed by the primary loading factor, cadmium in sediment (0.8449), and in PC2, the three components were lead in sandworms, cadmium in sandworms, and nickel in sediment. In PC3, arsenic in sediment was found to be an important factor. The PC distributions are presented in Figure 3.

3.6. Quality of Sediment around the Upper Gulf of Thailand
The results of the PCA found that location affects the level of heavy metals in sediment and macrofauna. This is also demonstrated by the PLI results, which indicate that the locations of this study are close to ideal. This shows that the coast of the upper zone of the Gulf of Thailand has good environmental quality. However, PLI analysis for surveying environmental quality is only one method for monitoring the environment [37], and it is important in long-term coastal environmental assessments and risk evaluations to use other indicators in conjunction with PLI analysis [38, 39]. Furthermore, in this study, the EF of all components and all areas was <2, which indicates a deficiency in mineral enrichment. Additionally, Igeo and CF signify no Al or Mn pollution, and Al is the major element of this soil, which is readily absorbed by bottom sediments in the form of metastable compounds [40]. Mn is a product of anthropogenic activity and street dust and is an influential indicator of coastal environmental quality [41, 42]. Its existence is an indicator of human activity on land, which affects the quantity and type of heavy metals in sediment [43, 44].
3.7. Performance of Sandworms as Bioindicators for Heavy Metals
The PCA and PC2 results showed that Pb and Cd had high loadings in macrofauna (0.9852 and 0.9654) and were important factors for this analysis. The sandworm (Perinereis spp.) is the optimal macrofauna to be collected and used to detect Zn, Cu, Pb, and Cd in the sediment [45, 46]. However, in the heavy metal concentration correlation analyses, only Mn was shown to be directly correlated between sediment and macrofauna (r = 0.976; ). Mn concentration in Perinereis spp. is indicative of human activity, which may affect sediment and cause chronic or acute pollution conditions [32]. The ability of polystyrene microplastics to contaminate seawater and sediment is also linked to bioaccumulation within Perinereis spp. [47].
4. Conclusion
The sediment heavy metal concentrations were found to decrease (Al > Mn > Pb > Ni > Cd > As), and the quantity of heavy metal contamination in the sediment was not over the emergency point defined by the Pollution Control Department in Thailand. Mn at the SH and SP collection sites had an Igeo, which is considered not polluted to moderately polluted and represents a gray area where there is the possibility of pollutants increasing to a crisis point if defensive measures are not taken in the future. The EF overall was deficient in mineral enrichment. The CF of the SH and SP sites had high Al and Mn values, indicating a low pollution level, but the PLI had an all-site perfection level. However, Al is a major component of parent soil sediment in the central region, although Mn levels are rising due to intense traffic and transit that includes industrial zones. This indicates that the present areas are close to ideal and not polluted, but the present result at a baseline level may increase in the future if no defensive measures or controls are taken. To measure heavy metals in macrofauna, P. quatrefagesi was assessed, and concentrations were found to decrease (Mn > Ni > Cd > Pb > As). The Mn between the sediment and the sandworm was directly correlated. The sandworm contributes as a bioindicator for the quality of coastal sediment, particularly when Mn is present. However, sediment and macrofauna are important instruments for measuring and monitoring the environmental conditions associated with pollution in the estuarine and Gulf of Thailand ecosystems around the upper zone of the Gulf of Thailand.
Data Availability
All the data generated and analyzed during the current study are included in this published article and are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Pakjirat Singhaboot and Kun Silprasit designed the research methodology. Arin Ngamniyom collected samples. Patarapong Kroeksakul planned the project and was the main writer of the manuscript.
Acknowledgments
This research was supported by the Srinakharinwirot University (Code. 562/2564).