Abstract
Calanus glacialis/marshallae is a dominant zooplankton species in the Pacific Arctic Ocean that is widely distributed in shelf areas, and it plays a vital role in connecting primary production to higher trophic levels. Its phenology is well adapted to hydrography, but there is little available information about regional and diel changes in population structure and grazing features. In this study, we investigated C. glacialis/marshallae during autumn 2019 in the Eastern and Northeastern Chukchi and Canadian basins to reveal geographic and diel variations in population structure, body size, grazing activity, and fatty acid composition. The abundance of C. glacialis/marshallae was found to be high on the slopes and low on the shelves. Body size (prosome length) was well described by the Bělehrádek equation combined with in-situ temperature throughout the sampling region. Cluster analyses based on hydrographic parameters were divided into four regions: southern shelf, northern shelf, slope, and basin. The southern shelf was dominated by copepodite stage five (C5) transported from the Bering Sea by Pacific waters. C4 and C5 were dominant on the northern shelf, suggesting that they grew slower than those on the southern shelf, and the populations also exhibited higher concentrations of fatty acids originating from dinoflagellates than those originating from the pan-Arctic Ocean, indicating low productivity in the region. The population on the slope had the highest abundance, C4 was dominant, and large amounts of diatom-derived eicosapentaenoic acid (EPA). These features are attributed to the upwelling of populations and nutrients that support diatom growth. In the basin, the early copepodite stages of composition were distinctly higher than those recorded in previous studies, because larger amounts of organisms flow into the region, resulting in more extended reproduction periods. In the basin, small and large forms of C5 were simultaneously found, and the small form exhibited a diel grazing activity pattern, but the large forms did not. These findings suggest their well adaptation in changing of the Pacific Arctic Ocean.
1 Introduction
Calanus Leach, 1816 is the dominant copepod genus in the zooplankton biomass within the Pacific Arctic Ocean (). Three Calanus species have been identified in different areas of the Pacific Arctic Ocean: Calanus hyperboreus in the basin (), Calanus glacialis on the shelf (; ), and Calanus marshallae in the southeastern Bering Sea (; Vidal and Smith, 1986). Of these, C. hyperboreus has the largest body size and can be easily distinguished from the other two species. Calanus glacialis and C. marshallae are morphologically similar and difficult to distinguish (), although they have been distinguished by genetic analysis (haplotype-based on 18S rRNA), and results have indicated that C. marshallae is transported to the shelf by the Pacific waters (Nelson et al., 2009; ). However, a genetic analysis was not conducted in the current study, and the C. glacialis/marshallae species complex is therefore referred to as C. glacialis in this work. Calanus glacialis has two populations (the Bering Sea and Arctic Basin populations that can also be distinguished using haplotype analyses (Nelson et al., 2009; ; Pinchuk and Eisner, 2017; ). Although these populations are believed to be morphologically indistinguishable, haplotype analysis indicates that the Bering Sea population is distributed on the Chukchi Sea shelf, and the Arctic Basin population is distributed north of the Beaufort Sea slope/Basin (Nelson et al., 2009). In addition, these populations appear to differ in their copepodid stage composition during the same season (; Pinchuk and Eisner, 2017), suggesting that the life cycles differ between populations.
Calanus glacialis is the dominant biomass species on the shelf of the Pacific Arctic Ocean (; ), and it plays an essential role in connecting primary production to higher trophic levels in fish, birds, and mammals (; ; ; ; ). The population peaks near the shelf border (the shelf break) and declines sharply in the basin (; Wassmann et al., 2015; ). Calanus glacialis has a generation length of 1–3 years (; ; ), and it grows from an egg to C3 and C4 within one year (Scott et al., 2000). The species has mixed reproductive strategies depending on food conditions (Plourde et al., 2005; ): capital breeding (based on internal energy storage before the onset of the spring bloom) and income breeding (utilizing ice algal or pelagic blooms) (Varpe et al., 2009). Due to this flexibility, their reproduction is initiated under sea ice and continues for one to two months using ice algae and pelagic phytoplankton blooms (Tourangeau and Runge, 1991; ; ). However, if individuals reproduce in the late season, they cannot grow to the diapause stage in a short period (; ). As the developmental time between the stages is negatively related to temperature under sufficient food conditions (; McLaren et al., 1988), it is acknowledged that the life cycle of this species exhibits plasticity in response to environmental conditions that are reflected in the biological features of the species, such as body size, grazing activity, and the accumulation of fatty acids.
The body size of copepods varies depending on the temperature and prey conditions (; Viitasalo et al., 1995; ). Low temperatures and high food availability stabilize the growth of copepods (Vidal, 1980), which increases the body size (Viitasalo et al., 1995). Geographic variations in the body size of this species have been reported in the Atlantic sector of the Arctic Ocean when conducting interspecific comparisons with Calanus finmarchicus (; ), but little information is available for the Pacific Arctic Ocean ().
Calanus glacialis is omnivorous, feeds mainly on diatoms and dinoflagellates (Søreide et al., 2008; Yeh et al., 2020), and adapts flexibly to changes in its prey composition (Stevens et al., 2004a; ; ). With respect to feeding, the fatty acid composition has been well analyzed, mainly in the Atlantic sector of the Arctic Ocean (Stevens et al., 2004a; ; ; Trudnowska et al., 2020). The fatty acid composition of zooplankton varies with diet (Stevens et al., 2004b). Zooplankton contains unique fatty acids that are taken up and preserved without chemical changes in the body; fatty acids can be analyzed to obtain an index of their feeding records (). For example, diatoms are rich in eicosapentaenoic acid (EPA), and dinoflagellates are rich in docosahexaenoic acid (DHA; ). The gut pigment is another index of the feeding activity of copepods, and it is useful for evaluating diel changes in grazing activity (); however, it has also not been well studied in the Pacific Arctic Ocean (Matsuno et al., 2015). Therefore, the regional differences in body size, feeding activity (fatty acid composition), and diel grazing patterns during autumn remain unknown. In this study, therefore, we investigated the population structure and ecological features (body size, feeding activity, and fatty acid composition) of C. glacialis to assess its current state in the changing Arctic Ocean.
2 Materials and methods
2.1 Field survey and onboard analysis
Zooplankton sampling was conducted via the R/V Mirai in the Pacific Arctic Ocean (66.50–78.04˚N, 145.00–168.75˚W) from October 8 to 27, 2019. Zooplankton samples were collected using vertical tows and a ring net (mouth diameter: 80 cm, mesh size: 335 µm) and a quad-NORPAC (North Pacific) net (mouth diameter: 45 cm, mesh sizes: 150 µm and 63 µm; ) at 10 m above the bottom or 150 m to the surface at 40 stations (Figure 1; Supplementary Table 1). The volume of the water filtered through the NORPAC net with the 150-µm mesh was estimated using a one-way flow meter (Rigosha CO., Ltd., Bunkyo-ku, Tokyo, Japan) mounted in the mouth of the net for quantitative collection of all copepodite stages in C. glacialis. Temperature, salinity, dissolved oxygen, and turbidity were simultaneously measured using vertical casts from a conductivity temperature depth (CTD) sensor (Sea-Bird Electronics Inc., SBE 911plus).
Figure 1
NORPAC net samples (mesh size: 150 µm) were immediately preserved in 5% v/v borax-buffered formalin to analyze the abundance and population structure of C. glacialis. The remaining samples collected by the ring net and the NORPAC net (mesh size: 63 µm) were used to conduct gut pigment and fatty acid composition analyses. Gut pigment analyses were conducted in samples from 11 stations on the Chukchi Plateau (77°N) obtained during October 21–24, and fatty acid analyses were conducted in samples from throughout the study area.
With respect to the gut pigment analysis, 10% v/v soda (saturated with CO2 in water) was immediately added to the fresh samples after sample collection to avoid copepod grazing, gut evacuation, and the decomposition of gut pigments (cf. Matsuno et al., 2015). Fresh specimens of C. glacialis copepodite stage five (C5) were sorted using a stereomicroscope (M165C-Ergo/MCI-HH/1; Leica Inc., Wetzlar, Germany). All specimens were sorted under low-temperature and dim-light conditions within 1 h, and the prosome length (PL) was measured using an eyepiece micrometer attached to a stereomicroscope with an accuracy of 0.001 mm. The sorted specimens (two to five specimens, average: three) were margined in small (2.5–3.3 mm) or large (3.4–3.8 mm) categories and transferred into cuvette tubes filled with 6 mL of N, N-dimethylformamide. One to four duplicates were obtained at each station, resulting in a total of 31 and 28 gut piments samples for the small and large categories, respectively. All samples were kept under cool and dark conditions for at least 24 h to extract chlorophyll and phaeopigments. After extraction, chlorophyll and phaeopigments were analyzed using a fluorometer (10-AU-005; Turner Designs Inc., California, United States), and the results were combined and expressed as the digestive tract pigment content (µg pigment ind.-1; cf. ).
The remaining sorted fresh samples were desalted using Milli-Q (Millipore Inc.), stored in 3-mL ointment bottles, and frozen at -80°C until the fatty acid composition was analyzed.
2.2 Sample analysis
In the laboratory, 19 zooplankton samples (Supplementary Table 1) collected by the NORPAC net with a 150-µm mesh size were split using a Motoda box splitter (Motoda, 1959). Calanus glacialis was identified at the copepodid stage level based on in the aliquots and counted under a dissecting microscope (SMZ1000, C-BD115, Nikon Corporation, Japan). To calculate the total species abundance (ind. m-3), the count at each stage was summarized and divided by the filtered volume (m3) of the net. Calanus marshallae is potentially distributed in the northern Bering Sea and Chukchi Sea (), but many zooplankton studies have treated this species as C. glacialis/marshallae owing to difficulties in their identification. Therefore, as mentioned in the Introduction, we refer to specimens from the C. glacialis/marshallae species complex as C. glacialis hereafter.
To determine geographical differences in the body size of C. glacialis, the PL of at least 10 C. glacialis C5 individuals at each station was measured using an eyepiece micrometer attached to a stereomicroscope with an accuracy of 0.001 mm. This measurement was performed for all zooplankton samples (Figure 1).
The length of lipid accumulation relative to the PL (lipid accumulation level) was measured at three levels: low, 0–4% of PL; medium, 4–40% of PL; and high, >40% of PL (). To compare C. glacialis data, the abundance (ind. m-2) was calculated by multiplying the population density (ind. m-3) with the sampling depth to eliminate the influence of sampling depth differences at each station. The mean C. glacialis copepodid stage (MCS) was calculated using Eq (1).
where i (1–6 indicates C1–C6) indicates the copepodid stage for C. glacialis, and Ai (ind. m-2) is the abundance of a copepodid stage (cf. ). Similarly, the mean lipid stage (MLS) of C. glacialis C5 was calculated using Eq (2).
where j (1, 2, and 3 indicates low, medium, and high, respectively) indicates the lipid accumulation level, and Bj (ind. m-2) is the ratio of the lipid accumulation level to the abundance.
Gonad maturation in adult females of the dominant copepod species was classified as stage I (immature), stage II (small oocytes in the ovary or oviduct), or stage III (large eggs or distended opacity in the oviduct) (cf. Miller et al., 1984; Niehoff, 1998).
2.3 Fatty acid analysis
Fatty acid analyses were performed using a simplified direct saponification/methylation method () after necessary modifications were made for wax ester-containing samples. Each C. glacialis sample was transferred to a pre-weighed screw-capped glass vial and dried in vacuo until a constant weight was achieved. Methyl tricosanoate and 1-tricosanol in toluene (2 mg mL-1 and 1 mg mL-1, respectively; 1 µL) and 1 M KOH in 95% ethanol (50 µL) were added to the vial, and after being crushed with a stainless-steel rod, the sample was directly saponified overnight in the dark at 23°C. The reaction was stopped by adding 3 M HCl in methanol (33 µL; Sigma-Aldrich), and the solvents and excess HCl were evaporated under a stream of nitrogen and then in vacuo. The fatty acids released into the vial were methylated in 3 M HCl in methanol (100 µL) at 80°C for 10 min, and the reagent was evaporated under a stream of nitrogen. The fatty acid methyl esters were purified via column chromatography using silica gel 60 (Merck) and hexane/diethyl ether (95:5, v/v), and fatty alcohols were recovered from the same column by elution with hexane/diethyl ether (50:50, v/v) and stored.
The purified fatty acid methyl esters were analyzed via gas chromatography using a GC-4000 gas chromatograph (GL Sciences, Tokyo, Japan) and an InertCap Pure-WAX column (30 m × 0.32 mm i.d., 0.25 µm film thickness; GL Sciences, Tokyo, Japan). The column temperature was programmed to increase from 90°C to 170°C at 20°C min-1, then to 240°C at 4°C min-1, and to be finally maintained at 240°C for 25 min. The injector and detector temperatures were set at 240°C. The carrier gas was helium at a linear velocity of 29.3 cm s-1 at 170°C. Methyl esters dissolved in hexane were injected into the column in the splitless injection mode, which was held for 1 min. A flame ionization detector (Shimadzu Corporation, Japan) and Shimadzu C-R3A integrator (Shimadzu Corporation, Japan) were used for detection and conducting peak area measurements. Peaks were identified by comparing their retention times with those of sardine oil fatty acids. The fatty acid content was calculated from the peak area ratios of methyl tricosanoate added to the initial C. glacialis samples.
2.4 Data analysis
To compare the biological parameters of C. glacialis C5 within the Pacific sector of the Arctic Ocean in 2019, a cluster analysis of environmental data was conducted. In this respect, water column (from the surface to 10 m above the bottom or 150 m above the bottom) averages of hydrographic data (temperature, salinity, dissolved oxygen, and turbidity) were calculated at each station. The water column data were then standardized, and similarities between stations were examined using the Euclidean distance. To group stations, similarity indices were coupled with hierarchical agglomerative clustering and the complete linkage method (unweighted pair group method with arithmetic mean, UPGMA) using the SIMPROF test. Based on the dendrogram, stations were separated to define regions with similar hydrographic conditions. A principal component analysis (PCA) was performed using standardized hydrographic data to evaluate hydrographic condition trends in each group. These analyses were conducted using Primer 7 software (PRIMER-E Ltd.).
Non-metric multi-dimensional scaling (NMDS) ordination was conducted using Primer 7 based on population structure, and Pearson’s regressions of standardized hydrographic data were then conducted to delineate the relationship between population structure and hydrography. The abundance of each copepodid stage (ind. m-2) was transformed into that of the 4th root stage (X-4). Similarities between samples were examined using the Bray-Curtis index according to differences in the stage composition.
To obtain the fatty acid composition, we calculated seven indices (docosahexaenoic acid (DHA)/eicosapentaenoic acid (EPA), 16 polyunsaturated fatty acids (PUFA)/18PUFA, 18:2n-6, 15:0 + 17:0, PUFA/saturated fatty acid [SFA], 18:1n-9/18:1n-7, and diatoms [D]/flagellates [F]) for the prey record of C. glacialis (; ; ; Stevens et al., 2004b; ).
Differences in biological parameters (abundance, MCS,PL, MLS, and fatty acid composition) of C. glacialis C5 among regions (identified by the cluster analysis based on hydrography) were tested using the max-t method with a heteroscedastic consistent covariance estimation (HC3) (). The tests were conducted using R software with the packages “multcomp” and “sandwich” (version 4.1.2, R Core Development Team, 2021). The relationship between the PL of C5 and temperature was regressed using the Bêlehrádek equation. A PCA based on normalized seven-index data was performed for fatty acids using Primer 7.
3 Results
3.1 Regions and water masses
The cluster analysis of biological parameters using hydrological data identified four regions in the Pacific sector of the Arctic Ocean (Figure 2A). As these regions were clearly divided geographically, they were named (from south to north) the southern shelf, northern shelf, slope, and basin (Figures 2A, C). The PCA results showed that temperature, salinity, and turbidity were higher on the southern shelf and lower in the basin (Figure 2B). Salinity was exceptionally high on the slope, and the dissolved oxygen level was high on the northern shelf and basin (Figure 2B).
Figure 2
defined water masses based on temperature and salinity as follows: warm Coastal Water (wCW), warm Shelf Water (wSW), Ice Melt Water, cool Coastal Water (IMW cCW), cool Shelf Water (cSW), Anadyr Water (AnW), Modified Winter Water (MWW), Winter Water (WW), and Atlantic water and bearing basin water (AtlW and BBW) (Figure 3A). Within the 20 m depth on the southern shelf, wSW was dominant on October 8 but cSW occupied this depth range on October 27 (Figure 3B; Supplementary Table 1). A comparison of the temperature and salinity on these two dates showed that salinity was similar (32–32.49 vs. 31.99–32.36) but temperature was lower on October 27 (decreased from 3.77–5.85°C to 2.36–2.84°C). In the northern shelf region, the water masses differed because the stations and currents were widely distributed. AnW was observed below a depth of 25 m at St. 54, and AtlW and BBW were found below depths of 100 m and 125 m, respectively, on the slope. In the basin, a similar vertical distribution was observed among the stations; the surface layer (0–40 m or 0–60 m) comprised IMW cCW, while cSW, MWW, or WW existed at depths below these water masses.
Figure 3
3.2 Geographical abundance and body size of C. glacialis
The abundance of C. glacialis ranged from 126 to 2,968 ind. m-2, with higher and lower abundances on the slope and shelf, respectively (Supplementary Figure 1). A comparison of the abundance and stage composition between the regions showed that the southern and northern shelves exhibited low abundance but had a high composition of C4 and C5 (Figure 4A). On the southern shelf, C5 was predominant in early October (St. 13 and 18), whereas adult females and males were predominant in late October (Sts. 61 and 70). The highest abundance of C3–C5 was observed on the slope. In the basin, the early stages (C1 or C2) occurred at 80% of stations (8/10), and the highest proportion of early stages was 44% (St. FP3; Figure 4). A few adult females were observed in each region, and almost all specimens showed stage I gonadal maturation (data not shown). NMDS based on the population structure showed relatively good separation with clustering by hydrography (Figure 4B). Among the hydrographic variables, temperature and turbidity had significant relationships with population structure. Of the basin group, Sts. 30, 34, and 35 were located slightly further from the main distribution of the basin community owing to the inflow of cSW at these stations (cf. Figure 3).
Figure 4
The body size (PL) of C. glacialis C5 ranged from 2.4–3.8 mm throughout the sampling region (Figure 5). The median PL value of each region showed an increasing trend from the Chukchi Sea toward the Arctic basin, in the order of the southern shelf, northern shelf, slope, and basin (Figure 5). Body size was well described by the Bêlehrádek equation (PL=23594 (T+76.17)-2.05; Figure 6).
Figure 5
Figure 6
3.3 Relationship between abundance, fatty acid index, gut pigment and geographical distribution
A comparison of the biological parameters of C. glacialis among regions showed significant differences in the abundance of C5, PL of C5, total fatty acid content per individual, and fatty acid index (16PUFA/18PUFA, PUFA/SFA, and D/F; Table 1). The abundance of C5 was significantly higher on the slope than on the other regions. (Figure 4; Table 1). Significant differences in the MCS were detected, and it decreased from south to north (Figure 4; Table 1). There was no significant difference in MLS among regions; however, a minimum value was observed on the southern shelf (Table 1; Supplementary Figure 2). The total fatty acid content per individual in the basin was significantly higher than that on the southern shelf (Table 1). With respect to the fatty acid index, 16PUFA/18PUFA, PUFA/SFA, and D/F were significantly lower on the northern shelf, but DHA/EPA, 18:2n-6, and 15:0 + 17:0 were higher, although not significantly (Figure 7; Table 1). The results of the PCA based on the fatty acid indices showed that the northern shelf was distant from the distribution of the other groups (Figure 8). As shown in the regional comparison using the max-t test, the northern shelf had high 15:0 + 17:0, DHA/EPA, and 18:2n-6 values (Figure 8), and the other groups had relatively high 16PUFA/18PUFA, D/F, and PUFA/SFA ratios. There was no significant difference between the three groups (excluding the northern shelf), and the variabilities within the southern shelf and slope were within the basin variability.
Table 1
| Parameters | Groups | Max-t test + HC3 | ||||||
|---|---|---|---|---|---|---|---|---|
| South_shelf | North_shelf | Slope | Basin | |||||
| Abundance of whole stages (ind. m-2) | 328 ± 128 | 266 ± 100 | 2412 ± 556 | 619 ± 302 | NS | |||
| Abundance of C5 (ind. m-2) | 294 ± 135 | 187 ± 99.0 | 590 ± 6.88 | 200 ± 127 | North_shelf2 | Basin2 | South_shelf2 | Slope1 |
| MCS (Mean copepodid stage) | 5.05 ± 0.05 | 4.68 ± 0.13 | 4.06 ± 0.04 | 3.67 ± 0.83 | Basin3 | Slope3 | North_shelf2 | South_shelf1 |
| Prosome length of C5 (mm) | 2.91 ± 0.21 | 3.11 ± 0.15 | 3.24 ± 0.26 | 3.34 ± 0.23 | South_shelf3 | North_shelf2 | Slope1 | Basin1 |
| MLS (Mean Lipid Stage) | 2.07 ± 0.53 | 2.54 ± 0.03 | 2.42 ± 0.08 | 2.48 ± 0.27 | NS | |||
| Fatty acid (µg ind.-1) | 138 ± 30.7 | 230 ± 97.3 | 233 ± 77.8 | 305 ± 149 | South_shelf2 | North_shelf1,2 | Slope1,2 | Basin1 |
| DHA/EPAa | 0.47 ± 0.12 | 0.92 ± 0.24 | 0.30 ± 0.02 | 0.33 ± 0.10 | NS | |||
| 16PUFA/18PUFAb | 0.65 ± 0.27 | 0.28 ± 0.02 | 0.89 ± 0.12 | 0.45 ± 0.23 | North_shelf2 | Basin1,2 | South_shelf1,2 | Slope1 |
| 18:2n-6c | 0.89 ± 0.14 | 2.30 ± 1.02 | 0.89 ± 0.20 | 0.92 ± 0.20 | NS | |||
| 15:0 + 17:0d | 1.50 ± 0.16 | 2.31 ± 0.35 | 1.66 ± 0.37 | 1.23 ± 0.20 | NS | |||
| PUFA/SFAe | 2.02 ± 0.29 | 0.92 ± 0.18 | 2.11 ± 0.53 | 1.92 ± 0.34 | North_shelf2 | Basin1 | South_shelf1 | Slope1 |
| 18:1n-9/18:1n-7f | 6.90 ± 1.82 | 6.66 ± 0.25 | 4.77 ± 1.63 | 6.02 ± 1.57 | NS | |||
| D/Fg | 1.98 ± 0.47 | 1.62 ± 0.23 | 3.06 ± 0.37 | 3.26 ± 0.95 | North_shelf2 | South_shelf2 | Slope1,2 | Basin1 |
Regional comparison of biological parameters of Calanus glacialis/marshallae in different areas of the Pacific sector of the Arctic Ocean during October 2019.
Different superscript numbers in the Max-t test + HC3 column indicate significant differences between regions. NS, not significant.
a: Dinoflagellates/diatoms, carnivory (); b: Diatoms/flagellates (; ); c: Terrestrial detritus or green algae (); d: Bacteria (); e: Carnivory (Stevens et al., 2004a); f: Carnivory or omnivory (Stevens et al., 2004b); g: Diatoms (16 polyunsaturated fatty acids (PUFA), 16:1n-7 + 16:1n-9, 20:5n-3)/flagellates (22:6n-3, 18PUFA, 18:2n-6; ).
Figure 7
Figure 8
In the basin, small (2.5–3.3 mm in PL) and large (3.4–3.8 mm in PL) forms occurred simultaneously (Figure 9). The cut off criteria between the forms is a length of 3.3 mm. Both small and large specimens were found in all 11 samples during October 21–24. The gut pigment showed a distinct diel pattern in the small forms, with high values during the night but low values during the daytime. However, the large forms did not exhibit a diel grazing activity pattern. The fatty acid indices were not significantly different between the small and large forms (data not shown).
Figure 9
4 Discussion
4.1 Geographical variation in body size
The body size of copepods varies with the environment (temperature and food concentration; ; Viitasalo et al., 1995; ). In the Chukchi Sea in 2004, 2009, and 2012, the body size of C. glacialis C5 was found to be negatively correlated with temperature (). Similar results were detected in our study for the Chukchi Sea and Arctic Basin.
In the pan-Arctic Ocean, the size range of the PL of species from each region has been well reported: 2.9–3.8 mm in Kandalaksha Bay and the White Sea (), 2.55–3.55 mm near the Svalbard Islands (), 2.8–4.0 mm in Disko Bay (Nielsen et al., 2014), 2.4–3.6 mm in Kongsfjorden (), and 2.4–3.7 mm in the Canadian Arctic and on the Atlantic coasts (Parent et al., 2011). Compared to those reports from other regions, the PL size range in our study (2.4–3.8 mm) was one of the largest, and this was because the border of the sampling region was situated in the Pacific sector of the Arctic Ocean. The body size of C5 in the Canada Basin in October 1997 was reported as 3.38 ± 0.03 mm (), which is similar to our results (basin: 3.34 ± 0.23 mm; slope: 3.24 ± 0.26 mm; Beaufort Sea: 3.46 ± 0.20 mm). However, smaller body sizes were observed on the southern and northern shelves in our study (2.91 ± 0.21 mm, and 3.11 ± 0.15 mm, respectively). It is considered that the co-occurrence of Calanus marshallae, which is slightly smaller than C. glacialis is (), may decrease the mean body size in both these regions.
Two C. glacialis haplotypes exist in the study area: Bering Sea and Arctic Basin populations (Nelson et al., 2009). The Bering Sea population extends from the Bering Strait to the northern Chukchi Sea owing to the inflow of Pacific Ocean water (Nelson et al., 2009; ). The inflow of Pacific water usually peaks from August to October (Woodgate, 2018) and is mainly composed of wSW (). As individuals in the Bering Sea population experience warmer conditions than do those in the Arctic Basin population, it is considered that the Bering Sea population would be smaller than would the Arctic Basin population be. However, a haplotype analysis could not be conducted in this study because ethanol samples were consumed for another scientific work, and the difference in the PL between the Bering Sea and Arctic Basin populations, therefore, remains unclear.
4.2 Reproduction and growth
Calanus glacialis is believed to undergo capital and income breeding depending on the food conditions (Plourde et al., 2005; ). Although our sampling was conducted in autumn, and our results do not, therefore, contain evidence of reproduction, the timing of reproduction in spring was back-calculated using the developmental time and embryonic duration () (see Supplementary Figure 3, Supplementary Table 3). This type of back-calculation helps to remove the sampling time gap among regions (). The estimation showed that almost all individuals had been spawned in July and August.
In the Pacific sector of the Arctic Ocean, peaks in phytoplankton blooms have been observed during May in the Chukchi Sea, and in June and July in the Siberian and Beaufort Seas (). In relation to this seasonality, early stages (C1 and C2) dominate within the southern Chukchi Sea (). However, there was a timing gap between the phytoplankton blooms (May–July) and our estimated reproduction period (July–August), which could have been caused by the shorter development time for the back-calculation of reproduction timing (Supplementary Table 3). The developmental time was assessed using the temperature from our sampling period (October); however, they survived in colder conditions (Spring and Summer). Since warmer temperature induce shorter development time (), therefore, the estimated development time was shorter than the in-situ value. Another hypothesis is that northward currents flushed out individuals during May–July, and individuals transported from the Bering Sea could not be established in the Chukchi Sea because of their short flushing time ().
Timing of phytoplankton blooms shifts from south to north in the Pacific Arctic Ocean (), subsequent reproduction of C. glacialis is expected to be associated with that. However, interestingly, it was suggested that the earliest reproduction period occurred in June in the basin region. This was potentially caused by the longer residence time of C5 than that of the other stages: the developmental time from nauplius three to copepodite five is almost similar for each molting, but the period from C5 to molt adults is longer (Peterson, 1986). Alternatively, it is possible that another previous generation survived because of the low metabolic rate at cold temperatures and the low grazing pressure by fish (Mundy et al., 2017; Skjoldal, 2022).
4.3 Southern shelf
In early October, the entire water column was occupied by wSW, but the surface layer (0–20 m) changed to cSW in late October, possibly due to cooling by low air temperature, as other shelf water masses (e.g., AnW and wCW) were not observed. Regarding the population structure, C5 was predominant in early October, whereas adults were predominant in late October. Based on these results, the population of C. glacialis appeared to grow, with some adults developing in late October. Incubation experiments revealed that 50% of C. glacialis individuals develop from C5 to adulthood in 20 days at 10°C (Peterson, 1986). In our development time estimation based on , it would take 10 days to undergo this process at 10°C and 20 days at 3°C (Supplementary Figure 3); therefore, it was not surprising that some adults were observed in October. A high abundance of adult females was observed in the Arctic Basin in February, but only a few such individuals were found throughout the year (). During its life cycle, C. glacialis is capable of diapause from C4 and C5 when it descends into the deep sea, and rich lipids are stored in its body during autumn and winter (; ; Pinchuk and Eisner, 2017). Because a sufficient depth (> 200 m) is required for the diapause of this species (; ), the shallow bottom depth (approximately 50 m; Nelson et al., 2009; ) and short residence time within the Chukchi Sea () imply that the population cannot be sustained/established.
The northward transportation of Pacific zooplankton has been reported in the Pacific Arctic Ocean during warm years (Matsuno et al., 2011; ). Calanus marshallae, which was originally distributed in the Northern California Current system to the coastal Gulf of Alaska and the Bering Sea (), should be included within the Pacific species; however, many studies have not distinguished it from C. glacialis owing to the morphological similarities between the species. Genetic approaches have shown that C. marshallae is present in the Chukchi Sea during the summer (Nelson et al., 2009 and ). In addition, the northward flow at the Bering Strait is increasing (0.01 Sv per year), and residence time in the Chukchi Sea has decreased by ~1.5 months (from 7.5 months to ~ 5 months) based on mooring results from 1990 to 2019 (Woodgate and Peralta-Ferriz, 2021), which subsequentially extend C. marshallae distribution. Therefore, the inflow of C. marshallae potentially distorts any changes identified in the C. glacialis population, especially in the shelf region during warm years. Genetic analyses are therefore required to avoid this distortion.
4.4 Northern shelf
The population structure of this area mainly comprises C4 and C5 C. glacialis. Assuming that this population on the shelf originates from water masses similar to those on the southern shelf (wSW and cSW) that occupy more than half of the water column, the population on the northern shelf (i.e., C4 and C5) is younger than that on the southern shelf (C5). This may be attributed to the shorter intermolt period at warmer temperatures (; Supplementary Table 3). The effect of temperature on copepod growth has also been observed during warm and cold years in the Chukchi Sea (Matsuno et al., 2011).
High variability in hydrography and water masses was observed in this study, and the fatty acid composition of the specimens collected on the northern shelf was lower in terms of 16 PUFA/18PUFA, PUFA/SFA, D, and F than that of specimens from other areas. The fatty acid composition is a valuable index for prey records because the composition of predators reflects that of their prey (). This result implies that individuals on the northern shelf do not actively utilize diatoms or protozoans (Stevens et al., 2004a; Stevens et al., 2004b). In contrast, higher values of DHA/EPA, 18:2n-6, and 15:0 + 17:0 were observed in specimens on the northern shelf than in those from other areas. Docosahexaenoic acid is abundant in dinoflagellates, and EPA is abundant in diatoms (), and the high DHA/EPA ratio, therefore, suggests that this species grazes mostly on dinoflagellates within this area. However, this feature is considered to be rare in the Arctic Ocean compared with that in the well-studied Atlantic Arctic Ocean. A survey in the southern shelf area west of Spitsbergen in July 2018 showed a DHA/EPA of 0.67 (Trudnowska et al., 2020), which is lower than that observed in our study (0.92). In addition, a DHA/EPA value of 0.21 was observed in Kongsfjord in the Atlantic Arctic Ocean during autumn in 1997 (), 0.15 in the Fram Strait during autumn in 2003 (Søreide et al., 2008), and 0.22 in Northwater during the autumn of 1999 (Stevens et al., 2004a). It is evident that this species can change its diet depending on the time and location (), and the results here show that the species on the northern shelf in the study area are highly dependent on dinoflagellates.
Considering the energy flow in lower ecosystems, eutrophic conditions induce the dominance of diatoms because their growth is rapid from ingesting rich nutrients; however, oligotrophic conditions stimulate microbial loop production (Pomeroy, 1974). In the oligotrophic state, bacteria act as primary producers using dissolved organic carbon, and several predator-prey interactions (bacteria vs. heterotrophic nanoflagellates, heterotrophic nanoflagellates vs. small ciliates) occur, reaching their production into heterotrophic dinoflagellates (; ). Therefore, the number of transfers between trophic levels is much greater under oligotrophic conditions, resulting in a low transfer efficiency. This suggests that the Calanus glacialis can survive in a low-productivity ecosystem. In addition, although a high 18:2n-6 value indicates that the species grazes on green algae of a terrestrial origin (), there have been no reports on the feeding on green algae by this species in the ocean; therefore, further research is required.
4.5 Slope
The highest abundance was observed on the slope, and these results are consistent with those obtained during studies conducted in August (; ) and September (). In this area, easterly winds stimulate upwelling along the shelf-break (Pickart et al., 2013), resulting in the occasional transport of diapause individuals (mainly C4 and C5) from the deep sea to the surface (; ). AtlW and BBW were only observed in this region below depths of 100 m or 125 m, which suggests the presence of upwelling (). A higher abundance was therefore observed on the slope because of the presence of individuals derived from the deep sea via upwelling (Pickart et al., 2013; Rutzen and Hopcroft, 2018). Another mechanism that increases abundance is anti-cyclonic, cold-core eddy formation from the boundary current flowing along the Chukchi Shelf (Pickart et al., 2005; ). Although the accumulation of zooplankton by eddies is not very strong (), the formation rate of eddies is very high, with up to 1–2 eddies per day from early spring to early autumn (Pickart et al., 2005). As previously mentioned, although shelf populations are flushed out into the basin (), a high abundance is frequently observed in this region, and this could potentially be related to the increased eddy formation rate.
The population structure derived in our study is considered to reflect population growth, and this has also been determined in other studies. In August, early stages (C1 and C2) are abundant (), C1–C4 are dominant in September (), and C3–C5 are dominant in October (this study). The abundance alters between months, and it decreases from August (mean 18,910 ind. m-2, ) and September (mean 8,000 ind. m-2, ) to October (1,855-2,968 ind. m-2). It is of note that all the abundance data were based on the samples collected using the 150 µm mesh. Assuming the mortality rate for broadcast spawner based on field measurements (), a value of 0.0335 d-1 was calculated at 0.27°C (Supplementary Table 3). Using this rate, the abundance in October was expected to be 2,448 ind. m-2 (=18,910 ind. m-2×(1-0.035)^(60 days)) from August and 2,878 ind. m-2 (=8,000 ind. m-2×(1-0.035)^(30 days)) from September. Predation mortality is independent of ambient temperature, and the ratio of predation mortality to total mortality increases to 3:4 in colder situations (). This indicates that the decrease in abundance on the slope was mainly caused by predator feeding. A large biomass of polar cod has been observed in Atlantic Water along the slope (). Therefore, the high biomass of C. glacialis is due to physical processes, and the rich prey is grazed on by organisms in the higher trophic levels, such as fish, birds, and marine mammals (; ; ).
With respect to the fatty acid composition, this region had the lowest DHA/EPA ratio and the highest 16PUFA/18PUFA ratio among the four regions, which indicates that individuals in this area consume more diatoms than do dinoflagellates and flagellates. This area is prone to upwelling (Pickart et al., 2013), which supplies rich nutrients from deeper layers to the euphotic layer (Pickart et al., 2013; Rutzen and Hopcroft, 2018). Consequently, silicate concentrations in the northern Chukchi Sea increase, providing a suitable environment for diatom growth (; Nishino et al., 2015). Therefore, geographical and physical features (i.e., upwelling) affect the C. glacialis population, and lipid accumulation relates to diatom grazing on the slope.
4.6 Basin
Based on the environmental data, Sts. 30, 34, and 35 were located in the basin, and the population structures at these stations differed considerably from those of the other stations on the Chukchi Plateau. This may be due to the transportation of shelf individuals by cSW inflow to depths of 50–100 m (; ). Previous studies found that the C5 stage dominates the basin population in autumn (; ). However, in our study, C5 was not abundant. The early copepodite stages (C1 or C2) occurred at 80% of the basin stations, and the maximum population percentage composition was 44%. Similar results, such as the occurrence of early stages on the surface of the basin, have been reported in late summer () and winter (). However, the composition of the early stage was lower than 10% because this species is capable of diapause after C4 (; ), and it is likely that C1 and C2 cannot overwinter in the basin because of their short growth periods ().
With respect to the population origin of the early stages, adult females were observed at some stations in the basin, but all specimens were immature. However, according to the developmental time (Supplementary Table 3), the early stages were spawned in August, and C5 was spawned in late June. This shows that long-term reproduction is supported by mixed reproductive strategies (Plourde et al., 2005; ). Primary production on the Chukchi Plateau is very low (), but the region is greatly affected by the inflow of Pacific Summer Water (PSW), which transports heat and organisms via eddies (Watanabe et al., 2014; Muramatsu et al., 2021). The composition of PSW in the basin has increased in recent years compared to that during the 1990s and the 2000s (; Muramatsu et al., 2021). Therefore, more organisms are transported by PSW with eddies (Watanabe et al., 2014), and this extends the reproduction period of the species. Potentially due to the extended reproduction period, the early stages eventually become more abundant than those in the results of and . Field measurements of longevity have shown that adults are quickly removed 16 days after reproduction (). These assumptions suggest that adult females reproduce over the long term using transported organisms (e.g., phytoplankton), then that are removed after reproduction. Similarly, overwintering individuals reproduce in the Chukchi Sea following sea ice retreat, and the early stages are transported to the basin (). However, no genetic evidence is available to support this hypothesis, and a haplotype analysis of the early stages is required to reveal the extended distribution and establishment of the Bering Sea population in the basin (; Skjoldal, 2022),
Small and large forms co-occur on the Chukchi Plateau. On the Atlantic side of the Arctic Ocean, it is difficult to distinguish between C. glacialis and the smaller Calanus finmarchicus (). investigated the possibility of identifying species based on PL, morphological characteristics (the 5th pair of swimming legs and mandible), antennules, genital somite pigmentation in live specimens, and genetic analyses. The represented PL helps to approximate the species composition, but no criteria for the species have been established; therefore, their size ranges overlap (). In our case, the small forms (2.5–3.3 mm) in C5 could be included as C. finmarchicus (), but this result was impossible to clarify because of a lack of genetic analysis.
An inter-form difference in the diel patterns of grazing activity was observed. Small forms show higher gut pigment levels at night than during the day, which is associated with diel vertical migration (DVM) (). The magnitude of DVM for Calanus spp. varies with season and copepodid stage, and the DVM intensity is greater during spring and autumn, when the diel changes in light penetration are large (). During spring, C. glacialis responds sensibly to phytoplankton dynamics via DVM, reaching sea ice and exhibiting diel patterns in grazing activity (Runge and Ingram, 1991). In September, in the northern Chukchi Sea, diel changes in gut pigments have been reported during small autumn phytoplankton blooms (Matsuno et al., 2015). In contrast, DVM and active grazing ceases during the night for Calanus spp., as determined by the lipids accumulated in C5 (; Matsuno et al., 2015). Based on this knowledge, it is hypothesized that the small forms perform DVM and graze at night, while the large forms cease DVM and active grazing and stay in the deep layer. Although our sampling design could not corroborate this hypothesis, both forms were distributed within depths of 0–150 m due to their co-occurrence during the day and night ().
5 Conclusions
This study describes geographic variations in the population structure, body size, fatty acid composition, and diel changes in the grazing activity of C. glacialis/marshallae in the Pacific region of the Arctic Ocean during autumn. The body size of this species was found to be negatively correlated with the water temperature from the Chukchi Sea to the Arctic basin. The population developed from C5 was dominant on the southern shelf, and it was mostly transported from the Bering Sea. C4 and C5 were dominant on the northern shelf, suggesting that their growth was slower than the growth of those on the southern shelf. The fatty acid composition revealed that the rate of prey dependence on dinoflagellates was exceptionally high in the Arctic Ocean, indicating that the species existed in a low-productivity ecosystem. The highest abundance and composition of diatom-derived EPA among the fatty acids were found on the slope. This could be explained by the upwelling accompanying the inflow of this species from the deep sea, and the nutrient supply for the growth of diatoms. In the basin, early-stage individuals were more abundant than those reported in previous studies because a larger number of organisms had flowed into the region, resulting in extended reproduction. Small and large C5 forms were found, showing the different diel patterns of grazing activity. It appears that C. glacialis/marshallae have adapted to the flexibility in their life cycle caused by changes in the physical and chemical environments of the Pacific Arctic Ocean. However, vertically stratified sampling and haplotype analyses are required to monitor their biological features and reveal population origins.
Statements
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
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
KM and AY designed the study. KT, KM, and KS performed field sampling. MI, KT, and YA analyzed the samples. MI and KM analyzed the data and MI and KM wrote the paper with contributions from all authors.
Funding
This work was funded by the Arctic Challenge for Sustainability (ArCS; Program Grant Number JPMXD1300000000) and Arctic Challenge for Sustainability II (ArCS II; program grant number JPMXD1420318865) projects. This study was also supported by the Environment Research and Technology Development Fund (JPMEERF20214002) of the Environmental Restoration and Conservation Agency of Japan and partially supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (grant numbers JP22H00374 (A), JP21H02263 (B), JP20K20573 (Pioneering), JP20H03054 (B), JP20J20410 (B), JP19H03037 (B), JP18K14506 (Early Career Scientists), and JPJSBP120238801 (Bilateral program)).
Acknowledgments
We thank the captain, officers, crew, and researchers onboard the R/V Mirai, the Canadian Coast Guard RV Amundsen, and IB S.W. Laurier for their outstanding efforts during field sampling.
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: https://www.frontiersin.org/articles/10.3389/fmars.2023.1168015/full#supplementary-material
References
1
AlfaroA. C.ThomasF.SergentL.DuxburyM. (2006). Identification of trophic interactions within an estuarine food web (northern New Zealand) using fatty acid biomarkers and stable isotopes. Estuar. Coast. Shelf. Sci.70, 271–286. doi: 10.1016/j.ecss.2006.06.017
2
ArrigoK. R.van DijkenG. L. (2011). Secular trends in Arctic Ocean net primary production. J. Geophys. Res.116, C09011. doi: 10.1029/2011JC007151
3
AshjianC. J.CampbellR. G.GelfmanC.AlataloP.ElliottS. M. (2017). Mesozooplankton abundance and distribution in association with hydrography in the Hanna Shoal on Hanna Shoal, NE Chukchi Sea, during August 2012 and 2013. Deep. Sea. Res. II.144, 21–36. doi: 10.1016/j.dsr2.2017.08.012
4
AshjianC. J.CampbellR. G.WelchH. E.ButlerM.Van KeurenD. (2003). Annual cycle in abundance, distribution, and size in relation to hydrography of important copepod species in the western Arctic Ocean. Deep. Sea. Res. I.50, 1235–1261. doi: 10.1016/S0967-0637(03)00129-8
5
AshjianC. J.PickartR. S.CampbellR. G.FengZ.GelfmanC.AlataloP.et al. (2021). Springtime renewal of zooplankton populations in the Chukchi Sea. Prog. Oceanogr.197, 102635. doi: 10.1016/j.pocean.2021.102635
6
BanasN. S.MøllerE. F.NielsenT. G.EisnerL. B. (2016). Copepod life strategy and population viability in response to prey timing and temperature: Testing a new model across latitude, time, and size spectrum. Front. Mar. Sci.3. doi: 10.3389/fmars.2016.00225
7
BengtsonJ. L.Hiruki-RaringL. M.SimpkinsM. A.BovengP. L. (2005). Ringed and bearded seal densities in the eastern Chukchi Sea (1999–2000. Polar. Biol.28, 833–845. doi: 10.1007/s00300-005-0009-1
8
BourgainP.GascardJ. C. (2012). The Atlantic and summer Pacific waters variability in the Arctic Ocean from 1997 to 2008. Geophys. Res. Lett.39, L05603. doi: 10.1029/2012GL051045
9
BrodskyK. A. (1967). Calanoida of the far eastern seas and polar basin of the USSR. Keys. to Fauna. USSR.35, 1–442.
10
BudgeS. M.ParrishC. C. (1998). Lipid biogeochemistry of plankton, settling matter and sediments in Trinity Bay, Newfoundland. II. Fatty acids. Org. Geochem.29, 1547–1559. doi: 10.1016/S0146-6380(98)00177-6
11
CampbellR. G.SherrE. B.AshjianC. J.PlourdeS.SherrB. F.HillV.et al. (2009). Mesozooplankton prey preferences and grazing impact in the western Arctic Ocean. Deep. Sea. Res. II.56, 1274–1289. doi: 10.1016/j.dsr2.2008.10.027
12
CarmackE. C.KulikovE. A. (1998). Wind-forced upwelling and internal Kelvin wave generation in Mackenzie Canyon, Beaufort Sea. J. Geophys. Res.103, 18447–18458. doi: 10.1029/98JC00113
13
ChoquetM.KosobokovaK. N.KwaśniewskiS.HatlebakkM.DhanasiriA. K. S.MelleW.et al. (2018). Can morphology reliably distinguish between the copepods Calanus finmarchicus and C. glacialis, or is DNA the only way? Limnol. Oceanogr. Methods16, 237–252. doi: 10.1002/lom3.10240
14
CodispotiL. A.KellyV.ThessenA.MatraiP.SuttlesS.HillV.et al. (2013). Synthesis of primary production in the Arctic Ocean: III. Nitrate and phosphate based estimates of the net community production. Prog. Oceanogr.110, 126–150. doi: 10.1016/j.pocean.2012.11.006
15
ConoverR. J. (1988). Comparative life histories in the genera Calanus and Neocalanus in high latitudes of the northern hemisphere. Hydrobiologia167–168, 127–142. doi: 10.1007/BF00026299
16
ConoverR. J.HuntleyM. (1991). Copepods in ice-covered seas—Distribution, adaptations to seasonally limited food, metabolism, growth patterns and life cycle strategies in polar seas. J. Mar. Syst.2, 1–41. doi: 10.1016/0924-7963(91)90011-I
17
CorkettC. J.McLarenI. A. (1978). The biology of Pseudocalanus. Adv. Mar. Biol.15, 1–231. doi: 10.1016/S0065-2881(08)60404-6
18
CorkettC. J.McLarenI. A.SevignyJ.-M. (1986). Rearing of the marine Calanoid Copepods Calanus finmarchicus (Gunnerus), C. glacialis Jaschnov, and C. hyperboreus Kroyer with comments on the equiproportional rule. Syllogeus58, 539–546.
19
CrawfordR. E.VagleS.CarmackE. C. (2012). Water mass and bathymetric characteristics of polar cod habitat along the continental shelf and slope of the Beaufort and Chukchi seas. Polar. Biol.35, 179–190. doi: 10.1007/s00300-011-1051-9
20
DaaseM.Falk-PetersenS.VarpeØ.DarnisG.SøreideJ. E.WoldA.et al. (2013). Timing of reproductive events in the marine copepod Calanus glacialis: A pan-Arctic perspective. Can. J. Fish. Aquat. Sci.70, 871–884. doi: 10.1139/cjfas-2012-0401
21
DalsgaardJ.St JohnM.KattnerG.Müller-NavarraD.HagenW. (2003). Fatty acid trophic markers in the pelagic marine environment. Adv. Mar. Biol.46, 225–340. doi: 10.1016/s0065-2881(03)46005-7
22
DanielsonS. L.AhkingaO.AshjianC.BasyukE.CooperL. W.EisnerL.et al. (2020). Manifestation and consequences of warming and altered heat fluxes over the Bering and Chukchi Sea continental shelves. Deep. Sea. Res. II.177, 104781. doi: 10.1016/j.dsr2.2020.104781
23
DeeveyG. B. (1960). Relative effects of temperature and food on seasonal variations in the lengths of marine copepods in eastern American and Western European waters. Bull. Bing. Oceanogr. Coll.17, 54–86.
24
DicksonD. L.GilchristH. G. (2002). Status of marine birds of the southeastern Beaufort Sea. Arctic55, 46–58. doi: 10.14430/arctic734
25
DigbyP. S. B. (1954). Biology of marine planktonic copepods of Scores by Sound, east Greenland. J. Anim. Ecol.23, 298–338. doi: 10.2307/1984
26
El-SabaawiR.DowerJ. F.KainzM.MazumderA. (2009). Characterizing dietary variability and trophic positions of coastal calanoid copepods: Insight from stable isotopes and fatty acids. Mar. Biol.156, 225–237. doi: 10.1007/s00227-008-1073-1
27
ErshovaE. A.HopcroftR. R.KosobokovaK. N. (2015). Inter-annual variability of summer mesozooplankton communities of the western Chukchi Sea: 2004–2012. Polar. Biol.38, 1461–1481. doi: 10.1007/s00300-015-1709-9
28
ErshovaE. A.KosobokovaK. N.BanasN. S.EllingsenI.NiehoffB.HildebrandtN.et al. (2021). Sea ice decline drives biogeographical shifts of key Calanus species in the central Arctic Ocean. Glob. Change Biol.27, 2128–2143. doi: 10.1111/gcb.15562
29
FalkenhaugT.TandeK. S.SemenovaT. (1997). Diel, seasonal and ontogenetic variations in the vertical distribution of four marine copepods. Mar. Ecol. Prog. Ser.149, 105–119. doi: 10.3354/meps149105
30
Falk-PetersenS.LeuE.BergeJ.KwasniewskiS.NygårdH.RøstadA.et al. (2008). Vertical migration in high Arctic waters during autumn 2004. Deep. Sea. Res. II.55, 2275–2284. doi: 10.1016/j.dsr2.2008.05.010
31
Falk-PetersenS.MayzaudP.KattnerG.SargentJ. R. (2009). Lipids and life strategy of Arctic Calanus. Mar. Biol. Res.5, 18–39. doi: 10.1080/17451000802512267
32
FreeseD.SøreideJ. E.NiehoffB. (2016). A year-round study on digestive enzymes in the Arctic copepod Calanus glacialis: Implications for its capability to adjust to changing environmental conditions. Polar. Biol.39, 2241–2252. doi: 10.1007/s00300-016-1891-4
33
FrostB. W. (1974). Calanus marshallae, a new species of calanoid copepod closely aligned to the sibling species C. finmarchicus and C. glacialis. Mar. Biol.26, 77–99. doi: 10.1007/BF00389089
34
GabrielsenT. M.MerkelB.SøreideJ. E.Johansson-KarlssonE.BaileyA.VogedesD.et al. (2012). Potential misidentifications of two climate indicator species of the marine arctic ecosystem: Calanus glacialis and C. finmarchicus. Polar. Biol.35, 1621–1628. doi: 10.1007/s00300-012-1202-7
35
GraingerE. H. (1961). The copepods Calanus glacial is Jaschnov and Calanus finmarchicus (Gunnerus) in Canadian Arctic-Subarctic Waters. J. Fish. Res. Board. Can.18, 663–678. doi: 10.1139/f61-051
36
GrebmeierJ. M.CooperL. W.FederH. M.SirenkoB. I. (2006). Ecosystem dynamics of the pacific-influenced northern bering and chukchi seas in the amerasian arctic. Prog. Oceanogr.71, 331–361. doi: 10.1016/j.pocean.2006.10.001
37
HamaN.AbeY.MatsunoK.YamaguchiA. (2019). Effect of net mesh size on filtering efficiency and zooplankton sampling efficiency using quad-NORPAC net. Bull. Fish. Sci. Hokkaido. Univ.69, 47–56. doi: 10.14943/bull.fish.69.1.47
38
HerberichE.SikorskiJ.HothornT. (2010). A Robust procedure for comparing multiple means under heteroscedasticity in unbalanced designs. PloS One5, e9788. doi: 10.1371/journal.pone.0009788
39
HillV.CotaG. (2005). Spatial patterns of primary production on the shelf, slope and basin of the western Arctic in 2002. Deep. Sea. Res. II.52, 3344–3354. doi: 10.1016/j.dsr2.2005.10.001
40
HircheH.-J. (1991). Distribution of dominant calanoid copepod species in the Greenland Sea during late fall. Polar. Biol.11, 351–362. doi: 10.1007/BF00239687
41
HirstA. G.KiørboeT. (2002). Mortality of marine planktonic copepods: Global rates and patterns. Mar. Ecol. Prog. Ser.230, 195–209. doi: 10.3354/meps230195
42
HopcroftR. R.ClarkeC.NelsonR. J.RaskoffK. A. (2005). Zooplankton communities of the Arctic Canada Basin: The contribution by smaller taxa. Polar. Biol.28, 198–206. doi: 10.1007/s00300-004-0680-7
43
JakobsenH. H.HansenP. J. (1997). Prey size selection, grazing and growth response of the small heterotrophic dinoflagellate Gymnodinium sp. and ciliate Balanion comatum: A comparative study. Mar. Ecol. Prog. Ser.158, 75–86. doi: 10.3354/meps158075
44
JiR.AshjianC. J.CampbellR. G.ChenC.GaoG.DavisC. S.et al. (2012). Life history and biogeography of Calanus copepods in the Arctic Ocean: An individual-based modeling study. Prog. Oceanogr.96, 40–56. doi: 10.1016/j.pocean.2011.10.001
45
KanedaT. (1991). Iso- and anteiso-fatty acids in bacteria: Biosynthesis, function, and taxonomic significance. Microbiol. Rev.55, 288–302. doi: 10.1128/mr.55.2.288-302.1991
46
KimuraF.MatsunoK.AbeY.YamaguchiA. (2022). Effects of early sea-ice reduction on zooplankton and copepod population structures in the Northern Bering Sea during the summers of 2017 and 2018. Front. Mar. Sci.9. doi: 10.3389/fmars.2022.808910
47
KobariT.TadokoroK.ShiomotoA.HashimotoS. (2003). Geographical variation in prosome length and body weight of Neocalanus copepods in the North Pacific. J. Oceanogr.59, 3–10. doi: 10.1023/A:1022895802468
48
KosobokovaK. N. (1999). The reproductive cycle and life history of the Arctic copepod Calanus glacialis in the White Sea. Polar. Biol.22, 254–263. doi: 10.1007/s003000050418
49
KosobokovaK. N.HircheH.-J. (2009). Biomass of zooplankton in the eastern Arctic Ocean: A baseline study. Prog. Oceanogr.82, 265–280. doi: 10.1016/j.pocean.2009.07.006
50
KosobokovaK. N.HopcroftR. R. (2010). Diversity and vertical distribution of mesozooplankton in the Arctic’s Canada Basin. Deep. Sea. Res. II.57, 96–110. doi: 10.1016/j.dsr2.2009.08.009
51
LaneP. V. Z.LlinásL.SmithS. L.PilzD. (2008). Zooplankton distribution in the western Arctic during summer 2002: Hydrographic habitats and implications for food chain dynamics. J. Mar. Syst.70, 97–133. doi: 10.1016/j.jmarsys.2007.04.001
52
LevinsenH.NielsenT. G. (2002). The trophic role of marine pelagic ciliates and heterotrophic dinoflagellates in arctic and temperate coastal ecosystems: A cross-latitude comparison. Limnol. Oceanogr.47, 427–439. doi: 10.4319/lo.2002.47.2.0427
53
LlinásL.PickartR. S.MathisJ. T.SmithS. L. (2009). Zooplankton inside an Arctic Ocean cold-core eddy: Probable origin and fate. Deep. Sea. Res. II.56, 1290–1304. doi: 10.1016/j.dsr2.2008.10.020
54
MackasD.BohrerR. (1976). Fluorescence analysis of zooplankton gut contents and an investigation of diel feeding patterns. J. Exp. Mar. Biol. Ecol.25, 77–85. doi: 10.1016/0022-0981(76)90077-0
55
MarinV. (1987). The oceanographic structure of the eastern Scotia sea—IV. Distribution of copepod species in relation to hydrography in 1981. Deep. Sea. Res. Part A. Oceanographic. Res. Papers.34, 105–121. doi: 10.1016/0198-0149(87)90125-7
56
MatsumotoY.AndoY.HiraokaY.TawaA.OhshimoS. (2018). A simplified gas chromatographic fatty-acid analysis by the direct saponification/methylation procedure and its application on wild tuna larvae. Lipids53, 919–929. doi: 10.1002/lipd.12098
57
MatsunoK.AbeY.YamaguchiA.KikuchiT. (2016). Regional patterns and controlling factors on summer population structure of Calanus glacialis in the western Arctic Ocean. Polar. Sci.10, 503–510. doi: 10.1016/j.polar.2016.09.001
58
MatsunoK.YamaguchiA.HirawakeT.ImaiI. (2011). Year-to-year changes of the mesozooplankton community in the Chukchi Sea during summers of 1991, 1992 and 2007, 2008. Polar. Biol.34, 1349–1360. doi: 10.1007/s00300-011-0988-z
59
MatsunoK.YamaguchiA.NishinoS.InoueJ.KikuchiT. (2015). Short-term changes in the mesozooplankton community and copepod gut pigment in the Chukchi Sea in autumn: Reflections of a strong wind event. Biogeosciences12, 4005–4015. doi: 10.5194/bg-12-4005-2015
60
McLarenI. A.SevignyJ.-M.CorkettJ. (1988). Body size, development rate, and genome size of Calanus species in Dev. Hydrobiol., eds G. A. Boxshall and H. K. Schminke. Biol. Copepods.47, 275–284. doi: 10.1007/978-94-009-3103-9_27
61
MillerC. B.FrostB. W.BatchelderH. P.ClemonsM. J.ConwayR. E. (1984). Life histories of large, grazing copepods in a subarctic ocean gyre Neocalanus plumchrus, Neocalanus cristatus, and Eucalanus bungii in the northeast Pacific. Prog. Oceanogr.13, 201–243. doi: 10.1016/0079-6611(84)90009-0
62
MotodaS. (1959). Device for a simple plankton apparatus. Mem. Fac. Fish. (Hokkaido. Univ.)7, 73–94.
63
MundyP. R.IngvaldsenR.SunnanåK. (2017) Synthesis of knowledge on fisheries science in the central Arctic Ocean and adjacent waters (ToR1) in Final report of the Fourth Meeting of Scientific Experts on Fish Stocks in the Central Arctic Ocean, FiSCAO, (NOAA). 34–79.
64
MuramatsuM.UenoH.WatanabeE.ItohM.OnoderaJ. (2021). Transport and heat loss of the Pacific Summer Water in the Arctic Chukchi Sea northern slope: Mooring data analysis. Polar. Sci.29, 100698. doi: 10.1016/j.polar.2021.100698
65
NelsonR. J.CarmackE. C.McLaughlinF. A.CooperG. A. (2009). Penetration of Pacific zooplankton into the western Arctic Ocean tracked with molecular population genetics. Mar. Ecol. Prog. Ser.381, 129–138. doi: 10.3354/meps07940
66
NiehoffB. (1998). The gonad morphology and maturation in Arctic Calanus species. J. Mar. Syst.15, 53–59. doi: 10.1016/S0924-7963(97)00048-1
67
NielsenT. G.KjellerupS.SmolinaI.HoarauG.LindequeP. (2014). Live discrimination of Calanus glacialis and C. finmarchicus females: Can we trust phenological differences? Mar. Biol.161, 1299–1306. doi: 10.1007/s00227-014-2419-5
68
NishinoS.KawaguchiY.InoueJ.HirawakeT.FujiwaraA.FutsukiR.et al. (2015). Nutrient supply and biological response to wind-induced mixing, inertial motion, internal waves, and currents in the northern Chukchi Sea. J. Geophys. Res. Oceans.120, 1975–1992. doi: 10.1002/2014JC010407
69
ParentG. J.PlourdeS.TurgeonJ. (2011). Overlapping size ranges of Calanus spp. off the Canadian Arctic and Atlantic Coasts: Impact on species’ abundances. J. Plankton. Res.33, 1654–1665. doi: 10.1093/plankt/fbr072
70
PetersonW. (1986). Development, growth, and survivorship of the copepod Calanus marshallae in the laboratory. Mar. Ecol. Prog. Ser.29, 61–72. doi: 10.3354/meps029061
71
PickartR. S.SchulzeL. M.MooreG. W. K.CharetteM. A.ArrigoK. R.van DijkenG.et al. (2013). Long-term trends of upwelling and their impacts on primary productivity in the Alaskan Beaufort Sea. Deep. Sea. Res. I.79, 106–121. doi: 10.1016/j.dsr.2013.05.003
72
PickartR. S.WeingartnerT. J.PrattL. J.ZimmermannS.TorresD. J. (2005). Flow of winter-transformed Pacific water into the western Arctic. Deep. Sea. Res. II.52, 3175–3198. doi: 10.1016/j.dsr2.2005.10.009
73
PinchukA. I.EisnerL. B. (2017). Spatial heterogeneity in zooplankton summer distribution in the eastern Chukchi Sea in 2012–2013 as a result of large-scale interactions of water masses. Deep. Sea. Res. II.135, 27–39. doi: 10.1016/j.dsr2.2016.11.003
74
PlourdeS.CampbellR. G.AshjianC. J.StockwellD. A. (2005). Seasonal and regional patterns in egg production of Calanus glacialis/marshallae in the Chukchi and Beaufort Seas during spring and summe. Deep. Sea. Res. II.52, 3411–3426. doi: 10.1016/j.dsr2.2005.10.013
75
PomeroyL. R. (1974). The ocean food web: A changing paradigm. BioScience24, 499–504. doi: 10.2307/1296885
76
R Core Team (2021). R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Available at: https://www.R-project.org/.
77
RungeJ. A.IngramR. G. (1991). Under-ice feeding and diel migration by the planktonic copepods Calanus glacialis and Pseudocalanus minutus in relation to the ice algal production cycle in southeastern Hudson Bay, Canada. Mar. Biol.108, 217–225. doi: 10.1007/BF01344336
78
RutzenI.HopcroftR. R. (2018). Abundance, biomass and community structure of epipelagic zooplankton in the Canada Basin. J. Plankton. Res.40, 486–499. doi: 10.1093/plankt/fby028
79
ScottC. L.KwasniewskiS.Falk-PetersenS.SargentJ. R. (2000). Lipids and life strategies of Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus in late autumn, Kongsfjorden, Svalbard. Polar. Biol.23, 510–516. doi: 10.1007/s003000000114
80
SkjoldalH. R. (Ed.) (2022). “Ecosystem assessment of the Central Arctic Ocean: Description of the ecosystem,” in ICES cooperative research reports. 355, 341. doi: 10.17895/ices.pub.20191787
81
SøreideJ. E.Falk-PetersenS.HegsethE. N.HopH.CarrollM. L.HobsonK. A.et al. (2008). Seasonal feeding strategies of Calanus in the high-Arctic Svalbard region. Deep. Sea. Res. II.55, 2225–2244. doi: 10.1016/j.dsr2.2008.05.024
82
StevensC. J.DeibelD.ParrishC. C. (2004a). Species-specific differences in lipid composition and omnivory indices in Arctic copepods collected in deep water during autumn (North water Polynya). Mar. Biol.144, 905–915. doi: 10.1007/s00227-003-1259-5
83
StevensC. J.DeibelD.ParrishC. C. (2004b). Incorporation of bacterial fatty acids and changes in a wax ester-based omnivory index during a long-term inCubation experiment with Calanus glacialis Jaschnov. J. Exp. Mar. Biol. Ecol.303, 135–156. doi: 10.1016/j.jembe.2003.11.008
84
TourangeauS.RungeJ. A. (1991). Reproduction of Calanus glacialis under ice in spring in southeastern Hudson Bay, Canada. Mar. Biol.108, 227–233. doi: 10.1007/BF01344337
85
TrudnowskaE.BalazyK.Stoń-EgiertJ.SmolinaI.BrownT.GluchowskaM. (2020). In a comfort zone and beyond—Ecological plasticity of key marine mediators. Ecol. Evol.10, 14067–14081. doi: 10.1002/ece3.6997
86
VarpeØ.JørgensenC.TarlingG. A.FiksenØ. (2009). The adaptive value of energy storage and capital breeding in seasonal environments. Oikos118, 363–370. doi: 10.1111/j.1600-0706.2008.17036.x
87
VidalJ. (1980). Physioecology of zooplankton. II. Effects of phytoplankton concentration, temperature, and body size on the development and molting rates of Calanus pacificus and Pseudocalanus sp. Mar. Biol.56, 135–146. doi: 10.1007/BF00397130
88
VidalJ.SmithS. L. (1986). Biomass, growth, and development of populations of herbivorous zooplankton in the southeastern Bering Sea during spring. Deep. Sea. Res. Part A. Oceanographic. Res. Papers.33, 523–556. doi: 10.1016/0198-0149(86)90129-9
89
ViitasaloM.KoskiM.PellikkaK.JohanssonS. (1995). Seasonal and long-term variations in the body size of planktonic copepods in the northern Baltic Sea. Mar. Biol.123, 241–250. doi: 10.1007/BF00353615
90
WassmannP.KosobokovaK. N.SlagstadD.DrinkwaterK. F.HopcroftR. R.MooreS. E.et al. (2015). The Contiguous domains of Arctic Ocean advection: Trails of life and death. Prog. Oceanogr.139, 42–65. doi: 10.1016/j.pocean.2015.06.011
91
WatanabeE.OnoderaJ.HaradaN.HondaM. C.KimotoK.KikuchiT.et al. (2014). Enhanced role of eddies in the Arctic marine biological pump. Nat. Commun.5, 3950. doi: 10.1038/ncomms4950
92
WoodgateR. A. (2018). Increases in the Pacific inflow to the Arctic from 1990 to 2015, and insights into seasonal trends and driving mechanisms from year-round Bering Strait mooring data. Prog. Oceanogr.160, 124–154. doi: 10.1016/j.pocean.2017.12.007
93
WoodgateR.Peralta-FerrizC. (2021). Warming and freshening of the Pacific inflow to the Arctic from 1990–2019 implying dramatic shoaling in Pacific winter water ventilation of the Arctic water column. Geophys. Res. Lett.48, e2021GL092528. doi: 10.1029/2021GL092528
94
YehH. D.QuestelJ. M.MaasK. R.BucklinA. (2020). Metabarcoding analysis of regional variation in the gut contents of the copepod Calanus finmarchicus in the North Atlantic Ocean. Deep. Sea. Res. II.180, 104738. doi: 10.1016/j.dsr2.2020.104738
Summary
Keywords
Calanus glacialis/marshallae, Pacific Arctic Ocean, population structure, body size, fatty acid, gut pigment
Citation
Ishihara M, Matsuno K, Tokuhiro K, Ando Y, Sato K and Yamaguchi A (2023) Geographic variation in population structure and grazing features of Calanus glacialis/marshallae in the Pacific Arctic Ocean. Front. Mar. Sci. 10:1168015. doi: 10.3389/fmars.2023.1168015
Received
17 February 2023
Accepted
23 August 2023
Published
20 September 2023
Volume
10 - 2023
Edited by
Guang Yang, Chinese Academy of Sciences (CAS), China
Reviewed by
Daria Martynova, Zoological Institute (RAS), Russia; Elizaveta Ershova, UiT The Arctic University of Norway, Norway; Vladimir G. Dvoretsky, Murmansk Marine Biological Institute, Russia; Emilia Trudnowska, Polish Academy of Sciences, Poland
Updates
Copyright
© 2023 Ishihara, Matsuno, Tokuhiro, Ando, Sato and Yamaguchi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Minami Ishihara, minami0525konakichi@eis.hokudai.ac.jp
Disclaimer
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.