ORIGINAL RESEARCH article

Front. Mar. Sci., 23 June 2020

Sec. Marine Ecosystem Ecology

Volume 7 - 2020 | https://doi.org/10.3389/fmars.2020.00391

Seasonal Changes in Zooplankton Community Structure and Distribution Pattern in the Yellow Sea, China

  • 1. Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Shandong Provincial Key Laboratory for Fishery Resources and Eco-Environment, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, China

  • 2. Laboratory for Marine Ecology and Environmental Science, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao, China

  • 3. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao, China

  • 4. Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, BC, Canada

  • 5. Institute for the Oceans and Fisheries, University of British Columbia, Vancouver, BC, Canada

  • 6. Hakai Institute, Heriot Bay, BC, Canada

Abstract

Changes in zooplankton community and distribution have significant influences on fishery resources due to their vital role in the marine food web. Seasonal variations in zooplankton community structure, abundance, and biomass of major taxa and distribution pattern in the Yellow Sea, China, were analyzed to obtain the current status of secondary production. A total of 73 taxa (mostly at the species level) were recorded with the highest species richness in summer and lowest in spring. The most abundant zooplankton species in all seasons were copepods Oithona similis, Paracalanus sp., and Calanus sinicus. The mean total zooplankton abundance was the highest in spring and followed a declining trend till winter. Cluster analysis identified two distinct zooplankton assemblages geographically found in the deep (>50 m depth) and shallow regions. Zooplankton in the shallow region had higher abundance and species richness with the exception of winter season. Small copepods exhibited higher biomass in the shallow region but few abundance differences between the two domains. Large copepods (LC) and giant crustaceans (GC) usually had higher biomass in the deep region with the exceptions of LC in spring and GC in summer. Both the abundance and biomass of chaetognaths were significantly higher in the shallow region, where cnidarian abundance was also higher from spring to autumn. Water currents contributed to the transport of zooplankton between the two domains. The high abundance of copepods in the shallow region satisfied the requirements for the larval fish survival from spring to autumn. However, the high abundance of carnivorous zooplankton in summer–autumn may compete with larval fish for prey, and may directly feed on fish larvae.

Introduction

Zooplankton are vital secondary producers in marine ecosystems, linking primary productivity and high trophic level species (Turner, 2004; Robert et al., 2014). Additionally, zooplankton are the main drivers of marine biological pump and have important contribution to the vertical transport of carbon (Steinberg and Landry, 2017). Therefore, changes in zooplankton community structure and distribution can influence the biogeochemical cycles and energy flows in aquatic ecosystems (Beaugrand et al., 2003; Mitra et al., 2014; Brun et al., 2019). Zooplankton are also good indicators of environmental variations caused by climate change or pollution (Hays et al., 2005; Batchelder et al., 2013). Understanding seasonal changes in zooplankton community structure and distribution pattern is of great value in supporting research on the ecosystem dynamics and the potential impacts of climate change.

The Yellow Sea is very productive and an important marine fishery production area in China. Affected by climate change and anthropogenic activities, the Yellow Sea ecosystem is undergoing structural and functional changes (Tang, 2009), including increases in seawater temperature and nutrient concentrations (Lin et al., 2005; Belkin, 2009), changes in plankton communities (Fu et al., 2012; Shi et al., 2016; Yang et al., 2016), increases in occurrence frequency and intensity of ecological disasters (Dong et al., 2010; Yu and Liu, 2016), and changes in major fishery stocks and fish community (Zhao et al., 2003; Zhang et al., 2007). Accordingly, the Chinese government has introduced a series of programs, such as closed seasons and zones, stock enhancement, and total allowable catch system, for sustainable fisheries and aquaculture in the Yellow Sea (Tang et al., 2016).

Evaluating the fishery resources scientifically and accurately is vital to build sustainable fisheries management, and it is important to know the secondary production dynamics of the supporting pelagic ecosystem (Sun, 2016). Many fishery species feed on zooplankton and undergo seasonal migration in the Yellow Sea, wintering in its deep region as well as spawning and nursing grounds in the shallow region (Jin et al., 2006; Xu and Chen, 2009). However, few studies focused on the seasonal changes in zooplankton community structure and distribution in the Yellow Sea.

Here we analyzed seasonal variations in zooplankton species composition, abundance, biomass, community structure, and distribution patterns in the Yellow Sea using most recent data. Our results will be beneficial for subsequent studies on the ecosystem ecological carrying capacity and long-term changes of zooplankton community, and contribute to the scientific evaluation and management of fishery resources in the Yellow Sea.

Materials and Methods

Sampling and Laboratory Processing

Four basin wide research surveys onboard the R.V. Beidou were conducted in the Yellow Sea from 21 December 2016 to 4 January 2017 (winter), from 10 to 22 May (spring), from 9 to 27 August (summer), and from 13 October to 1 November 2017 (autumn). Fifty-two stations, belong to nine transects (A–I), were sampled between 33.00°N and 37.00°N and between 120.50°E and 124.25°E (Figure 1).

Figure 1

Zooplankton was vertically sampled using a Juday–Bogorov type net (diameter: 50 cm, mesh size: 160 μm) from the near bottom to the surface at each station. Confined by available ship-time, the sampling times were not fixed at day or night. The net was retrieved at a speed of ~1 m s−1, and a calibrated flowmeter (Hydrobios) was used to measure the volume of filtered water. After collection, all samples were preserved in 5% neutral formalin seawater solution.

Vertical temperature and salinity profiles at each station were obtained by a Sea-Bird CTD instrument (SBE-19). Meanwhile, 500 mL seawater samples from each depth (0, 20, and 2–5 m from the bottom) were filtered onto a Whatman GF/F filter. The filters were extracted in 90% aqueous acetone for 24 h, then used to determine the content of chlorophyll a (Chl a) fluorometrically before and after acidification (Parsons et al., 1984).

In the laboratory, large zooplankton were firstly identified and counted. The remaining part was then subsampled until around 300–500 individuals were left in each subsample. All specimens in a subsample were identified to the species or the lowest taxonomic level and counted under a dissecting microscope. Zooplankton abundance (ind m−3) was calculated based on water volume filtered at each station. The abundance data of small copepods (SC), large copepods (LC), giant crustaceans (GC), and chaetognaths (CH) were transformed into dry weights (DW) according to length–dry weight relationships or references same as Sun et al. (2010). The dry weights of medusa (ME) and tunicates (TU) were not calculated due to their high water content.

Data Analysis

Multivariate analyses were performed using the PRIMER software V6.0 (Clarke and Gorley, 2006). Station clustering was performed for each survey dataset based on the Bray-Curtis similarity matrix of log (x+1) transformed zooplankton abundance and the average linkage group classification (Field et al., 1982). SIMPER (similarity percentages) analysis was used to assess the percent contributions of species to the similarity within groups and the dissimilarity between groups. ANOSIM (analysis of similarities) was used to test for the differences of zooplankton community structures between groups, with R-statistic value closer to 1 indicating greater difference.

The BIO-ENV procedure was used to evaluate the best sets of environmental factors (temperature, salinity, depth, Chl a) explaining the differences in zooplankton community structures. This process estimated the Spearman correlation coefficients (ρ) between the zooplankton and environmental factors similarity matrices, with ρ = 1 indicating a perfect match.

The zooplankton abundance, biomass, and environmental factors between groups were compared by one-way ANOVA using SPSS 19.0 software.

Results

Environmental Conditions

The distribution patterns of sea surface temperature (SST) and sea bottom temperature (SBT) were nearly consistent in winter, increasing southward and ranging from 7.5 to 13.0°C and from 7.9 to 12.8°C, respectively (Figure 2), as the water column was vertically well mixed (Figure 3). The SST also generally increased southward in the range of 11.8–20.2°C in spring, while the highest and lowest SBT occurred in the southwest and northwest coastal stations, respectively, ranging from 5.5 to 17.2°C (Figure 2). The SST ranges were 24.2–29.3°C and 18.5–21.0°C in summer and autumn, respectively. The Yellow Sea Cold Water Mass (YSCWM) begins to form in spring, well develops in summer, and decays in autumn. The YSCWM was located near the bottom of the central region of the Yellow Sea in summer and autumn, with <10°C SBT. However, the SBT did not exhibit the distribution pattern with lowest values in the deep area in spring, due to the presence of the Qingdao Cold Water Mass in the northwest coastal region (Figure 2). The highest SBT were 27.6 and 21.1°C in summer and autumn, respectively, and temperature fronts were found in the bottom layer at the edge of the YSCWM (Figure 2).

Figure 2

Figure 3

The sea surface salinity (SSS) and sea bottom salinity (SBS) showed small differences in winter like temperature, and a water tongue with high salinity extended westward along the 34.5°N transect (Figure 2). The SSS had similar distribution patterns in spring and summer, with lowest values in the farthest south stations. Generally, coastal stations had relatively low SSS and SBS in all four seasons, and YSCWM region stations had high SBS from spring to autumn (Figure 2).

The E (35.0°N) transect passed through the YSCWM region and was selected to depict the vertical profile of temperature and salinity (Figure 3). No thermocline was found in winter. The thermocline started to develop between 15 and 30 m water depth in spring. The SST at all stations was higher than 26°C, and a strong thermocline existed in summer. The YSCWM with <10°C SBT was well developed at stations E4–E8 in summer. The mixed layer depth expanded to about 0–38 m, and the YSCWM shrunk to stations E5–E8 in autumn. Both SSS and SBS increased from E1 to E8. Generally, bottom layer in the YSCWM region had the highest salinity.

In general, the chlorophyll a concentrations at surface, 20 m depth, and bottom layers were all higher at coastal shallow stations than those at deep stations (Figure 4). The bottom Chl a was very low in the YSCWM stations in summer and autumn (Figure 4).

Figure 4

Stations in the shallow region had depth < 52 m, while the depth range in the deep region was 52–84 m (Figure 1 and Table 1). With the exception in summer, SST showed no significant differences (P > 0.05) between the shallow and deep regions in other three seasons. SBT was significantly lower (P < 0.01) in the deep region than in the shallow region from spring to autumn. SSS and SBS had significantly higher values (P < 0.05) in the deep region in all four seasons with the exception of SSS in autumn. Average water column chlorophyll a concentration was significantly higher (P < 0.01) in the shallow region in all four seasons (Table 1).

Table 1

WinterSpringSummerAutumn
ShallowDeepShallowDeepShallowDeepShallowDeep
D (m)37.1 (19.0–50.0)72.9 (56.0–84.0)31.2 (9.0–48.0)70.2 (55.0–82.0)31.3 (12.0–52.0)68.5 (52.0–80.0)35.0 (12.0–52.0)72.8 (57.0–84.0)
SST (°C)11.0 (7.5–13.0)10.8 (9.1–12.2)16.8 (11.8–20.2)16.0 (13.4–19.0)27.2 (24.2–29.0)27.8 (26.3–29.3)19.6 (18.6–21.0)19.5 (18.5–20.4)
SBT (°C)11.1 (7.9–12.8)10.6 (8.8–12.0)11.9 (5.5–17.2)8.6 (5.9–11.0)16.6 (8.2–27.6)8.7 (7.2–10.5)19.0 (9.2–21.1)9.6 (7.9–15.0)
SSS31.7 (31.3–32.2)31.9 (31.3–32.4)31.7 (28.7–32.8)32.5 (32.1–33.5)31.1 (29.2–31.9)31.6 (30.9–32.0)31.9 (31.3–32.4)32.0 (31.5–32.7)
SBS31.9 (31.3–32.3)32.4 (31.8–33.0)32.4 (30.9–33.3)32.7 (32.2–33.6)32.2 (31.0–33.3)32.6 (32.3–33.1)31.9 (31.3–32.6)32.5 (32.1–32.9)
C (mg m−3)0.6 (0.4–1.0)0.4 (0.3–0.6)2.0 (0.6–6.4)0.5 (0.3–1.8)1.3 (0.4–3.8)0.5 (0.2–1.1)1.1 (0.6–2.4)0.5 (0.2–1.3)

Seasonal changes in environmental factors (mean and range) in the two domains in the Yellow Sea.

D, depth; SST, sea surface temperature; SBT, sea bottom temperature; SSS, sea surface salinity; SBS, sea bottom salinity; C, average water column chlorophyll a concentration.

Zooplankton Community Structure

Seventy-three taxa (mostly at the species level) were identified in our study. A detailed abundance description of 27 taxa (> 90% total zooplankton abundance at any station) is presented in Table 2. Copepods clearly dominated the zooplankton community, and the dominant species were Oithona similis, Paracalanus sp., and Calanus sinicus.

Table 2

WinterSpringSummerAutumn
ShallowDeepShallowDeepShallowDeepShallowDeep
Species number3731302055444533
Total abundance (ind m−3)2374.92380.316746.99435.810087.98281.57184.15657.6
Protozoa
Noctiluca scintillans46.8 (32)0.4 (4)3822.1 (54)17.3 (32)1234.0 (38)3.6 (4)104.6 (29)85.8 (11)
Cnidaria
Diphyes chamissonis2.5 (23)4.7 (73)0.5 (8)0.0 (0)13.6 (21)8 (33)33.7 (54)2.4 (19)
Muggiaea atlantica1.0 (23)1.6 (42)0.0 (0)0.0 (0)2.2 (8)2.1 (15)12.0 (38)0.1 (4)
Pteropoda
Creseis acicula0.0 (0)0.0 (0)0.0 (0)0.0 (0)0.0 (0)0.0 (0)13.7 (17)11.0 (30)
Cladocera
Penilia avirostris0.0 (0)0.0 (0)0.0 (0)0.0 (0)23.6 (42)36.0 (67)0.0 (0)0.0 (0)
Pseudevadne tergestina0.0 (0)0.0 (0)0.0 (0)0.0 (0)15.2 (33)174.1 (56)0.0 (0)0.0 (0)
Copepoda
Acartia hongi1.7 (27)9.1 (88)328.5 (96)856.8 (100)364.9 (75)379.3 (93)15.1 (38)5.8 (56)
Acartia pacifica0.0 (0)0.0 (0)0.0 (0)0.0 (0)142.6 (50)26.2 (19)62.7 (71)5.1 (41)
Calanus sinicus20.9 (86)77.7 (100)865.1 (96)411.4 (100)142.6 (67)249.6 (96)87.6 (79)123.1 (96)
Centropages abdominalis0.0 (0)0.1 (8)820.9 (79)174.6 (86)7.1 (21)0.9 (4)0.0 (0)0.0 (0)
Ditrichocorycaeus affinis19.4 (95)12.2 (100)321.1 (96)69 (93)254.9 (100)21.0 (85)195.4 (96)29.3 (81)
Euchaeta concinna3.7 (55)0.7 (35)0.2 (4)0.0 (0)10.2 (25)0.3 (7)28.7 (25)4.0 (19)
Labidocera euchaeta1.2 (36)0.0 (0)1.0 (8)0.0 (0)11.1 (21)0.1 (4)11.7 (38)0.0 (0)
Microsetella norvegica2.3 (36)0.0 (4)2.6 (8)0.1 (4)24.4 (13)0.5 (4)161.8 (21)0.2 (4)
Oithona similis1074.8 (100)1390.0 (100)5530.3 (96)6581.3 (100)2269.9 (100)3151.6 (100)1954.2 (100)2765.9 (100)
Paracalanus sp.1024.4 (100)777.9 (100)4205.9 (100)1128.4 (100)4696.2 (100)3665.4 (100)3411.5 (100)2139.1 (100)
Amphipoda
Themisto gracilipes0.6 (36)2.0 (92)10.0 (54)40 (100)4.9 (25)6.7 (63)0.8 (8)12.1 (59)
Euphausiacea
Euphausia pacifica1.6 (50)5.2 (88)14.7 (54)13.8 (79)3.1 (25)1.4 (52)1.4 (13)1.8 (48)
Pseudeuphausia sinica0.4 (27)0.5 (15)0.0 (0)0.0 (4)97.3 (25)0.0 (0)10.0 (50)0.4 (4)
Chaetognatha
Aidanosagitta crassa29.7 (100)9.5 (100)33.8 (63)1.1 (36)238.5 (79)98.9 (100)221.8 (96)93.5 (100)
Flaccisagitta enflata0.7 (32)0.0 (0)0.0 (0)0.0 (0)20.6 (29)3.7 (30)6.4 (25)9.6 (59)
Tunicata
Doliolum denticulatum0.0 (0)0.0 (0)0.0 (0)0.0 (0)0.2 (4)388.0 (78)29.1 (42)121.0 (59)
Oikopleura dioica12.5 (91)31.7 (100)163.8 (63)11.8 (14)68.4 (46)39.8 (67)477.3 (83)151.3 (100)
Pelagic larvae
Polychaeta larvae0.6 (32)0.2 (8)0.0 (0)0.0 (0)11.2 (46)0.5 (19)10.0 (42)0.1 (4)
Gastropoda larvae12.2 (77)9.4 (88)14.2 (25)0.5 (11)111.0 (67)5.8 (56)22.0 (38)12.6 (41)
Bivalvia larvae12.9 (59)4.3 (38)0.8 (4)0.0 (0)98.3 (75)1.4 (15)5.6 (25)0.4 (7)
Copepoda nauplii103.5 (95)42.3 (92)588 (71)128.5 (82)156.6 (58)9.2 (11)266.8 (96)80.8 (74)

Average abundances (ind m−3) of zooplankton species that contributed > 2% to community dissimilarities between domains (underlined numbers) and species number identified in the shallow and deep regions during four seasons in the Yellow Sea.

Values in parentheses represent the percentage of species occurrence within the domain (%).

In all four seasons, zooplankton assemblage in the deep stations was separated from that in the coastal shallow stations based on cluster analysis of zooplankton abundance (Figure 5). The two domains, named shallow and deep regions, were divided by the 50 m isobath (Figures 1, 5).

Figure 5

ANOSIM results showed that the zooplankton assemblages differed significantly (P < 0.001) between the shallow and deep regions in all four seasons (Table 3). SIMPER results indicated that the dissimilarity of zooplankton assemblages between the two domains was highest (47.3%) in summer and lowest (30.6%) in winter (Table 3). Zooplankton assemblage had higher similarities in the deep region than in the shallow region, and the similarities were highest in winter and lowest in summer in both domains (Table 3).

Table 3

SeasonANOSIMSIMPER average dissimilarity (%)SIMPER average similarity (%)
RP (%)ShallowDeep
Winter0.3920.130.669.678.2
Spring0.3420.131.168.078.2
Summer0.4730.147.353.568.9
Autumn0.3820.139.461.570.7

Comparison of the zooplankton community structures between the shallow and deep regions in the Yellow Sea according to ANOSIM (R value and significance level) and SIMPER.

Seasonal Variation in Zooplankton Community

The mean total zooplankton abundance in the Yellow Sea was the highest in spring and followed a declining trend till winter (Table 2). Zooplankton in the shallow region usually had higher abundance and higher species number with the exception of abundance in winter (Table 2). The species number was 37 and 31 in the shallow and deep regions, respectively, in winter, while the numbers decreased to 30 and 20 in the two domains in spring. Species richness peaked in summer with 55 in the shallow region and 44 in the deep region. Then 45 and 33 zooplankton species were identified in autumn in the shallow and deep regions, respectively (Table 2).

In winter, based on the SIMPER analysis, C. sinicus, Copepoda nauplii, Acartia hongi, Bivalvia larvae, Oikopleura dioica, Gastropoda larvae, Euphausia pacifica, Aidanosagitta crassa, and Diphyes chamissonis all contributed > 5% to the zooplankton community dissimilarity between the shallow and deep regions. In spring, SIMPER results showed that Noctiluca scintillans, Copepoda nauplii, O. dioica, Centropages abdominalis, Themisto gracilipes, Ditrichocorycaeus affinis, A. crassa, and A. hongi contributed the most to the dissimilarity between domains. In summer, Doliolum denticulatum, Copepoda nauplii, Bivalvia larvae, C. sinicus, A. hongi, N. scintillans, and Pseudevadne tergestina were the main contributors to the community dissimilarity. In autumn, the main contributors to the dissimilarity were Copepoda nauplii, D. affinis, Acartia pacifica, D. denticulatum, O. dioica, C. sinicus, N. scintillans, and D. chamissonis (Table 2). Additionally, SIMPER results showed that O. similis, Paracalanus sp., C. sinicus and D. affinis all contributed > 4% to the zooplankton community similarity within any domain in any season.

Species with > 10% occurrence of total stations and accounting for > 1% of the total abundance of each main zooplankton taxa were defined as dominant species (Table 4). No significant differences were found in SC abundances between the shallow and deep regions in all four seasons (Figure 6A). Only in winter LC and GC abundances were significantly higher (P < 0.01) in the deep region than in the shallow region. CH abundance in the shallow region was significantly higher (P < 0.05) in all four seasons. ME in the shallow region had significantly higher abundance (P < 0.01) in autumn, while TU showed significant abundance differences (P < 0.05) between the two domains in winter, spring, and summer (Figure 6A).

Table 4

Zooplankton taxonAcronymDominant species
Small copepodsSCOithona similis, Paracalanus sp., Acartia hongi, Ditrichocorycaeus affinis, Centropages abdominalis
Large copepodsLCCalanus sinicus, Euchaeta concinna, Labidocera euchaeta
Giant crustaceansGCEuphausia pacifica, Pseudeuphausia sinica, Themisto gracilipes
ChaetognathsCHAidanosagitta crassa, Flaccisagitta enflata
MedusaeMEDiphyes chamissonis, Muggiaea atlantica
TunicatesTUDoliolum denticulatum, Oikopleura dioica

Dominant species of six zooplankton taxa in the Yellow Sea.

Figure 6

The total biomass of the four zooplankton taxa (SC, LC, GC, and CH) was much lower in winter than that in other three seasons (Figure 6B). SC had higher biomass in the shallow region in all four seasons, with significant difference (P < 0.05) only in spring. LC had significantly higher biomass (P < 0.05) in the shallow region in spring and in the deep region in autumn. GC showed significant biomass difference (P < 0.01) between the two domains only in winter, with higher biomass in the deep region. CH biomass in the shallow region was significantly higher (P < 0.05) in all four seasons like abundance (Figure 6B).

GC biomass was usually higher in the deep stations, while CH usually had higher biomass in the shallow stations along the E transect (Figure 7). In spring, LC had high biomass in all the stations. LC and SC biomasses were much higher in station E7 in the YSCWM region in summer. The total biomass of SC, LC, and GC was relatively low in winter and autumn (Figure 7).

Figure 7

Factors Affecting Community Structure

According to the BIO-ENV analysis (Table 5), a combination of depth, SST, and Chl a were the best environment variables to explain the variance in the community structure based on abundance in winter (ρ = 0.483). A combination of depth, SST, SBT, and SSS explained the most in spring (ρ = 0.516), while SBT was the best indicator in summer (ρ = 0.679), and depth and SBS were the best combination to indicate differences in the zooplankton community in autumn (ρ = 0.492) (Table 5).

Table 5

K*WinterSpringSummerAutumn
1D (0.377), SST (0.286)D (0.417), SBT (0.391)SBT (0.679), D (0.508)D (0.421), SBT (0.347)
2D, SST (0.478)D, SST (0.478)SBT, C (0.644)D, SBS (0.492)
3D, SST, C (0.483)D, SST, SBT (0.507)D, SBT, SBS (0.648)D, SBT, SBS (0.477)
4D, SST, SBT, C (0.470)D, SST, SBT, SSS (0.516)D, SBT, SSS, SBS (0.659)D, SBT, SSS, SBS (0.470)
5D, SST, SBT, SBS, C (0.464)D, SST, SBT, SSS, SBS (0.487)D, SBT, SSS, SBS, C (0.664)D, SBT, SSS, SBS, C (0.468)
6D, SST, SBT, SSS, SBS, C (0.456)D, SST, SBT, SSS, SBS, C (0.461)D, SST, SBT, SSS, SBS, C (0.595)D, SST, SBT, SSS, SBS, C (0.457)

BIO-ENV analysis of zooplankton community structure in the Yellow Sea during four seasons to depth (D), sea surface temperature (SST), sea bottom temperature (SBT), sea surface salinity (SSS), sea bottom salinity (SBS) and average water column chlorophyll a concentration (C).

*

means a combination of K variables. ρ is given in parentheses. Best combinations explaining clustering for zooplankton community structure are in bold.

Discussion

Community Structure and Distribution

Zooplankton species assemblages were associated with water mass distributions (Zuo et al., 2006; Domínguez et al., 2017). Based on the community structure, two distinct zooplankton assemblages, namely the Yellow Sea neritic and basin assemblages, spatially coincided with the shallow and deep regions of the study area (Zheng, 1965; Chen et al., 1980; Zuo et al., 2006). Zooplankton in the shallow region showed relatively low similarities (Table 3), due to the wide latitude range and seasonal variations in boundaries between zooplankton assemblages (Zheng, 1965; Chen et al., 2011). However, since two domains were separated by the 50 m isobath that demarcates ecologically important spawning-nursing and wintering grounds for many fishery species in the Yellow Sea, we assumed that the boundaries between two zooplankton assemblages were stable across four seasons.

Small copepods O. similis and Paracalanus sp. were the most abundant species in the present study with different distribution patterns. Oithona similis and A. hongi had higher abundances in the deep region (Table 2), which was consistent with previous observations in the Yellow Sea (Zuo et al., 2006; Shi et al., 2018); while the other small copepods, Paracalanus sp., D. affinis, C. abdominalis, A. pacifica, and Microsetella norvegica, were mostly distributed in the shallow region (Table 2). The comparable abundances between O. similis and Paracalanus sp. resulted in the small differences of SC abundances in the two domains. When transformed into dry weights, SC biomass was higher in the shallow region benefited from the higher DW of Paracalanus sp. (Figure 6).

Calanus sinicus dominated LC group, and most contributed to the LC abundance and biomass in the Yellow Sea (Sun et al., 2010; Huo et al., 2012). Higher LC biomass and C. sinicus abundance were found in the shallow region only in spring, confirming the C. sinicus life strategy that the population mainly propagated in the shallow region in spring, then moved toward the YSCWM region to survive through the hot summer (Sun, 2005).

The dominant species of GC were E. pacifica, T. gracilipes, and Pseudeuphausia sinica. Previous studies showed that high GC biomass was mainly located in the deep region in the Yellow Sea (Sun et al., 2010), as well as the E. pacifica and T. gracilipes populations (Wang and Zuo, 2004; Sun et al., 2011). Euphausia pacifica preferred to stay in the deep water (Wang and Zuo, 2004), making it seldom be sampled which might result in the low abundance in the present study. Compared with that higher abundance of T. gracilipes occurred in the deep region, E. pacifica did not show clear distribution preference from the view of abundance; however, both the two species exhibited much higher occurrence percentage in the deep region (Table 2). The high GC biomass in the shallow region in summer was contributed by the high P. sinica abundance.

Significantly higher CH abundance and biomass were found in the shallow region during all four seasons in the Yellow Sea, and ME had higher abundance in the shallow region from spring to autumn, suggesting their coastal distribution preference (Sun et al., 2012; Wang et al., 2013, 2016). TU dominant species D. denticulatum and O. dioica showed different distribution patterns, that D. denticulatum was distributed in the deep region while O. dioica was mainly located in the shallow region.

In the Japan Sea, a clear seasonal succession of zooplankton community structure occurred from dominance by cold-water copepods in winter–spring to prevalence of gelatinous and carnivorous plankton and small warm-water copepods in summer–autumn (Chiba and Saino, 2003). In the present study, the abundances of CH, ME, and tropical species Euchaeta concinna and D. denticulatum also increased in summer–autumn.

Environmental Factors Affecting Zooplankton Community

Depth was the most important single factor that affected the zooplankton community in the Yellow Sea (Table 5). Shallow water may enhance the mortality rate of zooplankton (Uye, 2000). Additionally, depth was the main contributor to the differences of environmental factors between the two domains in the present study. In the shallow region, the water column was relatively vertically well mixed due to tidal mixing (Liu et al., 2003), and an upwelling occurred in the frontal zone of the YSCWM in summer (Wei et al., 2016, 2018), both benefiting nutrient supply and phytoplankton reproduction. However, the more stratified conditions in the deep region in summer–autumn reduced nutrient concentrations at surface and maintained the cold and low Chl a YSCWM at bottom (Wei et al., 2016).

Temperature influenced the zooplankton community in two ways. Based on the view of seasonal succession, gelatinous plankton and warm-water species peaked in summer–autumn due to higher temperatures. Meanwhile, the YSCWM in the deep region provided a refuge for some low-temperature species, like C. sinicus and E. pacifica, to survive in summer (Sun, 2005; Sun et al., 2011). Higher abundances of low-saline species Labidocera euchaeta, A. pacifica, and Bivalvia larvae occurred in the shallow region suggesting the effect of salinity on the zooplankton distribution (Zuo et al., 2006; Shi et al., 2018). Rich nutrients and high Chl a in the shallow region benefited heterotrophic dinoflagellate N. scintillans population development (Glibert et al., 2018), where was also an ideal area for copepod recruitment (Sun, 2005; Shi et al., 2015), resulting in higher abundances of many small copepods and Copepoda nauplii during all four seasons as well as higher C. sinicus abundance in spring in the shallow region during the present study.

Water currents affected the zooplankton community by advective transport (Lü et al., 2013; Smeti et al., 2015; Peterson et al., 2017). The Yellow Sea Warm Current (YSWC), which is a prominent feature in winter in the Yellow Sea (Lin et al., 2011), brought tropical zooplankton species into the Yellow Sea and affected zooplankton biomass and migration patterns in the intrusion area (Wang and Zuo, 2004; Lü et al., 2013; Chen et al., 2020). Cladocerans P. tergestina and Penilia avirostris are low-saline species, with distribution associated with coastal water masses (Xu et al., 2007; Domínguez et al., 2017). High abundance of cladocerans that occurred in the deep region in summer in the present study suggested their transport by coastal current and river front.

Ecological Function of Zooplankton Taxa and Its Implications

The crustaceans SC, LC, and GC are the main food sources of many fish. As larval fish grow and the body length increases, the prey selectivity shifts from SC to LC, and then to GC (Meng, 2003). Importantly, small copepods are preferred prey during the critical fish larval stage due to their high abundance, appropriate size, and good nutritional value (Llopiz, 2013; Robert et al., 2014), affecting fish larvae survival and a subsequent recruitment to fish stocks (Hjort, 1914; Cushing, 1990).

The deep region of the Yellow Sea served as wintering grounds for many fishery species, while the shallow region is known to provide spawning, nursery, and foraging habitats (Jin et al., 2006; Xu and Chen, 2009). Higher abundances of LC and GC in the deep region in winter may support the wintering fish community, while higher SC and LC biomass and Copepoda nauplii abundance in the shallow region in spring could provide stable feeding grounds for larval fish. Additionally, the sampling net used in the present study (160 μm mesh) may have led to an underestimation of small copepods and copepod nauplii abundances (Atkinson et al., 2012; Chen et al., 2016).

Although SC biomass was relatively low when compared with LC and GC, SC had the highest production rate (Huo et al., 2012). Prey concentrations for maximum growth of larval marine fish range between 5 and 10 copepodites l−1 (Munk, 1995; Peck and Daewel, 2007). In the present study, the abundance of small copepods in the shallow region satisfied the requirements of > 5 prey l−1 for the larval fish survival from spring to autumn (Peck and Daewel, 2007). However, the high abundance of carnivorous CH and ME in summer–autumn may compete with larval fish for prey, and may directly feed on fish larvae. This study provides a foundation for future research to understand the inter-annual variability in seasonal development of plankton communities and main factors driving this variability.

Statements

Data availability statement

The datasets generated for this study are available on request to the corresponding author.

Author contributions

JW, XS, and XJ contributed to the conception of the study. YS, JS, and WY contributed to sample collection and data collation. YS and TZ analyzed the samples and conducted the data analyses. YS and EP wrote the manuscript and contributed to the definition of aims. All authors contributed to the improvement of the manuscript before approving the submission.

Funding

This study was supported by the National Key R&D Program of China (2018YFD0900901), the Creative Team Project of QNLM LMEES (LMEES-CTSP-2018-4), the National Natural Science Foundation of China (31872692), the Central Public-interest Scientific Institution Basal Research Fund, YSFRI, CAFS (20603022020009), the Open Fund of QNLM LMEES (LMEES201807), Central Public-interest Scientific Institution Basal Research Fund, CAFS (2020TD01), and the China Scholarship Council (201803260001).

Acknowledgments

We are grateful to the captain and crew of the R.V. Beidou for their efforts in the field and to colleagues who provided support during our sampling efforts. We are grateful to the University of British Columbia's Department of Earth, Ocean and Atmospheric Sciences for accommodating YS during his research visit.

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.

References

  • 1

    AtkinsonA.WardP.HuntB. P. V.PakhomovE. A.HosieG. W. (2012). An overview of Southern Ocean zooplankton data: abundance, biomass, feeding and functional relationships. CCAMLR Sci.19, 171218.

  • 2

    BatchelderH. P.DalyK. L.DavisC. S.JiR.OhmanM. D.PetersonW. T.et al. (2013). Climate impacts on zooplankton population dynamics in coastal marine ecosystems. Oceanography26, 3451. 10.5670/oceanog.2013.74

  • 3

    BeaugrandG.BranderK. M.LindleyJ. A.SouissiS.ReidP. C. (2003). Plankton effect on cod recruitment in the North Sea. Nature426, 661664. 10.1038/nature02164

  • 4

    BelkinI. M. (2009). Rapid warming of Large Marine Ecosystems. Prog. Oceanogr.81, 207213. 10.1016/j.pocean.2009.04.011

  • 5

    BrunP.StamieszkinK.VisserA. W.LicandroP.PayneM. R.KiørboeT. (2019). Climate change has altered zooplankton-fuelled carbon export in the North Atlantic. Nat. Ecol. Evol.3, 416423. 10.1038/s41559-018-0780-3

  • 6

    ChenH. J.QiY. P.LiuG. X. (2011). Spatial and temporal variations of macro- and mesozooplankton community in the Huanghai Sea (Yellow Sea) and East China Sea in summer and winter. Acta Oceanol. Sin.30, 8495. 10.1007/s13131-011-0108-5

  • 7

    ChenH. J.YuH.LiuG. X. (2016). Comparison of copepod collection efficiencies by three commonly used plankton nets: A case study in Bohai Sea, China. J. Ocean Univ. China15, 10071013. 10.1007/s11802-016-3122-6

  • 8

    ChenH. J.ZhuY. Z.LiuG. X. (2020). Characteristics of zooplankton community in North Yellow Sea unveiled an indicator species for the Yellow Sea Warm Current in winter: Euchaeta plana. J. Oceanol. Limnol. 10.1007/s00343-019-9165-y. [Epub ahead of print].

  • 9

    ChenQ. C.ChenY. Q.HuY. Z. (1980). Preliminary study on the plankton communities in the southern Yellow Sea and the East China Sea [In Chinese, with English Abstract]. Acta Oceanol. Sin.2, 149157.

  • 10

    ChibaS.SainoT. (2003). Variation in mesozooplankton community structure in the Japan/East Sea (1991–1999) with possible influence of the ENSO scale climatic variability. Prog. Oceanogr.57, 317339. 10.1016/S0079-6611(03)00104-6

  • 11

    ClarkeK. R.GorleyR. N. (2006). PRIMER v6: User Manual/Tutorial.Plymouth: PRIMER-E.

  • 12

    CushingD. (1990). Plankton production and year-class strength in fish populations: an update of the match/mismatch hypothesis. Adv. Mar. Biol.26, 249293. 10.1016/S0065-2881(08)60202-3

  • 13

    DomínguezR.GarridoS.SantosA. M. P.dos SantosA. (2017). Spatial patterns of mesozooplankton communities in the Northwestern Iberian shelf during autumn shaped by key environmental factors. Estuar. Coast. Shelf Sci.198, 257268. 10.1016/j.ecss.2017.09.008

  • 14

    DongZ. J.LiuD. Y.KeesingJ. K. (2010). Jellyfish blooms in China: Dominant species, causes and consequences. Mar. Pollut. Bull.60, 954963. 10.1016/j.marpolbul.2010.04.022

  • 15

    FieldJ. G.ClarkeK. R.WarwickR. M. (1982). A practical strategy for analysing multispecies distribution patterns. Mar. Ecol. Prog. Ser.8, 3752. 10.3354/meps008037

  • 16

    FuM. Z.WangZ. L.PuX. M.XuZ. J.ZhuM. Y. (2012). Changes of nutrient concentrations and N:P:Si ratios and their possible impacts on the Huanghai Sea ecosystem. Acta Oceanol. Sin.31, 101112. 10.1007/s13131-012-0224-x

  • 17

    GlibertP. M.BerdaletE.BurfordM. A.PitcherG. C.ZhouM. J. (2018). Global Ecology and Oceanography of Harmful Algal Blooms. Ecological Studies 232. Springer. 10.1007/978-3-319-70069-4

  • 18

    HaysG. C.RichardsonA. J.RobinsonC. (2005). Climate change and marine plankton. Trends Ecol. Evol.20, 337344. 10.1016/j.tree.2005.03.004

  • 19

    HjortJ. (1914). Fluctuations in the great fisheries of northern Europe viewed in the light of biological research. Cons. Int. Explor. Mer.20, 1228.

  • 20

    HuoY. Z.SunS.ZhangF.WangM. X.LiC. L.YangB. (2012). Biomass and estimated production properties of size-fractionated zooplankton in the Yellow Sea, China. J. Mar. Syst.94, 18. 10.1016/j.jmarsys.2011.09.013

  • 21

    JinX. S.ChengJ. S.QiuS. Y.LiP. J.CuiY.DongJ. (2006). Comprehensive Research and Evaluation of Fishery Resources in the Yellow Sea and Bohai Sea.Beijing: China Ocean Press (In Chinese).

  • 22

    LinC.NingX.SuJ.LinY.XuB. (2005). Environmental changes and the responses of the ecosystems of the Yellow Sea during 1976–2000. J. Mar. Syst.55, 223234. 10.1016/j.jmarsys.2004.08.001

  • 23

    LinX.YangJ.GuoZ.ZhangZ.YinY.SongX.et al. (2011). An asymmetric upwind flow, Yellow Sea Warm Current: 1. New observations in the western Yellow Sea. J. Geophys. Res.116, C04026. 10.1029/2010JC006513

  • 24

    LiuG. M.SunS.WangH.ZhangY.YangB.JiP. (2003). Abundance of Calanus sinicus across the tidal front in the Yellow Sea, China. Fish. Oceanogr.12, 291298. 10.1046/j.1365-2419.2003.00253.x

  • 25

    LlopizJ. K. (2013). Latitudinal and taxonomic patterns in the feeding ecologies of fish larvae: a literature synthesis. J. Mar. Syst.109–110, 6977. 10.1016/j.jmarsys.2012.05.002

  • 26

    L. G.WangX.WangH. W.LiL.YangG. B. (2013). The variations of zooplankton biomass and their migration associated with the Yellow Sea Warm Current. Cont. Shelf Res.64, 1019. 10.1016/j.csr.2013.05.007

  • 27

    MengT. X. (2003). Studies on the feeding of anchovy (Engraulis japonicus) at different life stages on zooplankton in the Middle and Southern Waters of the Yellow Sea [In Chinese, with English Abstract]. Mar. Fish. Res.24, 19.

  • 28

    MitraA.CastellaniC.GentlemanW. C.JónasdóttirS. H.FlynnK. J.BodeA.et al. (2014). Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link. Prog. Oceanogr.129, 176199. 10.1016/j.pocean.2014.04.025

  • 29

    MunkP. (1995). Foraging behaviour of larval cod (Gadus morhua) influenced by prey density and hunger. Mar. Biol.122, 205212.

  • 30

    ParsonsT. R.MaitaY.LalliC. M. (1984). A Manual of Chemical and Biological Methods for Seawater Analysis.Oxford: Pergamon Press.

  • 31

    PeckM. A.DaewelU. (2007). Physiologically based limits to food consumption, and individual-based modeling of foraging and growth of larval fishes. Mar. Ecol. Prog. Ser.347, 171183. 10.3354/meps06976

  • 32

    PetersonW. T.FisherJ. L.StrubP. T.DuX.RisienC.PetersonJ.et al. (2017). The pelagic ecosystem in the Northern California Current off Oregon during the 2014–2016 warm anomalies within the context of the past 20 years. J. Geophys. Res. Oceans122, 72677290. 10.1002/2017JC012952

  • 33

    RobertD.MurphyH. M.JenkinsG. P.FortierL. (2014). Poor taxonomical knowledge of larval fish prey preference is impeding our ability to assess the existence of a “critical period” driving year-class strength. ICES J. Mar. Sci.71, 20422052. 10.1093/icesjms/fst198

  • 34

    ShiY. Q.SunS.LiC. L.ZhangG. T. (2016). Interannual changes in the abundance of zooplankton functional groups in the southern Yellow Sea in early summer [In Chinese, with English Abstract]. Oceanol. Limnol. Sin.47, 18.

  • 35

    ShiY. Q.SunS.ZhangG. T.WangS. W.LiC. L. (2015). Distribution pattern of zooplankton functional groups in the Yellow Sea in June: a possible cause for geographical separation of giant jellyfish species. Hydrobiologia754, 4358. 10.1007/s10750-014-2070-7

  • 36

    ShiY. Q.ZuoT.YuanW.SunJ. Q.WangJ. (2018). Spatial variation in zooplankton communities in relation to key environmental factors in the Yellow Sea and East China Sea during winter. Cont. Shelf Res.170, 3341. 10.1016/j.csr.2018.10.004

  • 37

    SmetiH.PaganoM.MenkesC.Lebourges-DhaussyA.HuntB. P. V.AllainV.et al. (2015). Spatial and temporal variability of zooplankton off New Caledonia (Southwestern Pacific) from acoustics and net measurements. J. Geophys. Res. Oceans120, 26762700. 10.1002/2014JC010441

  • 38

    SteinbergD. K.LandryM. R. (2017). Zooplankton and the ocean carbon cycle. Annu. Rev. Mar. Sci.9, 413444. 10.1146/annurev-marine-010814-015924

  • 39

    SunS. (2005). Over-summering strategy of Calanus sinicus. GLOBEC Int. Newsl.11, 34.

  • 40

    SunS. (2016). Fishery Version 3.0. [In Chinese, with English Abstract]. Bull. Chin. Acad. Sci.31, 13321338. 10.16418/j.issn.1000-3045.2016.12.007

  • 41

    SunS.HuoY. Z.YangB. (2010). Zooplankton functional groups on the continental shelf of the yellow sea. Deep Sea Rea. II57, 10061016. 10.1016/j.dsr2.2010.02.002

  • 42

    SunS.TaoZ. C.LiC. L.LiuH. L. (2011). Spatial distribution and population structure of Euphausia pacifica in the Yellow Sea (2006–2007). J. Plankton Res.33, 873889. 10.1093/plankt/fbq160

  • 43

    SunS.ZhangF.LiC. L.YangB.JiP. (2012). The distribution pattern of small medusae in the Yellow Sea [In Chinese, with English Abstract]. Oceanol. Limnol. Sin.43, 429437.

  • 44

    TangQ. S. (2009). Changing states of the Yellow Sea Large Marine Ecosystem: anthropogenic forcing and climate impacts, in Sustaining the World's Large Marine Ecosystems, eds ShermanK.AquaroneM.C.AdamsS. (Gland: IUCN), 7788.

  • 45

    TangQ. S.YingY. P.WuQ. (2016). The biomass yields and management challenges for the Yellow sea large marine ecosystem. Environ. Dev.17, 175181. 10.1016/j.envdev.2015.06.012

  • 46

    TurnerJ. T. (2004). The importance of small planktonic copepods and their roles in pelagic marine food webs. Zool. Stud.43, 255266.

  • 47

    UyeS. I. (2000). Why does Calanus sinicus prosper in the shelf ecosystem of the Northwest Pacific Ocean?ICES J. Mar. Sci.57, 18501855. 10.1006/jmsc.2000.0965

  • 48

    WangR.ZuoT. (2004). The Yellow Sea Warm Current and the Yellow Sea Cold Bottom Water, their impact on the distribution of zooplankton in the southern Yellow Sea. J. Korean Soc. Oceanogr.39, 113.

  • 49

    WangX.JiangM. J.LiuP.ZhangX. L.WangY.WangZ. L. (2016). Distribution pattern of zooplankton and its influencing factors in the South Yellow Sea in autumn [In Chinese, with English Abstract]. Haiyang Xuebao38, 125134. 10.3969/j.issn.0253-4193.2016.10.013

  • 50

    WangX.WangZ. L.PuX. M.LiuP. (2013). Analysis of the distribution of zooplankton in the Southern Yellow Sea in summer and its influencing factors [In Chinese, with English Abstract]. Acta Oceanol. Sin.35, 147155. 10.3969/j.issn.0253-4193.2013.05.016

  • 51

    WeiQ. S.LiX. S.WangB. D.FuM. Z.GeR. F.YuZ. G. (2016). Seasonally chemical hydrology and ecological responses in frontal zone of the central southern Yellow Sea. J. Sea Res.112, 112. 10.1016/j.seares.2016.02.004

  • 52

    WeiQ. S.WangB. D.YaoQ. Z.FuM. Z.SunJ. C.XuB. C.et al. (2018). Hydro-biogeochemical processes and their implications for Ulva prolifera blooms and expansion in the world's largest green tide occurrence region (Yellow Sea, China). Sci. Total Environ.645, 257266. 10.1016/j.scitotenv.2018.07.067

  • 53

    XuZ. L.ChenJ. J. (2009). Analysis on migratory routine of Larimichthy polyactis [In Chinese, with English Abstract]. J. Fish. Sci. China16, 931940. 10.3321/j.issn:1005-8737.2009.06.014

  • 54

    XuZ. L.GaoQ.ChenH.ChenJ. J.CaiM. (2007). Ecological adaptation of pelagic Cladocera and Cumacea in East China Sea [In Chinese, with English Abstract]. Chin. J. Ecol.26, 17821787. 10.13292/j.1000-4890.2007.0319

  • 55

    YangY.SunX. X.ZhengS.ZhaoY. F. (2016). Distribution and long-term variation of net phytoplankton in south Yellow Sea in early summer [In Chinese, with English Abstract]. Oceanol. Limnol. Sin.47, 755763. 10.11693/hyhz20160200028

  • 56

    YuR. C.LiuD. Y. (2016). Harmful algal blooms in the coastal waters of China: current situation, long-term changes and prevention strategies [In Chinese, with English Abstract]. Bull. Chin. Acad. Sci.31, 11671174. 10.16418/j.issn.1000-3045.2016.10.005

  • 57

    ZhangB.TangQ. S.JinX. S. (2007). Decadal-scale variations of trophic levels at high trophic levels in the Yellow Sea and the Bohai Sea ecosystem. J. Mar. Syst.67, 304311. 10.1016/j.jmarsys.2006.04.015

  • 58

    ZhaoX. Y.HamreJ.LiF. G.JinX. S.TangQ. S. (2003). Recruitment, sustainable yield and possible ecological consequences of the sharp decline of the anchovy (Engraulis japonicus) stock in the Yellow Sea in the 1990s. Fish. Oceanogr.12, 495501. 10.1046/j.1365-2419.2003.00262.x

  • 59

    ZhengZ. Z. (1965). The structure of zooplankton communities and its seasonal variation in the Yellow Sea and in the western East China Sea [In Chinese, with English Abstract]. Oceanol. Limnol. Sin.7, 199204.

  • 60

    ZuoT.WangR.ChenY. Q.GaoS. W.WangK. (2006). Autumn net copepod abundance and assemblages in relation to water masses on the continental shelf of the Yellow Sea and East China Sea. J. Mar. Syst.59, 159172. 10.1016/j.jmarsys.2005.09.008

Summary

Keywords

zooplankton, community structure, distribution, seasonal change, Yellow sea, ecological function

Citation

Shi Y, Wang J, Zuo T, Shan X, Jin X, Sun J, Yuan W and Pakhomov EA (2020) Seasonal Changes in Zooplankton Community Structure and Distribution Pattern in the Yellow Sea, China. Front. Mar. Sci. 7:391. doi: 10.3389/fmars.2020.00391

Received

03 March 2020

Accepted

07 May 2020

Published

23 June 2020

Volume

7 - 2020

Edited by

Michael Arthur St. John, Technical University of Denmark, Denmark

Reviewed by

Song Sun, Institute of Oceanology (CAS), China; Hongju Chen, Ocean University of China, China

Updates

Copyright

*Correspondence: Jun Wang Tao Zuo

This article was submitted to Marine Ecosystem Ecology, a section of the journal Frontiers in Marine Science

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics