Two sets of SQGs developed for marine and estuarine ecosystems (
MacDonald et al., 1996;
Long and MacDonald, 1998) were applied in this study to assess the ecotoxicological risk assessment of metals in sediments: (1) the effect range low (ERL)/effect range median (ERM) and (2) the threshold effect level (TEL)/probable effect level (PEL) values. Low range values (i.e., ERLs or TELs) are concentrations below which adverse effects upon sediment dwelling fauna would infrequent be expected. In contrast, the ERMs and PELs represent chemical concentrations above which adverse effects are likely to occur (
Long and MacDonald, 1998). At all sampling sites, concentrations of Pb and Zn were lower than the relevant values of effects range low (ERL). A decrease in As, Cr and Cu contents were in the range between ERL and ERM, proposing that these metals would be anticipated rarely to cause adverse biological effects on biota (
Long et al., 1995) in the studied areas. Furthermore, Ni levels were higher than ERM at some sampling sites, implicating that negative eco-risk effects often occur at Stas B and F. All elements, except for Cr and V, showed high positive correlations with mud, depth and TOM. The results presented here, suggested the impact of sediment grain size, depth and TOM content on metal distributions (Pb, Zn, As, Cu, Co and Ni) in the studied area. Generally, sediment grain size becomes finer as depth increases. Sandy sediments are made in coastal areas due to their much weight, and the amount of large-sized sediments decreases off shore, while fine-sized sediments increase when coming on shore. The fine grains, representing the higher rate of surface to volume and ionic absorption power, are more capable in the absorption of contaminated organic and inorganic materials (
McCave, 1984;
Horowitz and Elrick, 1987). Generally, fine-grained sediments carrying lots of organic matter are more contaminated than coarse-grained sediments. By contrast, the presence of a week correlation between Cr and V with sediment grain size, depth and TOM might indicate the possible effect of other parameters on their concentrations. Sediment grain size and organic matters has been proposed as dominant factors elucidating a significant part of the variance in species composition of nematode assemblages (
Heip et al., 1985). Additionally, water depths seemed to be potentially important which influence the structure of nematode assemblage, probably by determining other factors like as TOM and the stability of physico-chemical factors. Both the amount of organic matter and oxygen content are also determined by the sediment grain size, as the highest oxygen penetration depth with the lowest organic matter is generally observed in coarser sediment (
Castro and Huber, 2003). Several investigators have reported the increased density and diversity of marine nematodes with increased sediment grain size (
Vanaverbeke et al., 2011;
Fonseca et al., 2014), though contradictory result has also been published elsewhere (
Maria et al., 2013), implying that their response is species-specific. In our survey, the abundance and diversity indices of nematodes were negatively affected by depth, mud and TOM while they showed a positive correlation with sand. This is probably due to the organic matter enrichment in muddy station which can decrease total density and species number of nematode by creating a hypoxia inside the sediment; thus change their community structure (
Armenteros et al., 2009,
2010). But a higher density and species number is commonly observed in coarser sediments, where a deep penetration of oxygen is occurred by active flows through the sediments (
de Beer et al., 2005).