We test a current assumption that Atlantic cod (Gadus morhua) off Newfoundland and Labrador, Canada, do not feed during the protracted spawning season (March-September). Stomach contents were analysed from 10 473 cod from four Northwest Atlantic Fisheries Organization regions (2J, 3K, 3L, and 3Ps) over 9 years from which gonads were also analysed to determine sex and maturity status. Adult cod in spawning condition did feed in all regions, usually at rates equivalent to or even greater than non-spawning fish and juveniles. Both sexes fed during spawning, though females consumed lesser amounts. Regional differences were evident. The total fullness index was greater in the northern (2J) than the southern (3Ps) region, with no consistent differences between spawners and non-spawners. The most southerly region (3Ps) exhibited the greatest prey diversity, the northern region (2J) the least. Shrimp was the major diet item in the northern regions. Capelin, zooplankton, crab, and other fish increased in importance to the south. Differences in prey items between non-spawning and spawning individuals of both sexes were possibly related to spawning behaviour. Models using consumption rates should not assume that cod do not feed during the protracted spawning season in these waters.

Species- and size-selective overexploitation often have ecosystem-wide impacts that are evident in community size spectra. To both derive potential ecosystem targets for community rebuilding and assess contemporary indicators relative to these targets, we constructed theoretical size spectra to predict pristine biomass densities using a combination of species- and size-specific nitrogen stable isotope signatures and a range of trophic efficiencies and primary productivity estimates within and among 3 sub-regions of the Newfoundland and Labrador shelf. Theoretical size spectra were compared to empirically derived size spectra using trawl survey data. The descending slopes of the empirically determined size spectra were between 1.25 and 2.42 times steeper than the theoretical slopes. The percentage of the theoretical distribution represented by the empirical size structure ranged between 1.1 and 29.4%, with the closest and furthest estimates associated with the smallest- and largest-sized fishes, respectively, and strongly influenced by estimates of trophic efficiency and primary productivity. Regional variation was also observed, with southern regions reaching 1.3-32.3% of the theoretical biomass density and the northernmost region reaching only 0-8.3%. Importantly, the descending slopes varied depending on fish size, with biomass density of larger sizes decreasing faster than that of smaller sizes. Variations among sub-regions and fish guilds were also observed. These analyses provide a means to derive potential ecosystem targets and indicators through which recovery of fish communities can be monitored and assessed.

Exploitation and changing ocean conditions have resulted in altered species interactions and varied population dynamics within marine fish communities off northeast Newfoundland and southern Labrador, Canada. To understand contemporary species interactions, we quantified the isotopic niches, niche overlap, and ontogenetic niche change among seven dominant fish species using stable isotope analyses. Analyses used fishes from three regions differing in fish and prey diversities. Differences in fish and diet composition diversity among regions were found using Simpson’s inverse diversity index. The regions of lowest diversities had higher instances of niche overlap and higher percentage of niche overlap area. The region of highest diversity had the widest spread of niches with greater distances from the community centroid. Ontogenetic shifts were observed such that larger individuals shifted towards the community centroid with the exception of Atlantic cod. Atlantic cod in particular was found to consistently be the top predator of the analyzed species. Our results reveal: (a) overlap in isotopic niches and spread within niche space was correlated with fish and diet diversity; (b) ontogenetic shifts are important when considering a species’ niche and quantifying spatial variation in community niche profiles.

Spatial heterogeneity in food web structure and interactions may reconcile spatial variation in population and community dynamics in large marine ecosystems. In order to assess food web contributions to the different community recovery dynamics along the Newfoundland and Labrador shelf ecosystem, we quantified species interactions using stable isotope mixing models and food web metrics within three sub-regions. Representative samples of each species caught in trawls and plankton tows were analyzed for stomach contents and stable isotope ratios (δ15N and δ13C) to parameterize isotope mixing models. Regional variation, highlighted by the diets of three economically important species, was observed such that the southern region demonstrated a variety of trophic pathways of nutrient flow into the higher food web while the diets of fish in the northern regions were typically dominated by one or two pathways via dominant prey species, specifically shrimp (Pandalus sp.) and hyperiids. Food web metrics indicated that the low-diversity northern regions had higher connectance and shorter food chain lengths. This observed regional variation contributes to our understanding of the role of specific forage species to the ecosystem which is an essential contribution towards ecosystem-based management decisions.