Have climate induced changes in zooplankton communities led to poor conditions for larval fish growth? A test on cod larvae on the UK shelf

Climate effects have been shown to be at least partly responsible for the reorganisation in the plankton ecosystem on the shelf seas of NW Europe over the last 50 years. Most fish larvae feed primarily on zooplankton, so changes in zooplankton quantity, quality and seasonal timing have been hypothesized to be a key factor affecting their survival. To investigate this we have implemented a 1-dimensional trophodynamic growth model of cod larvae for the waters around the UK covering the period 1960 to 2003. Larval growth is modelled as the difference between the amount of food absorbed by the larva and its various metabolic costs. Prey availability is based upon the biomass and size of available preys (i.e. adults and nauplii copepods and cladocerans) taken from the Continuous Plankton Recorder dataset. Temperature and wind forcing are also taken into account. Results suggest that observed changes in plankton community structure may have had less impact than previously suggested. This is because changes in prey availability may be compensated for by increased temperatures resulting in little overall impact on potential larval growth. Stock recovery, at least in the short term is likely to be more dependent upon conserving the year classes recruited to allow spawning stock biomass to rebuild. If as our model suggests, the larvae are still able to survive in the changing environment, reduction in fishing on the adults is needed to allow the stock to recover.

An individual-based model of the early life history of mackerel (Scomber scombrus) in the eastern North Atlantic, simulating transport, growth and mortality

The main purpose of this paper is to provide the core description of the modelling exercise within the Shelf Edge Advection Mortality And Recruitment (SEAMAR) programme. An individual-based model (IBM) was developed for the prediction of year-to-year survival of the early life-history stages of mackerel (Scomber scombrus) in the eastern North Atlantic. The IBM is one of two components of the model system. The first component is a circulation model to provide physical input data for the IBM. The circulation model is a geographical variant of the HAMburg Shelf Ocean Model (HAMSOM). The second component is the IBM, which is an i-space configuration model in which large numbers of individuals are followed as discrete entities to simulate the transport, growth and mortality of mackerel eggs, larvae and post-larvae. Larval and post-larval growth is modelled as a function of length, temperature and food distribution; mortality is modelled as a function of length and absolute growth rate. Each particle is considered as a super-individual representing 10 super(6) eggs at the outset of the simulation, and then declining according to the mortality function. Simulations were carried out for the years 1998-2000. Results showed concentrations of particles at Porcupine Bank and the adjacent Irish shelf, along the Celtic Sea shelf-edge, and in the southern Bay of Biscay. High survival was observed only at Porcupine and the adjacent shelf areas, and, more patchily, around the coastal margin of Biscay. The low survival along the shelf-edge of the Celtic Sea was due to the consistently low estimates of food availability in that area.

Mechanisms of long-term decline in size of lesser sandeels in the North Sea explored using a growth and phenology model

The lesser sandeel Ammodytes marinus is a key species in the North Sea ecosystem, transferring energy from planktonic producers to top predators. Previous studies have shown a long-term decline in the size of 0-group sandeels in the western North Sea, but they were unable to pinpoint the mechanism (later hatching, slower growth or changes in size-dependent mortality) or cause. To investigate the first 2 possibilities we combined 2 independent time series of sandeel size, namely data from chick-feeding Atlantic puffins Fratercula arctica and from the Continuous Plankton Recorder (CPR), in a novel statistical model implemented using Markov Chain Monte Carlo (MCMC). The model estimated annual mean length on 1 July, as well as hatching date and growth rate for sandeels from 1973 to 2006. Mean length-at-date declined by 22% over this period, corresponding to a 60% decrease in energy content, with a sharper decline since 2002. Up to the mid-1990s, the decline was associated with a trend towards later hatching. Subsequently, hatching became earlier again, and the continued trend towards smaller size appears to have been driven by lower growth rates, particularly in the most recent years, although we could not rule out changes in size-dependent mortality. Our findings point to major changes in key aspects of sandeel life history, which we consider are most likely due to direct and indirect temperature-related changes over a range of biotic factors, including the seasonal distribution of copepods and intra- and inter-specific competition with planktivorous fish. The results have implications both for the many predators of sandeels and for age and size of maturation in this aggregation of North Sea sandeels.

Simulation of mackerel (Scomber scombrus) recruitment with an individual-based model and comparison with field data

An individual-based model (IBM) for the simulation of year-to-year survival during the early life-history stages of the north-east Atlantic stock of mackerel (Scomber scombrus) was developed within the EU funded Shelf-Edge Advection, Mortality and Recruitment (SEAMAR) programme. The IBM included transport, growth and survival and was used to track the passive movement of mackerel eggs, larvae and post-larvae and determine their distribution and abundance after approximately 2 months of drift. One of the main outputs from the IBM, namely distributions and numbers of surviving post-larvae, are compared with field data as recruit (age-0/age-1 juveniles) distribution and abundance for the years 1998, 1999 and 2000. The juvenile distributions show more inter-annual and spatial variability than the modelled distributions of survivors; this may be due to the restriction of using the same initial egg distribution for all 3 yr of simulation. The IBM simulations indicate two main recruitment areas for the north-east Atlantic stock of mackerel, these being Porcupine Bank and the south-eastern Bay of Biscay. These areas correspond to areas of high juvenile catches, although the juveniles generally have a more widespread distribution than the model simulations. The best agreement between modelled data and field data for distribution (juveniles and model survivors) is for the year 1998. The juvenile catches in different representative nursery areas are totalled to give a field abundance index (FAI). This index is compared with a model survivor index (MSI) which is calculated from the total of survivors for the whole spawning season. The MSI compares favourably with the FAI for 1998 and 1999 but not for 2000; in this year, juvenile catches dropped sharply compared with the previous years but there was no equivalent drop in modelled survivors.

Empirically modelling the potential effects of changes in temperature and prey availability on the growth of cod larvae in UK shelf seas

It has been hypothesized that changes in zooplankton community structure over the past four decades led to reduced growth and survival of prerecruit Atlantic cod (Gadus morhua) and that this was a key factor underlying poor year classes, contributing to stock collapse, and inhibiting the recovery of stocks around the UK. To evaluate whether observed changes in plankton abundance, species composition and temperature could have led to periods of poorer growth of cod larvae, we explored the effect of prey availability and temperature on early larval growth using an empirical trophodynamic model. Prey availability was parameterized using species abundance data from the Continuous Plankton Recorder. Our model suggests that the observed changes in plankton community structure in the North Sea may have had less impact on cod larval growth, at least for the first 40 days following hatching, than previously suggested. At least in the short term, environmental and prey conditions should be able to sustain growth of cod larvae and environmental changes acting on this early life stage should not limit stock recovery.

Growth and production of Pacific ocean perch (Sebastes alutus) in nursery habitats of the Gulf of Alaska

Nursery areas for juvenile fishes are often important for determining recruitment in marine populations by providing habitats that can maximize growth and thereby minimize mortality. Pacific ocean perch (POP, Sebastes alutus) have an extended juvenile period where they inhabit rocky nursery habitats. We examined POP nursery areas to link growth potential to recruitment. Juvenile POP were captured from nursery areas in 2004 and 2008, and estimated growth rates ranged from −0.19 to 0.60 g day−1 based on differences in size between June and August. Predicted growth rates from a bioenergetics model ranged from 0.05 to 0.49 g day−1 and were not significantly different than observed. Substrate preferences and the distribution of their preferred habitats were utilized to predict the extent of juvenile POP nursery habitat in the Gulf of Alaska. Based on densities of fish observed on underwater video transects and the spatial extent of nursery areas, we predicted 278 and 290 million juvenile POP were produced in 2004 and 2008. Growth potential for juvenile POP was reconstructed using the bioenergetics model, spring zooplankton bloom timing and duration and bottom water temperature for 1982–2008. When a single outlying recruitment year in 1986 was removed, growth potential experienced by juvenile POP in nursery areas was significantly correlated to the recruitment time-series from the stock assessment, explaining ∼30% of the variability. This research highlights the potential to predict recruitment using habitat-based methods and provides a potential mechanism for explaining some of the POP recruitment variability observed for this population.

Vulnerability of coastal ecosystems to changes in harmful algal bloom distribution in response to climate change: projections based on model analysis

Harmful algal blooms (HABs), those proliferations of algae that can cause fish kills, contaminate seafood with toxins, form unsightly scums, or detrimentally alter ecosystem function have been increasing in frequency, magnitude, and duration worldwide. Here, using a global modeling approach, we show, for three regions of the globe, the potential effects of nutrient loading and climate change for two HAB genera, pelagic Prorocentrum and Karenia, each with differing physiological characteristics for growth. The projections (end of century, 2090-2100) are based on climate change resulting from the A1B scenario of the Intergovernmental Panel on Climate Change Institut Pierre Simon Laplace Climate Model (IPCC, IPSL-CM4), applied in a coupled oceanographic-biogeochemical model, combined with a suite of assumed physiological 'rules' for genera-specific bloom development. Based on these models, an expansion in area and/or number of months annually conducive to development of these HABs along the NW European Shelf-Baltic Sea system and NE Asia was projected for both HAB genera, but no expansion (Prorocentrum spp.), or actual contraction in area and months conducive for blooms (Karenia spp.), was projected in the SE Asian domain. The implications of these projections, especially for Northern Europe, are shifts in vulnerability of coastal systems to HAB events, increased regional HAB impacts to aquaculture, increased risks to human health and ecosystems, and economic consequences of these events due to losses to fisheries and ecosystem services.

Multiple environmental controls on phytoplankton growth strategies determine adaptive responses of the N : P ratio

The controls on the 'Redfield' N:P stoichiometry of marine phytoplankton and hence the N:P ratio of the deep ocean remain incompletely understood. Here, we use a model for phytoplankton ecophysiology and growth, based on functional traits and resource-allocation trade-offs, to show how environmental filtering, biotic interactions, and element cycling in a global ecosystem model determine phytoplankton biogeography, growth strategies and macromolecular composition. Emergent growth strategies capture major observed patterns in marine biomes. Using a new synthesis of experimental RNA and protein measurements to constrain per-ribosome translation rates, we determine a spatially variable lower limit on adaptive rRNA:protein allocation and hence on the relationship between the largest cellular P and N pools. Comparison with the lowest observed phytoplankton N:P ratios and N:P export fluxes in the Southern Ocean suggests that additional contributions from phospholipid and phosphorus storage compounds play a fundamental role in determining the marine biogeochemical cycling of these elements.

Trait-based representation of diatom functional diversity in a plankton functional type model of the eutrophied southern North Sea

We introduce a trait-based description of diatom functional diversity to an existing plankton functional type (PFT) model, implemented for the eutrophied coastal ecosystem in the Southern Bight of the North Sea. The trait-based description represents a continuum of diatom species, each characterized by a distinct cell volume, and includes size dependence of four diatom traits: the maximum growth rate, the half-saturation constants for nutrient uptake, the photosynthetic efficiency, and the relative affinity of copepods for diatoms. Through competition under seasonally varying forcing, the fitness of each diatom varies throughout time, and the outcome of competition results in a changing community structure. The predicted seasonal change in mean cell volume of the community is supported by field observations: smaller diatoms, which are more competitive in terms of resource acquisition, prevail during the first spring bloom, whereas the summer bloom is dominated by larger species which better resist grazing. The size-based model is used to determine the ecological niche of diatoms in the area and identifies a range of viable sizes that matches observations. The general trade-off between small, competitive diatoms and large, grazing-resistant species is a convenient framework to study patterns in diatom functional diversity. PFT models and trait-based approaches constitute promising complementary tools to study community structure in marine ecosystems.

Fundamental shift in vitamin B12 eco-physiology of a model alga demonstrated by experimental evolution.

A widespread and complex distribution of vitamin requirements exists over the entire tree of life, with many species having evolved vitamin dependence, both within and between different lineages. Vitamin availability has been proposed to drive selection for vitamin dependence, in a process that links an organism's metabolism to the environment, but this has never been demonstrated directly. Moreover, understanding the physiological processes and evolutionary dynamics that influence metabolic demand for these important micronutrients has significant implications in terms of nutrient acquisition and, in microbial organisms, can affect community composition and metabolic exchange between coexisting species. Here we investigate the origins of vitamin dependence, using an experimental evolution approach with the vitamin B(12)-independent model green alga Chlamydomonas reinhardtii. In fewer than 500 generations of growth in the presence of vitamin B(12), we observe the evolution of a B(12)-dependent clone that rapidly displaces its ancestor. Genetic characterization of this line reveals a type-II Gulliver-related transposable element integrated into the B(12)-independent methionine synthase gene (METE), knocking out gene function and fundamentally altering the physiology of the alga.