Description: The Humboldt Current Upwelling System (HCS) is the most productive eastern boundary upwelling system (EBUS) in terms of fishery yield on the planet. EBUSs are considered hotspots of climate change with predicted expansion of mesopelagic oxygen minimum zones (OMZs) and related changes in the frequency and intensity of upwelling of nutrient-rich, low-oxygen deep water. To increase our mechanistic understanding of how upwelling impacts plankton communities and trophic links, we investigated mesozooplankton community succession and gut fluorescence, fatty acid and elemental compositions (C, N, O, P), and stable isotope (δ13C, δ15N) ratios of dominant mesozooplankton and microzooplankton representatives in a mesocosm setup off Callao (Peru) after simulated upwelling with OMZ water from two different locations and different N:P signatures (moderate and extreme treatments). An oxycline between 5 and 15 m with hypoxic conditions (〈50 µmol L−1) below ∼10 m persisted in the mesocosms throughout the experiment. No treatment effects were determined for the measured parameters, but differences in nutrient concentrations established through OMZ water additions were only minor. Copepods and polychaete larvae dominated in terms of abundance and biomass. Development and reproduction of the dominant copepod genera Paracalanus sp., Hemicyclops sp., Acartia sp., and Oncaea sp. were hindered as evident from accumulation of adult copepodids but largely missing nauplii. Failed hatching of nauplii in the hypoxic bottom layer of the mesocosms and poor nutritional condition of copepods suggested from very low gut fluorescence and fatty acid compositions most likely explain the retarded copepod development. Correlation analysis revealed no particular trophic relations between dominant copepods and phytoplankton groups. Possibly, particulate organic matter with a relatively high C:N ratio was a major diet of copepods. C:N ratios of copepods and polychaetes ranged 4.8–5.8 and 4.2–4.3, respectively. δ15N was comparatively high (∼13 ‰–17 ‰), potentially because the injected OMZ source water was enriched in δ15N as a result of anoxic conditions. Elemental ratios of dinoflagellates deviated strongly from the Redfield ratio. We conclude that opportunistic feeding of copepods may have played an important role in the pelagic food web. Overall, projected changes in the frequency and intensity of upwelling hypoxic waters may make a huge difference for copepod reproduction and may be further enhanced by varying N:P ratios of upwelled OMZ water masses.

Description: Zooplankton plays a notable role in ocean biogeochemical cycles. However, it is often simulated as one generic group and top closure term in ocean biogeochemical models. This study presents the description of three zooplankton functional types (zPFTs, micro-, meso- and macrozooplankton) in the ocean biogeochemical model FESOM-REcoM. In the presented model, microzooplankton is a fast-growing herbivore group, mesozooplankton is another major consumer of phytoplankton, and macrozooplankton is a slow-growing group with a low temperature optimum. Meso- and macrozooplankton produce fast-sinking fecal pellets. With three zPFTs, the annual mean zooplankton biomass increases threefold to 210 Tg C. The new food web structure leads to a 25% increase in net primary production and a 10% decrease in export production globally. Consequently, the export ratio decreases from 17% to 12% in the model. The description of three zPFTs reduces model mismatches with observed dissolved inorganic nitrogen and chlorophyll concentrations in the South Pacific and the Arctic Ocean, respectively. Representation of three zPFTs also strongly affects phytoplankton phenology: Fast nutrient recycling by zooplankton sustains higher chlorophyll concentrations in summer and autumn. Additional zooplankton grazing delays the start of the phytoplankton bloom by 3 weeks and controls the magnitude of the bloom peak in the Southern Ocean. As a result, the system switches from a light-controlled Sverdrup system to a dilution-controlled Behrenfeld system. Overall, the results suggest that representation of multiple zPFTs is important to capture underlying processes that may shape the response of ecosystems and ecosystem services to on-going and future environmental change in model projections.