We quantified the distribution, abundance and assemblage structure of macrofauna at 22 stations in Gorgan bay, seasonally in 2012-2013. Also, depth, temperature, salinity, DO, TOM and sediment particle size were measured in each station. The highest concentration of TOM was measured near the western littoral zone (10.22) while the mouth part and north-eastern area was characterized by the lowest values (2.65% and 4.69). A total of 31658 individuals belonging to 12 families and 14 species were identified. Polychaeta with 3 species was the most dominant group in terms of abundance. The four most abundant taxa making up 85% of all specimens wereStreblospio gynobranchiata, Tubificidae,Hediste diversicolor andAbra segmentum. The maximum density (7,893 ind/m2) was obtained at station 1 while the minimum (1,777 ind/m2) was observed at station 16. The western area was characterized by the highest species diversity (H', 1.94) and the stations 10, 8 and 7 were characterized by the lowest diversity indices (H', 0.72, 0.77 and 0.87, respectively). The PCA showed that water parameters with more temporary variations had a greater significance in explaining the system variability, and a not marked but evident difference between the two parts of Gorgan Bay was observed and supported by nmMDS test. So Gorgan Bay presents transitional macrobenthic assemblages that are spatially distributed along substrate gradients but it seems that the coastal ecosystem of the south Caspian Sea and mouth-eastern part of Gorgan Bay is very dynamic and some species are forming a metapopulation toward western sites.

Keywords:Spatial, Temporal, Benthic macrofauna, Metapopulation, Gorgan Bay