Climate-induced effects on the meroplankton and the benthic-pelagic ecology of the North Sea

Analyses of long-term time series of North Sea plankton and sea surface temperature (SST) data reveal that the annual planktonic larval abundance of three benthic phyla, Echinodermata, Arthropoda, and Mollusca, responds positively and immediately to SST. Long-term outcomes for the planktonic abundance of these three phyla are different, however. The planktonic larvae of echinoderms and decapod crustaceans have increased in abundance from 1958 to 2005, and especially since the mid-1980s, as North Sea SST has increased. In contrast, the abundance of bivalve mollusc larvae has declined, despite the positive year-to-year relationship between temperature and bivalve larval abundance continuing to hold. We argue that the changes in meroplankton abundance, coincident with increased phytoplankton and declining holoplankton, reflect the synchronous effect of rising SST and related changes in the pelagic community on the reproduction and recruitment of many benthic marine invertebrates. Under this scenario, the long-term decline in bivalve mollusc larvae will reflect increased predation on the settled larvae and adults by benthic decapods. These alterations in the zooplankton may therefore describe an ecosystem-wide restructuring of North Sea trophic interactions.

Influence of climate change and trophic coupling across four trophic levels in the Celtic Sea

Climate change has had profound effects upon marine ecosystems, impacting across all trophic levels from plankton to apex predators. Determining the impacts of climate change on marine ecosystems requires understanding the direct effects on all trophic levels as well as indirect effects mediated by trophic coupling. The aim of this study was to investigate the effects of climate change on the pelagic food web in the Celtic Sea, a productive shelf region in the Northeast Atlantic. Using long-term data, we examined possible direct and indirect ‘bottom-up’ climate effects across four trophic levels: phytoplankton, zooplankton, mid-trophic level fish and seabirds. During the period 1986–2007, although there was no temporal trend in the North Atlantic Oscillation index (NAO), the decadal mean Sea Surface Temperature (SST) in the Celtic Sea increased by 0.66±0.02°C. Despite this, there was only a weak signal of climate change in the Celtic Sea food web. Changes in plankton community structure were found, however this was not related to SST or NAO. A negative relationship occurred between herring abundance (0- and 1-group) and spring SST (0-group: p = 0.02, slope = −0.305±0.125; 1-group: p = 0.04, slope = −0.410±0.193). Seabird demographics showed complex species–specific responses. There was evidence of direct effects of spring NAO (on black-legged kittiwake population growth rate: p = 0.03, slope = 0.0314±0.014) as well as indirect bottom-up effects of lagged spring SST (on razorbill breeding success: p = 0.01, slope = −0.144±0.05). Negative relationships between breeding success and population growth rate of razorbills and common guillemots may be explained by interactions between mid-trophic level fish. Our findings show that the impacts of climate change on the Celtic Sea ecosystem is not as marked as in nearby regions (e.g. the North Sea), emphasizing the need for more research at regional scales.

Drivers of euphausiid species abundance and numerical abundance in the Atlantic Ocean

Mid-ocean ridges are common features of the world’s oceans but there is a lack of understanding as to how their presence affects overlying pelagic biota. The Mid-Atlantic Ridge (MAR) is a dominant feature of the Atlantic Ocean. Here, we examined data on euphausiid distribution and abundance arising from several international research programmes and from the continuous plankton recorder. We used a generalized additive model (GAM) framework to explore spatial patterns of variability in euphausiid distribution on, and at either side of, the MAR from 60°N to 55°S in conjunction with variability in a suite of biological, physical and environmental parameters. Euphausiid species abundance peaked in mid-latitudes and was significantly higher on the ridge than in adjacent waters, but the ridge did not influence numerical abundance significantly. Sea surface temperature (SST) was the most important single factor influencing both euphausiid numerical abundance and species abundance. Increases in sea surface height variance, a proxy for mixing, increased the numerical abundance of euphausiids. GAM predictions of variability in species abundance as a function of SST and depth of the mixed layer were consistent with present theories, which suggest that pelagic niche availability is related to the thermal structure of the near surface water: more deeply-mixed water contained higher euphausiid biodiversity. In addition to exposing present distributional patterns, the GAM framework enables responses to potential future and past environmental variability including temperature change to be explored.

Effects of environmental conditions on the seasonal distribution of phytoplankton biomass in the North Sea

We study the spatial and seasonal variability of phytoplankton biomass (as phytoplankton color) in relation to the environmental conditions in the North Sea using data from the Continuous Plankton Recorder survey. By using only environmental fields and location as predictor variables we developed a nonparametric model (generalized additive model) to empirically explore how key environmental factors modulate the spatio-temporal patterns of the seasonal cycle of algal biomass as well as how these relate to the ,1988 North Sea regime shift. Solar radiation, as manifest through changes of sea surface temperature (SST), was a key factor not only in the seasonal cycle but also as a driver of the shift. The pronounced increase in SST and in wind speed after the 1980s resulted in an extension of the season favorable for phytoplankton growth. Nutrients appeared to be unimportant as explanatory variables for the observed spatio-temporal pattern, implying that they were not generally limiting factors. Under the new climatic regime the carrying capacity of the whole system has been increased and the southern North Sea, where the environmental changes have been more pronounced, reached a new maximum.

Variation among northeast Atlantic regions in the responses of zooplankton to climate change: Not all areas follow the same path

Broad scale climate forcing can interact with local environmental processes to affect the observed ecological phenomena. This causes potential problems of over-extrapolation for results from a limited number of sites or the averaging out of region-specific responses if data from too wide an area are combined. In this study, an area similar in extent to the Celtic Biscay Large Marine Ecosystem, but including off-shelf areas, was partitioned using clustering of satellite chlorophyll (chl-a) measurements. The resulting clusters were used to define areas over which to combine copepod data from the Continuous Plankton Recorder. Following filtering due to data limitations, nine regions were defined with sufficient records for analysis. These regions were consistent with known oceanographic structure in the study area. Off-shelf regions showed a progressively later timing in the seasonal peak of chl-a measurements moving northwards. Generalised additive models were used to estimate seasonal and multiannual signals in the adult and juvenile stages of Calanus finmarchicus, C. helgolandicus and the Paracalanus–Pseudocalanus group. Associations between variables (sea surface temperature (SST), phenology and annual abundance) differed among taxonomic groups, but even within taxonomic groups, relationships were not consistent across regions. For example, in the deep waters off Spain and Portugal the annual abundance of Calanus finmarchicus has a weak positive association with SST, in contrast to the pattern in most other regions. The regions defined in this study provide an objective basis for investigations into the long term dynamics of plankton populations and suggest suitable sub regions for deriving pelagic system indicators.

Spatio-Temporal Variability of the North Sea Cod Recruitment in Relation to Temperature and Zooplankton

The North Sea cod (Gadus morhua, L.) stock has continuously declined over the past four decades linked with overfishing and climate change. Changes in stock structure due to overfishing have made the stock largely dependent on its recruitment success, which greatly relies on environmental conditions. Here we focus on the spatio-temporal variability of cod recruitment in an effort to detect changes during the critical early life stages. Using International Bottom Trawl Survey (IBTS) data from 1974 to 2011, a major spatio-temporal change in the distribution of cod recruits was identified in the late 1990s, characterized by a pronounced decrease in the central and southeastern North Sea stock. Other minor spatial changes were also recorded in the mid-1980s and early 1990s. We tested whether the observed changes in recruits distribution could be related with direct (i.e. temperature) and/or indirect (i.e. changes in the quantity and quality of zooplankton prey) effects of climate variability. The analyses were based on spatially-resolved time series, i.e. sea surface temperature (SST) from the Hadley Center and zooplankton records from the Continuous Plankton Recorder Survey. We showed that spring SST increase was the main driver for the most recent decrease in cod recruitment. The late 1990s were also characterized by relatively low total zooplankton biomass, particularly of energy-rich zooplankton such as the copepod Calanus finmarchicus, which have further contributed to the decline of North Sea cod recruitment. Long-term spatially-resolved observations were used to produce regional distribution models that could further be used to predict the abundance of North Sea cod recruits based on temperature and zooplankton food availability.

The influence of phytoplankton productivity, temperature and environmental stability on the control of copepod diversity in the North East Atlantic

The patterns of copepod species richness (S) and their relationship with phytoplankton productivity, temperature and environmental stability were investigated at climatological, seasonal and year-to-year time scales as well as scales along latitudinal and oceanic–neritic gradients using monthly time series of the Continuous Plankton Recorder (CPR) Survey collected in the North East Atlantic between 1958 and 2006. Time series analyses confirmed previously described geographic patterns. Equatorward and towards neritic environments, the climatological average of S increases and the variance explained by the seasonal cycle decreases. The bi-modal character of seasonality increases equatorward and the timing of the seasonal cycle takes place progressive earlier equatorward and towards neritic environments. In the long-term, the climatological average of S decreased significantly (p < 0.001) between 1958 and 2006 in the Bay of Biscay and North Iberian shelf at a rate of ca. 0.04 year−1, and increased at the same rate between 1991 and 2006 in the northernmost oceanic location. The climatological averages of S correlate positively with those of the index of seasonality of phytoplankton productivity (ratio between the minimum and maximum monthly values of surface chlorophyll) and sea surface temperature, and negatively with those of the proxy for environmental stability (monthly frequency of occurrence of daily averaged wind speed exceeding 10 m s−1). The seasonal cycles of S and phytoplankton productivity (surface chlorophyll as proxy) exhibit similar features in terms of shape, timing and explained variance, but the relationship between the climatological averages of both variables is non-significant. From year-to-year, the annual averages of S correlate negatively with those of phytoplankton productivity and positively with those of sea surface temperature along the latitudinal gradient, and negatively with those of environmental stability along the oceanic–neritic gradient. The annual anomalies of S (i.e. factoring out geographic variation) show a unimodal relationship with those of sea surface temperature and environmental stability, with S peaking at intermediate values of the anomalies of these variables. The results evidence the role of seasonality of phytoplankton productivity on the control of copepod species richness at seasonal and climatological scales, giving support to the species richness–productivity hypothesis. Although sea surface temperature (SST) is indeed a good predictor of richness along the latitudinal gradient, it is unable to predict the increase of richness form oceanic to neritic environments, thus lessening the generality of the species richness–energy hypothesis. Meteo-hydrographic disturbances (i.e. SST and wind speed anomalies as proxies), presumably through its role on mixed layer depth dynamics and turbulence and hence productivity, maximise local diversity when occurring at intermediate frequency and or intensity, thus providing support to the intermediate disturbance hypothesis on the control of copepod diversity.

Macroscale factors affecting diatom abundance: a synergistic use of Continuous Plankton Recorder and satellite remote sensing data

Diatoms exist in almost every aquatic regime; they are responsible for 20% of global carbon fixation and 25% of global primary production, and are regarded as a key food for copepods, which are subsequently consumed by larger predators such as fish and marine mammals. A decreasing abundance and a vulnerability to climatic change in the North Atlantic Ocean have been reported in the literature. In the present work, a data matrix composed of concurrent satellite remote sensing and Continuous Plankton Recorder (CPR) in situ measurements was collated for the same spatial and temporal coverage in the Northeast Atlantic. Artificial neural networks (ANNs) were applied to recognize and learn the complex non-monotonic and non-linear relationships between diatom abundance and spatiotemporal environmental factors. Because of their ability to mimic non-linear systems, ANNs proved far more effective in modelling the diatom distribution in the marine ecosystem. The results of this study reveal that diatoms have a regular seasonal cycle, with their abundance most strongly influenced by sea surface temperature (SST) and light intensity. The models indicate that extreme positive SSTs decrease diatom abundances regardless of other climatic conditions. These results provide information on the ecology of diatoms that may advance our understanding of the potential response of diatoms to climatic change.

Enhancement of primary production in the North Atlantic outside of the spring bloom, identified by remote sensing of ocean colour and temperature

The heterogeneity in phytoplankton production in the North Atlantic after the spring bloom is poorly understood. We analysed merged microwave and infrared satellite sea surface temperature (SST) data and ocean colour phytoplankton size class biomass, primary production (PP) and new production (ExP) derived from SeaWiFS data, to assess the spatial and temporal frequency of surface thermal fronts and areas of enhanced PP and ExP. Strong and persistent surface thermal fronts occurred at the Reykjanes Ridge (RR) and sub-polar front (SPF), which sustain high PP and ExP and, outside of the spring bloom, account for 9% and 15% of the total production in the North Atlantic. When normalised by area, PP at the SPF is four times higher than the RR. Analysis of 13 years of satellite ocean colour data from SeaWiFS, and compared with MODIS-Aqua and MERIS, showed that there was no increase in Chla from 1998 to 2002, which then decreased in all areas from 2002 to 2007 and was most pronounced in the RR. These time series also illustrated that the SPF exhibited the highest PP and the lowest variation in Chla over the ocean colour record. This implies that the SPF provides a high and consistent supply of carbon to the benthos irrespective of fluctuations in the North Atlantic Oscillation.

Quantifying the impact of environmental variables upon catch-per-unit-effort of the Blue Shark Prionace glauca in the Western English Channel

The effect of environmental variables on blue shark Prionace glauca catch per unit effort (CPUE) in a recreational fishery in the western English Channel, between June and September 1998–2011, was quantified using generalized additive models (GAMs). Sea surface temperature (SST) explained 1·4% of GAM deviance, and highest CPUE occurred at 16·7° C, reflecting the optimal thermal preferences of this species. Surface chlorophyll a concentration (CHL) significantly affected CPUE and caused 27·5% of GAM deviance. Additionally, increasing CHL led to rising CPUE, probably due to higher productivity supporting greater prey biomass. The density of shelf-sea tidal mixing fronts explained 5% of GAM deviance, but was non-significant, with increasing front density negatively affecting CPUE. Time-lagged frontal density significantly affected CPUE, however, causing 12·6% of the deviance in a second GAM and displayed a positive correlation. This outcome suggested a delay between the evolution of frontal features and the subsequent accumulation of productivity and attraction of higher trophic level predators, such as P. glauca.

Seasonal shelf-sea front mapping using satellite ocean colour and temperature to support development of a marine protected area network

Front detection and aggregation techniques were applied to 300m resolution MERIS satellite ocean colour data for the first time, to describe frequently occurring shelf-sea fronts near to the Scottish coast. Medium resolution (1km) thermal and colour data have previously been used to analyse the distribution of surface fronts, though these cannot capture smaller frontal zones or those in close proximity to the coast, particularly where the coastline is convoluted. Seasonal frequent front maps, derived from both chlorophyll and SST data, revealed a number of key frontal zones, a subset of which were based on new insights into the sediment and plankton dynamics provided exclusively by the higher-resolution chlorophyll fronts. The methodology is described for applying colour and thermal front data to the task of identifying zones of ecological importance that could assist the process of defining marine protected areas. Each key frontal zone is analysed to describe its spatial and temporal extent and variability, and possible mechanisms. It is hoped that these tools can provide guidance on the dynamic habitats of marine fauna towards aspects of marine spatial planning and conservation.

Composite front maps for improved visibility of dynamic sea-surface features on cloudy SeaWiFS and AVHRR data

Novel techniques have been developed for increasing the value of cloud-affected sequences of Advanced Very High Resolution Radiometer (AVHRR) sea-surface temperature (SST) data and Sea-viewing Wide Field-of-view Sensor (SeaWiFS) ocean colour data for visualising dynamic physical and biological oceanic processes such as fronts, eddies and blooms. The proposed composite front map approach is to combine the location, strength and persistence of all fronts observed over several days into a single map, which allows intuitive interpretation of mesoscale structures. This method achieves a synoptic view without blurring dynamic features, an inherent problem with conventional time-averaging compositing methods. Objective validation confirms a significant improvement in feature visibility on composite maps compared to individual front maps. A further novel aspect is the automated detection of ocean colour fronts, correctly locating 96% of chlorophyll fronts in a test data set. A sizeable data set of 13,000 AVHRR and 1200 SeaWiFS scenes automatically processed using this technique is applied to the study of dynamic processes off the Iberian Peninsula such as mesoscale eddy generation, and many additional applications are identified. Front map animations provide a unique insight into the evolution of upwelling and eddies.

Quantifying the impact of environmental variables upon catch-per-unit-effort of the Blue Shark Prionace glauca in the Western English Channel

The effect of environmental variables on blue shark Prionace glauca catch per unit effort (CPUE) in a recreational fishery in the western English Channel, between June and September 1998–2011, was quantified using generalized additive models (GAMs). Sea surface temperature (SST) explained 1·4% of GAM deviance, and highest CPUE occurred at 16·7° C, reflecting the optimal thermal preferences of this species. Surface chlorophyll a concentration (CHL) significantly affected CPUE and caused 27·5% of GAM deviance. Additionally, increasing CHL led to rising CPUE, probably due to higher productivity supporting greater prey biomass. The density of shelf-sea tidal mixing fronts explained 5% of GAM deviance, but was non-significant, with increasing front density negatively affecting CPUE. Time-lagged frontal density significantly affected CPUE, however, causing 12·6% of the deviance in a second GAM and displayed a positive correlation. This outcome suggested a delay between the evolution of frontal features and the subsequent accumulation of productivity and attraction of higher trophic level predators, such as P. glauca.

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.

Satellite estimates of net community production indicate predominance of net autotrophy in the Atlantic Ocean

There is ongoing debate as to whether the oligotrophic ocean is predominantly net autotrophic and acts as a CO2 sink, or net heterotrophic and therefore acts as a CO2 source to the atmosphere. This quantification is challenging, both spatially and temporally, due to the sparseness of measurements. There has been a concerted effort to derive accurate estimates of phytoplankton photosynthesis and primary production from satellite data to fill these gaps; however there have been few satellite estimates of net community production (NCP). In this paper, we compare a number of empirical approaches to estimate NCP from satellite data with in vitro measurements of changes in dissolved O2 concentration at 295 stations in the N and S Atlantic Ocean (including the Antarctic), Greenland and Mediterranean Seas. Algorithms based on power laws between NCP and particulate organic carbon production (POC) derived from 14C uptake tend to overestimate NCP at negative values and underestimate at positive values. An algorithm that includes sea surface temperature (SST) in the power function of NCP and 14C POC has the lowest bias and root-mean square error compared with in vitro measured NCP and is the most accurate algorithm for the Atlantic Ocean. Nearly a 13 year time series of NCP was generated using this algorithm with SeaWiFS data to assess changes over time in different regions and in relation to climate variability. The North Atlantic subtropical and tropical Gyres (NATL) were predominantly net autotrophic from 1998 to 2010 except for boreal autumn/winter, suggesting that the northern hemisphere has remained a net sink for CO2 during this period. The South Atlantic subtropical Gyre (SATL) fluctuated from being net autotrophic in austral spring-summer, to net heterotrophic in austral autumn–winter. Recent decadal trends suggest that the SATL is becoming more of a CO2 source. Over the Atlantic basin, the percentage of satellite pixels with negative NCP was 27%, with the largest contributions from the NATL and SATL during boreal and austral autumn–winter, respectively. Variations in NCP in the northern and southern hemispheres were correlated with climate indices. Negative correlations between NCP and the multivariate ENSO index (MEI) occurred in the SATL, which explained up to 60% of the variability in NCP. Similarly there was a negative correlation between NCP and the North Atlantic Oscillation (NAO) in the Southern Sub-Tropical Convergence Zone (SSTC),which explained 90% of the variability. There were also positive correlations with NAO in the Canary Current Coastal Upwelling (CNRY) and Western Tropical Atlantic (WTRA)which explained 80% and 60% of the variability in each province, respectively. MEI and NAO seem to play a role in modifying phases of net autotrophy and heterotrophy in the Atlantic Ocean.