Difference in Macrobenthic Community Structures at Thermal Effluent Discharge Areas of Two Nuclear Power Plants in Korea
Abstract
This study investigated and compared sediment properties and macrobenthic community structures within heated effluent plumes at the discharge areas of Kori Nuclear Power Plant (KNPP) and Sinkori Nuclear Power Plant (SNPP) in Korea, which have different thermal effluent discharge systems. There were significant differences in sand, clay and organic carbon contents between sediments at the two discharge areas. Species richness and abundance of macrobenthos were higher at the SNPP discharge area than at that of the KNPP, although the values at both areas were comparable to previous studies in coastal areas of eastern Korea.
초록
본 연구는 고리원자력발전소와 신고리원자력발전소로부터 배출되는 온배수의 영향을 받는 해역에서 퇴적물 특성과 대형 저서동물의 군집구조를 조사하고 비교하였다. 두 배출해역의 퇴적물 특성들 중에서 모래와 펄, 유기물의 함량이 통계적으로 유의한 차이를 보였다. 대형 저서동물의 출현종과 개체 수는 고리원자력발전소보다 신고리원자력발전소의 배출해역에서 더 많았으며, 그 수는 이전의 다른 연구결과들과 유사한 수준이었다. 두 원자력발전소의 온배수 배출해역은 퇴적물 특성과 대형 저서동물의 출현종 수에서 유의한 차이를 보였다.
Keywords:
Discharge system, Macrobenthos, Nuclear power plant, Sediment properties, Thermal effluent키워드:
배출방식, 대형 저서동물, 원자력발전소, 퇴적물 특성, 온배수1. Introduction
Benthic organisms are capable of actively responding to physicochemical variables in aquatic environments, and their spatial and temporal distributions are largely influenced by sediment properties such as sediment particle sizes and organic carbon contents (Gray[1974]; Pearson and Rosenberg[1978]; Weston[1990]). Their distribution and community structure are highly affected by various human activities including logging, mining, reclamation, dredging, dumping, water pollution, oil spills, fishing, fish farming and so forth (Gray[1997]; Pearson and Rosenberg[1978]; Thrush and Dayton[2002]; Weston[1990]). In addition, altered temperatures via thermal pollution by fossil- fuel and nuclear power plants have potentials to affect physiological processes of all aquatic life forms and thus have caused great concerns in residents and fishermen of affected areas (Krishnakumar et al.[1991]).
Kori Nuclear Power Plant (KNPP), constructed in 1978, was the first large-scale nuclear power plant in Korea (Fig. 1). Composed of four turbines, it uses 5.1 billion tons of coastal seawater per year as condenser coolant and directly discharges heated effluent onto sea surface (KHNP[2013]). This discharge system has been widely adopted by most thermal power plants in Korea because of its economic feasibility (GESAMP[1984]); however, it has caused considerable public concerns and debates on probable adverse effects on coastal ecosystems. Sinkori Nuclear Power Plant (SNPP), constructed in 2012, has four turbine generators and uses 2.8 billion tons of coastal seawater per year (Fig. 1). In an effort to minimize environmental impacts associated with heated effluent, it adopted a novel outfall design where heated effluent goes through pipe lines submerged into seabed and discharges through diffusers at 15-20 m depth into overlying water column (KHNP[2013]). This discharge system is expected to mitigate ecological disturbance by rapid dilution and heat dispersal of effluent through mixing with underwater currents. As a result, the extent of heated effluent plume at the SNPP discharge area is much narrower than that at the KNPP (Fig. 1).
In this study, we compared sediment properties and macrobenthic community structures at the thermal effluent discharge areas of KNPP and SNPP. This study will provide baseline data that will improve not only our understanding of the ecological effects of thermal effluent, but also our ability to take effective conservation measures in the impacted areas.
2. Materials and methods
2.1 Sediment sampling
Study areas for KNPP and SNPP were limited to areas with a delta temperature (Δt) of 3 °C, which indicates a 3 °C increase in seawater temperature caused directly by heated effluent plumes (Fig. 1). Sediment samples were collected from each of 21 stations in soft bottoms at the two discharge areas between June 11-13, 2013 (Fig. 1), after considering the tide which had considerable effects on thermal plumes. Roughly 30 g of surface sediment was collected with a van Veen grab sampler (surface area: 0.1 m2). To analyze macrobenthic community structure, tworeplicate samples were collected and pooled together (total area: 0.2 m2). Fine sediment particles were removed with a 1-mm metal sieve, and the remaining sample was fixed with 10% neutralized formalin. Fixed samples were kept in watertight plastic bags and transported to the laboratory on ice in a cooler.
2.2 Sediment properties
To analyze sediment texture and type, sediment samples were air-dried and pretreated with boiling sodium acetate (NaCH3CO2) to remove carbonates and with hydrogenperoxide (H2O2) to oxidize organic matters. Pretreated samples were sieved using 0.062-mm metal mesh (phi scale, ϕ = 4). Sandy particles retained on the sieve were separated according to the phi scale using a Roe-tap sieve shaker for 15 min (Ingram[1971]), and the resulting filtrate was analyzed using the pipette method of McBride[1971]. Sediment texture and type were classified according to Folk[1968] and Folk and Ward[1957], respectively.
Organic carbon content was determined according to Santisteban et al.[2004]. About 5 g of each sediment sample was placed in a preweighed porcelain crucible and dried at 105 °C for 24 h in a drying oven. After cooling to room temperature, sample weight was measured using an analytical balance (Adventurer™, Ohaus Co., Pine Brook, NJ, USA). Dried samples were then ignited in a muffle furnace at 550 °C for 30 min and weighed after cooling to room temperature in a desiccator. Percent organic carbon was calculated as the sediment weight lost after ignition compared with before.
2.3 Macrobenthic community structure
Macrobenthic community structure at each discharge area was examined by delicately separating each fixed sediment sample into aquatic organisms and sediment particles under a magnifying lens or stereomicroscope (SZ51, Olympus, Japan). The organisms were identified to the lowest possible taxonomic level. Macrobenthic species richness of each station was determined as the number of species residing on and in the sediment and expressed by the sampling volume (species 0.2 m-2).
Abundances were determined by the number of individuals of each species from each station. Total pooled biomass of all organisms was measured after removing as much water as possible by blotting with a paper towel. Abundance and biomass were converted to the number of individuals per square meter (ind. m-2) and wet weight gram per square meter (ww g m-2), respectively.
In addition, to evaluate the health status of these macrobenthic assemblages, the Benthic Pollution Index (BPI) (Choi et al. [2003]) and Azti’s Marine Biotic Index (AMBI) (Borja et al. [2003]) were calculated based on species abundance.
2.4 Statistical analysis
Macrobenthic communities at the KNPP and SNPP discharge areas were compared with a nonparametric Mann Whitney’s U-test using SPSS software for Windows ver. 12.0 (SPSS Inc., Chicago, IL, USA). For multivariate analysis, square root transformation of the abundance data of macrobenthic fauna was applied. Cluster analysis and non-metric multidimensional scaling (nMDS) ordination based on the Bray-Curtis similarity matrix were performed using the PRIMER software ver. 5.0 (Clark and Warwick [1994]).
3. Results
3.1 Sediment properties
Sediment textures at the KNPP and SNPP discharge areas are shown in Table 1. At the discharge area of KNPP, gravel content ranged from 0.0-41.1% [6.5 ± 9.5% (mean ± SD)], while sand was 11.8-93.6% (47.6 ± 29.9%), silt was 0.1-37.3% (19.4 ± 15.3%) and clay was 0.1-48.2% (26.5 ± 19.5%). Their mean particle size (ϕ) at KNPP was 4.3 ± 3.2. At the SNPP discharge area, gravel content ranged from 0.0-15.1% (4.0 ± 3.8%), sand was 33.4- 94.9% (70.0 ± 17.0%), silt was 0.1-27.1% (11.4 ± 8.1%) and clay was 0.1-39.5% (14.7 ± 10.9%). Their mean particle size (ϕ) at SNPP was 3.1 ± 1.7. Overall, sediment composition was similar between the two discharge areas, with sand comprising the majority of each area, followed by clay, silt and gravel. However, the SNPP discharge area had significantly higher sand content (70.0% vs. 47.6%) and lower clay content (14.7% vs. 26.5%) than the KNPP (both P < 0.05).
Spatial distribution patterns of the sediment types (Folk[1968]) at the KNPP and SNPP discharge areas are shown in Fig. 2. Sediment samples were divided into four types: sand, muddy sand, sandy mud and clayey sand. Samples from KNPP (n = 21) were predominantly sandy mud (n = 12), followed by sand (n = 7), muddy sand (n = 1) and clayey sand (n = 1). In contrast, sediment samples from SNPP (n = 21) were mostly muddy sand (n = 15), followed by sand (n = 5) and sandy mud (n = 1). At both KNPP and SNPP, sandy sediment was typically found at stations proximal to the outfalls under the direct physical influence of discharged effluent, and both sandy mud and muddy sand were found distal to the outfall.
Organic carbon content was significantly higher at the KNPP discharge area (range: 2.6-8.7%; mean ± SD: 4.9 ± 1.5%) than that at the SNPP (2.2-6.2%; 3.9 ± 1.0%) (Table 1; P < 0.05).
3.2 Macrobenthic community structure
Species richness, abundance and biomass of macrobenthos at the thermal effluent discharge areas of KNPP and SNPP are compared in Table 2. Species richness at the KNPP discharge area ranged from 5 to 52 species (St. 3 and 21, respectively), with a mean of 27 ± 12 species. A total of 153 species were found in the KNPP discharge area, of which abundances ranged from 50 to 1,517 (St. 3 and 2, respectively) with a mean of 743 ± 380 ind. m-2. Total biomass ranged from 0.30 to 1,833.50 ww g m-2 (St. 3 and 10, respectively) with mean of 156.14 ± 388.38 ww g m-2. Meanwhile, at the SNPP discharge area, species richness ranged from 15 to 52 (St. 6 and 4, respectively) with a mean of 34 ± 10 species. A total of 195 species were found. Abundances ranged from 338 to 2,088 m-2 (St. 7 and 15, and St. 6, respectively) with a mean of 884 ± 394 ind. m-2, and biomass ranged from 16.75 to 636.80 ww g m-2 (St. 16 and 18, and St. 20, respectively) with a mean of 127.39 ± 154.15 ww g m-2. Species richness was significantly higher at the SNPP discharge area compared with KNPP (P < 0.05), but there were no statistical differences in abundance and biomass.
Spatial distributions of macrobenthic species richness at the KNPP and SNPP discharge areas (Fig. 3) were similar in pattern to the distributions of sediment types (Fig. 2). Species richness was low at the stations close to the discharge outfalls and high at the ones on the outer edges of both discharge areas. Fig. 4 shows macrobenthic community compositions in terms of species richness, abundance and biomass of the KNPP and SNPP discharge areas. Among the macrobenthos, Polychaeta was the most dominant taxon at both discharge areas, with the largest number of species (50.4% and 46.2% of species richness at KNPP and SNPP, respectively), followed by Arthropoda (17.6% and 23.1%), Mollusca (12.4% and 12.3%) and Echinodermata (7.2% and 6.7%). It was also the most abundant taxon at both KNPP and SNPP (62.8% and 65.9% of samples, respectively). At KNPP, the second most abundant taxon was Arthropoda (12.1%), followed by Echinodermata (9.8%) and Mollusca (5.8%), whereas Echinodermata was the second most abundant at SNPP (10.9%), followed by Arthropoda (7.3%) and Mollusca (6.4%). Composition ratios by biomass were distorted by the occurrences of large-bodied invertebrates, which included the cnidarian Cerianthus filiformis at the KNPP discharge area and the molluscan Siphonalia fusoides at SNPP.
BPI and AMBI values were calculated to evaluate the health status of the benthic habitat at the two thermal effluent discharge areas of KNPP and SNPP based on the macrobenthic community structures (Table 3). BPI values ranged from 47 to 92 with a mean of 66 ± 13 at KNPP and from 18 to 76 with a mean of 61 ± 11 at SNPP. According to Choi et al.[2003]), both discharge areas resembled offshore environments without organic pollution. Similarly, AMBI values at KNPP (0.3-3.0; mean ± SD: 1.9 ± 0.7) and SNPP (1.2-2.7; 1.9 ± 0.4) indicated that the ecological status of both discharge areas was ‘good’ (Borja et al.[2003]; five classes: high, good, moderate, poor and bad).
3.3 Cluster analysis and nMDS ordination
From the results of cluster analysis and nMDS ordination, the KNPP and SNPP discharge areas were divided into four groups with similar species compositions at 18% similarity level. Group A was composed of only one station (St. 6) of the SNPP discharge area and had no prominent dominant species. Group B was also composed of only one station (St. 3) of the KNPP discharge area and showed the lowest values of the species richness and abundance. Both stations were characterized by close locations to discharge outfalls, high ratio of sand and shallow water depth. Group C consisted of four stations (St. 4, 5, 13 and 17) of the KNPP discharge area. Diogenes edwardsii (Arthropoda) dominantly appeared as a dominant species in this group. The sediments of the three stations (St. 4, 5 and 17) were characterized by high ratio of sand and the other one (St. 13) did by muddy sand. Finally, Group D included all the other stations of both KNPP and SNPP discharge areas and were characterized by small sediment particle sizes and high species richness and abundance. Overall, our results showed that there were no prominent differences between the two discharge areas of KNPP and SNPP in terms of macrobenthic fauna besides several stations close to discharge outfalls.
4. Discussion
This study investigated and compared the sediment properties and macrobenthic community structures at the thermal effluent discharge areas of two nuclear power plants with differing discharge systems. Sediments of the discharge areas of KNPP and SNPP were both composed primarily of sand, followed by clay, silt and gravel. However, there was significantly less sand and more clay at the KNPP discharge area (47.6% and 26.5%, respectively) than those in the SNPP (70.0% and 14.7%, respectively). These differences in the ratio of sand and clay contents were also reflected in a larger mean particle size at the KNPP discharge area (ϕ = 4.3) than that at the SNPP (ϕ = 3.1). Furthermore, sediment texture differences were well aligned with those of sediment types; the KNPP discharge area was mostly composed of sandy mud, whereas SNPP was dominated by coarser muddy sand. Correlation analysis confirmed that most sediment properties (sand, clay and organic carbon content, and mean particle size) were significantly correlated with each other at both discharge areas (data not shown).
In agreement with previous studies (Kim et al.[2011]; Loi and Wilson[1979]), species richness and abundance of macrobenthos were higher at the SNPP discharge area where the sediments have higher sand and lower organic carbon contents than those at the KNPP, although only species richness was significantly higher. Sediment properties play a pivotal role in driving spatial distributions of macrobenthos (Gray[1974]; Pearson and Rosenberg[1978]; Weston[1990]). Though our findings cannot be directly compared with previous studies because of differences in seasonality, sampling volume and expertise in taxonomy, the species richness of the KNPP and SNPP discharge areas (153 and 195 species, respectively) were quite similar to those found previously at the KNPP discharge area (157 species, Kim et al.[2011]), as well as in unaffected eastern coastal areas (136 species, Choi et al.[2000]; 117 species, Shin et al.[2001]). They were also quite similar to those of Wolseong Nuclear Power Plant on the eastern coast of Korea (163 species, Seo et al.[2009]), but substantially lower than the species richness in uncontaminated coastal areas in Korea (319 species, Paik et al.[2007]; 284 species, Yoon et al.[2009]). These comparisons suggest that species richness was not significantly decreased (e.g., Krishnakumar et al.[1991]) by thermal effluent at the KNPP and SNPP discharge areas.
Abundance of macrobenthos at the KNPP and SNPP discharge areas (743 and 884 ind. m-2, respectively) were comparable but slightly higher than abundances previously found at the KNPP discharge area (552 ind. m-2, Kim et al.[2011]). They were also comparable to abundances in other eastern coastal areas (1,168 ind. m-2, Choi et al.[2000]; 535 ind. m-2, Shin et al.[2001]) and at the Wolseong Nuclear Power Plant discharge area (1,005 ind. m-2, Seo et al.[2009]), but less than half those of uncontaminated coastal areas in Korea (1,972 ind. m-2, Paik et al.[2007]; 2,202 ind. m-2, Yoon et al.[2009]).
Compositions by taxonomic group of macrobenthos were similar between the KNPP and SNPP discharge areas both in terms of species richness and abundance. Nearly half the species found at both discharge areas belonged to Polychaeta. It was also the most abundant taxon, comprising slightly more than 60% of the total abundances at both discharge areas. Previous studies of other coastal areas of Korea similarly showed that Polychaeta dominated in terms of species richness and abundance (e.g., Kim et al.[2011]; Seo et al.[2009]; Yoon et al.[2009]). Although highly contaminated areas typically exhibit abnormally high biomass of opportunistic species (e.g., Bamber and Spencer[1984]), the KNPP and SNPP discharge areas showed no obvious detrimental effects of heated effluents in terms of macrobenthic community structures.
Because our study did not include replicated control sites that experienced only natural temperatures, it is not possible to determine the direct effects of heated effluent on macrobenthic Fig. 5. Result of cluster analysis and nMDS ordination plot based on square-root transformation of abundance data of macrobenthos at the thermal effluent discharge areas of Kori Nuclear Power Plant (KNPP) and Sinkori Nuclear Power Plant (SNPP). All stations were divided into four groups at 18% similarity level. communities (e.g., Kim et al.[2013]; Lardicci et al.[1999]; Riera et al.[2011]). However, agreements between our study and previous ones at nearby areas (Kim et al.[2011]; Seo et al.[2009]) to our study areas and foreign areas (Lardicci et al.[1999]; Loi and Wilson[1979]; Riera et al.[2011]) suggest no detrimental effects of heated effluent from either KNPP or SNPP on macrobenthic community structure. However, a few previous studies have reported detrimental effects of heated effluent (Bamber and Spencer[1984]; Kailasam and Sivakami[2004]). Kim et al.[2013] also reported that species richness and abundance of meiobenthic nematodes were significantly different between impacted and control sampling locations at the KNPP discharge area.
Spatial distributions of sediment types within the heated effluent plumes were similar between the KNPP and SNPP discharge areas. At both areas, stations close to the discharge outfalls were composed primarily of sand, while offshore stations were mostly sandy mud at the KNPP discharge areas and coarser muddy sand at the SNPP. In addition, spatial distribution patterns of macrobenthic species richness were aligned with the sediment types. At both discharge areas, species richness was low at the proximal stations to the discharge outfalls of which sediments were relatively coarse, and high at the outer stations of which sediments were characterized by fine sediment types. These distribution patterns are consistent with those of previous studies (Kailasam and Sivakami[2004]; Kim et al.[2011]; Riera et al.[2011]) and partially explained by strong flow intensity removing fine sediments near the discharge outfalls.
Acknowledgments
This study was conducted by Korea Environment Institute (KEI) grant support. We have to disclose that this study was based on the partial extract of ‘Studies on Marine Environment Impact Prediction and Minimizing Measures according to Thermal Effluent Discharge from Power Plants’ of KEI.
References
- Bamber, R.N., and Spencer, J.F., (1984), “The benthos of a coastal power station thermal discharge canal”, J. Mar. Biol. Ass. UK, 64(3), p603-623. [https://doi.org/10.1017/S0025315400030290]
- Borja, A., Franco, J., and Muxika, I., (2003), “Classification tools for marine ecological quality assessment: the usefulness of macrobenthic communities in an area affected by a submarine outfall”, ICES CM, 2, p1-10.
- Choi, J.W., Je, J.G., Lee, J.H., and Lim, H.S., (2000), “Distributional pattern of macrobenthic invertebrates on the shallow subtidal sandy bottoms near Kangrung, East Coast of Korea”, J. Korean Soc. Oceanogr, 5(4), p346-356, (in Korean).
- Choi, J.W., Hyun, S., and Chang, M., (2003), “The summer benthic environmental conditions assessed by the functional groups of macrobenthic fauna in Gwangyang Bay, southern coast of Korea”, Korean J. Environ. Biol, 21(2), p101-113, (in Korean).
- Clarke, K.R., and Warwick, R.M., (1994), Changes in Marine Communities: An Approach to Statistical Analysis and Interpretation, UK, Plymouth Marine Laboratory, p144.
- Folk, R.L., (1968), Petrology of Sedimentary Rock, Hemphill, Austin.
- Folk, R.L., and Ward, W.C., (1957), “Brazos river bar: A study in the significance of grain size parameters”, J. Sed. Petol, 27(1), p3-26. [https://doi.org/10.1306/74D70646-2B21-11D7-8648000102C1865D]
- GESAMP, (1984), Thermal Discharges in the Marine Environment, UNEP Regional Seas Reports and Studies No. 45, United Nations Environment Programme.
- Gray, J.S., (1974), “Animal-sediment relationships”, Oceanogr. Mar. Biol. Rev, 12, p223-261.
- Gray, J.S., (1997), “Marine biodiversity: patterns, threats and conservation needs”, Int. J. Biodivers. Conserv, 6(1), p153-175. [https://doi.org/10.1023/A:1018335901847]
- Ingram, R.L., (1971), Procedures in Sedimentary Petrology, Wiley, New York, p49-68.
- Kailasam, M., and Sivakami, S., (2004), “Effect of thermal effluent discharge on benthic fauna off Tuticorin bay, south east coast of India”, Indian J. Mar. Sci, 33(2), p194-201.
- KHNP, (2013), A White Paper on Nuclear Power Generation, Korea Hydro and Nuclear Power Co. Ltd, (in Korean).
- Kim, D.I., Seo, I.S., Moon, C.H., Choi, B.M., Jung, R.H., and Son, M.H., (2011), “Community structure of macrobenthic assemblages around Gijang Province, East Sea of Korea”, J. Korean Soc. Oceanogr, 16(2), p97-105, (in Korean). [https://doi.org/10.7850/jkso.2011.16.2.097]
- Kim, H.G., Rho, H.S., and Oh, C.W., (2013), “Seasonal and spatial variations in nematode assemblages affected by thermal influence of nuclear power plant in Korea (East Sea, Pacific Ocean)”, Mar. Biol. Res, 9(8), p725-738. [https://doi.org/10.1080/17451000.2013.766346]
- Krishnakumar, V., Sastry, J.S., and Swamy, G.N., (1991), “Implication of thermal discharges into the sea-a review”, Indian J. Environ. Prot, 11(7), p525-527.
- Lardicci, C., Rossi, F., and Maltagliati, F., (1999), “Detection of thermal pollution: variability of benthic communities at two different spatial scales in an area influenced by a coastal power station”, Mar. Pollut. Bull, 38(4), p296-303. [https://doi.org/10.1016/S0025-326X(98)00149-0]
- Loi, T.N., and Wilson, B.J., (1979), “Macroinfaunal structure and effects of thermal discharges in a mesohaline habitat of Chesapeake Bay, near a nuclear power plant”, Mar. Biol, 55(1), p3-16. [https://doi.org/10.1007/BF00391711]
- McBride, E.F., (1971.), Procedures in Sedimentary Petrology, Wiley, New York, p109-127.
- Paik, S.G., Kang, R.S., Jeon, J.O., Lee, J.H., and Yun, S.G., (2007), “Distribution patterns of sandy bottom macrobenthic community on the Hupo coastal area, in the East Sea of Korea”, Ocean Polar Res, 29(2), p123-134, (in Korean). [https://doi.org/10.4217/OPR.2007.29.2.123]
- Pearson, T.H., and Rosenberg, R., (1978), “Macrobenthic succession in relation to organic enrichment and pollution of the marine environment”, Oceanogr. Mar. Biol. Ann. Rev, 16, p229-311.
- Riera, R., Núñez, J., and Martín, D., (2011), “Effects of thermal pollution on the soft-bottoms surrounding a power station in the Canary Islands (NE Atlantic Ocean)”, Oceanology, 51(6), p1040-1046. [https://doi.org/10.1134/S0001437011060142]
- Santisteban, J.I., Mediavilla, R., López-Pamo, E., Dabrio, C.J., Zapata, M.B.R., García, M.J.G., Castaño, S., and Martínez-Alfaro, P.E., (2004), “Loss on ignition: a qualitative or quantitative method for organic matter and carbonate mineral content in sediments?”, J. Paleolimnol, 32(3), p287-299. [https://doi.org/10.1023/B:JOPL.0000042999.30131.5b]
- Seo, I.S., Moon, H.T., Choi, B.M., Kim, M.H., Kim, D.I., Yun, J.S., Byun, J.Y., Choi, H.C., and Son, M.H., (2009), “Community structure of macrobenthic assemblages around the Wolseong Nuclear Power Plant, East Sea of Korea”, Korean J. Environ. Biol, 27(4), p341-352, (in Korean).
- Shin, H.C., Yoon, S.M., and Koh, C.H., (2001), “Spatial distribution of benthic macrofaunal community in Ulsan Bay and Onsan Bay, eastern coast of Korea”, J. Korean Soc. Oceanogr, 6(3), p180-189, (in Korean).
- Thrush, S.F., and Dayton, P.K., (2002), “Disturbance to marine benthic habitats by trawling and dredging: implications for marine biodiversity”, Annu. Rev. Ecol. Syst, 33, p449-473. [https://doi.org/10.1146/annurev.ecolsys.33.010802.150515]
- Weston, D.P., (1990), “Quantitative examination of macrobenthic community changes along an organic enrichment gradient”, Mar. Ecol. Prog. Ser, 61(3), p233-244. [https://doi.org/10.3354/meps061233]
- Yoon, K.T., Seo, I.S., Kim, K.B., Choi, B.M., and Son, M.H., (2009), “Community structure of macrobenthic fauna in the Hallyeohaesang National Park from Korea Strait, Korea”, Korean J. Environ. Biol, 27(1), p124-133, (in Korean).