Bottom-up effects of climate on fish populations: data from the Continuous Plankton Recorder

The Continuous Plankton Recorder (CPR) dataset on fish larvae has an extensive spatio-temporal coverage that allows the responses of fish populations to past changes in climate variability, including abrupt changes such as regime shifts, to be investigated. The newly available dataset offers a unique opportunity to investigate long-term changes over decadal scales in the abundance and distribution of fish larvae in relation to physical and biological factors. A principal component analysis (PCA) using 7 biotic and abiotic parameters is applied to investigate the impact of environmental changes in the North Sea on 5 selected taxa of fish larvae during the period 1960 to 2004. The analysis revealed 4 periods of time (1960–1976; 1977–1982; 1983–1996; 1997–2004) reflecting 3 different ecosystem states. The larvae of clupeids, sandeels, dab and gadoids seemed to be affected mainly by changes in the plankton ecosystem, while the larvae of migratory species such as Atlantic mackerel responded more to hydrographic changes. Climate variability seems more likely to influence fish populations through bottom-up control via a cascading effect from changes in the North Atlantic Oscillation (NAO) impacting on the hydro dynamic features of the North Sea, in turn impacting on the plankton available as prey for fish larvae. The responses and adaptability of fish larvae to changing environmental conditions, parti cularly to changes in prey availability, are complex and species-specific. This complexity is enhanced with fishing effects interacting with climate effects and this study supports furthering our under - standing of such interactions before attempting to predict how fish populations respond to climate variability

Scale-dependent foraging ecology of a marine top predator modelled using passive acoustic data

1.Understanding which environmental factors drive foraging preferences is critical for the development of effective management measures, but resource use patterns may emerge from processes that occur at different spatial and temporal scales. Direct observations of foraging are also especially challenging in marine predators, but passive acoustic techniques provide opportunities to study the behaviour of echolocating species over a range of scales. 2.We used an extensive passive acoustic data set to investigate the distribution and temporal dynamics of foraging in bottlenose dolphins using the Moray Firth (Scotland, UK). Echolocation buzzes were identified with a mixture model of detected echolocation inter-click intervals and used as a proxy of foraging activity. A robust modelling approach accounting for autocorrelation in the data was then used to evaluate which environmental factors were associated with the observed dynamics at two different spatial and temporal scales. 3.At a broad scale, foraging varied seasonally and was also affected by seabed slope and shelf-sea fronts. At a finer scale, we identified variation in seasonal use and local interactions with tidal processes. Foraging was best predicted at a daily scale, accounting for site specificity in the shape of the estimated relationships. 4.This study demonstrates how passive acoustic data can be used to understand foraging ecology in echolocating species and provides a robust analytical procedure for describing spatio-temporal patterns. Associations between foraging and environmental characteristics varied according to spatial and temporal scale, highlighting the need for a multi-scale approach. Our results indicate that dolphins respond to coarser scale temporal dynamics, but have a detailed understanding of finer-scale spatial distribution of resources.

Case history and persistence of the non-indigenous diatom Coscinodiscus wailesii in the north-east Atlantic

The introduction of non-indigenous marine plankton species can have a considerable ecological and economic effect on regional systems. Their presence, however, can go unnoticed until they reach nuisance status and as a consequence few case histories exist containing information on their initial appearance and their spatio-temporal patterns. Here we report on the occurrence of the non-indigenous diatom Coscinodiscus wailesii in 1977 in the English Channel, its subsequent geographical spread into European shelf seas, and its persistence as a significant member of the diatom community in the north-east Atlantic from 1977-1995.

Relating sardine recruitment in the Northern Benguela to satellite-derived sea surface height using a neural network pattern recognition approach

Processes of enrichment, concentration and retention are thought to be important for the successful recruitment of small pelagic fish in upwelling areas, but are difficult to measure. In this study, a novel approach is used to examine the role of spatio-temporal oceanographic variability on recruitment success of the Northern Benguela sardine Sardinops sagax. This approach applies a neural network pattern recognition technique, called a self-organising map (SOM), to a seven-year time series of satellite-derived sea level data. The Northern Benguela is characterised by quasi-perennial upwelling of cold, nutrient-rich water and is influenced by intrusions of warm, nutrient-poor Angola Current water from the north. In this paper, these processes are categorised in terms of their influence on recruitment success through the key ocean triad mechanisms of enrichment, concentration and retention. Moderate upwelling is seen as favourable for recruitment, whereas strong upwelling, weak upwelling and Angola Current intrusion appear detrimental to recruitment success. The SOM was used to identify characteristic patterns from sea level difference data and these were interpreted with the aid of sea surface temperature data. We found that the major oceanographic processes of upwelling and Angola Current intrusion dominated these patterns, allowing them to be partitioned into those representing recruitment favourable conditions and those representing adverse conditions for recruitment. A marginally significant relationship was found between the index of sardine recruitment and the frequency of recruitment favourable conditions (r super(2) = 0.61, p = 0.068, n = 6). Because larvae are vulnerable to environmental influences for a period of at least 50 days after spawning, the SOM was then used to identify windows of persistent favourable conditions lasting longer than 50 days, termed recruitment favourable periods (RFPs). The occurrence of RFPs was compared with back-calculated spawning dates for each cohort. Finally, a comparison of RFPs with the time of spawning and the index of recruitment showed that in years where there were 50 or more days of favourable conditions following spawning, good recruitment followed (Mann-Whitney U-test: p = 0.064, n = 6). These results show the value of the SOM technique for describing spatio-temporal variability in oceanographic processes. Variability in these processes appears to be an important factor influencing recruitment in the Northern Benguela sardine, although the available data time series is currently too short to be conclusive. Nonetheless, the analysis of satellite data, using a neural network pattern-recognition approach, provides a useful framework for investigating fisheries recruitment problems.

The meso-scale response of subarctic North Pacific seabird community structure to lower trophic level abundance and diversity

Understanding the mechanisms that structure communities and influence biodiversity are fundamental goals of ecology. To test the hypothesis that the abundance and diversity of upper-trophic level predators (seabirds) is related to the underlying abundance and diversity of their prey (zooplankton) and ecosystem-wide energy availability (primary production), we initiated a monitoring program in 2002 that jointly and repeatedly surveys seabird and zooplankton populations across a 7,500 km British Columbia-Bering Sea-Japan transect. Seabird distributions were recorded by a single observer (MH) using a strip-width technique, mesozooplankton samples were collected with a Continuous Plankton Recorder, and primary production levels were derived using the appropriate satellite parameters and the Vertically Generalized Production Model (Behrenfeld and Falkowski 1997). Each trophic level showed clear spatio-temporal patterns over the course of the study. The strongest relationship between seabird abundance and diversity and the lower trophic levels was observed in March/April ('spring') and significant relationships were also found through June/July ('summer'). No discernable relationships were observed during the September/October ('fall') months. Overall, mesozooplankton abundance and biomass explained the dominant portion of seabird abundance and diversity indices (richness, Simpson's Index, and evenness), while primary production was only related to seabird richness. These findings underscore the notion that perturbations of ocean productivity and lower trophic level ecosystem constituents influenced by climate change, such as shifts in timing (phenology) and synchronicity (match-mismatch), could impart far-reaching consequences throughout the marine food web.

Seasonal plankton dynamics along a cross-shelf gradient off northern Spain

The development of population models able to reproduce the dynamics of zooplankton is a central issue when trying to understand how a changing environment would affect zooplankton in the future. Using 10 years of monthly data on phytoplankton and zooplankton abundance in the Bay of Biscay from the IEO's RADIALES time-series programme, we built non-parametric Generalized Additive Models (GAMs) able to reproduce the dynamics of plankton on the basis of environmental factors (nutrients, temperature, upwelling and photoperiod). We found that the interaction between these two plankton components is approximately linear, whereas the effects of environmental factors are non-linear. With the inclusion of the environmental variability, the main seasonal and inter-annual dynamic patterns observed within the studied plankton assemblage indicate the prevalence of bottom-up regulatory control. The statistically deduced models were used to simulate the dynamics of the phytoplankton and zooplankton. A good agreement between observations and simulations was obtained, especially for zooplankton. We are presently developing spatio-temporal GAM models for the North Sea based on the Continuous Plankton Recorder database.

Non-linearities, regime shifts and recovery: The recent influence of climate on Black Sea chlorophyll

The Black Sea ecosystem experienced severe eutrophication-related degradation during the 1970s and 1980s. However, in recent years the Black Sea has shown some signs of recovery which are often attributed to a reduction in nutrient loading. Here, SeaWiFS chlorophyll a (chl a), a proxy for phytoplankton biomass, is used to investigate spatio-temporal patterns in Black Sea phytoplankton dynamics and to explore the potential role of climate in the Black Sea's recovery. Maps of chl a anomalies, calculated relative to the 8 year mean, emphasize spatial and temporal variability of phytoplankton biomass in the Black Sea, particularly between the riverine-influenced Northwest Shelf and the open Black Sea. Evolution of phytoplankton biomass has shown significant spatial variability of persistence of optimal bloom conditions between three major regions of the Black Sea. With the exception of 2001, chl a has generally decreased during our 8 year time-series. However, the winter of 2000–2001 was anomalously warm with low wind stress, resulting in reduced vertical mixing of the water column and retention of nutrients in the photic zone. These conditions were associated with anomalously high levels of chl a throughout much of the open Black Sea during the following spring and summer. The unusual climatic conditions occurring in 2001 may have triggered a shift in the Black Sea's chl a regime. The long-term significance of this recent shift is still uncertain but illustrates a non-linear response to climate forcing that makes future ecosystem changes in the pelagic Black Sea ecosystem difficult to predict.

Identifying four phytoplankton functional types from space: An ecological approach

Deriving maps of phytoplankton taxa based on remote sensing data using bio-optical properties of phytoplankton alone is challenging. A more holistic approach was developed using artificial neural networks, incorporating ecological and geographical knowledge together with ocean color, bio-optical characteristics, and remotely sensed physical parameters. Results show that the combined remote sensing approach could discriminate four major phytoplankton functional types (diatoms, dinoflagellates, coccolithophores, and silicoflagellates) with an accuracy of more than 70%. Models indicate that the most important information for phytoplankton functional type discrimination is spatio-temporal information and sea surface temperature. This approach can supply data for large-scale maps of predicted phytoplankton functional types, and an example is shown.

Basin-scale phenology and effects of climate variability on global timing of initial seaward migration of Atlantic salmon (Salmo salar)

Migrations between different habitats are key events in the lives of many organisms. Such movements involve annually recurring travel over long distances usually triggered by seasonal changes in the environment. Often, the migration is associated with travel to or from reproduction areas to regions of growth. Young anadromous Atlantic salmon (Salmo salar) emigrate from freshwater nursery areas during spring and early summer to feed and grow in the North Atlantic Ocean. The transition from the freshwater (parr') stage to the migratory stage where they descend streams and enter salt water (smolt') is characterized by morphological, physiological and behavioural changes where the timing of this parr-smolt transition is cued by photoperiod and water temperature. Environmental conditions in the freshwater habitat control the downstream migration and contribute to within- and among-river variation in migratory timing. Moreover, the timing of the freshwater emigration has likely evolved to meet environmental conditions in the ocean as these affect growth and survival of the post-smolts. Using generalized additive mixed-effects modelling, we analysed spatio-temporal variations in the dates of downstream smolt migration in 67 rivers throughout the North Atlantic during the last five decades and found that migrations were earlier in populations in the east than the west. After accounting for this spatial effect, the initiation of the downstream migration among rivers was positively associated with freshwater temperatures, up to about 10 degrees C and levelling off at higher values, and with sea-surface temperatures. Earlier migration occurred when river discharge levels were low but increasing. On average, the initiation of the smolt seaward migration has occurred 2.5days earlier per decade throughout the basin of the North Atlantic. This shift in phenology matches changes in air, river, and ocean temperatures, suggesting that Atlantic salmon emigration is responding to the current global climate changes.

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.

Basking sharks and oceanographic fronts: quantifying associations in the north‐east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.

Basking sharks and oceanographic fronts: quantifying associations in the north-east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.

Long-term variability of the siphonophores Muggiaea atlantica and M. kochi in the Western English Channel

We investigated long-term variability of the calycophoran siphonophores Muggiaea atlantica and Muggiaea kochi in the Western English Channel (WEC) between 1930 and 2011. Our aims were to describe long-term changes in abundance and temporal distribution in relation to local environmental dynamics. In order to better understand mechanisms that regulate the species’ populations, we identified periods that were characteristic of in situ population growth and the environmental optima associated with these events. Our results show that between 1930 and the 1960s both M. atlantica and M. kochi were transient components of the WEC ecosystem. In the late 1960s M. atlantica, successfully established a resident population in the WEC, while the occurrence of M. kochi became increasingly sporadic. Once established as a resident species, the seasonal abundance and distribution of M. atlantica increased. Analysis of environmental conditions associated with in situ population growth revealed that temperature and prey were key determinants of the seasonal distribution and abundance of M. atlantica. Salinity was shown to have an indirect effect, likely representing a proxy for water circulation in the WEC. Anomalies in the seasonal cycle of salinity, indicating deviation from the usual circulation pattern in the WEC, were negatively associated with in situ growth, suggesting dispersal of the locally developing M. atlantica population. However, our findings identified complexity in the relationship between characteristics of the environment and M. atlantica variability. The transition from a period of transiency (1930–1968) to residency (1969–2011) was tentatively attributed to structural changes in the WEC ecosystem that occurred under the forcing of wider-scale hydroclimatic changes.

Basking sharks and oceanographic fronts: quantifying associations in the north-east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.