Journal cover Journal topic
Biogeosciences An interactive open-access journal of the European Geosciences Union
Journal topic
Biogeosciences, 15, 2055-2073, 2018
https://doi.org/10.5194/bg-15-2055-2018
Biogeosciences, 15, 2055-2073, 2018
https://doi.org/10.5194/bg-15-2055-2018

Research article 09 Apr 2018

Research article | 09 Apr 2018

# Variability in copepod trophic levels and feeding selectivity based on stable isotope analysis in Gwangyang Bay of the southern coast of the Korean Peninsula

Variability in copepod trophic levels and feeding selectivity
Mianrun Chen1,2, Dongyoung Kim2, Hongbin Liu3, and Chang-Keun Kang2 Mianrun Chen et al.
• 1South China Sea Institute of Planning and Environmental Research, South China Sea Branch, SOA, Guangzhou, P.R. China
• 2School of Earth Sciences & Environmental Engineering, Gwangju Institute of Science and Technology, Gwangju 61005, Republic of Korea
• 3Division of Life Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong SAR, China
Abstract

Trophic preference (i.e., food resources and trophic levels) of different copepod groups was assessed along a salinity gradient in the temperate estuarine Gwangyang Bay of Korea, based on seasonal investigation of taxonomic results in 2015 and stable isotope analysis incorporating multiple linear regression models. The δ13C and δ15N values of copepods in the bay displayed significant spatial heterogeneity as well as seasonal variations, which were indicated by their significant relationships with salinity and temperature, respectively. Both spatial and temporal variations reflected those in isotopic values of food sources. The major calanoid groups (marine calanoids and brackish water calanoids) had a mean trophic level of 2.2 relative to nanoplankton as the basal food source, similar to the bulk copepod assemblage; however, they had dissimilar food sources based on the different δ13C values. Calanoid isotopic values indicated a mixture of different genera including species with high δ15N values (e.g., Labidocera, Sinocalanus, and Tortanus), moderate values (Calanus sinicus, Centropages, Paracalanus, and Acartia), and relatively low δ15N values (Eurytemora pacifica and Pseudodiaptomus). Feeding preferences of different copepods probably explain these seasonal and spatial patterns of the community trophic niche. Bayesian mixing model calculations based on source materials of two size fractions of particulate organic matter (nanoplankton at < 20 µm vs. microplankton at 20–200 µm) indicated that Acartia and Centropages preferred large particles; Paracalanus, Calanus, Eurytemora, and Pseudodiaptomus apparently preferred small particles. Tortanus was typically carnivorous with low selectivity on different copepods. Labidocera preferred marine calanoids Acartia, Centropages, and harpacticoids; on the other hand, Sinocalanus and Corycaeus preferred brackish calanoids Paracalanus and Pseudodiaptomus. Overall, our results depict a simple energy flow of the planktonic food web of Gwangyang Bay: from primary producers (nanoplankton) and a mixture of primary producers and herbivores (microplankton) through omnivores (Acartia, Calanus, Centropages, and Paracalanus) and detritivores (Pseudodiaptomus, Eurytemora, and harpacticoids) to carnivores (Corycaeus, Tortanus, Labidocera, and Sinocalanus).

1 Introduction

Mesozooplankton constitute essential trophic mediators of marine food webs in transferring energy and materials by linking the microbial food web to higher trophic levels. Copepods are a diverse assemblage-dominating mesozooplankton communities. With broad feeding spectra and flexible feeding strategies, the bulk copepod assemblage displays varying degree of herbivory, omnivory, or carnivory, depending on dominant species or groups (Graeve et al., 1994; Sell et al., 2001; Turner, 2004; Vadstein et al., 2004; Gifford et al., 2007; Chen et al., 2017). The role of copepods in planktonic food webs can be determined by their overall trophic levels (TLs) relative to primary producers. In turn, the TL of a diverse copepod assemblage is balanced from different groups with different feeding preferences and is ultimately determined by species composition. Because copepods rely significantly on phytoplankton as prey, the seasonal and spatial changes in the composition and availability of phytoplankton determine the abundance and feeding behavior of the copepod assemblages.

The most dominant copepod species, such as Neocalanus, Calanus, Temora, and Paracalanus, are filter feeders that perform a size-selective feeding behavior depending on particles effectively retained by the feeding appendages of copepods. Large phytoplankton (> 20 µm; mainly diatoms and dinoflagellates) are generally grazed at high rates by copepods, as shown by many field studies in coastal and estuarine waters (e.g., Liu et al., 2005a, b; Chen et al., 2017). Many other field studies have reported that omnivorous species dominate copepod assemblages because of high feeding selectivity on larger microzooplankton that are considered to have higher nutritional quality (e.g., Berk et al., 1977; Fessenden and Cowles, 1994; Calbet and Saiz, 2005; Gifford et al., 2007; Chen et al., 2013). These omnivorous copepods might induce increases in phytoplankton levels indirectly through trophic cascades as they graze intensely on microzooplankton (e.g., ciliates and heterotrophic dinoflagellates) (Nejstgaard et al., 2001; Stibor et al., 2004; Sommer and Sommer, 2006; Zöllner et al., 2009; Chen and Liu, 2011; Chen et al., 2013). Therefore, the assessment of the trophic position (herbivores, omnivores, or carnivores) of copepods within a complex planktonic food web is critical to understanding the ecological relationships between predators and prey.

Stable isotope analysis (SIA) is a reliable technique providing insight into the trophic positions of copepods relative to basal food sources (Grey et al., 2001; Sommer et al., 2005; Hannides et al., 2009; Kürten et al., 2013). Isotopic comparisons with food sources enable us to analyze prey selectivity during predators' feeding history as well as within food web structures (Fry, 2007; Layman et al., 2012). In general, the carbon stable isotope ratio (δ13C) can be useful for tracing food sources because of small fractionation (0.5–1 ‰ per TL) during trophic transfer, particularly when different food sources at a given period in a specific system have distinct δ13C values. By contrast, the nitrogen stable isotope ratio (δ15N) can be useful for estimating relative TLs because δ15N values of consumers generally increase with TL (an average 3.2 ‰ of enrichment per TL; Post, 2002; Michener and Kaufman, 2007). The development of linear mixing models and the Bayesian mixing model has allowed researchers to predict the proportions of different food sources in the diets assimilated by grazers (Phillips and Koch, 2002; Phillips and Gregg, 2003; Moore and Semmens, 2008; Ward et al., 2010; Parnell et al., 2010, 2013).

Coastal and estuarine environments often experience rapid fluctuations of inorganic carbon and nitrogen inputs in response to diverse oceanographic processes (e.g., coastal currents, upwelling, tidal mixing, and river discharges), which drive spatial and seasonal heterogeneities in biogeochemical dynamics and isotopic signatures (Rolff, 2000). Indeed, the δ13C values of suspended particulate organic matter (POM) in estuarine systems increase progressively from the head to the mouth of each estuary because of the lower δ13C values in terrestrial carbon or sewage materials through river discharge (Cifuentes et al., 1988). In contrast, the δ15N values of primary producers increase from being nutrient-sufficient (high fractionation) to nutrient-limiting (low fractionation) and are especially high in anthropogenic wastewater nitrogen inputs (McClelland et al., 1997). In addition, different phytoplankton groups utilize different nitrogen sources with different enrichment factors, possibly offering different isotopic pools to grazers (Gearing et al., 1984; Rolff, 2000; Montoya et al., 2002). For example, diatoms primarily utilize nitrate with varying fractionation factor on 15N (0.7–6.2 ‰) depending on species (Waser et al., 1998; Needoba et al., 2003), while flagellates primarily utilize ammonia with an enrichment factor of 6.5–8 ‰ (Montoya et al., 1991).

Given that isotopic values of copepods vary in association with copepods' food source by one or two increases in TL values, seasonal and spatial patterns generally follow the trends of their food sources or dominant prey (Grey et al., 2001; Montoya et al., 2002; Kürten et al., 2013). Higher δ15N values of copepods caused by fractionation rather than food source or by averaging from mixed food sources are evident considering the lowered isotopic values of fecal pellets (Checkley and Entzeroth, 1985; Checkley and Miller, 1989; Tamelander et al., 2006). Furthermore, the effect of the microbial food web on the elevated δ15N values of copepods cannot be ignored (Rolff, 2000; Kürten et al., 2013). Therefore, variations in isotope signatures of both copepods and POM (including phytoplankton, bacteria, ciliates, and detritus) help to depict the biogeochemical cycles of specific systems (Grey et al., 2001; Montoya et al., 2002; Francis et al., 2011). Nevertheless, because copepods graze preferentially on larger phytoplankton (diatoms and dinoflagellates) and microzooplankton (ciliates and heterotrophic dinoflagellates), we hypothesize that isotopic values of the copepod assemblage will be much closer to those of larger rather than smaller food source plankton.

However, high diversity of the assemblage and size overlap among different species make it hard to determine the relative trophic positions of different subgroups or species. Isotope analysis for different subgroups requires great expertise in isolating species from highly complex mixtures. Moreover, the number of individuals of a specific genus is often insufficient for analysis because of limited instrument sensitivity. Thus, to our knowledge, direct comparisons of different mesozooplankton groups or copepod species are seldom found in the literature (Schmidt et al., 2003; Sommer et al., 2005; Hannides et al., 2009). Here, we estimated isotope values of different copepods by mass balancing linear mixing models from values of bulk samples and taxonomic data of copepods. The allocated masses of calanoids and cyclopoids were achieved from the literature and empirical formulas.

Overall, we aimed (1) to understand the seasonal variations and spatial heterogeneity of copepod δ13C and δ15N values in a temperate estuarine system (Gwangyang Bay, Korea); (2) to compare the trophic positions of different copepods; and finally (3) to elucidate the compositions of two major size classes (< 20 and > 20 µm) of POM in grazer diets. The dietary composition (nano- vs. microplankton) of copepods was estimated using Bayesian isotopic mixing models (Parnell et al., 2010, 2013). The results of this study will provide insights into trophic preference information (i.e., food resources and trophic levels) for different copepod groups and help in understanding the biogeochemistry of this estuarine system.

2 Materials and methods

## 2.1 Study area

Gwangyang Bay is a semi-enclosed bay system, located on the southern coast of the Korean Peninsula, and is one of the most industrialized coastal areas exposed to anthropogenic pressure. It starts from Seomjin River and goes through Yeosu Channel (between Yeosu Peninsula and Namhae Island) to the open ocean (the East China Sea). The bay area covers approximately 145 km2, and water depth is generally shallow at 2.4–8.0 m in the northern upper-middle Seomjin River estuarine channel compared with 10–30 m in the deep-bay channel (Kim et al., 2014). The annual freshwater discharges of Seomjin River are 10.7–39.3 × 108 tons. The seasonality of nutrient input from the catchment area (ca. 5 × 103 km2), including agricultural and forested land, is profound (Kwon et al., 2002). The wet season starts from late spring and the discharge peaks during the summer monsoon period.

Figure 1Map showing the location of Gwangyang Bay (a), the appearance of the bay before the reclamation of tidal flats in 1982 (b), and the sampling stations in the bay (c) and in the estuarine channel (d). The broken line represents the lowest water line (b); the dotted areas show intertidal beds and the dark gray areas Zostera beds (c); and the darker areas show Phragmites beds (d).

Accordingly, the maximum median river discharge varies from 30–95 m3 s−1 in the dry season to 300–400 m3 s−1 in the summer monsoon, with an annual mean of ca. 120 m3 s−1 (Kim et al., 2014). The tidal cycle of the bay is semidiurnal with maximum ranges of 3.40 m during spring tides and 1.10 m during neap tides. Tidal currents from the Yeosu Channel also strongly influence the system, and approximately 82 % of the Seomjin River flux is discharged toward this channel. Overall, increasing industrial pollution facilitates eutrophic conditions in the estuarine and related bay waters. Diatoms dominate in the phytoplankton community, and density is high in the middle part of Gwangyang Bay (Kim et al., 2009; data from our parallel study not shown). The distribution patterns of copepods in the Seomjin River during summer were represented by three main salinity zones: an oligohaline zone (dominated by Pseudodiaptomas koreanus, Sinocalanus tenellus, and Tortanus dextrilobatus), a mesohaline zone (dominated by Acartia ohtsukai and Acartia forticrusa), and a polyhaline zone (dominated by Acartia erythraea, Calanus sinicus, Centropages dorsispinatus, Labidocera rotunda, and Paracalanus parvus) (Park et al., 2015).

## 2.2 Sampling and processing

Surface water and net-tow samples were collected seasonally (February, May, August, and November) at nine stations from the head to the mouth of Gwangyang Bay in 2015 (Fig. 1). The stations were chosen based on salinity regime and different geographic characteristics. Stations 1–3 were located in the Seomjin River, stations 4 and 5 were in Gwangyang Bay (the middle part of the estuary), and stations 6–9 were located from the offshore deep-bay channel to the southern mouth of the estuary. On each sampling occasion, water temperature and salinity were determined in situ using an YSI Model 85 probe (YSI Inc., Yellow Springs, OH, USA).

Zooplankton taxonomic samples were collected during daytime by net towing using a plankton net (45 cm diameter, 200 µm mesh size) equipped with a flowmeter (Model 2030R Mechanical Flowmeter, General Oceanics Inc., Miami, FL) and gently hauled horizontally at a subsurface depth of 0.5–1 m with the ship speed at about 1 knot (0.5 m s−1). The average volume filtered per tow was 16.7 ± 5.1 m3 (mean ± se). Samples were fixed in formalin solution with a final concentration of 5 % and then identified and enumerated under a stereomicroscope (SMZ 645; Nikon, Tokyo, Japan) in the laboratory. At each station, one additional net tow was collected for isotope analysis. After collection, specimens were transferred immediately into plastic bottles and preserved in a refrigerator (4 C) until analysis. In the laboratory, subsamples were picked out from the mixed zooplankton samples under a dissecting microscope. Easily distinguishable zooplankton groups such as harpacticoids were separated from a mixture of calanoids and cyclopoids. Copepodites were counted but grouped together with adults. All subsamples were lyophilized and then homogenized by pulverizing them with a mortar and pestle before isotope analysis.

POM in surface water (0.5–1 m depth) was collected using a 5 L Niskin bottle at a midday high tide at the same time as zooplankton collection. The approximately 20 L seawater collected was first screened through a 200 µm Nitex mesh to remove zooplankton and large-sized particles. The prescreened water samples were transported to the laboratory as soon as possible within 1–2 h. In the laboratory, water samples were filtered again through a 20 µm Nitex mesh and then filtered onto precombusted (450 C for 4 h) Whatman GF/F glass fiber filters to determine isotope ratios of fine POM (< 20 µm) representing pico- and nano-sized plankton. To obtain enough plankton cells for isotope analysis of coarse POM ( 20 µm), we collected POM samples by net towing with a plankton net of 50 cm diameter and 20 µm mesh size. After collection, each sample was prefiltered through a 200 µm Nitex mesh to remove large particles and zooplankton. Both size fractions of samples were prepared in duplicate. Samples for δ13C measurements were acidified by fuming for about 5 h over concentrated HCl in a vacuum desiccator to remove carbonates, while the samples for δ15N measurements were not acidified. All the samples were lyophilized and pulverized with a mortar and pestle before isotope analysis.

For chlorophyll a (Chl a) determination, 1 L subsamples of surface water were filtered through Waterman GF/F glass fiber filters. The filters for Chl a (including other photosynthetic pigments) were extracted with 95 % methanol (5 mL) for 12 h in the dark at 20 C and sonicated for 5 min to foster cell disruption. Aliquots of 1 mL of the supernatants were mixed with 300 µL of water; 100 µL of this solution was analyzed by reverse-phase high-performance liquid chromatography (HPLC, LC-20A HPLC system, Shimadzu Co., Kyoto, Japan) using a Water Symmetry C8 (4.6 × 150 mm, particle size: 3.5 µm, 100 Å pore size) column (Waters, Milford, MA, USA) and a method derived from Zapata et al. (2000). Quantification of standard pigments was calculated by a spectrophotometer with the known specific extinction coefficients after Jeffrey (1997). Sample peaks were identified based on their retention time compared with those of pure standards. Further details on analysis, calibration, and quantification have been given elsewhere (Lee et al., 2011; Kwak et al., 2017).

## 2.3 Isotope analysis

For measurements of carbon and nitrogen stable isotope ratios, all pretreated samples were analyzed using a continuous-flow isotope ratio mass spectrometer (CF-IRMS; Isoprime100, Cheadle, UK) connected to an elemental analyzer (vario Micro cube, Hanau, Germany) following the procedure described by Park et al. (2016). Briefly, powdered samples were sealed in tin combustion cups and filter samples were wrapped with a tin plate. All prepared samples were put into the elemental analyzer to oxidize at high temperature (1030 C). CO2 and N2 gases were introduced into the CF-IRMS with the carrier being helium gas. Data of isotope values are shown in terms of δX, indicating the relative differences between isotope ratios of the sample and conventional standard reference materials (Vienna Pee Dee Belemnite for carbon, and atmospheric N2 for nitrogen), which were calculated by the following equations: $\mathit{\delta }X=\left[\left({R}_{\mathrm{sample}}/{R}_{\mathrm{standard}}\right)-\mathrm{1}\right]×{\mathrm{10}}^{\mathrm{3}}$, where X is 13C or 15N and Rsample and Rstandard are the ratios of heavy to light isotope for samples and standards, respectively. International standards of sucrose (ANU C12H22O11; National Institute of Standards and Technology (NIST), Gaithersburg, MD, USA) for carbon and of ammonium sulfate ([NH4]2SO4; NIST) for nitrogen were used for calibration after analyzing every 5–10 samples. The analytical precision for 20 replicates of urea were approximately  0.1 and  0.05 ‰ for δ13C and δ15N, respectively.

## 2.4 Data analysis

We used multiple linear regression to assign the isotopic value to each species from a mixture sample. As we did not measure the dry weights of each species directly, we calculated them using published empirical equations of the relationship between body length and dry weight of each species, and we searched the literature for the range of body length and their living environment of each species in the “World of Copepods Database” and related references (Table S1 in the Supplement). Those calanoid species primarily living in a marine environment based on the definition of the “World of Copepods Database” and related references were grouped as marine calanoids in this paper and included Acartia hudsonica, Acartia omorii, Bestiolina coreana, Calanus sinicus, Clausocalanus furcatus, Centropages abdominalis, Centropages dorsispinatus, Paracalanus aculeatus, Paraeuchaeta plana, Labidocera rotunda, Labidocera euchaeta, Tortanus dextrilobatus, and Tortanus forcipatus. The group of brackish calanoids included Acartia ohtsukai, Acartia erythraea, Eurytemora Pacifica, Paracalanus parvus, Pseudodiaptomus koreanus, Pseudodiaptomus marinus, and Sinocalanus tenellus.

Assuming that the proportion we picked to do isotope analysis was the same as the proportion in samples for composition analysis, once the weights of different groups or different genera/species were assigned, we computed the isotope ratio for each group by multiple linear regression models as in the following equations:

$\begin{array}{ll}& m×{\mathit{\delta }}^{\mathrm{13}}\mathrm{C}=\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{1}}+\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{2}}+\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{3}}\\ \text{(1)}& & \phantom{\rule{1em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{m}_{X}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{X}+\mathrm{error},& m×{\mathit{\delta }}^{\mathrm{15}}\mathrm{N}=\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{1}}+\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{2}}+\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{3}}\\ \text{(2)}& & \phantom{\rule{1em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{m}_{X}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{X}+\text{error},\end{array}$

where m is the weight of the total community and m1mX are the weight of different groups or genera of each group. δ13C1δ13CX and δ15N1δ15NX are the δ13C and δ15N values of each group or genus, respectively. We used R software to do the estimation using the whole sampling data set. Insignificant results for sparse species, such as Bestiolina coreana, Clausocalanus furcatus, Paraeuchaeta plana, Oithona davisae, and Oncaea venella, found in Gwangyang Bay are not tested separately, while they were incorporated to respective groups based on their living environment.

Given that the isotopic values of consumers come from their diets and thereby from mixed proportions of different sources, the proportions of each source could be simulated by linear mixing models with a fractionation factor (also called a trophic enrichment factor). For instance, a mass balance mixing model is given by

$\begin{array}{ll}& {\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{consumer}}={f}_{\mathrm{1}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}\mathrm{1}}+{f}_{\mathrm{2}}{\mathit{\delta }}^{\mathrm{13}}{\mathrm{C}}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}\mathrm{2}}+\mathrm{\dots }\\ \text{(3)}& & \phantom{\rule{1em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{f}_{n}{\mathit{\delta }}^{\mathrm{13}}{C}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}n}+{\mathit{\alpha }}_{\mathrm{Carbon}},& {\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{consumer}}={f}_{\mathrm{1}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}\mathrm{1}}+{f}_{\mathrm{2}}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}\mathrm{2}}+\mathrm{\dots }\\ \text{(4)}& & \phantom{\rule{1em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{f}_{n}{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{source}\phantom{\rule{0.25em}{0ex}}n}+{\mathit{\alpha }}_{\mathrm{Nitrogen}},\text{(5)}& & {f}_{\mathrm{1}}+{f}_{\mathrm{2}}+\mathrm{\dots }+\phantom{\rule{0.125em}{0ex}}{f}_{n}=\mathrm{1},\end{array}$

where f1fn are the proportion of different sources and αCarbon and αNitrogen are trophic enrichment factors for δ13C and δ15N values, respectively.

Here, a Bayesian isotopic mixing model (available as an open source Stable Isotope Analysis package in R: SIAR) was performed to estimate the relative contribution of nanoplankton (defined by fine POM in the present study) and microplankton (coarse POM) to the copepod diets, as well as copepods to the carnivore diets (Parnell et al., 2010, 2013). The model assumes that each isotopic ratio of consumers follows the pattern of a Gaussian distribution with an unknown mean and standard deviation. The structure of mean values of consumers is a weighted combination of the food sources' isotopic values. The weights make up dietary proportions (given by a Dirichlet prior distribution). The standard deviation is divided up between the uncertainty around the fractionation corrections and the natural variability between all individuals within a defined group (Parnell et al., 2010, 2013). Because the values of consumers calculated from bulk copepod samples using the previous multiple linear regression models were only means and standard errors, we generated a vector consisting of 250 numbers for each group by a random normal distribution function. We then used the default iteration numbers (iterations = 500 000, burin = 50 000) provided by the SIAR package to perform our analysis. Fractionation factors used in the model estimation were estimated by difference of TLs multiplied by 3 ‰ per TL for δ15N and by 0.5 ‰ per TL for δ13C, respectively. TLs were calculated from the δ15N difference between consumer and source as follows (Post, 2002): (TL =$\mathrm{1}+\left({\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{consumer}}-{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{source}}\right)/\mathrm{3}$). Concentrations of isotope per mass among different diets (nanoplankton, microplankton, and major copepod genera) were not considered in this study. Model fitting was done via a Markov Chain Monte Carlo (MCMC) protocol that produces simulations of plausible values of the dietary proportions of each source. More details on model simulation can be found elsewhere (Parnell et al., 2010, 2013).

Table 1Seasonal variations in basic environmental factors, including temperature (T), salinity (S), and chlorophyll a (Chl a), copepod abundance, and dominant species and the percentage (%) of dominant species in total copepods at the nine stations in the Gwangyang Bay system.

All statistical analyses were performed using R 3.4.0 software (https://cran.r-project.org/bin/windows/base/). Regression analyses of copepod isotopic values were performed by generalized additive models (GAMs) using the mgcv library (Wood and Wood, 2015). Data were smoothed by cubic regression splines and fitted by the family of Gaussian functions. One-way analysis of variance (ANOVA) was adopted to test seasonal differences in environmental factors and copepod abundances, and Student's t tests were used to test for significant differences in mean δ13C and δ15N values between nano- and microplankton. Before applying ANOVA and t tests, the data were tested for normality of distribution and equal variance; significance was assumed at P=0.05. For the coefficients of variation, we calculated them by dividing the standard error with mean value.

3 Results

## 3.1 Environmental variability and zooplankton abundances

Environmental factors including temperature, salinity, Chl a levels, copepod abundance, dominant species, and percentages of total copepods are shown in Table 1. Water temperature was significantly higher in summer and lower in winter (ANOVA, P< 0.001). The spatial variability in salinity was significant, with extremely low values at stations 1 and 2 (the river mouth), and then the values gradually increased to station 5 (the middle of the bay). Chl a concentrations ranged from 0.1 to 6.8 µg L−1, and they were significantly higher in spring and summer than in winter and autumn (ANOVA, P< 0.01). The highest Chl a concentration occurred in the middle of the bay during the spring, while the lowest concentration was found at the river station during the autumn. The seasonal variability in copepod abundance was significant (ANOVA, P< 0.01), with higher abundances in winter when temperatures and Chl a concentrations were low.

Figure 2Temporal and spatial variations in plankton δ13C (‰) and δ15N (‰).

The detailed abundance composition of copepods is shown in Table S2. Seasonal and spatial variations in dominant species (> 10 % of total abundance) of copepods were apparent (Table 1). The marine calanoid Acartia dominated at the river mouth to the middle part of the bay, while Paracalanus dominated at the mouth of the bay during winter. Acartia also dominated at the most highly saline stations in summer, except for station 7, where the community was dominated by Labidocera rotunda. A brackish water-preferring calanoid species, Pseudodiaptomus, dominated stations 1 and 2 at the river mouth in spring, and another brackish calanoid species, Sinocalanus, dominated station 1 in autumn. At the river mouth stations in summer, copepods were unexpectedly dominated by the marine calanoid species Tortanus dextrilobatus. The cyclopoid species Corycaeus affinis mainly dominated the most highly saline stations in spring and autumn.

## 3.2 Variability in plankton δ13C and δ15N values

The δ13C values of size-fractionated plankton (< 20 and 20–200 µm) and mixed copepod samples showed distinct spatial variations in each season (Fig. 2a–c). The δ13C values of nanoplankton (< 20 µm POM) ranged from 27.6 to 19.4 ‰, with a mean of 22.7 ‰ (Fig. 2a). The lowest δ13C value of nanoplankton was found at station 1 (the upper stream station of Seomjin River) in spring and the highest at station 9 (the mouth of the estuary) in summer. The δ13C values of microplankton (20–200 µm POM) ranged from 26.3 to 17.8 ‰, with a mean of 20.8 ‰ (Fig. 2b), being significantly higher than those of nanoplankton (paired t test, t=7.6, P< 0.001). Its lowest δ13C value was found at station 2 in spring and the highest at station 8 in winter. Overall, similar to nanoplankton, the microplankton δ13C values were more negative at the river portion (stations 1–3) and less negative at the mouth of the estuary (stations 7–9). The δ13C values of mixed copepods ranged from 25.9 to 16.4 ‰, with a mean of 20.1 ‰ (Fig. 2c). The lowest was found at station 1 in autumn and the highest at station 8 in summer. The spatial variability in copepod δ13C values followed the pattern of POM δ13C values. However, the copepod δ13C values were significantly higher than those of nanoplankton (paired t test, t=8.6, P< 0.001) and microplankton (t=3.1, P=0.004). Their δ13C values were higher in summer and winter than in spring and autumn. At stations 1–3, river input lowered the δ13C values of nanoplankton during the wet season (spring to summer). At stations 4–9, significantly lower δ13C values were observed in autumn than in other seasons (ANOVA, F=13.4, P< 0.001). For copepods, the autumn values were significantly lower than those in other seasons (ANOVA, F=5.9, P=0.004).

Table 2Coefficients (F) and significance levels (p) of the effects of variable predictors on δ15N of nanoplankton (plankton < 20 µm) and microplankton (plankton > 20 µm) using a generalized additive model test. P value < 0.05 indicates significance. The abbreviation “n.s.” stands for “no significance”.

The δ15N values exhibited wider fluctuations than δ13C values (coefficients of variation = 29.3 % vs. 9.0 %, 21.5 % vs. 11.8 %, and 18.8 % vs. 13.1 % for nanoplankton, microplankton, and copepods, respectively). The δ15N values of nanoplankton ranged from 3.2 ‰ (station 4 in summer) to 8.8 ‰ (station 1 in winter), with a mean of 5.6 ‰ (Fig. 2d). There were distinct patterns in the three locations of the bay. The δ15N values tended to decline with distance from the river mouth, then increased in the middle of the bay, and decreased again toward the mouth of the estuary. The nanoplankton δ15N values were higher in winter than in other seasons (paired t test, t=5.4, P=0.001 for spring; t=3.0, P=0.017 for summer; t=4.1, P=0.004 for autumn). As indicated by regression analyses between the distribution of nanoplankton δ15N and environmental factors (Table 2), significant increases in the nanoplankton δ15N values depend on ammonia (GAM, F=4.1, P=0.029) and Chl a (GAM, F=3.8, P=0.044). In addition, the seasonal distribution of nanoplankton δ15N values was well indicated by their relationship with temperature among different seasons (GAM, $F=-\mathrm{5.5}$, P=0.013), with decreasing values in summer to autumn.

Mean δ15N value of microplankton (7.6 ‰), ranging from 4.8 ‰ (station 2 in spring) to 10.2 ‰ (station 6 in spring), was significantly higher than that of nanoplankton (paired t test, t=4.9, P< 0.001). The microplankton δ15N values were higher in summer than in other seasons (ANOVA, F=4.6, P=0.009), with the spatial trend vanishing in summer. Indeed, spatial trends differed between seasons, increasing progressively from the river mouth to the bay mouth in spring and autumn and decreasing in winter. As tested by GAM analysis, the microplankton δ15N values were elevated stepwise by environmental factors including temperature (GAM, F=5.0, P=0.015), salinity (GAM, F=5.0, P=0.031), ammonia (GAM, F=4.5, P=0.031), and nitrate (GAM, F=7.8, P=0.010). Similar to nanoplankton, regression analysis also showed that the microplankton δ15N values increased significantly with increasing Chl a concentrations (GAM, F=4.2, P=0.043).

Copepod δ15N values ranged from 6.6 to 12.3 ‰ and were higher in summer than in other seasons (ANOVA, F=15.6, P< 0.001). Generalized additive model analysis showed that the deviances of copepod δ13C and δ15N values explained by the GAMs were 92.7 and 76.9 %, respectively (Table 3). Copepod δ13C values changed significantly toward increasing salinity (GAM, F=7.9, P=0.005), for both the nanoplankton (GAM, F=6.2, P=0.008) and the microplankton (GAM, F=16.4, P< 0.001; Fig. S1a–c in the Supplement). In contrast, temperature was the most important factor to explain the variability in copepod δ15N values (GAM, F=13.6, P< 0.001; Fig. S1d). The microplankton δ15N value was another important contributor to the variability in copepod δ15N values (GAM, F=3.5, P=0.034; Fig. S1e), while nanoplankton δ15N was not (GAM, P> 0.05). The Chl a concentration influenced the variability in copepod δ15N values significantly (GAM, F=3.3, P=0.047; Fig. S1f).

Figure 3Bi-plots of major plankton groups and genus isotopes in Gwangyang Bay.

Table 3Coefficients (F) and significance levels (p) of the effects of variable predictors on δ13C and δ15N copepods using a generalized additive model test. P value < 0.05 indicates significance. The abbreviation “n.s.” stands for “no significance”.

## 3.3 Trophic positions of major groups

Multiple linear regression analyses to estimate mean isotopic values of different copepod groups (i.e., brackish calanoids, marine calanoids, and cyclopoids) from mixed copepod values (excluding harpacticoids) were all significant (R2=0.94, P< 0.001 for δ13C; R2=0.78, P< 0.001 for δ15N). The intercepts of the model indicating the errors of Eqs. (1) and (2) were 1.1 ± 0.8 and 0.3 ± 0.8 ‰ for δ13C and δ15N, respectively. The two calanoid groups displayed a close mean δ15N value (around 9 ‰) but significantly different δ13C values (26.3 ± 1.5, 20.6 ± 1.6, and 20.1 ± 1.2 ‰ for brackish water calanoids and marine calanoids, respectively; Fig. 3a). Cyclopoids occupied a relatively broad trophic niche based on a big coefficient of variation (41.6 %), and the mean δ13C value of cyclopoids (31.0 ± 6.0 ‰) was even lower than those of brackish calanoids. The δ15N values of the three major groups were all higher than the basal food resource (nanoplankton) and relatively higher than microplankton with some overlap of error bar. The values of harpacticoids isolated from the winter (stations 2–9) and spring samples (stations 1) were measured directly. The mean δ15N values of harpacticoids (6.9 ± 0.6 ‰) were lower than those of other copepods and microplankton, but the δ13C values (16.7 ± 1.4 ‰) of Harpacticoids were relatively higher than other groups.

The multiple linear regression analysis performed for major copepod genera/species was also significant (R2=0.99, P< 0.001 for δ13C; R2=0.88, P< 0.001 for δ15N). The intercepts were 0.4 ± 0.7 and 0.8 ± 0.3 ‰ for δ13C and δ15N, respectively. There were two to three different trophic positions in the mixed copepod assemblages based on different δ15N values (see patterns in Fig. 3b). The mean δ15N values of Corycaeus affinis (13.7 ± 5.6 ‰) were the highest among the taxa, followed by those of Labidocera (12.3 ± 3.9 ‰), Tortanus (12.3 ± 1.9 ‰), Sinocalanus (11.3 ± 1.7 ‰), and Paracalanus (11.0 ± 3.6 ‰). The mean values of Acartia (9.4 ± 2.6 ‰), Calanus Sinicus (9.0 ± 1.1 ‰), Centropages (8.1 ± 1.7 ‰), and Eurytemora pacifica (7.4 ± 2.5 ‰) indicated their trophic positions were relatively lower than those carnivorous species, while they were higher than nanoplankton. Pseudodiaptomus had the lowest value of δ15N (5.3 ± 2.5 ‰), and it was even lower than microplankton and not much different to nanoplankton. Compared to of the two putative food resources (nanoplankton and microplankton), those brackish calanoid genera/species (Pseudodiaptomus, E. pacifica, Paracalanus and Sinocalanus) and Corycaeus had lower δ13C than nanoplankton, while marine calanoid genera/species (Tortanus, Labidocera, Calanus, Acartia and Centropages) had higher values (Fig. 3b).

The trophic level in this study was defined as trophic position relative to nanoplankton, which was considered as the trophic baseline. Considering that the TL of nanoplankton is 1, the TL value of microplankton was calculated to be 0.7 times higher than that of nanoplankton (Fig. 4). As a whole assemblage balanced from different feeding behaviors, as indicated by the standard errors, copepods occupied a 1.2 level higher TL than that of nanoplankton, indicating herbivory (here, herbivory means a trophic level of 2) on nanoplankton with slight omnivory (TL = 2–3) on other dietary sources. The TLs of two major calanoid groups (marine calanoids and brackish calanoids) were similar to the bulk copepod assemblage with mean levels slightly higher than 2. The mean TL value of cyclopoids apparently had a higher trophic level than calanoids. In contrast, the mean TL of harpacticoids was very low, reflected by their low δ15N values.

Figure 4Trophic levels (TLs) of different groups. Nanoplankton were set as TL = 1, while consumers' trophic levels were calculated as TL =$\mathrm{1}+\left({\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{consumer}}-{\mathit{\delta }}^{\mathrm{15}}{\mathrm{N}}_{\mathrm{Nanoplakton}}\right)/\mathrm{3}$. The reference line indicates the herbivores relative to nanoplankton. However, nanoplankton here might not be truly primary producers as the bulk samples might include heterotrophic flagellates, dinoflagellates, and ciliates, which we could not separate from the collected samples.

Figure 5Trophic enrichment (or fractionation factor) from basal food items (nanoplankton: a, c; microplankton: b, d), based on the difference of each sample's δ15N between higher trophic level and lower trophic level; a 0.5 ‰ per one trophic level was used to calculate the δ13C enrichment for each group. Reference lines indicate mean values from all groups.

Among calanoids, the mean TLs of Eurytemora pacifica (1 ± 0.5) and Pseudodiaptomus (1 ± 0.5) are indicative of their herbivorous and/or detrivorous characteristics. The mean TLs of Calanus sinicus (2.0 ± 0.6) and Centropages (1.6 ± 0.6) indicated they were primarily herbivorous; and those of Acartia and Paracalanus were slightly higher than 2, indicating they were primarily omnivorous, feeding on both nanoplankton and microplankton. The levels of Sinocalanus (3.0 ± 0.7), Tortanus, and Labidocera were higher than 3.

Based on TLs, the mean 15N enrichments of the copepod assemblage were estimated to be 3.4 and 1.7 ‰ for nanoplankton and microplankton, respectively (Fig. 5a, b). The enrichment values on nanoplankton for both marine and brackish calanoids, as well as for the genera of Acartia, were close to the average value of the copepod assemblage. The trophic enrichment of Calanus sinicus on nanoplankton (3.0 ± 1.7 ‰) was slightly lower than averaged marine calanoids and total copepods, while the enrichment on microplankton (1.7 ± 1.5 ‰) was similar. Centropages (1.9 ± 1.8 ‰) had much lower enrichment on nanoplankton compared to total copepods and marine calanoids. Four high-TL genera Tortanus, Labidocera, Sinocalanus, and Corycaeus had high enrichments > 6 on nanoplankton and > 3 on microplankton.

Figure 6Comparison of the dietary compositions of major copepod groups: (a) Marine calanoids, including Acartia hudsonica, Acartia omorii, Bestiolina coreana, Calanus sinicus, Clausocalanus furcatus, Centropages abdominalis, Centropages dorsispinatus, Paracalanus aculeatus, Paraeuchaeta plana, Labidocera rotunda, Labidocera euchaeta, Tortanus dextrilobatus, and Tortanus forcipatus; (b) brackish calanoids, including Acartia ohtsukai, Acartia erythraea, Eurytemora Pacifica, Paracalanus parvus, Pseudodiaptomus koreanus, Pseudodiaptomus marinus, and Sinocalanus tenellus; and (c) harpacticoids.

Figure 7Comparison of the dietary compositions of the major omnivorous copepod genera.

The 13C enrichments for total copepods were on average 0.6 and 0.3 ‰ when feeding on nanoplankton and microplankton, respectively (Fig. 5c, d). Patterns were same with enrichments on 15N. The four high-TL genera Tortanus, Labidocera, Sinocalanus, and Corycaeus had high enrichments > 1.2 on nanoplankton and > 0.6 on microplankton. The enrichments of four herbivorous/omnivorous species increased from Centropages, Calanus, and Acartia to Paracalanus. In contrast, the brackish calanoid genera Eurytemora and Pseudodiaptomus had extremely low enrichments of both 15N and 13C.

## 3.4 Contribution of size-fractionated POM to copepod diets

The Bayesian mixing model calculations showed that the contributions of different sizes of POM to copepod diets varied significantly with season (Student's t test, P< 0.001 for all cases; Fig. 7). Size-selective feeding phenomena were particularly apparent in winter (Fig. S2a) and summer (Fig. S2c). Mean contributions of microplankton accounted for about two-thirds of their assimilated diets at all stations in winter and summer and were almost equal to that of nanoplankton in spring and autumn (Fig. S2b, d). The proportions of the two size fractions of POM averaged from all four seasons contributing to copepod diets were also distinctly different at stations 2–5 and station 9 (Fig. S2). The mean contributions of microplankton to the copepod diets increased gradually from the river mouth up to a peak (0.81 ± 0.11) in the middle part of the bay. Then, the proportion declined gradually to a trough (0.31 ± 0.18) at the deep-bay channel. The proportion then rebounded to a high level again at the bay mouth station.

Figure 8Comparison of the dietary composition of major carnivorous genera. Credibility intervals of 95 % (dots), 75 % (whiskers), and 25 % (boxes) and mean values (lines in the boxes) are shown in box plots for each source.

Three major groups of copepods showed contrasting size-selective feeding behaviors (Fig. 6). Marine calanoids typically preferred feeding on larger particles, with a contributing proportion of 0.63 ± 0.03 (range: 0.50–0.71) for microplankton (Fig. 6a). Harpacticoids had a more apparent size-selective feeding behavior, preferring microplankton (0.89 ± 0.10), and merely fed on nanoplankton (extremely low reliance of < 0.11; Fig. 6c). In contrast, brackish calanoids preferred feeding on nanoplankton (0.68 ± 0.05) to on microplankton (0.32 ± 0.05). For herbivorous/omnivorous copepod genera/species (e.g., Acartia, C. sinicus, Centropages, E. pacifica, Paracalanus, and Pseudodiaptomus), we tested their feeding preference on the two size fractions of POM using the Bayesian mixing model (Fig. 7). Acartia and Centropages significantly preferred large to small particles (Student's t test, P< 0.001) with a reliance of > 0.9 on microplankton (Fig. 7a, c); on the other hand, C. sinicus and Paracalanus apparently preferred small particles (Student's t test, P< 0.001) (Fig. 7b, e). The two “low trophic level” species/genera Eurytemora pacifica and Pseudodiaptomus selected nanoplankton more while there were large proportions of both potential food sources in diets (Fig. 7d, f). Based on the TL shown above, we considered that the four genera (Labidocera, Sinocalanus, Tortanus, and Corycaeus) are carnivorous and able to predate on other copepods. Thus, we tested the two fractions of POM and some herbivorous/omnivorous species/genera as potential food sources for them. The dietary compositions estimated by Bayesian mixing models showed a slight difference among them (Fig. 8). Except for Sinocalanus, the other three carnivorous species all showed only negligible reliance on the two size fractions of POM. The marine calanoid genus Labidocera preferred predating on Centropages (0.35 ± 0.08), followed by harpacticoids (0.15 ± 0.10), Acartia (0.13 ± 0.09), and Calanus (0.10 ± 0.07) (Fig. 8a). By contrast, besides feeding on nanoplankton (0.23 ± 0.03), the brackish calanoid genus Sinocalanus also preferred predating on other brackish calanoids Paracalanus (0.37 ± 0.02), Pseudodiaptomus (0.15 ± 0.03), and cyclopoids (0.09 ± 0.02) (Fig. 8b). Tortanus frequently co-occurred with many other copepods, and its diet was composed of many different copepod species/genera without apparent selectivity (Fig. 8c). The cyclopoid species Corycaeus affinis primarily predated on Paracalanus (0.67 ± 0.05) and Pseudodiaptomus (0.22 ± 0.06).

4 Discussion

## 4.1 Variability in δ13C and δ15N values of plankton with time and space

We found that seasonal variations and the spatial heterogeneity of copepod δ13C and δ15N values in Gwangyang Bay followed those of nanoplankton (POM < 20 µm) and microplankton (POM > 20 µm) (Figs. 2 and 3). Based on the results of regression analyses, we found that the variability in copepod isotopic values was influenced by salinity (spatial variations), temperature (temporal variations), and isotopic values of food sources (both spatial and temporal variations; Fig. 3). In general, spatial variations were much more pronounced because of the effect of river input and thereby riverine carbon in different salinity regimes. More negative values of three plankton groups (nanoplankton, microplankton, and copepods) were measured near the river mouth, and then the values increased progressively to the mouth of the estuary, indicating an apparently decreasing effect of river runoff and thus the uptake of carbon derived from river-borne terrestrial organic matter. These results are consistent with other studies in estuarine environments (Cifuentes et al., 1988; Matson and Brinson, 1990; Thornton and McManus, 1994; Deegan and Garritt, 1997; Fry, 2002). Such spatial distribution patterns have also been found for other primary producers such as seagrasses (reviewed by Hemminga and Mateo, 1996), macroalgae (Lee, 2000), as well as benthic microalgae (Kang et al., 2003), and the pattern will further propagate to consumers such as fish (Melville and Connolly, 2003; Herzka, 2005), oysters (Fry, 2002), mollusks (Antonio et al., 2010), and other benthic macroinvertebrates (Choy et al., 2008).

Seasonal successions of δ13C values were also apparent, probably because of high river input in spring, elevated productivity in summer, and species successions in autumn. When the wet season started in spring and phytoplankton started to bloom, river input lowered the δ13C values. The values increased again in summer because of the persistence of phytoplankton bloom (a low fractionation effect because of source limitation) and elevated productivity in summer. Although both river discharge and the input of light carbon were low in autumn, the observed δ13C values were low. This was probably because of the lack of a heavy carbon pool (a dissolved inorganic carbon pool in which the carbon was primarily composed of heavy 13C) due to microbial respiration and species succession (Rau et al., 1990). During the post-bloom period in autumn, phytoplankton show low productivity and Chl a concentrations, with low abundance of diatoms but a dominance of flagellates (Baek et al., 2015). Flagellates are known to have more negative δ13C values than those of diatoms arising from different fractionation effects (Gearing et al., 1984; Cifuentes et al., 1988; Rolff, 2000).

The δ15N variability in the three major plankton groups was relatively complex spatially. The seasonal pattern of nanoplankton δ15N values, with decreasing values in summer to autumn, can be well explained by the relationship between the δ15N values and temperature, indicated by significant GAM analysis. Spatial trends in the nanoplankton δ15N values can be explained by three distinct distribution patterns. The first pattern found in the river mouth area exemplifies a declining trend expected from the mixing of freshwater planktonic materials, which grew up in water with high levels of dissolved inorganic nitrogen. The second pattern found in the middle part of the bay, in which nitrate inputs were low while the concentrations of ammonia increased (Kwon et al., 2004; our own data, not shown), characterizes an increase in δ15N values in association with high Chl a concentrations. Fractionations by autotrophic assimilation and bacterial utilization were the most likely source of the 15N-enriched ammonia in the nutrient pools of the middle of the bay (Cifuentes et al., 1988). The elevated POM δ15N values in the middle of the bay may be explained by 15N-enriched ammonia remaining after algal uptake in the river mouth channel (Sato et al., 2006). The input of sewage-derived 15N-enriched ammonia (domestic sewage and livestock waste) could contribute substantially and overwhelm the other subtle processes to increase δ15N values of nanoplankton. The third distribution pattern represents declining δ15N values toward the offshore bay mouth in association with a reduction in the supply of 15N-enriched nutrients from terrestrial sewage. Furthermore, the fractionation effect of phytoplankton will be reduced when nutrients substantially decrease, and phytoplankton would take up nitrogen with little fractionation and store relatively light nitrogen isotope (Cifuentes et al., 1988; Fogel and Cifuentes, 1993; Granger et al., 2004). As indicated by GAM regression analyses between the distribution of nanoplankton δ15N and environmental factors, significant increases in the nanoplankton δ15N values depend on ammonia and Chl a, further supporting our explanation. The microplankton δ15N values were elevated stepwise by environmental factors including temperature, salinity, ammonia, and nitrate, indicated by significant results of GAM analysis (Table 2). As with nanoplankton, regression analysis also showed that the microplankton δ15N values increased significantly with increasing Chl a concentrations (Table 2). One possible mechanism for this pattern is that higher phytoplankton abundance will result in a 15N-enriched nutrient pool because of fractionation during nutrient assimilation (Kang et al., 2009). Nitrate was important for microplankton, indicative of the role of diatoms (preferring nitrate) in controlling the variation in microplankton δ15N values, whereas nanoflagellates (preferring ammonia) probably controlled the variation in nanoplankton δ15N values.

As indicated by GAM analysis, the seasonality of copepod δ15N values was primarily enhanced by temperature, which probably caused an elevated fractionation effect during the rapid assimilation of copepods. The lower explained deviance of δ15N indicated reflects the food web processes that affect δ15N disproportionally and were not included in the GAM analysis. The regression relationship between larger plankton and copepods was significant, whereas the patterns were somewhat decoupled, as they were primarily observed in spring and autumn. This kind of decoupling has also been reported in the open ocean (Montoya et al., 2002), where the transfer of nitrogen from primary producers to zooplankton is weak. A time lag in zooplankton development might cause the mismatch of zooplankton to 15N-enriched POM at the initial stage of nutrient supplies. Indeed, here we found that the high δ15N values of copepods were primarily observed in summer, while the corresponding δ15N values of POM started to increase from the winter, and phytoplankton blooming occurred in the spring.

## 4.2 Trophodynamics and trophic enrichments of copepods

The variability in copepod isotopic values in Gwangyang Bay suggests that the TLs of the copepod assemblage were highly dynamic. Because of different feeding behaviors and fractionation effects of copepods, the variability in the average community trophic position depends on the overall composition of species and is determined by the dominant species. Direct measurements of copepod isotopic values for species levels have been poorly conducted in the literature, although there are still clear patterns in existing reports. In the Southern Ocean copepods, the known carnivores Euchaeta and Heterorhabdus had high δ15N values while the acknowledged omnivores Calanoides and Metridia were intermediate in position and Rhincalanus had the lowest values (Schmidt et al., 2003). A mesocosm study found that the δ15N values were increasingly higher in the order Temora<Pseudocalanus<Centropages, suggesting an increase of carnivory in the same manner (Sommer et al., 2005). The trophic positions of primary consumers (Oithona and Neocalanus) and secondary consumers (Pleuromamma and Euchaeta) in the North Pacific Subtropical Gyre are estimated to be 2.1 and 2.9, respectively (Hannides et al., 2009). Furthermore, Kürten et al. (2013) reported that the relative trophic positions of zooplankton in the North Sea were high when the assemblage was mainly composed of Sagitta and Calanus but low when the assemblage was dominated by Pseudocalanus and zoea larvae.

Our study has demonstrated the trophodynamics of estuarine copepods by using multiple linear mixing model analysis based on the values of bulk samples and percentages in total biomass, by which the results of estimated δ13C and δ15N values were both significant (P< 0.01). The estimated δ13C values varied greatly by up to 10.9 ‰ among groups (from the lowest for cyclopoids to the highest for harpacticoids) and 14.2 ‰ among genera (from the lowest for Paracalanus to the highest for Centropages) (Fig. 4). The estimated δ15N values also varied somewhat by up to 4.9 ‰ among groups, indicating one TL difference among groups if we consider 3–4 ‰ trophic enrichment of δ15N between two adjacent TLs (Post, 2002; Michener and Kaufman, 2007). The trophic enrichments calculated by the differences from the basal food sources indicate that the overall enrichments of copepods were around 3.4 and 1.7 ‰ from nanoplankton and microplankton, respectively (Fig. 6). Our estimation shows that both marine calanoids and brackish water calanoids overall occupy a similar trophic niche (i.e., similar mean δ15N values) but have contrasting food sources (i.e., different δ13C values; Fig. 4a). The isotopic values of calanoid copepods indicate a mixture of different genera including both high and low δ15N values. For example, marine calanoids were mixed by high δ15N genera (Tortanus and Labidocera) and low δ15N genera (Acartia, Calanus, Centropages) and moderate δ15N genera (Paracalanus and Acartia). Brackish types were mixed by high δ15N genera (Sinocalanus) and low δ15N genera (Eurytemora and Pseudodiaptomus). Consequently, the two major calanoid groups (marine and brackish water types) as well as the mixture of all copepod groups were estimated to be on average one TL higher than the nanoplankton base (Fig. 4). However, we do not necessarily conclude that they are herbivores because the nanoplankton studied here represent POM with a size range of 2–20 µm, which may include ciliates, heterotrophic nanoflagellates, and heterotrophic dinoflagellates. Instead, all assemblages mentioned above might be omnivorous with varying levels of relative trophic positions depending on dominant species.

Among calanoids, brackish water species had significantly lower δ13C values than marine species, indicative of an apparent effect of riverine carbon sources on brackish species through the food web (Fig. 3). These results were consistent with the distribution pattern that most brackish species occurred or were dominant in the upper part of Gwangyang Bay. Sinocalanus tenellus had a contrasting TL value to those of Pseudodiaptomus and Eurytemora (Fig. 4b), indicating a mixture of brackish water calanoids being close to omnivory with a broad feeding size spectrum. However, as Sinocalanus and Pseudodiaptomus were dominant in different seasons (Table 1), the TLs of the copepod assemblage at a specific condition will become relatively carnivorous (Sinocalanus dominating) or omnivorous (Pseudodiaptomus and Paracalanus parvus dominating). P. parvus is an important small species (body length  1 mm) that is widely distributed in coastal and estuarine waters worldwide (Turner, 2004). Our results showed that, similar to other brackish water calanoids, this species was greatly influenced by 13C-depleted dietary sources, dominating both brackish stations in autumn and saline stations in winter (Table 1). This result indicates that Paracalanus was well adapted to fluctuating estuarine environments by feeding on prey originating from freshwater or prey that depends on riverine carbon sources. Acartia, one of the commonest genera in Gwangyang Bay throughout the year (Table 1), included both marine types (A. hudsonica and A. Omorii) and brackish types (A. ohtsukai, and A. erythraea), based on the literature (Table S1). Marine types of Acartia primarily occurred or dominated during the winter and spring, while brackish types dominated during the summer. Overall the genus Acartia had higher δ13C values than those of Paracalanus, suggesting different food sources. Acartia is also a widely distributed genus, with a switching feeding behavior in response to the status of food composition (Kiøboe et al., 1996; Rollwagen-Bollens and Penry, 2003; Chen et al., 2013). The isotopic values of Acartia were similar to those of the assemblage of marine calanoids, indicating that this genus is omnivorous, as typical of marine calanoids. The marine calanoids species C. sinicus had a close trophic niche to total marine calanoids based on similar δ15N value and trophic level (Fig. 4) that was 1 TL higher than nanoplankton, suggesting that this species is a typical marine herbivore relative to nanoplankton. Conversely, two other marine calanoid genera, Tortanus (T. dextrilobatus and T. forcipatus) and Labidocera (L. euchaeta and L. rotunda), were primarily carnivorous as indicated by their δ15N values (Figs. 3b and 4). These estimated results are consistent with the former experimental tests and field investigations (Ambler and Frost, 1974; Landry, 1978; Conley and Turner, 1985; Hooff and Bollens, 2004). Cyclopoids (primarily Corycaeus affinis) dominated copepod assemblages in spring and autumn at the middle part and deep-bay channel of the bay (Table 1). Our data reveal that Corycaeus was primarily carnivorous, being 2 TLs higher than nanoplankton (Fig. 4), and prefers feeding on 13C-depleted dietary sources (Fig. 3). This result is consistent with previous reports that the Corycaeus genus is carnivorous (Gophen and Harris, 1981; Landry et al., 1985; Turner, 1986).

Isotopic values of microplankton indicate that they are roughly half a TL value higher than nanoplankton. Considering that the sizes of most ciliates and heterotrophic dinoflagellates primarily fall within this size spectrum (20–200 µm), this result suggests an omnivorous trend among the mixed microplankton groups. Similarly, although measured only in winter samples, the benthic copepod group harpacticoids, represented by the species Euterpina sp., also differed from calanoids and cyclopoids with low δ15N values. The TL of harpacticoids estimated from this approach was somewhat misleading because of their unexpectedly low δ15N values, which probably reflect feeding on detritus or dead organisms that are depleted in 15N (Sautour and Castel, 1993).

## 4.3 Selective feeding of copepods

Feeding preferences of different groups or genera on two size fractions of POM are of particular importance to explain seasonal and spatial patterns of community trophic niches and in turn will predict the impacts of the grazer community on lower TLs including phytoplankton and microzooplankton. Because not all possible food sources, such as bacteria, picoplankton, fecal pellets, and dead detritus, were investigated, our Bayesian mixing model calculations might have led to some biased results. Nevertheless, the model results might provide an estimation of what size fractions of dietary sources the grazer community ingest and assimilate. In general, our results highlight that the copepod assemblages have size-selective feeding behaviors and that these vary with season (Fig. S2) and space (Fig. S3). The feeding selectivity of the bulk copepod assemblage was a balance of ingestion among different groups. The whole copepod assemblage assimilated two-thirds of its food requirement from microplankton in winter and summer, but it fed nearly equally on both size fractions of POM in spring and autumn (Fig. S2). Our results suggest that groups that preferred large-sized prey played a more important role in the total assemblage of copepods in winter and summer, during which the assemblages were primarily dominated by marine calanoids (Table 1). On the other hand, in the spring and autumn when the assemblage was primarily dominated by carnivorous species (Corycaeus, Tortanus) and dominated by both Paracalanus (our result suggested this genus was an omnivorous species and the size selectivity was less pronounced preferring small-sized particle) and Acartia (preferring large-sized particle), the assemblage overall did not show an apparent size selectivity.

Based on the model results for major copepod groups and genera, marine calanoids preferred feeding on large particles (Figs. 6a and 7a) Such a size-selective feeding preference has been widely reported in many field investigations (Liu et al., 2005a; Jang et al., 2010; Chen et al., 2017). By contrast, brackish calanoids as a group preferred feeding on small-sized plankton. These results were not well reported in the literature. Among calanoids, Acartia and Paracalanus both contain marine species and brackish species.Our model results showed that Acartia preferred feeding on large-sized plankton, while Paracalanus was contrastingly different. Acartia, well studied in the literature, is reported to prefer feeding on phytoplankton larger than 20 µm in the coastal water adjacent to the present study area (Jang et al., 2010). Besides, due to nutritional requirements, Acartia was reported to prefer predating on microzooplankton such as ciliates and heterotrophic dinoflagellates (Wiadnyana and Rassoulzadegan, 1989; Chen et al., 2013). On the other hand, another dominant calanoid genus, Paracalanus, preferred feeding on nanoplankton (Fig. 7e), suggesting that the feeding efficiency of the grazers' feeding appendage to retain particles differs between these two genera. Because of such different feeding preferences, they dominated in different stations or seasons with differing food conditions (Table 1) and frequently co-occurred with little overlap of preferred food particle sizes. Besides, a low δ13C value of P. parvus indicated that this wildly distributed species preferred feeding on particles from terrestrial carbon.

Although C. sinicus was not dominant in abundance, this species contributes a great amount in biomass due to its large body size (Table S1). Similar to Paracalanus, C. Sinicus preferred feeding on small-sized particles, while the difference is that the relatively high δ13C value of C. Sinicus indicated that this species preferred prey originating from offshore areas. The result is consistent with reports that C. Sinicus is a dominant species in the Yellow Sea and the East China Sea (the offshore shelf waters of this study estuary) (Uye, 2000; Liu et al., 2003). Another marine calanoid genus similar to Acartia, Centropages (C. abdominalis and dorsispinatus), apparently preferred feeding on large-sized plankton (this study and Wiadnyana and Rassoulzadegan, 1989). For the other two marine calanoid genera – Labidocera (L. rotunda and L. euchaeta) and Tortanus (T. dextrilobatus and T. forcipatus) – feeding on two size fractions of POM did not occur, based on no contribution of the two size fractions of POM to their diets (Fig. 8a, c). Model results indicated that Labidocera and Tortanus were carnivorous genera, consistent with many reports (Landry, 1978; Mullin, 1979; Turner, 1984; Uye and Kayano, 1994; Hooff and Bollens, 2004). Our result also showed that Labidocera had an apparent selectivity in predating on their prey, preferring marine calanoids (Acartia, Calanus and Centropages) and harpacticoids. This feeding selectivity was consistent with the distribution of this species. It was never distributed at river stations (stations 1–3) (Table S2) and even dominated at station 7 during the summer (Table 1). By contrast, Tortanus did not show selectivity among different prey (Fig. 8c), and thus this genus occurred frequently in different stations (Table S2).

The model results indicated that the size selectivity of brackish water calanoids such as Pseudodiaptomus and Eurytemora was also apparently similar to Paracalanus for nanoplankton (Fig. 7d, f). To our knowledge, the feeding habit of Pseudodiaptomus is unknown in the current literature, whereas some field studies suggest that estuarine Pseudodiaptomus flourishes by feeding on small phytoplankton cells (< 20 µm) (Froneman, 2004), consistent with the present results. Incorporating the results of low δ13C values of Pseudodiaptomus and Eurytemora, we believe that these two genera were able to feed on that prey with small sizes (2–20 µm) and low δ13C from terrestrial sources, whereas a low trophic position and trophic enrichments indicated that these two genera may ingest prey with extremely low δ15N compared to nanoplankton (e.g., detritus). The detritus with low δ15N may contribute to the balance between the δ15N in Pseudodiaptomus and in Eurytemora, and, thus, these may act as detritivores as well as herbivores. Another brackish calanoid species – Sinocalanus tenellus – was more diverse in feeding selectivity (Fig. 8b). It can act as a suspension feeder, preferring small-sized plankton, as well as a raptorial feeder, preferring to predate on other brackish calanoids and cyclopoids. This species was even reported as a cannibalistic feeder in an estuary (Hada and Uye, 1991).

Figure 9A simplified energy flow figure of Gwangyang Bay plankton and copepods. Arrows indicate feeding relationships. Thick arrows indicate feeding preference.

In contrast to both marine calanoids and brackish water calanoids, the cyclopoid species Corycaeus affinis was a carnivorous species preferring to predate on two brackish species – Paracalanus and Pseudodiaptomus (Fig. 8d). This species is widely distributed in all the world's oceans (Turner, 2004) and also frequently dominates in copepod assemblages in spring and autumn in Gwangyang Bay (Table 1). Although the isotopic data for harpacticoid copepods were limited to only one season, we still obtained a clear feeding selectivity pattern based on the Bayesian mixing model (Fig. 7c). Harpacticoids in winter preferred microplankton to nanoplankton if benthic food sources were not considered. Their feeding selectivity contributed to the overall feeding preference of total copepods in this season.

5 Conclusions

Here we have demonstrated the temporal and spatial variability in stable isotope ratios of copepods, which was determined by the isotopic values of two size fractions of POM and strongly influenced by salinity (spatiality) and temperature (temporality). Such a characteristic is a key in understanding the biogeochemical cycles of carbon and nitrogen in Gwangyang Bay. We further used a simple linear mixing model and a Bayesian mixing model to extrapolate from the information derived from the isotopic analysis of bulk copepod samples. The model results were robust and allowed the estimation of the relative TLs and trophic enrichment (fractionation effect) of different groups and dominant genera of copepods as well as their diet compositions. Temporal and spatial patterns of copepod isotopic traits were further explained by size selectivity on plankton size fractions, as well as the feeding preference of dominant species. Based on such relative trophic positions and feeding preference, we can depict a simple energy flow diagram of the Gwangyang Bay planktonic food web: from primary producers (nanoplankton) and a mixture of primary producers and herbivores (microplankton) through omnivores (represented by Calanus, Centropages, Acartia, and Paracalanus) and detrivores (represented by Pseudodiaptomus and Eurytemora) to carnivores (represented by Corycaeus, Tortanus, Labidocera, and Sinocalanus) (see a simplified energy flow diagram, Fig. 9).

Data availability
Data availability.

Data are available upon request.

Supplement
Supplement.

Author contributions
Author contributions.

MC and CKK designed the experiments. DK and CKK conducted field work, sample collection, and laboratory analyses. MC, HL, and CKK developed the growth model and performed the simulations. MC, HL, and CKK prepared the manuscript with contributions from all coauthors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This research was supported by “Long-term change of structure and function in marine ecosystems of Korea” funded by the Ministry of Oceans and Fisheries, Korea. Mianrun Chen is also supported by the China Scholarship Council (CSC: 201604180018) and the National Natural Science Foundation of China (NSFC: 41306168). We kindly thank Young-Jae Lee, Hee-Yoon Kang, Changseong Kim, and Jenny Seong for sampling assistance and sample analysis.

Edited by: Emilio Marañón
Reviewed by: four anonymous referees

References

Ambler, J. W. and Frost, B. W.: The feeding behavior of a predatory planktonic copepod, Torlanus discaudatus, Limnol. Oceanogr., 19, 446–451, 1974.

Antonio, E. S., Kasai, A., Ueno, M., Kurikawa, Y., Tsuchiya, K., Toyohara, H., Inshihi, Y., Yokoyama, H., and Yamashita, Y.: Consumption of terrestrial organic matter by estuarine molluscs determined by analysis of their stable isotopes and cellulase activity, Estuar. Coast. Shelf S., 86, 401–407, 2010.

Baek, S. H., Kim, D., Son, M., Yun, S. M., and Kim, Y. O.: Seasonal distribution of phytoplankton assemblages and nutrient-enriched bioassays as indicators of nutrient limitation of phytoplankton growth in Gwangyang Bay, Korea, Estuar. Coast. Shelf S., 163, 265–278, 2015.

Berk, S. G., Brownlee, D. C., Heinle, D. R., Kling, H. J., and Colwell, R. R.: Ciliates as a food source for marine planktonic copepods, Microb. Ecol., 4, 27–40, 1977.

Calbet, A. and Saiz, E.: The ciliate-copepod link in marine ecosystems, Aquat. Microb. Ecol., 38, 157–167, 2005.

Checkley, D. M. and Entzeroth, L. C.: Elemental and isotopic fractionation of carbon and nitrogen by marine, planktonic copepods and implications to the marine nitrogen cycle, J. Plankton Res., 7, 553–568, 1985.

Checkley, D. M. and Miller, C. A.: Nitrogen isotope fractionation by oceanic zooplankton, Deep-Sea Res. Pt. I, 36, 1449–1456, 1989.

Chen, M. and Liu, H.: Experimental simulation of trophic interactions among omnivorous copepods, heterotrophic dinoflagellates and diatoms, J. Exp. Mar. Biol. Ecol., 403, 65–74, 2011.

Chen, M., Liu, H., and Li, H.: Effect of mesozooplankton feeding selectivity on the dynamics of algae in the presence of intermediate grazers a laboratory simulation, Mar. Ecol.-Prog. Ser., 486, 47–58, 2013.

Chen, M., Liu, H., and Chen, B.: Seasonal variability of mesozooplankton feeding rates on phytoplankton in subtropical coastal and estuarine waters, Front. Mar. Sci., 4, 186, https://doi.org/10.3389/fmars.2017.00186, 2017.

Choy, E. J., An, S., and Kang, C. K.: Pathways of organic matter through food webs of diverse habitats in the regulated Nakdong River estuary (Korea), Estuar. Coast. Shelf S., 78, 215–226, 2008.

Cifuentes, L. A., Sharp, J. H., and Fogel, M. L.: Stable carbon and nitrogen isotope biogeochemistry in the Delaware estuary, Limnol. Oceanogr., 33, 1102–1115, 1988.

Conley, W. J. and Turner, J. T.: Omnivory by the coastal marine copepods Centropages hamatus and Labidocera aestiva, Mar. Ecol.-Prog. Ser., 241, 113–120, 1985.

Deegan, L. A. and Garritt, R. H.: Evidence for spatial variability in estuarine food webs, Mar. Ecol.-Prog. Ser., 147, 31–47, 1997.

Fessenden, L. and Cowles, T. J.: Copepod predation on phagotrophic ciliates in Oregon coastal waters, Mar. Ecol.-Prog. Ser., 107, 103–111, 1994.

Fogel, M. L. and Cifuentes, L. A.: Isotope fractionation during primary production, in: Organic Geochemistry, edited by: Engel, M. H. and Macko, S. A., Plenum Press, New York, 73–98, 1993.

Francis, T. B., Schindler, D. E., Holtgrieve, G. W., Larson, E. R., Scheuerell, M. D., Semmens, B. X., and Ward, E. J.: Habitat structure determines resource use by zooplankton in temperate lakes, Ecol. Lett., 14, 364–372, 2011.

Froneman, P. W.: In situ feeding rates of the copepods, Pseudodiaptomus hessei and Acartia longipatella, in a temperate, temporarily open/closed Eastern Cape estuary, S. Afr. J. Sci., 100, 577–583, 2004.

Fry, B.: Conservative mixing of stable isotopes across estuarine salinity gradients: a conceptual framework for monitoring watershed influences on downstream fisheries production, Estuar. Coast., 25, 264–271, 2002.

Fry, B.: Stable Isotope Ecology, Springer Science & Business Media, New York, USA, 2007.

Gearing, J. N., Gearing, P. J., Rudnick, D. T., Requejo, A. G., and Hutchins, M. J.: Isotopic variability of organic carbon in a phytoplankton-based, temperate estuary, Geochim. Cosmochim. Ac., 48, 1089–1098, 1984.

Gifford, S. M., Rollwagen-Bollens, G., and Bollens, S. M.: Mesozooplankton omnivory in the upper San Francisco Estuary, Mar. Ecol.-Prog. Ser., 348, 33–46, 2007.

Gophen, M. and Harris, R. P.: Visual predation by a marine cyclopoid copepod, Corycaeus anglicus, J. Mar. Biol. Assoc. UK, 61, 391–399, 1981.

Graeve, M., Hagen, W., and Kattner, G.: Herbivorous or omnivorous? On the significance of lipid compositions as trophic markers in Antarctic copepods, Deep-Sea Res. Pt. I, 41, 915–924, 1994.

Granger, J., Sigman, D. M., Needoba, J. A., and Harrison, P. J.: Coupled nitrogen and oxygen isotope fractionation of nitrate during assimilation by cultures of marine phytoplankton, Limnol. Oceanogr., 49, 1763–1773, 2004.

Grey, J., Jones, R. I., and Sleep, D.: Seasonal changes in the importance of the source of organic matter to the diet of zooplankton in Loch Ness, as indicated by stable isotope analysis, Limnol. Oceanogr., 46, 505–513, 2001.

Hada, A. and Uye, S.: Cannibalistic feeding behavior of the brackish-water copepod Sinocalanus tenellus, J. Plankton. Res., 13, 155–166, 1991.

Hannides, C., Popp, B. N., Landry, M. R., and Graham, B. S.: Quantification of zooplankton trophic position in the North Pacific Subtropical Gyre using stable nitrogen isotopes, Limnol. Oceanogr., 54, 50–61, 2009.

Hemminga, M. A. and Mateo, M. A.: Stable carbon isotopes in seagrasses: variability in ratios and use in ecological studies, Mar. Ecol.-Prog. Ser., 140, 285–298, 1996.

Herzka, S. Z.: Assessing connectivity of estuarine fishes based on stable isotope ratio analysis, Estuar. Coast. Shelf S., 64, 58–69, 2005.

Hooff, R. C. and Bollens, S. M.: Functional response and potential predatory impact of Tortanus dextrilobatus, a carnivorous copepod recently introduced to the San Francisco Estuary, Mar. Ecol.-Prog. Ser., 277, 167–179, 2004.

Jang, M.-C., Shin, K., Lee, T., and Noh, I.: Feeding selectivity of calanoid copepods on phytoplankton in Jangmok Bay, South Coast of Korea, Ocean Sci. J., 45, 101–111, 2010.

Jeffrey, S. W.: Application of pigment methods to oceanography, in: Phytoplankton Pigments in Oceanography: Guidelines to Modern Methods, edited by: Jeffery, S. W., Mantoura, R. F. C., and Wright, S. W., UNESCO Publishing, Paris, France, 127–166, 1997.

Kang, C. K., Kim, J. B., Lee, K. S., Kim, J. B., Lee, P. Y., and Hong, J. S.: Trophic importance of benthic microalgae to macrozoobenthos in coastal bay systems in Korea: dual stable C and N isotope analyses, Mar. Ecol.-Prog. Ser., 259, 79–92, 2003.

Kang, C. K., Choy, E. J., Hur, Y. B., and Myeong, J. I.: Isotopic evidence of particle size-dependent food partitioning in cocultured sea squirt Halocynthia roretzi and Pacific oyster Crassostrea gigas, Aquat. Biol., 6, 289–302, 2009.

Kim, B. J., Ro, Y. J., Jung, K. Y., and Park, K. S.: Numerical modeling of circulation characteristics in the Kwangyang Bay estuarine system, J. Korean Soc. Coast. Ocean Eng., 26, 253–266, 2014.

Kim, S. Y., Moon, C. H., Cho, H. J., and Lim, D. I.: Dinoflagellate cysts in coastal sediments as indicators of eutrophication: a case of Gwangyang Bay, South Sea of Korea, Estuar. Coast., 32, 1225–233, 2009.

Kiørboe, T., Saiz, E., and Viitasalo, M.: Prey switching behaviour in the planktonic copepod Acartia tonsa, Mar. Ecol.-Prog. Ser., 143, 65–75, 1996.

Kürten, B., Painting, S. J., Struck, U., Polunin, N. V., and Middelburg, J. J.: Tracking seasonal changes in North Sea zooplankton trophic dynamics using stable isotopes, Biogeochemistry, 113, 1–21, 2013.

Kwak, J. H., Han, E., Lee, S. H., Park, H. J., Kim, K. R., and Kang, C. K.: A consistent structure of phytoplankton communities across the warm–cold regions of the water mass on a meridional transect in the East/Japan Sea, Deep-Sea Res. Pt. II, 143, 36–44, 2017.

Kwon, K. Y., Moon, C. H., Kang, C. K., and Kim, Y. N.: Distribution of particulate organic matters along the salinity gradients in the Seomjin River estuary, Korean J. Fish. Aquat. Sci., 35, 86–96, 2002.

Kwon, K. Y., Moon, C. H., Lee, J. S., Yang, S. R., Park, M. O., and Lee, P. Y.: Estuarine behavior and flux of nutrients in the Seomjin River estuary, The Sea J. Korean Soc. Oceanogr., 9, 153–163, 2004.

Landry, M. R.: Predatory feeding behavior of a marine copepod, Labidocera trispinosa, Limnol. Oceanogr., 23, 1103–1113, 1978.

Landry, M. R., Lehner-Fournier, J. M., and Fagerness, V. L.: Predatory feeding behavior of the marine cyclopoid copepod Corycaeus anglicus, Mar. Biol., 85, 163–169, 1985.

Layman, C. A., Araujo, M. S., Boucek, R., Hammerschlag-Peyer, C. M., Harrison, E., Jud, Z. R., and Post, D. M.: Applying stable isotopes to examine food-web structure: an overview of analytical tools, Biol. Rev., 87, 545–562, 2012.

Lee, S. Y.: Carbon dynamics of Deep Bay, eastern Pearl River estuary, China. II: Trophic relationship based on carbon-and nitrogen-stable isotopes, Mar. Ecol.-Prog. Ser., 205, 1–10, 2000.

Lee, Y. W., Park, M. O., Kim, Y. S., Kim, S. S., and Kang, C. K.: Application of photosynthetic pigment analysis using a HPLC and CHEMTAX program to studies of/ phytoplankton community composition, The Sea J. Korean Soc. Oceanogr., 16, 117–124, 2011.

Liu, G. M., Sun, J., Wang, H., Zhang, Y., Yang, B., and Ji, B.: Abundance of Calanus sinicus across the tidal front in the Yellow Sea, China, Fish. Oceanogr., 12, 291–298, 2003.

Liu, H., Dagg, M. J., and Strom, S.: Grazing by the calanoid copepod Neocalanus cristatus on the microbial food web in the coastal Gulf of Alaska, J. Plankton Res., 27, 647–662, 2005a.

Liu, H., Dagg, M., Wu, C., and Chiang, K.: Mesozooplankton consumption of microplankton in the Mississippi River plume, with special emphasis on planktonic ciliates, Mar. Ecol.-Prog. Ser., 286, 133–144, 2005b.

Matson, E. A. and Brinson, M. M.: Stable carbon isotopes and the C : N ratio in the estuaries of the Pamlico and Neuse Rivers. North Carolina, Limnol. Oceanogr., 35, 1290–1300, 1990.

McClelland, J. W., Valiela, I., Michener, R. H.: Nitrogen-stable isotope signatures in estuarine food webs: A record of increasing urbanization in coastal watersheds, Limnol. Oceanogr., 42, 930–937, 1997.

Melville, A. J. and Connolly, R. M.: Spatial analysis of stable isotope data to determine primary sources of nutrition for fish, Oecologia, 136, 499–507, 2003.

Michener, R. H. and Kaufman, L.: Stable isotope ratios as tracers in marine food webs: an update, Stable Isotopes in Ecology and Environmental Science, 2, 238–282, 2007.

Montoya, J. P., Horrigan, S. G., and McCarthy, J. J.: Rapid, storm-induced changes in the natural abundance of 15N in a planktonic ecosystem, Chesapeake Bay, USA, Geochim. Cosmochim. Ac., 55, 3627–3638, 1991.

Montoya, J. P., Carpenter, E. J., and Capone, D. G.: Nitrogen fixation and nitrogen isotope abundances in zooplankton of the oligotrophic North Atlantic, Limnol. Oceanogr., 47, 1617–1628, 2002.

Moore, J. W. and Semmens, B. X.: Incorporating uncertainty and prior information into stable isotope mixing models, Ecol. Lett., 11, 470–480, 2008.

Mullin, M. M.: Differential predation by the carnivorous marine copepod Tortanus discautus, Limnol. Oceanogr., 24, 774–777, 1979.

Needoba, J. A., Waser, N. A., Harrison, P. J., and Calvert, S. E.: Nitrogen isotope fractionation in 12 species of marine phytoplankton during growth on nitrate, Mar. Ecol.-Prog. Ser., 255, 81–91, 2003.

Nejstgaard, J. C., Naustvoll, L. J., and Sazhin, A.: Correcting for underestimation of microzooplankton grazing in bottle incubation experiments with mesozooplankton, Mar. Ecol.-Prog. Ser., 221, 59–75, 2001.

Park, E. O., Suh, H. L., and Soh, H. Y.: Spatio-temporal distribution of Acartia (Copepoda: Calanoida) species along a salinity gradient in the Seomjin River estuary, South Korea, J. Nat. Hist., 49, 2799–2812, 2015.

Park, H. J., Han, E., Lee, Y. J., and Kang, C. K.: Trophic linkage of a temperate intertidal macrobenthic food web under opportunistic macroalgal blooms: A stable isotope approach, Mar. Pollut. Bull., 111, 86–94, 2016.

Parnell, A. C., Inger, R., Bearhop, S., and Jackson, A. L.: Source partitioning using stable isotopes: coping with too much variation, PloS One, 5, e9672, https://doi.org/10.1371/journal.pone.0009672, 2010.

Parnell, A. C., Phillips, D. L., Bearhop, S., Semmens, B. X., Ward, E. J., Moore, J. W., Jackson, A. L., Grey, J., Kelly, D. J., and Inger, R.: Bayesian stable isotope mixing models, Environmetrics, 24, 387–399, 2013.

Phillips, D. L. and Gregg, J. W.: Source partitioning using stable isotopes: coping with too many sources, Oecologia, 136, 261–269, 2003.

Phillips, D. L. and Koch, P. L.: Incorporating concentration dependence in stable isotope mixing models, Oecologia, 130, 114–125, 2002.

Post, D. M.: Using stable isotopes to estimate trophic position: models, methods, and assumptions, Ecology, 83, 703–718, 2002.

Rau, G. H., Teyssie, J. L., Rassoulzadegan, F., and Fowler, S. W.: 13C/12C and 15N/14N variations among size-fractionated marine particles: implications for their origin and trophic relationships, Mar. Ecol.-Prog. Ser., 59, 33–38, 1990.

Rolff, C.: Seasonal variation in δ13C and δ15N of size-fractionated plankton at a coastal station in the northern Baltic proper, Mar. Ecol.-Prog. Ser., 203, 47–65, 2000.

Rollwagen-Bollens, G. C. and Penry, D. L.: Feeding dynamics of Acartia spp. copepods in a large, temperate estuary (San Francisco Bay, CA), Mar. Ecol.-Prog. Ser., 257, 139–158, 2003.

Sato, T., Miyajima, T., Ogawa, H., Umezawa, Y., and Koike, I.: Temporal variability of stable carbon and nitrogen isotopic composition of size-fractionated particulate organic matter in the hypertrophic Sumida River Estuary of Tokyo Bay, Japan, Estuar. Coast. Shelf S., 68, 245–258, 2006.

Sautour, B. and Castel, J.: Feeding behaviour of the coastal copepod Euterpina acutifrons on small particles, Cah. Biol. Mar., 34, 239–251, 1993.

Schmidt, K., Atkinson, A., Stübing, D., McClelland, J. W., Montoya, J. P., and Voss, M.: Trophic relationships among Southern Ocean copepods and krill: some uses and limitations of a stable isotope approach, Limnol. Oceanogr., 48, 277–289, 2003.

Sell, A. F., van Keuren, D., and Madin, L. P.: Predation by omnivorous copepods on early developmental stages of Calanus finmarchicus and Pseudocalanus spp., Limnol. Oceanogr., 46, 953–959, 2001.

Sommer, F., Saage, A., Santer, B., Hansen, T., and Sommer, U.: Linking foraging strategies of marine calanoid copepods to patterns of nitrogen stable isotope signatures in a mesocosm study. Mar. Ecol.-Prog. Ser., 286, 99–106, 2005.

Sommer, U. and Sommer, F.: Cladocerans versus copepods: the cause of contrasting top–down controls on freshwater and marine phytoplankton, Oecologia, 147, 183–194, 2006.

Stibor, H., Vadstein, O., Diehl, S., Gelzleichter, A., Hansen, T., Hantzsche, F., Katechakis, A., Lippert, B., Løseth, K., Peters, C., Roederer, W., Sandow, M., Sundt-Hansen, L., and Olsen, Y.: Copepods act as a switch between alternative trophic cascades in marine pelagic food webs, Ecol. Lett., 7, 321–328, 2004.

Tamelander, T., Søreide, J. E., Hop, H., and Carroll, M. L.: Fractionation of stable isotopes in the Arctic marine copepod Calanus glacialis: effects on the isotopic composition of marine particulate organic matter, J. Exp. Mar. Biol. Ecol., 333, 231–240, 2006.

Thornton, S. F. and McManus, J.: Application of organic carbon and nitrogen stable isotope and C  N ratios as source indicators of organic matter provenance in estuarine systems: evidence from the Tay Estuary, Scotland, Estuar. Coast. Shelf S., 38, 219–233, 1994.

Turner, J. T.: Zooplankton feeding ecology: contents of fecal pellets of the copepods Acartia tonsa and Labidocera aestiva from continental shelf waters near the mouth of the Mississippi River, Mar. Ecol., 5, 265–282, 1984.

Turner, J. T.: Zooplankton feeding ecology: contents of fecal pellets of the cyclopoid copepods Oncaea venusta, Corycaeus amazonicus, Oithona plumifera, and O. simplex from the northern Gulf of Mexico, Mar. Ecol., 7, 289–302, 1986.

Turner, J. T.: The importance of small planktonic copepods and their roles in pelagic marine food webs, Zool. Stud., 43, 255–266, 2004.

Uye, S.: Why does Calanus sinicus prosper in the shelf ecosystem of the Northwest Pacific Ocean?, ICES J. Mar. Sci, 57, 1850–1855, 2000.

Uye, S. and Kayano, Y.: Predatory feeding behavior of Tortanus (Copepoda: Calanoida): life-stage differences and the predation impact on small planktonic crustaceans, J. Crustacean Biol., 14, 473–483, 1994.

Vadstein, O., Stibor, H., Lippert, B., Løseth, K., Roederer, W., Sundt-Hansen, L., and Olsen, Y.: Moderate increase in the biomass of omnivorous copepods may ease grazing control of planktonic algae, Mar. Ecol.-Prog. Ser., 270, 199–207, 2004.

Ward, E. J., Semmens, B. X., and Schindler, D. E.: Including source uncertainty and prior information in the analysis of stable isotope mixing models, Environ. Sci. Technol., 44, 4645–4650, 2010.

Waser, N. A. D., Harrison, P. J., Nielsen, B., Calvert, S. E., and Turpin, D. H.: Nitrogen isotope fractionation during the uptake and assimilation of nitrate, nitrite, ammonium, and urea by a marine diatom, Limnol. Oceanogr., 43, 215–224, 1998.

Wiadnyana, N. N. and Rassoulzadegan, F.: Selective feeding of Acartia clause and Centropages typicus on microzooplankton, Mar. Ecol.-Prog. Ser., 53, 37–45, 1989.

Wood, S. and Wood, M. S.: Package “mgcv”, R Package Version 1.7–29, 2015.

Zapata, M., Rodríguez, F., and Garrido, J. L.: Separation of chlorophylls and carotenoids from marine phytoplankton: a new HPLC method using a reversed phase C8 column and pyridine-containing mobile phases, Mar. Ecol.-Prog. Ser., 195, 29–45, 2000.

Zöllner, E., Hoppe, H. G., Sommer, U., and Jürgens, K.: Effect of zooplankton-mediated trophic cascades on marine microbial food web components (bacteria, nanoflagellates, ciliates), Limnol. Oceanogr., 54, 262–275, 2009.