Assessing marine plankton community structure from long-term monitoring data with multivariate autoregressive (MAR) models: a comparison of fixed station versus spatially distributed sampling data

We examined how marine plankton interaction networks, as inferred by multivariate autoregressive (MAR) analysis of time-series, differ based on data collected at a fixed sampling location (L4 station in the Western English Channel) and four similar time-series prepared by averaging Continuous Plankton Recorder (CPR) datapoints in the region surrounding the fixed station. None of the plankton community structures suggested by the MAR models generated from the CPR datasets were well correlated with the MAR model for L4, but of the four CPR models, the one most closely resembling the L4 model was that for the CPR region nearest to L4. We infer that observation error and spatial variation in plankton community dynamics influenced the model performance for the CPR datasets. A modified MAR framework in which observation error and spatial variation are explicitly incorporated could allow the analysis to better handle the diverse time-series data collected in marine environments.

PCR from the CPR offers a historical perspective on marine population ecology.

The Continuous Plankton Recorder (CPR) survey has collected plankton samples from regular tracks across the world's oceans for almost 70 y. Over 299,000 spatially extensive CPR samples are archived and stored in buffered formalin. This CPR archive offers huge potential to study changes in marine communities using molecular data from a period when marine pollution, exploitation and global anthropogenic impact were much less pronounced. However, to harness the amount of data available within the CPR archive fully, it is necessary to improve techniques of larval identification, to genus and species preferably, and to obtain genetic information for historical studies of population ecology. To increase the potential of the CPR database this paper describes the first extraction, amplification by the polymerase chain reaction and utilization of a DNA sequence (mitochondrial 16S rDNA) from a CPR sample, a formalin fixed larval sandeel.

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species-based Dynamic Bioclimate Envelope Model (DBEM) with a size-based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness-of-fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter-specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.

Synergistic effects of climate and fishing in a marine ecosystem

Current climate change and overfishing are affecting the productivity and structure of marine ecosystems. This situation is unprecedented for the marine biosphere and it is essential to understand the mechanisms and pathways by which ecosystems respond. We report that climate change and overfishing are likely to be responsible for a rapid restructuring of a highly productive marine ecosystem with effects throughout the pelagos and the benthos. In the mid-1980s, climate change, consequent modifications in the North Sea plankton, and fishing, all reduced North Sea cod recruitment. In this region, production of many benthic species respond positively and immediately to temperature. Analysis of a long-term, spatially extensive biological (plankton and cod) and physical (sea surface temperature) dataset suggests that synchronous changes in cod numbers and sea temperature have established an extensive trophic cascade favoring lower trophic level groups over economic fisheries. A proliferation of jellyfish that we detect may signal the climax of these changes. This modified North Sea ecology may provide a clear indication of the synergistic consequences of coincident climate change and overfishing. The extent of the ecosystem restructuring that has occurred in the North Sea suggests we are unlikely to reverse current climate and human-induced effects through ecosystem resource management in the short term. Rather, we should understand and adapt to new ecological regimes. This implies that fisheries management policies will have to be fully integrated with the ecological consequences of climate change to prevent a similar collapse in an exploited marine ecosystem elsewhere.

Spatial variability of the plankton trophic interaction in the North Sea: a new feature after the early 1970s

Traditionally, marine ecosystem structure was thought to be bottom-up controlled. In recent years, a number of studies have highlighted the importance of top-down regulation. Evidence is accumulating that the type of trophic forcing varies temporally and spatially, and an integrated view – considering the interplay of both types of control – is emerging. Correlations between time series spanning several decades of the abundances of adjacent trophic levels are conventionally used to assess the type of control: bottom-up if positive or top-down if this is negative. This approach implies averaging periods which might show time-varying dynamics and therefore can hide part of this temporal variability. Using spatially referenced plankton information extracted from the Continuous Plankton Recorder, this study addresses the potential dynamic character of the trophic structure at the planktonic level in the North Sea by assessing its variation over both temporal and spatial scales. Our results show that until the early-1970s a bottom-up control characterized the base of the food web across the whole North Sea, with diatoms having a positive and homogeneous effect on zooplankton filter-feeders. Afterwards, different regional trophic dynamics were observed, in particular a negative relationship between total phytoplankton and zooplankton was detected off the west coast of Norway and the Skagerrak as opposed to a positive one in the southern reaches. Our results suggest that after the early 1970s diatoms remained the main food source for zooplankton filter-feeders east of Orkney–Shetland and off Scotland, while in the east, from the Norwegian Trench to the German Bight, filter-feeders were mainly sustained by dinoflagellates.

Spatial patterns and trends in abundance of larval sandeels in the North Sea: 1950-2005

Early recruitment indices based on larval fish data from the Continuous Plankton Recorder (CPR) have the potential to inform stock assessments of Ammodytes marinus in the North Sea. We evaluate whether the CPR data are reliable for sandeel larvae. Spatially, CPR larval data were comparable with catches by dedicated larval samplers (Gulf and bongo nets) during ICES coordinated surveys in 2004 and 2009. ICES data are also used to explore environmental influences on sandeel distributions. Temporally, CPR data correlate with larval data from plankton surveys off Stonehaven (1999–2005), with sandeel 0-group trawl data at the east Fair Isle ground (since 1984), and with recruitment data (since 1983) for the Dogger Banks stock assessment area. Therefore, CPR data may provide an early recruit index of relative abundance for the Dogger Banks assessment area, where the majority of the commercial catch of A. marinus is taken, and the Wee Bankie area that is particularly important for seabird foraging. While warm conditions may stimulate the production of sandeel larvae, their natural mortality is typically greater, in the Dogger Banks and Wadden Sea areas, when the larvae are hatched in warm years and/or with abundant 1-year-old sandeel that are likely to be cannibalistic.

Covariance among North Sea ecosystem state indicators during the past 50 years - contrasts between coastal and open waters

Regime shift and principal component analysis of a spatially disaggregated database capturing time-series of climatic, nutrient and plankton variables in the North Sea revealed considerable covariance between groups of ecosystem indicators. Plankton and climate time-series span the period 1958–2003, those of nutrients start in 1980. In both regions, the period from 1989 to 2001 identified in principal component 1 had warmer surface waters, higher Atlantic inflow and stronger winds, than the periods before or after. However, it was preceded by a regime shift in both open (PC2) and coastal (PC3) waters during 1977 towards more hours of sunlight and higher water temperature, which lasted until 1997. The relative influence of nutrient availability and climatic forcing differed between open and coastal North Sea regions. Inter-annual variability in phytoplankton dynamics of the open North Sea was primarily regulated by climatic forcing, specifically by sea surface temperature, Atlantic inflow and co-varying wind stress and NAO. Coastal phytoplankton variability, however, was regulated by insolation and sea surface temperature, as well as Si availability, but not by N or P. Regime shifts in principal components of hydrographic and climatic variables (explaining 55 and 61% of the variance in coastal and open water variables) were detected using Rodionov's sequential t-test. These shifts in hydroclimatic variables which occurred around 1977, 1989, 1997 and 2001, were synchronized in open and coastal waters, and were tracked by open water chlorophyll and copepods, but not by coastal plankton. North–central–south or open-coastal spatial breakdowns of the North Sea explained similar amounts of variability in most ecosystem indicators with the exception of diatom abundance and chlorophyll concentration, which were clearly better explained using the open-coastal configuration.

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.

Plankton indicators and ocean observing systems: support to the marine ecosystem state assessment

Ecological indicators are used extensively as tools to manage environmental resources. In the oceans, indicators of plankton can be measured using a variety of observing systems including: mooring stations, ships, autonomous floats and ocean colour remote sensing. Given the broad range of temporal and spatial sampling resolutions of these different observing systems, as well as discrepancies in measurements obtained from different sensors, the estimation and interpretation of plankton indicators can present significant challenges. To provide support to the assessment of the state of the marine ecosystem, we propose a suite of plankton indicators and subsequently classify them in an ecological framework that characterizes key attributes of the ecosystem. We present two case studies dealing with plankton indicators of biomass, size structure and phenology, estimated using the most spatially extensive and longest in situ and remote-sensing observations. Discussion of these studies illustrates how some of the challenges in estimating and interpreting plankton indicators may be addressed by using for example relative measurement thresholds, interpolation procedures and delineation of biogeochemical provinces. We demonstrate that one of the benefits attained, when analyzing a suite of plankton indicators classified in an ecological framework, is the elucidation of non-trivial changes in composition, structure and functioning of the marine ecosystem.

Quantifying ecosystem-level consequences of ocean acidification and warming by scaling up individual level responses of a predator snail and its trophic interactions

Understanding long‐term, ecosystem‐level impacts of climate change is challenging because experimental research frequently focuses on short‐term, individual‐level impacts in isolation. We address this shortcoming first through an interdisciplinary ensemble of novel experimental techniques to investigate the impacts of 14‐month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterize a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual‐level responses, while acidification had a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual‐level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large‐scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local‐environmental conditions and resource availability. Such changes in macro‐scale distributions cannot be predicted by investigating individual‐level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long‐term, multiscale responses to multiple stressors, in an ecosystem context.

Scaling up experimental ocean acidification and warming research: from individuals to the ecosystem

Understanding long-term, ecosystem-level impacts of climate change is challenging because experimental research frequently focuses on short-term, individual-level impacts in isolation. We address this shortcoming first through an inter-disciplinary ensemble of novel experimental techniques to investigate the impacts of 14-month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterise a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual level responses, while acidification has a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large-scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local environmental conditions and resource availability. Such changes in macro-scale distributions cannot be predicted by investigating individual level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long-term, multi-scale responses to multiple stressors, in an ecosystem context.

Microsporidia-nematode associations in methane seeps reveal basal fungal parasitism in the deep sea

The deep sea is Earth’s largest habitat but little is known about the nature of deep-sea parasitism. In contrast to a few characterized cases of bacterial and protistan parasites, the existence and biological significance of deep-sea parasitic fungi is yet to be understood. Here we report the discovery of a fungus-related parasitic microsporidium, Nematocenator marisprofundi n. gen. n. sp. that infects benthic nematodes at Pacific Ocean methane seeps on the Pacific Ocean floor. This infection is species-specific and has been temporally and spatially stable over two years of sampling, indicating an ecologically consistent host-parasite interaction. A high distribution of spores in the reproductive tracts of infected males and females and their absence from host nematodes’ intestines suggests a sexual transmission strategy in contrast to the fecal-oral transmission of most microsporidia. N. marisprofundi targets the host’s body wall muscles causing cell lysis, and in severe infection even muscle filament degradation. Phylogenetic analyses placed N. marisprofundi in a novel and basal clade not closely related to any described microsporidia clade, suggesting either that microsporidia-nematode parasitism occurred early in microsporidia evolution or that host specialization occurred late in an ancient deep-sea microsporidian lineage. Our findings reveal that methane seeps support complex ecosystems involving interkingdom interactions between bacteria, nematodes, and parasitic fungi and that microsporidia parasitism exists also in the deep sea biosphere.

Foraging success of biological Levy flights recorded in situ

It is an open question how animals find food in dynamic natural environments where they possess little or no knowledge of where resources are located. Foraging theory predicts that in environments with sparsely distributed target resources, where forager knowledge about resources’ locations is incomplete, Lévy flight movements optimize the success of random searches. However, the putative success of Lévy foraging has been demonstrated only in model simulations. Here, we use high-temporal-resolution Global Positioning System (GPS) tracking of wandering (Diomedea exulans) and black-browed albatrosses (Thalassarche melanophrys) with simultaneous recording of prey captures, to show that both species exhibit Lévy and Brownian movement patterns. We find that total prey masses captured by wandering albatrosses during Lévy movements exceed daily energy requirements by nearly fourfold, and approached yields by Brownian movements in other habitats. These results, together with our reanalysis of previously published albatross data, overturn the notion that albatrosses do not exhibit Lévy patterns during foraging, and demonstrate that Lévy flights of predators in dynamic natural environments present a beneficial alternative strategy to simple, spatially intensive behaviors. Our findings add support to the possibility that biological Lévy flight may have naturally evolved as a search strategy in response to sparse resources and scant information.

Simple shade plots aid better long-term choices of data pre-treatment in multivariate assemblage studies

Shade plots, simple visual representations of abundance matrices from multivariate species assemblage studies, are shown to be an effective aid in choosing an overall transformation (or other pre-treatment) of quantitative data for long-term use, striking an appropriate balance between dominant and less abundant taxa in ensuing resemblance-based multivariate analyses. Though the exposition is entirely general and applicable to all community studies, detailed illustrations of the comparative power and interpretative possibilities of shade plots are given in the case of two estuarine assemblage studies in south-western Australia: (a) macrobenthos in the upper Swan Estuary over a two-year period covering a highly significant precipitation event for the Perth area; and (b) a wide-scale spatial study of the nearshore fish fauna from five divergent estuaries. The utility of transformations of intermediate severity is again demonstrated and, with greater novelty, the potential importance seen of further mild transformation of all data after differential down-weighting (dispersion weighting) of spatially clumped' or schooled' species. Among the new techniques utilized is a two-way form of the RELATE test, which demonstrates linking of assemblage structure (fish) to continuous environmental variables (water quality), having removed a categorical factor (estuary differences). Re-orderings of sample and species axes in the associated shade plots are seen to provide transparent explanations at the species level for such continuous multivariate patterns.