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.

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.

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.

Photoacclimation of natural phytoplankton communities

Phytoplankton regulate internal pigment concentrations in response to light and nutrient availability. Chlorophyll to carbon ratios (Chl:Cphyto) are commonly reported as a function of growth irradiance (Eg) for evaluating the photoacclimation response of phytoplankton. In contrast to most culture experiments, natural phytoplankton communities experience fluctuating environmental conditions making it difficult to compare field and lab observations. Observing and understanding photoacclimation in nature is important for deciphering changes in Chl:Cphyto resulting from environmental forcings and for accurately estimating net primary production (NPP) in models which rely on a parameterized description of photoacclimation. Here we employ direct analytical measurements of Cphyto and parallel high-resolution biomass estimates from particulate backscattering (bbp) and flow cytometry to investigate Chl:Cphyto in natural phytoplankton communities. Chl:Cphyto observed over a wide range of Eg in the field was consistent with photoacclimation responses inferred from satellite observations. Field-based photoacclimation observations for a mixed natural community contrast with laboratory results for single species grown in continuous light and nutrient replete conditions. Applying a carbon-based net primary production (NPP) model to our field data for a north-south transect in the Atlantic Ocean results in estimates that closely match 14C depth-integrated NPP for the same cruise and with historical records for the distinct biogeographic regions of the Atlantic Ocean. Our results are consistent with previous satellite and model observations of cells growing in natural or fluctuating light and showcase how direct measurements of Cphyto can be applied to explore phytoplankton photophysiology, growth rates, and production at high spatial resolution in-situ.

Individual-based modelling of the development and transport of a Karenia mikimotoi bloom on the North-West European continental shelf

In 2006, a large and prolonged bloom of the dinoflagellate Karenia mikimotoi occurred in Scottish coastal waters, causing extensive mortalities of benthic organisms including annelids and molluscs and some species of fish ( Davidson et al., 2009). A coupled hydrodynamic-algal transport model was developed to track the progression of the bloom around the Scottish coast during June–September 2006 and hence investigate the processes controlling the bloom dynamics. Within this individual-based model, cells were capable of growth, mortality and phototaxis and were transported by physical processes of advection and turbulent diffusion, using current velocities extracted from operational simulations of the MRCS ocean circulation model of the North-west European continental shelf. Vertical and horizontal turbulent diffusion of cells are treated using a random walk approach. Comparison of model output with remotely sensed chlorophyll concentrations and cell counts from coastal monitoring stations indicated that it was necessary to include multiple spatially distinct seed populations of K. mikimotoi at separate locations on the shelf edge to capture the qualitative pattern of bloom transport and development. We interpret this as indicating that the source population was being transported northwards by the Hebridean slope current from where colonies of K. mikimotoi were injected onto the continental shelf by eddies or other transient exchange processes. The model was used to investigate the effects on simulated K. mikimotoi transport and dispersal of: (1) the distribution of the initial seed population; (2) algal growth and mortality; (3) water temperature; (4) the vertical movement of particles by diurnal migration and eddy diffusion; (5) the relative role of the shelf edge and coastal currents; (6) the role of wind forcing. The numerical experiments emphasized the requirement for a physiologically based biological model and indicated that improved modelling of future blooms will potentially benefit from better parameterisation of temperature dependence of both growth and mortality and finer spatial and temporal hydrodynamic resolution.

Individual-based modelling of the development and transport of a Karenia mikimotoi bloom on the North-West European continental shelf

In 2006, a large and prolonged bloom of the dinoflagellate Karenia mikimotoi occurred in Scottish coastal waters, causing extensive mortalities of benthic organisms including annelids and molluscs and some species of fish ( Davidson et al., 2009). A coupled hydrodynamic-algal transport model was developed to track the progression of the bloom around the Scottish coast during June–September 2006 and hence investigate the processes controlling the bloom dynamics. Within this individual-based model, cells were capable of growth, mortality and phototaxis and were transported by physical processes of advection and turbulent diffusion, using current velocities extracted from operational simulations of the MRCS ocean circulation model of the North-west European continental shelf. Vertical and horizontal turbulent diffusion of cells are treated using a random walk approach. Comparison of model output with remotely sensed chlorophyll concentrations and cell counts from coastal monitoring stations indicated that it was necessary to include multiple spatially distinct seed populations of K. mikimotoi at separate locations on the shelf edge to capture the qualitative pattern of bloom transport and development. We interpret this as indicating that the source population was being transported northwards by the Hebridean slope current from where colonies of K. mikimotoi were injected onto the continental shelf by eddies or other transient exchange processes. The model was used to investigate the effects on simulated K. mikimotoi transport and dispersal of: (1) the distribution of the initial seed population; (2) algal growth and mortality; (3) water temperature; (4) the vertical movement of particles by diurnal migration and eddy diffusion; (5) the relative role of the shelf edge and coastal currents; (6) the role of wind forcing. The numerical experiments emphasized the requirement for a physiologically based biological model and indicated that improved modelling of future blooms will potentially benefit from better parameterisation of temperature dependence of both growth and mortality and finer spatial and temporal hydrodynamic resolution.

Environmental selection and resource allocation determine spatial patterns in picophytoplankton cell size

Here we describe a new trait-based model for cellular resource allocation that we use to investigate the relative importance of different drivers for small cell size in phytoplankton. Using the model, we show that increased investment in nonscalable structural components with decreasing cell size leads to a trade-off between cell size, nutrient and light affinity, and growth rate. Within the most extreme nutrient-limited, stratified environments, resource competition theory then predicts a trend toward larger minimum cell size with increasing depth. We demonstrate that this explains observed trends using a marine ecosystem model that represents selection and adaptation of a diverse community defined by traits for cell size and subcellular resource allocation. This framework for linking cellular physiology to environmental selection can be used to investigate the adaptive response of the marine microbial community to environmental conditions and the adaptive value of variations in cellular physiology.

An updated climatology of surface dimethlysulfide concentrations and emission fluxes in the global ocean

The potentially significant role of the biogenic trace gas dimethylsulfide (DMS) in determining the Earth's radiation budget makes it necessary to accurately reproduce seawater DMS distribution and quantify its global flux across the sea/air interface. Following a threefold increase of data (from 15,000 to over 47,000) in the global surface ocean DMS database over the last decade, new global monthly climatologies of surface ocean DMS concentration and sea-to-air emission flux are presented as updates of those constructed 10 years ago. Interpolation/extrapolation techniques were applied to project the discrete concentration data onto a first guess field based on Longhurst's biogeographic provinces. Further objective analysis allowed us to obtain the final monthly maps. The new climatology projects DMS concentrations typically in the range of 1–7 nM, with higher levels occurring in the high latitudes, and with a general trend toward increasing concentration in summer. The increased size and distribution of the observations in the DMS database have produced in the new climatology substantially lower DMS concentrations in the polar latitudes and generally higher DMS concentrations in regions that were severely undersampled 10 years ago, such as the southern Indian Ocean. Using the new DMS concentration climatology in conjunction with state-of-the-art parameterizations for the sea/air gas transfer velocity and climatological wind fields, we estimate that 28.1 (17.6–34.4) Tg of sulfur are transferred from the oceans into the atmosphere annually in the form of DMS. This represents a global emission increase of 17% with respect to the equivalent calculation using the previous climatology. This new DMS climatology represents a valuable tool for atmospheric chemistry, climate, and Earth System models.

Effect of phytoplankton size classes on bio-optical properties of phytoplankton in the Western Iberian coast: Application of models

Chlorophyll-a satellite products are routinely used in oceanography, providing a synoptic and global view of phytoplankton abundance. However, these products lack information on the community structure of the phytoplankton, which is crucial for ecological modelling and ecosystem studies. To assess the usefulness of existing methods to differentiate phytoplankton functional types (PFT) or phytoplankton size classes from satellite data, in-situ phytoplankton samples collected in the Western Iberian coast, on the North-East Atlantic, were analysed for pigments and absorption spectra. Water samples were collected in five different locations, four of which were located near the shore and another in an open-ocean, seamount region. Three different modelling approaches for deriving phytoplankton size classes were applied to the in situ data. Approaches tested provide phytoplankton size class information based on the input of pigments data (Brewin et al., 2010), absorption spectra data (Ciotti et al., 2002) or both (Uitz et al., 2008). Following Uitz et al. (2008), results revealed high variability in microphytoplankton chlorophyll-specific absorption coefficients, ranging from 0.01 to 0.09 m2 (mg chl)− 1 between 400 and 500 nm. This spectral analysis suggested, in one of the regions, the existence of small cells (< 20 μm) in the fraction of phytoplankton presumed to be microphytoplankton (based on diagnostic pigments). Ciotti et al. (2002) approach yielded the highest differences between modelled and measured absorption spectra for the locations where samples had high variability in community structure and cell size. The Brewin et al. (2010) pigment-based model was adjusted and a set of model coefficients are presented and recommended for future studies in offshore water of the Western Iberian coast.

Mechanisms shaping size structure and functional diversity of phytoplankton communities in the ocean

The factors regulating phytoplankton community composition play a crucial role in structuring aquatic food webs. However, consensus is still lacking about the mechanisms underlying the observed biogeographical differences in cell size composition of phytoplankton communities. Here we use a trait-based model to disentangle these mechanisms in two contrasting regions of the Atlantic Ocean. In our model, the phytoplankton community can self-assemble based on a trade-off emerging from relationships between cell size and (1) nutrient uptake, (2) zooplankton grazing, and (3) phytoplankton sinking. Grazing 'pushes' the community towards larger cell sizes, whereas nutrient uptake and sinking 'pull' the community towards smaller cell sizes. We find that the stable environmental conditions of the tropics strongly balance these forces leading to persistently small cell sizes and reduced size diversity. In contrast, the seasonality of the temperate region causes the community to regularly reorganize via shifts in species composition and to exhibit, on average, bigger cell sizes and higher size diversity than in the tropics. Our results raise the importance of environmental variability as a key structuring mechanism of plankton communities in the ocean and call for a reassessment of the current understanding of phytoplankton diversity patterns across latitudinal gradients.