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.

How should sparse marine in situ measurements be compared to a continuous model: an example

This work demonstrates an example of the importance of an adequate method to sub-sample model results when comparing with in situ measurements. A test of model skill was performed by employing a point-to-point method to compare a multi-decadal hindcast against a sparse, unevenly distributed historic in situ dataset. The point-to-point method masked out all hindcast cells that did not have a corresponding in situ measurement in order to match each in situ measurement against its most similar cell from the model. The application of the point-to-point method showed that the model was successful at reproducing the inter-annual variability of the in situ datasets. Furthermore, this success was not immediately apparent when the measurements were aggregated to regional averages. Time series, data density and target diagrams were employed to illustrate the impact of switching from the regional average method to the point-to-point method. The comparison based on regional averages gave significantly different and sometimes contradicting results that could lead to erroneous conclusions on the model performance. Furthermore, the point-to-point technique is a more correct method to exploit sparse uneven in situ data while compensating for the variability of its sampling. We therefore recommend that researchers take into account for the limitations of the in situ datasets and process the model to resemble the data as much as possible.

Continuous Plankton Records - Persistence In Time-Series Of Annual Means Of Abundance Of Zooplankton

Time-series of annual means of abundance of zooplankton of the north-east Atlantic Ocean and the North Sea, for the period 1948 to 1977, show considerable associations between successive years. The seasonal dynamics of the stocks appear to be consistent with at least a proportion of this being due to inherent persistence from year-to-year. Experiments with a simple model suggest that the observed properties of the time-series cannot be reproduced as a response to simple random forcing. The extent of trends and long wavelength variations can be simulated by introducing fairly extensive persistence into the perturbations, but this underestimates the extent of shorter wavelength variability in the observed time-series. The effect of persistence is to increase the proportion of trend and long wavelength variability in time-series of annual means, but stocks can respond to short wavelength perturbations provided these have a clearly defined frequency.

Oceanic controls on the primary production of the northwest European continental shelf: model experiments under recent past conditions and a potential future scenario

In this paper we clearly demonstrate that changes in oceanic nutrients are a first order factor in determining changes in the primary production of the northwest European continental shelf on time scales of 5–10 yr. We present a series of coupled hydrodynamic ecosystem modelling simulations, using the POLCOMS-ERSEM system. These are forced by both reanalysis data and a single example of a coupled ocean-atmosphere general circulation model (OA-GCM) representative of possible conditions in 2080–2100 under an SRES A1B emissions scenario, along with the corresponding present day control. The OA-GCM forced simulations show a substantial reduction in surface nutrients in the open-ocean regions of the model domain, comparing future and present day time-slices. This arises from a large increase in oceanic stratification. Tracer transport experiments identify a substantial fraction of on-shelf water originates from the open-ocean region to the south of the domain, where this increase is largest, and indeed the on-shelf nutrient and primary production are reduced as this water is transported on-shelf. This relationship is confirmed quantitatively by comparing changes in winter nitrate with total annual nitrate uptake. The reduction in primary production by the reduced nutrient transport is mitigated by on-shelf processes relating to temperature, stratification (length of growing season) and recycling. Regions less exposed to ocean-shelf exchange in this model (Celtic Sea, Irish Sea, English Channel, and Southern North Sea) show a modest increase in primary production (of 5–10%) compared with a decrease of 0–20% in the outer shelf, Central and Northern North Sea. These findings are backed up by a boundary condition perturbation experiment and a simple mixing model.

Carbon fluxes and pelagic ecosystem dynamics near two western Antarctic Peninsula Adélie penguin colonies: an inverse model approach

An inverse food-web model for the western Antarctic Peninsula (WAP) pelagic food web was constrained with data from Palmer Long Term Ecological Research (PAL-LTER) project annual austral summer sampling cruises. Model solutions were generated for 2 regions with Adelie penguin Pygoscelis adeliae colonies presenting different population trends (a northern and a southern colony) for a 12 yr period (1995-2006). Counter to the standard paradigm, comparisons of carbon flow through bacteria, microzooplankton, and krill showed that the diatom-krill-top predator food chain is not the dominant pathway for organic carbon exchanges. The food web is more complex, including significant contributions by microzooplankton and the microbial loop. Using both inverse model results and network indices, it appears that in the northern WAP the food web is dominated by the microbial food web, with a temporal trend toward its increasing importance. The dominant pathway for the southern WAP food web varies from year to year, with no detectable temporal trend toward dominance of microzooplankton versus krill. In addition, sensitivity analyses indicated that the northern colony of Adelie penguins, whose population size has been declining over the past 35 yr, appears to have sufficient krill during summer to sustain its basic metabolic needs and rear chicks, suggesting the importance of other processes in regulating the Adelie population decline.

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.

Regional and annual variation in black-legged kittiwake breeding productivity is related to sea surface temperature

Spatiotemporal variation in seabird demographic parameters is often pronounced and may be an important source of information on the state of marine ecosystems. Black-legged kittiwakes Rissa tridactyla in Britain and Ireland show strong regional structure in breeding productivity, and both temporal and spatial variation are probably related to abundance of the principal prey of breeding kittiwakes, the lesser sandeel Ammodytes marinus. Annual regional estimates of sandeel abundance do not exist, prohibiting direct tests of this hypothesis. We examined relationships between kittiwake breeding productivity and 2 potential proxies of sandeel abundance, winter sea surface temperature (SST) and abundance of Calanus copepods, within and among 6 regions in Britain and Ireland from 1986 to 2004. Means and trends in winter SST differed among regions, with higher means and less pronounced increasing trends in western (Atlantic) regions than in eastern (North Sea) regions. A negative relationship between breeding productivity and winter SST in the previous year was found within 2 regions (East Scotland and Orkney), as well as in a cross-regional analysis. Results were inconclusive for Calanus abundance, with a positive relationship in East Scotland and negative in Orkney. These results demonstrate that although a single environmental driver (SST) is related to both within- and between-region variation in a key demographic parameter, regional heterogeneity in SST trends as well as the importance of other factors may lead to highly variable responses. Understanding this heterogeneity is critical for predicting long-term effects of climate change or other anthropogenic drivers on marine ecosystems.

Exploring the Red Sea seasonal ecosystem functioning using a three-dimensional biophysical model

The Red Sea exhibits complex hydrodynamic and biogeochemical dynamics, which vary both in time and space. These dynamics have been explored through the development and application of a 3-D ecosystem model. The simulation system comprises two off-line coupled submodels: the MIT General Circulation Model (MITgcm) and the European Regional Seas Ecosystem Model (ERSEM), both adapted for the Red Sea. The results from an annual simulation under climatological forcing are presented. Simulation results are in good agreement with satellite and in situ data illustrating the role of the physical processes in determining the evolution and variability of the Red Sea ecosystem. The model was able to reproduce the main features of the Red Sea ecosystem functioning, including the exchange with the Gulf of Aden, which is a major driving mechanism for the whole Red Sea ecosystem and the winter overturning taking place in the north. Some model limitations, mainly related to the dynamics of the extended reef system located in the southern part of the Red Sea, which is not currently represented in the model, still need to be addressed.

Comparison of new and primary production models using SeaWiFS data in contrasting hydrographic zones of the northern North Atlantic.

The accuracy of two satellite models of marine primary (PP) and new production (NP) were assessed against 14C and 15N uptake measurements taken during six research cruises in the northern North Atlantic. The wavelength resolving model (WRM) was more accurate than the Vertical General Production Model (VGPM) for computation of both PP and NP. Mean monthly satellite maps of PP and NP for both models were generated from 1997 to 2010 using SeaWiFS data for the Irminger basin and North Atlantic. Intra- and inter-annual variability of the two models was compared in six hydrographic zones. Both models exhibited similar spatio-temporal patterns: PP and NP increased from April to June and decreased by August. Higher values were associated with the East Greenland Current (EGC), Iceland Basin (ICB) and the Reykjanes Ridge (RKR) and lower values occurred in the Central Irminger Current (CIC), North Irminger Current (NIC) and Southern Irminger Current (SIC). The annual PP and NP over the SeaWiFS record was 258 and 82 gC m-2 yr-1 respectively for the VGPM and 190 and 41 gC m-2 yr-1 for the WRM. Average annual cumulative sum in the anomalies of NP for the VGPM were positively correlated with the North Atlantic Oscillation (NAO) in the EGC, CIC and SIC and negatively correlated with the multivariate ENSO index (MEI) in the ICB. By contrast, cumulative sum of the anomalies of NP for the WRM were significantly correlated with NAO only in the EGC and CIC. NP from both VGPM and WRM exhibited significant negative correlations with Arctic Oscillation (AO) in all hydrographic zones. The differences in estimates of PP and NP in these hydrographic zones arise principally from the parameterisation of the euphotic depth and the SST dependence of photo-physiological term in the VGPM, which has a greater sensitivity to variations in temperature than the WRM. In waters of 0 to 5C PP using the VGPM was 43% higher than WRM, from 5 to 10C the VGPM was 29% higher and from 10 to 15C the VGPM was 27% higher.