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.

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.

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.

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.

Assessing the sensitivity of subtidal sedimentary habitats to pressures associated with marine activities. Phase 1 Report: Rationale and proposed ecological groupings for Level 5 biotopes against which sensitivity assessments would be best undertaken.

The purpose of this study is to produce a series of Conceptual Ecological Models (CEMs) that represent sublittoral rock habitats in the UK. CEMs are diagrammatic representations of the influences and processes that occur within an ecosystem. They can be used to identify critical aspects of an ecosystem that may be studied further, or serve as the basis for the selection of indicators for environmental monitoring purposes. The models produced by this project are control diagrams, representing the unimpacted state of the environment free from anthropogenic pressures. It is intended that the models produced by this project will be used to guide indicator selection for the monitoring of this habitat in UK waters. CEMs may eventually be produced for a range of habitat types defined under the UK Marine Biodiversity Monitoring R&D Programme (UKMBMP), which, along with stressor models, are designed to show the interactions within impacted habitats, would form the basis of a robust method for indicator selection. This project builds on the work to develop CEMs for shallow sublittoral coarse sediment habitats (Alexander et al 2014). The project scope included those habitats defined as ‘sublittoral rock’. This definition includes those habitats that fall into the EUNIS Level 3 classifications A3.1 Atlantic and Mediterranean high energy infralittoral rock, A3.2 Atlantic and Mediterranean moderate energy infralittoral rock, A3.3 Atlantic and Mediterranean low energy infralittoral rock, A4.1 Atlantic and Mediterranean high energy circalittoral rock, A4.2 Atlantic and Mediterranean moderate energy circalittoral rock, and A4.3 Atlantic and Mediterranean low energy circalittoral rock as well as the constituent Level 4 and 5 biotopes that are relevant to UK waters. A species list of characterising fauna to be included within the scope of the models was identified using an iterative process to refine the full list of species found within the relevant Level 5 biotopes. A literature review was conducted using a pragmatic and iterative approach to gather evidence regarding species traits and information that would be used to inform the models and characterise the interactions that occur within the sublittoral rock habitat. All information gathered during the literature review was entered into a data logging pro-forma spreadsheet that accompanies this report. Wherever possible, attempts were made to collect information from UK-specific peer-reviewed studies, although other sources were used where necessary. All data gathered was subject to a detailed confidence assessment. Expert judgement by the project team was utilised to provide information for aspects of the models for which references could not be sourced within the project timeframe. A multivariate analysis approach was adopted to assess ecologically similar groups (based on ecological and life history traits) of fauna from the identified species to form the basis of the models. A model hierarchy was developed based on these ecological groups. One general control model was produced that indicated the high-level drivers, inputs, biological assemblages, ecosystem processes and outputs that occur in sublittoral rock habitats. In addition to this, seven detailed sub-models were produced, which each focussed on a particular ecological group of fauna within the habitat: ‘macroalgae’, ‘temporarily or permanently attached active filter feeders’, ‘temporarily or permanently attached passive filter feeders’, ‘bivalves, brachiopods and other encrusting filter feeders’, ‘tube building fauna’, ‘scavengers and predatory fauna’, and ‘non-predatory mobile fauna’. Each sub-model is accompanied by an associated confidence model that presents confidence in the links between each model component. The models are split into seven levels and take spatial and temporal scale into account through their design, as well as magnitude and direction of influence. The seven levels include regional to global drivers, water column processes, local inputs/processes at the seabed, habitat and biological assemblage, output processes, local ecosystem functions, and regional to global ecosystem functions. The models indicate that whilst the high level drivers that affect each ecological group are largely similar, the output processes performed by the biota and the resulting ecosystem functions vary both in number and importance between groups. Confidence within the models as a whole is generally high, reflecting the level of information gathered during the literature review. Physical drivers which influence the ecosystem were found to be of high importance for the sublittoral rock habitat, with factors such as wave exposure, water depth and water currents noted to be crucial in defining the biological assemblages. Other important factors such as recruitment/propagule supply, and those which affect primary production, such as suspended sediments, light attenuation and water chemistry and temperature, were also noted to be key and act to influence the food sources consumed by the biological assemblages of the habitat, and the biological assemblages themselves. Output processes performed by the biological assemblages are variable between ecological groups depending on the specific flora and fauna present and the role they perform within the ecosystem. Of particular importance are the outputs performed by the macroalgae group, which are diverse in nature and exert influence over other ecological groups in the habitat. Important output processes from the habitat as a whole include primary and secondary production, bioengineering, biodeposition (in mixed sediment habitats) and the supply of propagules; these in turn influence ecosystem functions at the local scale such as nutrient and biogeochemical cycling, supply of food resources, sediment stability (in mixed sediment habitats), habitat provision and population and algae control. The export of biodiversity and organic matter, biodiversity enhancement and biotope stability are the resulting ecosystem functions that occur at the regional to global scale. Features within the models that are most useful for monitoring habitat status and change due to natural variation have been identified, as have those that may be useful for monitoring to identify anthropogenic causes of change within the ecosystem. Biological, physical and chemical features of the ecosystem have been identified as potential indicators to monitor natural variation, whereas biological factors and those physical /chemical factors most likely to affect primary production have predominantly been identified as most likely to indicate change due to anthropogenic pressures.

Heavy metal accumulation reflecting natural sedimentary processes and anthropogenic activities in two contrasting coastal wetland ecosystems, eastern China

Due to the impacts of natural processes and anthropogenic activities, different coastal wetlands are faced with variable patterns of heavy metal contamination. It is important to quantify the contributions of pollutant sources, in order to adopt appropriate protection measures for local ecosystems. The aim of this research was to compare the heavy metal contamination patterns of two contrasting coastal wetlands in eastern China. In addition, the contributions from various metal sources were identified and quantified, and influencing factors, such as the role of the plant Spartina alterniflora, were evaluated. Materials and methods Sediment samples were taken from two coastal wetlands (plain-type tidal flat at the Rudong (RD) wetland vs embayment-type tidal flat at Luoyuan Bay (LY)) to measure the content of Al, Fe, Co, Cr, Cu, Mn, Mo, Ni, Sr, Zn, Pb, Cd, and As. Inductively coupled plasma atomic emission spectrometry, flame atomic absorption spectrometry, and atomic fluorescence spectrometry methods were used for metal detection. Meanwhile, the enrichment factor and geoaccumulation index were applied to assess the pollution level. Principle component analysis and receptor modeling were used to quantify the sources of heavy metals. Results and discussion Marked differences in metal distribution patterns between the two systems were present. Metal contents in LY were higher than those in RD, except for Sr and Mo. The growth status of S. alterniflora influenced metal accumulations in RD, i.e., heavy metals were more easily adsorbed in the sediment in the following sequence: Cu > Cd > Zn > Cr > Al > Pb ≥ Ni ≥ Co > Fe > Sr ≥ Mn > As > Mo as a result of the presence and size of the vegetation. However, this phenomenon was not observed in LY. A higher potential ecological risk was associated with LY, compared with RD, except for Mo. Based on a receptor model output, sedimentary heavy metal contents at RD were jointly influenced by natural sedimentary processes and anthropogenic activities, whereas they were dominated by anthropogenic activities at LY. Conclusions A combination of geochemical analysis and modeling approaches was used to quantify the different types of natural and anthropogenic contributions to heavy metal contamination, which is useful for pollution assessments. The application of this approach reveals that natural and anthropogenic processes have different influences on the delivery and retention of metals at the two contrasting coastal wetlands. In addition, the presence and size of S. alterniflora can influence the level of metal contamination in sedimentary environments.