Trophic structure of micronekton in the Northwest Pacific Ocean

Introduction

The Northwest Pacific is the largest marine ecosystem on Earth (Brodeur et al., 1999) and is located where two strong western boundary currents (the Kuroshio Current and Oyashio Current) converge; as a result, according to the FAO (Food and Agriculture Organization of the United Nations), this is the location of the richest fishery in the world. Since reaching over 20 million tons in 1984, the total output of the fishery in the Northwest Pacific Ocean has been stable and was about 19.4 million tons in 2020; this accounted for 21.2% of the total global marine catch (FAO, 2022)。Many economically important species are found here in abundance, including Pacific sardine (Sardinops sagax), skipjack tuna (Katsuwonus pelamis), Pacific saury (Cololabis saira), Japanese flying squid (Todarodes pacificus) and neon flying squid (Ommastrephes bartramii) (Yatsu et al., 2013;Chang et al., 2019)。Smaller but more abundant species are found in mesopelagic waters, typically at depths of 100-1000m (Irigoien et al., 2014)。The majority of mesopelagic organisms have a diel migratory feeding habit, migrating upward to the surface region at night and later migrating downward a few hundred meters to the depths where they remain during the day (Salvanes and Kristofersen, 2001;Catul et al., 2011)。In recent years, many scholars have conducted research on biological resources and fishing grounds, including research on species such as the neon flying squid and Pacific saury (Chen et al., 2010;Huang and Huang, 2015)。In previous studies on the trophic ecology of organisms in the Northwest Pacific Ocean,Lee et al. (2012)focused on the trophic structure and biomass of plankton communities, and found that there was a positive correlation between phytoplankton biomass and both ciliates and mesozooplankton biomass and the increase of mesozooplankton biomass was closely related with ciliate and dinoflagellates biomass.Ohshimo et al. (2019)researched different oceanographic processes as well as trophic ecology and movement of large fishery organisms in the Northwest Pacific Ocean based on stable isotope analysis to map spatial isopleths of top predator prey organisms including fishes and cephalopods and found that the differences in δ¹³C and δ¹⁵N values in the Northwest Pacific might be influenced by production dynamics, ocean currents and upwelling, and the degree of carbon and nitrogen fixation in the atmosphere.Van der Lingen et al. (2009)investigated the trophic ecology and planktonic interactions of small pelagic fishes such as anchovy and Pacific sardine in the Northwest Pacific Ocean, as well as the possible effects of climate on some species, and suggested that changes in zooplankton community structure may cause species alternations between anchovy and sardine, while climate change has more indirect effects.

The trophic structure is one of the most important features of an ecosystem. The production and transfer of biomass by trophic relationships is the core function of an ecosystem, and the relationship between ecosystem productivity and niche size is an important factor in determining ecosystem function. Moreover, an understanding of species distribution and trophic structure gives a better insight into the ecology of a community (Andrade et al., 2016;Lesser et al., 2020)。In the study of community trophic structure, stomach content analysis (SCA) and stable isotope analysis (SIA) are often used. SCA can only provide information about how an organism was feeding within the short period of time before it was caught (Ibáñez et al., 2021), whereas SIA can provide longer-term information about its feeding (Li et al., 2014)。In addition, different tissues tend to be synthesized and replaced at different rates, which means that SIA can also provide information about temporal variations in individuals’ feeding and TLs. For example, muscle (Malpica-Cruz et al., 2012) can provide isotope information covering periods of a few months to a year, whereas metabolically inert tissues such as fish vertebrae (Carlisle et al., 2015) can usually record biochemical information over the entire period of an individual’s development, from which information about changes in feeding habits and habitat can be derived (Inger and Bearhop, 2008)。

Additionally, SIA can be used to study biological carbon sources, energy sources and feeding habits, as well as temporal and spatial trophic relationships and migration, thus revealing how substances cycle through ecosystems and details of the trophic relationship between consumers and producers in the food web (Hansson et al., 1997;Zeng et al., 2008)。The carbon stable isotope value (δ¹³C ) increases negligibly from a consumer’s diet (0-1‰); as it reflects the source of primary production, this value can be used to discriminate between inshore and offshore species or between pelagic and benthic feeders (凯利et al ., 2006;Liu et al., 2019)。The stable nitrogen isotope value (δ¹⁵N ) is usually used to estimate the trophic level (TL) of consumers due to¹⁵N enrichment with trophic transfer. Compared with the complexity of previous methods such as SCA, SIA, by virtue of its versatility and simplicity, can improve the understanding of elemental and material-energy cycles in ecosystems (Zhang et al., 2019)。

The description of a niche derived using the stable carbon and nitrogen isotope values of biological tissues not only gives information about the spatial position occupied by a species within a certain environmental range at a certain time, it also reflects the impact of environmental and other factors on the species as well as the material circulation and energy flow within the ecosystem or community. The niche breadth reflects the diversity of the food resources available to a species: the greater the breadth, the stronger the environmental adaptability of the species. The niche overlap also reflects the degree of competition and the similarity of resource utilization between different species (Hurlbert, 1978;Bearhop et al., 2004)。The isotopic niche can reflect variability associated with TL, primary production sources, and possibly isotope baselines in the case of migratory species (Layman et al., 2007;Newsome et al., 2007)。

Carbon and nitrogen stable isotope analysis has been widely used in research into trophic structure, trophic levels and seasonal variations in relation to biological communities in fisheries (Jennings et al., 2002)。KoichiYoshii et al. (1999)used carbon and nitrogen SIA to analyze the structural features of the pelagic food web in Lake Baikal and studied the composition of the diets and changes in foraging habits of individual organisms.Huang et al. (2019a)established a continuous trophic spectrum for the main fishery in the South China Sea central and western fishing grounds. However, the number of studies on the trophic structure of micronekton in the Northwest Pacific Ocean has been limited (but seeOhshimo et al., 2019)。

Our study was based on micronekton that were collected during a fishery resource survey in the Kuroshio−Oyashio mixing zone in March and September 2019. The stable carbon and nitrogen isotope values of these organisms were measured and their trophic levels estimated. These results were then used to establish a continuous trophic spectrum and combined with isotopic niche to research interspecies differences in community trophic structure and potential impacts on ecosystems.

Material and methods

Sample collection

The sampling vessel used in this study was Shanghai Ocean University’s oceanic fisheries resources survey vessel “Song Hang”. This vessel has a total length of 85m, a moulded depth of 8.7 m and a displacement of 3271.4 tonnes when fully loaded. A four piece mesopelagic trawl net with a single-capsule structure was used for all collections; the overall dimensions of the net were 434 m × 97.1m (44.98 m). The net port part had large mesh and the net body part had a machine-woven mesh. Double-leaf mesh board had single-hand steel connection. Zooplankton were collected using a Bongo net following the same collection protocol for each survey. Specific sampling information on the trawl and sampled zooplankton is available inTables 2and3.

Fish, cephalopod and shrimp samples were collected by trawling within the Kuroshio−Oyashio mixing zone in the Northwest Pacific from March 21st to March 23rd and from September 5th to September 9th, 2019. The individual sampling sites were located at 148°E, 33°N; 148°E, 31°59′N; 148°5′E, 38°45′N; 150°13′E, 35°56′N and 148°2′E, 34°43′N (Table 1andFigure 1) with tow depths of 150m, 100m, 409m, 604m and 243m, respectively (Table 2)。After classification and identification at sea, the samples were stored in a freezer at −20°C.

The zooplankton samples were obtained by vertical trawling using a large plankton net at a depth of 100 m. The trawling speed was 0.5−0.8 m/s. Trawling was carried out for 30 minutes at each of the five sites listed above (Table 3andFigure 1)。Trawl-collected samples were washed into 1000ml jars using seawater, and then filtered through a 100 um mesh sieve to end up with approximately 450ml of remaining zooplankton, which were also stored in a freezer at -20°C.

Sample processing

After using a microscope to pick out copepods in the laboratory, the copepod samples were filtered using GF/C filter membranes. They were then wrapped in tinfoil and dried to a constant mass in a drying oven at 60°C before being ground to a fine, homogeneous powder using a RETSCH MM200 automatic ball mill in the laboratory. After sieving, between 1 and 2 mg of this powder was then put into a 0.3-mg tin capsule for preservation before further testing. Before being dried, the GF/C filter membranes used for concentrating the samples were heated in a muffle furnace at 450°C for 4 hours to remove any organic matter present.

A total of 90 fishes, 162 cephalopods and 5 shrimp were collected (Table 1)。After thawing and cleaning the trawl-caught organisms in the laboratory, the determination of biological data and muscle sample collection were carried out in accordance with the Specifications for Oceanographic Survey——Part 6: Marine biological survey (GB/T12763.6-2007) (China National Standardization Management Committee, 2007)。Standard (fishes) and mantle (cephalopods) length ( ± 0.01 mm) and weight ( ± 0.01 g) were measured and recorded for each organism. Approximately 10 g of soft tissue (dorsal muscle for fish, abdomen muscle for shrimp and mantle muscle for squid) were obtained from selected individuals to determine the δ¹³C and δ¹⁵N values. All of the selected samples were rinsed with Milli-Q ultrapure water for 5 minutes to remove possible contaminants; they were then freeze-dried in German Christ Alpha 1-4 LSG freeze-dryers at −50°C for at least 24 hours before being ground to a fine, homogeneous powder using an automatic RETSCH MM200 ball mill. Between 1 and 2 mg of this powder was then put into a 0.3-mg tin capsule for preservation before further testing.

Stable carbon and nitrogen isotope analysis

The stable isotope values were determined at the Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai, China, using an ISOPRIME 100 isotope-ratio mass spectrometer (Isoprime Corporation, Cheadle, UK) and a vario ISOTOPE cube elemental analyzer (Elementar Analysensysteme GmbH, Hanau, Germany). The standard substance used for determining δ¹³C was Vienna Pee Dee Belemnite (V-PDB) and the standard substance used for determining δ¹⁵N was atmospheric nitrogen (N₂)。During the sample testing, in order to ensure the stability of the instrument and the accuracy of the test results, three standard isotope samples were added to every tenth sample for testing with an analytical precision of ±0.06‰. The values of δ¹³C , δ¹⁵N , %C and %N were measured and the δ¹³C and δ¹⁵N measurements were made with a precision of 0.2‰. Stable isotope values were calculated using the following formula:

where X is¹³C or¹⁵N and δ is the ratio of the heavy isotope to the light isotope in the sample. R_sampleand R_standardare the atomic ratios of¹³C/¹²C and¹⁵N/¹⁴N在样本和标准物质,分别ly.

考虑到弱泳生物通常是脂质丰富,we chose the correction model developed byHoffman and Sutton (2010)for deep-sea fishes to correct the δ¹³C values of micronekton with C: N >3.76 in this study.

Trophic level calculation

The trophic level of each sample was estimated and a continuous trophic spectrum was constructed. The formula used to estimate the trophic level was

where TL_baserepresents the trophic level of the primary consumer; its value was set as 2 (Du et al., 2020)。The nitrogen stable isotope values of the consumer and copepod baseline are represented by δ¹⁵N_consumerand δ¹⁵N_base, respectively. In this study, copepods obtained near the same site location and during the same season as the consumer species were selected as the baseline organisms (Table 4)。我们选择使用桡足类作为基线形式sm to reduce variability associated with other zooplankton components that could affect baseline value estimates (e.g., size, trophic position, biochemical composition). TEF is the trophic enrichment factor. In this study, the TEF was assumed to be 3.4‰ (Minagawa and Wada, 1984;Post, 2002)。

Isotopic niche estimation

The SIBER package (Jackson et al., 2011) in R 4.0.4 (R Core Team, 2022) was used to draw the isotopic niches of the micronekton. The niche areas were then calculated to determine the niche breadths, and the SIAR package in R 4.0.4 was used to calculate the isotopic niche overlap. The standard ellipse area corrected for small sample sizes (SEAc) contained about 40% of the stable carbon and nitrogen isotope data, which meant that it could be used to represent the core niche of each individual and also to calibrate the sample sizes (Layman et al., 2007;Layman et al., 2012)。It could also be used to estimate the breadth of the isotopic niches and the overlap between the isotopic niches of different species. The overlap values ranged from 0 to 1 using a scale proposed byLangton (1982)in which values from 0 to 0.29 indicate a low degree of overlap, values from 0.30 to 0.60 indicate medium overlap and values higher than 0.60 indicate high overlap.

Data analysis

Excel, Sigmaplot 12.5 and R were used for the data analysis, drawing and calculations. The difference in δ¹³C and δ¹⁵N弱泳生物的不同组织之间tested using analysis of variance (ANOVA). Prior to data analysis, normality and homogeneity of variance tests were run. If the results of the aforementioned tests were significant, nonparametric tests were run to determine the levels of significance of the differences. Multiple comparisons were then carried out for fishes, cephalopods and shrimp including Bonferroni corrections. The correlation between isotope values and standard length (SL; fishes) or mantle length (ML; cephalopods) were tested using linear regression.

Results and analysis

Stable carbon and nitrogen isotope values

A total of 43 micronekton species from 16 families and 33 genera were collected, including 30 fish species, 12 cephalopod species and 1 shrimp species. Enoploteuthidae were the most abundant, accounting for 58.37% of all samples and the micronekton in Stomiidae were the most diverse, accounting for 30.23% of all micronekton. The most common species found in the samples were the fishesBenthodesmus tenuis,Chauliodus pammelas,Chauliodus sloani,Lestidium prolixum, the cephalopodsEnoploteuthis chunii,Abralia similisand the shrimpOplophorus gracilirostris(Table 1)。

The δ¹³C值介于−21.9‰−17.9‰(−19.8‰±0。6‰) and the δ¹⁵N values from 7.0‰ to 12.4‰ (10.0‰ ± 1.0‰). For the fish samples, δ¹³C ranged from −21.9‰ to −18.7‰ (−20.0‰ ± 0.5‰) and δ¹⁵N ranged from 7.0‰ to 12.4‰ (9.6‰ ± 1.2‰); for the cephalopods, δ¹³C ranged from −21.3‰ to −17.7‰ (−19.7‰ ± 0.7‰) and δ¹⁵N ranged from 8.2‰ to 12.2‰ (10.3‰ ± 0.8‰); and for the shrimp, δ¹³C ranged from −20.4‰ to −19.5‰ (−20.0‰ ± 0.4‰) and δ¹⁵N ranged from 7.6‰ to 10.6‰ (9.3‰ ± 1.1‰) (Table 5)。In addition, the copepod δ¹³C values ranged from −23.4‰ to −21.0‰ (−23.2‰ ± 0.2‰) and the δ¹⁵N values from 2.4‰ to 4.5‰ (3.5‰ ± 1.1‰). Detailed information of each micronekton is presented inTable 5.

Population variance of δ¹³C and δ¹⁵N values violated the assumptions of homogeneity of variance (P_C=0.013<0.05;P_N=0.001<0.05), and the results of the non-parametric tests showed that the difference in δ¹³C and δ¹⁵N values between the fish, cephalopod and shrimp groups were significant (F_C=10.076,P_C=0.006<0.05;F_N=32.472,P_N=0.000<0.05). Multiple comparisons showed significant differences in δ¹³C and δ¹⁵N values between fish and cephalopod groups (P_C=0.006<0.05,P_N=0.000<0.05), while δ¹³C and δ¹⁵N values between shrimp and fish groups (P_C=1.000>0.05,P_N=1.000>0.05) and between shrimp and cephalopod groups (P_C=0.922>0.05,P_N=0.121>0.05) were not significant. In addition, for cephalopods, there was a significant correlation between ML and δ¹³C values (F=22.652,P=0.000<0.05), while the correlation between δ¹⁵N values and ML was not significant (P=0.585>0.05); for fishes, neither δ¹³C or δ¹⁵N values were significantly correlated with SL (P_C=0.074>0.05;P_N=0.849>0.05).

Trophic levels and continuous trophic spectrum

The continuous trophic spectrum that was constructed (Figure 2) showed that, for the samples obtained, the TL values ranged from 2.67 (forGrimalditeuthis bonplandi) to 4.80 (forEchiostoma barbatum); the average TL value for cephalopods was 3.30 ± 0.34, for fishes it was 3.66 ± 0.59 and for shrimp it was 4.04 ± 0.32. A total of twenty-one species (48.8%) had TL values ranging from 3.0 to 4.0. The TL values of 19 species were greater than 4.0 − these species included seventeen species of fish (among which wereScombrolabrax heterolepisandChauliodus sloani), one species of cephalopod (Histioteuthis inermis), as well as the shrimpO. gracilirostris. In addition, there were three species with TL values greater than 4.5:Photostomias guernei,Diaphus gigasandEchiostoma barbatum.

Isotopic niches

Among different groups of micronekton, the fish group had the highest SEAc value − 1.80‰². Cephalopod and shrimp groups had relatively low values of 1.59‰²and 1.73‰², respectively (Figure 3A)。There was a moderate amount of trophic overlap between two of the three groups, with the trophic overlap between cephalopod and fish being the greatest (cepholopod-fish: 0.55; cephalopod-shrimp: 0.33; fish-shrimp: 0.45).

我们选择七个物种相对较大的sample sizes (n≥5) among all sampled species, which included two cephalopods, four fishes and a shrimp, to compare their isotopic niche. Among different species,O. gracilirostrishad the largest SEAc value of 1.73‰², whereasB. tenuisandLestidium prolixumagain had low values of 0.45‰²and 0.46‰², respectively. In this case,Enoploteuthis chuniiandChauliodus pammelas(0.53) orAbralia similis(0.51) had a greater degree of overlap. Trophic overlap was moderate betweenChauliodus pammelasandE. chuniiorA. similisorO. gracilirostris, and betweenE. chuniiandO. gracilirostris, respectively. In addition, there was relatively low trophic overlap betweenB. tenuisandA. similisorC. pammelasorL. prolixum(Table 6andFigure3B)。

Discussion

Stable carbon and nitroge n isotope values

δ¹³C can be used to gain information about the sources of carbon and energy used by organisms (Peterson and Fry, 1987;Fry, 2006)。The range of δ¹³C values (CR) can be used to estimate the diversity of food sources in the food web; the higher its value, the greater the diversity (Madurell et al., 2008;Jackson et al., 2011)。In this study, the range of δ¹³C values was 4.1‰, meaning that a wide range of food sources was available. The confluence of the Kuroshio Current and Oyashio Current brings a rich composition of nutrients and species (Nishikawa et al., 2020), thus the range of food sources was complex and the degree of feeding specialization by the different types of fishes, cephalopods and shrimp was relatively low. As a result, the differences in δ¹³C among fishes, cephalopods and shrimp were significant. In addition, in this study, the difference in δ¹³C among fishes, cephalopods and shrimp may be related to different vertical and horizontal ranges of some species (Sherwood and Rose, 2005)。

The range of δ¹⁵N values (NR) mainly reflects the trophic length of the community: the greater NR is, the greater the trophic length of the community and the more TLs there will be in the community (Jackson et al., 2011)。Ohshimo et al. (2019)investigated a total of 84 species of pelagic forage fish and squid in the Northwest Pacific Ocean. CR had a range of 6.3‰ (from −21.5‰ to −15.2‰) and NR had a range of 10.5‰ (from 4.3‰ to 14.8‰). In comparison, the ranges of CR values in our study were relatively low, which may be related to the smaller benthic sample size and small sampling range that was used. Benthic organisms tend to have higher δ¹³C values than pelagic organisms (Le Loc'h et al., 2008;Kopp et al., 2015)。The range of NR values was also quite small due to the lack of large predators.

Trophic levels

When stable isotope techniques are used to investigate micronekton food webs and TLs, the selection of the δ¹⁵N baseline and the trophic enrichment factor (TEF) affects the determination of the trophic positions. In recent years, primary consumers, including zooplankton, have mainly been chosen as baseline species (Post, 2002)。Chiba et al. (2012)estimated the average copepod (Neocalanus cristatus,Neocalanus flemingeriandNeocalanus plumchrus) δ¹⁵N values in Oyashio and the results gave values of 7‰−10‰ with δ¹⁵N variation often associated with shifts in feeding strategies (i.e., between herbivorous and omnivorous/carnivorous). In our study, the average δ¹⁵N values of copepods obtained in March (3.5‰) and September (6.7‰) were lower thanChiba et al. (2012)estimates. Copepods sampled byChiba et al. (2012)may have included more omnivorous individuals that would be feeding at a higher TL compared to the copepods in our study. In addition, the low sea surface temperature (SST) in March compared with September (SST_Mar: 18.00°C; SST_Sep: 26.37°C) (https://resources.marine.copernicus.eu/products), could affect photosynthesis and nitrogen absorption as well as the fractionation of phytoplankton (Takai et al., 2000)。As zooplankton feed on various kinds of phytoplankton, this will ultimately be reflected in the difference of copepod δ¹⁵N values between March and September in this study (Xu et al., 2010)。

在这项研究中,采样弱泳生物它d to three TLs with the average TL value being 3.44. Other studies have shown that, on average, there are between four and six TLs between the primary predators and top predators in marine ecosystems (Arreguín-Sánchez et al., 1993;Bulman et al., 2002), and the average TL value of the global marine ecosystem based on δ¹⁵N has been estimated as 4.0 ± 0.5 (Vander Zanden and Fetzer, 2007)。Compared with these results, the samples in this study corresponded to a lower number of TLs and the average TL value was also lower. This can be explained by the fact that the samples mainly consisted of micronekton and did not include large fishes from high TLs.

In addition, some of the listed TL estimates are a bit surprising given the physiology of the species and their relative TL estimates, especially some of the lowest estimated TLs. For example, the deep-sea squidGrimalditeuthis bonplandihas the lowest estimated TL of< 3.0 whereasHoving et al. (2013)noted that this species was believed to prey on crustaceans and cephalopods. Relatedly,Ommastrephes bartramiihad the second lowest TL estimate of the dataset but others have observed high TL for this species (Sthenoeuthis oualaniensisoccupied trophic level about 2.72~3.08 and adultO. bartramiiwere 1.3TLs higher thanS. oualaniensis) (Parry, 2008;Huang et al., 2019b)。Combined with the sampling location and sampling date of theO. bartramiiin this study, as well as its migratory route throughout the life cycle (Ichii et al., 2009), we believe that theO. bartramiiin this study is in the juvenile stage. In this stage the species tends to prey not only on zooplankton, but also small fishes and cephalopods (Watanbe et al., 2004)。As the individual grows, its ML and beak size increase, which means that they can prey on fishes and cephalopods with higher TLs (Jin and Chen, 2014)。Size explains the lower TL ofO. bartramiiin this study compared with other studies, but the TL estimated in this study forO. bartramiiis still lower than what would be expected with their presumed diet.O. bartramiimigrations cross different habitats throughout its life cycle (Yu et al., 2015), which could create a mismatch between the baseline isotope value used for our TL calculation and the baseline reflected inO. bartramiimantle tissue. The use of sampling site baseline values introduces some spatial bias if the subjects migrate to the sampling sites within a short period of time before being caught, which also introduces some error in the TL measurements (Young et al., 2015)。In addition, Ammonical cephalopod species (e.g.,G. bonplandiandO. bartramii)也可能有不同的组织生化混合涂料sition that could bias TL estimates for these species as has been observed with urea retention in elasmobranchs (Kim and Koch, 2012)。

Moreover, in the trophic level results (Figure 2), we also found that the TLs of some mesopelagic fishes (3.45 to 4.80) were almost the same as those of large oceanic predators like tunas (3.0~4.8) (Sarà and Sarà, 2007) or sharks (3.8~4.1) (Cortés, 1999), which implies that the TLs may be overestimated in this study.Richards et al. (2020)indicated that the δ¹⁵N values in POM increase with depth and we assume that baseline δ¹⁵N values at depth may be higher than surface waters. Considering the copepods used in our study largely reflect the δ¹⁵N values of surface waters, using copepods as the TL baseline could result in artificially high TL estimates for mesopelaic species.

Myctophidae and Stomiidae are considered byChoy et al. (2012)to be the two dominant groups of meso-pelagic fish consumers.Kenaley (2012)showed that Stomiidae is capable of consuming prey whose SLs are greater than 50% of their SLs due to their large jaw and wide gape, which also explains that in our study, compared with other species of the same size such asAphanopus carbo,B. tenuisandLestidium prolixum, Stomiidae occupied higher trophic positions. Myctophids undertake daily vertical migrations and have a diverse prey composition including copepods, seston, fishes, and chaetognaths (Eduardo et al., 2021)。Olivar et al. (2019)estimated the TL of lanternfishes (Myctophidae) in the tropical and equatorial Atlantic and estimated lanternfish TL to be between 2 and 4. The higher TL estimates of myctophids in our study suggest that they fed on a higher proportion of prey like fishes rather than zooplankton.

Community isotopic niches

In this study, among different groups of micronekton, fishes had a larger SEAc than cephalopods and shrimp. In this study, there is no significant correlation between fish δ¹³C values and SL, indicating that the carbon source did not change significantly throughout the ontogenetic process. It is therefore considered that the baseline values of their feeding environment did not change significantly (Chen et al., 2021)。However, the majority of the fishes included in this study are meso- and bathypelagic species. According to Brierley (2014), most marine pelagic communities exhibit diel vertical migration. Compared with cephalopods and shrimp, the fishes in this study had a greater SEAc, which could be due to their broader range of SLs among species and corresponding foraging diversity. Additionally, the medium overlap among fish, cephalopod and shrimp groups suggests a more competitive relationship between them regarding the use of resources as well as a certain degree of similarity.

Among different species of micronekton,O. gracilirostrishad the largest SEAc. According toBurdett et al. (2017), the Oplophorid shrimp species are extensively dispersed throughout the meso-pelagic zone of every ocean. Shrimp are omnivorous and their main sources of food include algae, organic detritus, phytoplankton, zooplankton and zoobenthos as well as some types of fish and other shrimp (Wilcox and Jeffries, 1974) and also serve as important prey organisms for top predators such as tunas and squid (Hopkins et al., 1994)。O. gracilirostrisis a vertically migrating deep-sea species which means that they can easily adapt to the environment of different habitats and have a large area of activity (Chan and Yu, 1986)。B. tenuisis a benthopelagic species whose low SEAc value may be explained by a relatively constrained habitat and diet (Nakamura and Parin, 1993)。The SEAc ofC. sloaniwas small, which is consistent with the conclusion drawn byEduardo et al. (2020)that the value ofC. sloani’s SEAc is limited by the diversity of the species’ prey: they mainly feed on zooplanktivorous fishes, of which myctophids account for more than 50%.C. sloaniare meso- and bathypelagic fish. Temperature affects the vertical migration mode ofC. sloani, limiting its food sources and habitat to a certain extent.L. prolixumare small to medium-sized deep-sea fish whose small SEAc value may be related to its relatively narrow SL range.

The values of the niche overlap betweenB. tenuisandL. prolixumwere low; the same was also true for the overlap betweenB. tenuisandC. pammelasorA. similis. These patterns indicate that there is a certain degree of feeding separation betweenB. tenuisand the other three species (L. prolixum,C. pammelasandA. similis) in the ecosystem, thus the degree of competition for resources and habitats between these species is low (Wang et al., 2018)。However, the higher niche overlap betweenE. chuniiandA. similisorC. pammelasreflect that there is a greater degree of interspecific resource competition betweenE. chuniiand these two other species.

Conclusion and caveats

In this study, finding an appropriate baseline species for the purposes of TL estimates was challenging and is a major potential source of error in TL calculations. Both temporal and spatial (vertical and horizontal) variation in the baseline used can have an impact on the calculated TL values. The copepods used in this study may artificially overestimate the TL of mesopelagic species due to vertical baseline differences. In addition, ammonium accumulation in deep water cephalopod muscle tissues may affect δ¹⁵N values and then influence TL estimates.

Due to the limits on the trawling frequency, the trawling equipment used and the number of sampling sites this study did not sample all of the existing biological communities present. Sample sizes of the species that were collected were often small, so even the species represented in this study were not always sampled comprehensively. Samples in this study were collected during different months before pooling. In the future we will expand the temporal and spatial range of the sampling, as well as the number of trawling sites, to collect larger sample sizes over a larger regional scale. Considering that individual copepod size was not measured at each site, and that only one stable isotope value of copepods was measured at each site, in future studies we will collect a larger number of zooplankton, and then conduct species identification of zooplankton at each site. Given the potential errors associated with TL estimates based on bulk δ¹⁵N and surface zooplankton baseline samples, future studies would benefit from the inclusion of deep water particulate organic matter (POM) as an alternative baseline or the use of amino acid compound specific stable isotope analysis (AA-CSIA) of micronekton tissues as an alternative approach to TL calculation (Popp et al., 2007;McMahon and Newsome, 2019)。

Despite these caveats, our study provides general trophic information on the micronekton community of the Northwest Pacific Ocean. Our results have shown that the degree of feeding specialization by the different types of fishes, cephalopods and shrimp was relatively low. As the two dominant groups of meso-pelagic fish consumers, the species of Myctophidae and Stomiidae occupied higher TLs whileOplophorus gracilirostris,Enoploteuthis chuniiandAbralia similisalso played important roles in the transfer of energy and materials in the meso-pelagic ecosystem.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics statement

本月的动物研究是进行审核和批准itutional animal care and use Committee of Shanghai Ocean University.

Author contributions

JZ contributed to the idea and experimental design. BL provided samples. JZ completed the experiments and statistical analyses and data interpretation with the help of BL and wrote the manuscript with the help of BL, SH. YG provided the data of zooplankton stable carbon and nitrogen values. All authors contributed to the article and approved the submitted version.

Funding

This work was sponsored by Program on the Survey, Monitoring and Assessment of Global Fishery Resources (Comprehensive scientific survey of fisheries resources at the high seas) sponsored by the Ministry of Agriculture and Rural Affairs; the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning under contract 0810000243.

Acknowledgments

We thank those people for helping to collect, process samples, and analyze the isotopic signatures. We thank Project on the Survey and Monitor-Evaluation of Global Fishery Resources sponsored by Ministry of Agriculture and Rural Affairs.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at://www.gosselinpr.com/articles/10.3389/fmars.2022.971649/full#supplementary-material

References