Abstract
Understanding niche differentiation and interspecific coexistence patterns of macrobenthos provides critical references for ecological assessment of bivalve aquaculture. Seasonal field surveys (Autumn 2023–Summer 2024) across 12 sampling sites were conducted in an intensive Meretrix meretrix tidal flat farm in Rudong, Jiangsu. Synchronous water depth, sediment grain size and total organic carbon (TOC), and salinity were measured at each station to characterize habitat gradients; no unfarmed reference tidal flats were included, so all community patterns only reflect within-farm variation. The study area has sustained 10-year continuous bottom culture with baseline juvenile stocking density of 250 ind./m2 and twice-yearly manual harvesting. A total of 58 macrobenthic taxa were identified, and 13 dominant species were screened with a dominance index threshold Y > 0.02. All abundance data were fourth-root-transformed before Bray–Curtis community analysis to reduce bias of hyper-abundant taxa. Standardized Shannon niche breadth (0–1 normalized, r = 48 site–season composite units) and permutation-tested (999 replicates) Pianka niche overlap/similarity coefficients were calculated to quantify resource utilization similarity. Niche breadth ranged from 0.018 to 0.971, with farmed M. meretrix showing the maximum value, indicating broader seasonal resource adaptation. High niche overlap (>0.50) was detected between M. meretrix and four wild dominant species (Nephtys californiensis, Moerella iridescens, Lingula anatina, Umbonium thomasi). Notably, niche metrics derived solely from spatial abundance data cannot directly confirm interspecific competition or trophic segregation; the stable coexistence observed in our assemblage is tentatively hypothesized to be driven by vertical sediment partitioning and divergent feeding guilds based on published functional studies, without direct in situ trophic or stratified sediment evidence from this work. Hierarchical UPGMA clustering (0.35 Bray–Curtis distance threshold) and corrected NMDS ordination (48 site–season units, dominant-species fitted vectors) illustrated distinct functional groups. Based on the observed niche convergence among filter-feeding bivalves under local farming conditions, the range of 230–280 ind./m2 is presented as an observational reference. Controlled manipulative experiments with gradient stocking densities are needed to test whether this range can balance aquaculture yield and benthic biodiversity.