Scale-dependent response diversity of seabirds to prey in the North Sea

Functional response diversity is defined as the diversity of responses to environmental change among species that contribute to the same ecosystem function. Because different ecological processes dominate on different spatial and temporal scales, response diversity is likely to be scale dependent. Using three extensive data sets on seabirds, pelagic fish, and zooplankton, we investigate the strength and diversity in the response of seabirds to prey in the North Sea over three scales of ecological organization. Two-stage analyses were used to partition the variance in the abundance of predators and prey among the different scales of investigation: variation from year to year, variation among habitats, and variation on the local patch scale. On the year-to-year scale, we found a strong and synchronous response of seabirds to the abundance of prey, resulting in low response diversity. Conversely, as different seabird species were found in habitats dominated by different prey species, we found a high diversity in the response of seabirds to prey on the habitat scale. Finally, on the local patch scale, seabirds were organized in multispecies patches. These patches were weakly associated with patches of prey, resulting in a weak response strength and a low response diversity. We suggest that ecological similarities among seabird species resulted in low response diversity on the year-to-year scale. On the habitat scale, we suggest that high response diversity was due to interspecific competition and niche segregation among seabird species. On the local patch scale, we suggest that facilitation with respect to the detection and accessibility of prey patches resulted in overlapping distribution of seabirds but weak associations with prey. The observed scale dependencies in response strength and diversity have implications for how the seabird community will respond to different environmental disturbances.

Rapid biogeographical plankton shifts in the North Atlantic Ocean

Large-scale biogeographical changes in the biodiversity of a key zooplankton group (calanoid copepods) were detected in the north-eastern part of the North Atlantic Ocean and its adjacent seas over the period 1960–1999. These findings provided key empirical evidence for climate change impacts on marine ecosystems at the regional to oceanic scale. Since 1999, global temperatures have continued to rise in the region. Here, we extend the analysis to the period 1958–2005 using all calanoid copepod species assemblages (nine species assemblages based on an analysis including a total of 108 calanoid species or taxa) and show that this phenomenon has been reinforced in all regions. Our study reveals that the biodiversity of calanoid copepods are responding quickly to sea surface temperature (SST) rise by moving geographically northward at a rapid rate up to about 23.16 km yr−1. Our analysis suggests that nearly half of the increase in sea temperature in the northeast Atlantic and adjacent seas is related to global temperature rises (46.35% of the total variance of temperature) while changes in both natural modes of atmospheric and oceanic circulation explain 26.45% of the total variance of temperature. Although some SST isotherms have moved northwards by an average rate of up to 21.75 km yr−1 (e.g. the North Sea), their movement cannot fully quantify all species assemblage shifts. Furthermore, the observed rates of biogeographical movements are far greater than those observed in the terrestrial realm. Here, we discuss the processes that may explain such a discrepancy and suggest that the differences are mainly explained by the fluid nature of the pelagic domain, the life cycle of the zooplankton and the lesser anthropogenic influence (e.g. exploitation, habitat fragmentation) on these organisms. We also hypothesize that despite changes in the path and intensity of the oceanic currents that may modify quickly and greatly pelagic zooplankton species, these organisms may reflect better the current impact of climate warming on ecosystems as terrestrial organisms are likely to significantly lag the current impact of climate change.

Changes in marine dinoflagellate and diatom abundance under climate change

Marine diatoms and dinoflagellates play a variety of key ecosystem roles as important primary producers (diatoms and some dinoflagellates) and grazers (some dinoflagellates). Additionally some are harmful algal bloom (HAB) species and there is widespread concern that HAB species may be increasing accompanied by major negative socio-economic impacts, including threats to human health and marine harvesting1, 2. Using 92,263 samples from the Continuous Plankton Recorder survey, we generated a 50-year (1960–2009) time series of diatom and dinoflagellate occurrence in the northeast Atlantic and North Sea. Dinoflagellates, including both HAB taxa (for example, Prorocentrum spp.) and non-HAB taxa (for example, Ceratium furca), have declined in abundance, particularly since 2006. In contrast, diatom abundance has not shown this decline with some common diatoms, including both HAB (for example, Pseudo-nitzschia spp.) and non-HAB (for example, Thalassiosira spp.) taxa, increasing in abundance. Overall these changes have led to a marked increase in the relative abundance of diatoms versus dinoflagellates. Our analyses, including Granger tests to identify criteria of causality, indicate that this switch is driven by an interaction effect of both increasing sea surface temperatures combined with increasingly windy conditions in summer.

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.

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.

The influence of ocean acidification on nitrogen regeneration and nitrous oxide production in the North-West European shelf sea

The assimilation and regeneration of dissolved inorganic nitrogen, and the concentration of N2O, was investigated at stations located in the NW European shelf sea during June/July 2011. These observational measurements within the photic zone demonstrated the simultaneous regeneration and assimilation of NH4+, NO2− and NO3−. NH4+ was assimilated at 1.82–49.12 nmol N L−1 h−1 and regenerated at 3.46–14.60 nmol N L−1 h−1; NO2− was assimilated at 0–2.08 nmol N L−1 h−1 and regenerated at 0.01–1.85 nmol N L−1 h−1; NO3− was assimilated at 0.67–18.75 nmol N L−1 h−1 and regenerated at 0.05–28.97 nmol N L−1 h−1. Observations implied that these processes were closely coupled at the regional scale and nitrogen recycling played an important role in sustaining phytoplankton growth during the summer. The [N2O], measured in water column profiles, was 10.13 ± 1.11 nmol L−1 and did not strongly diverge from atmospheric equilibrium indicating that sampled marine regions where neither a strong source nor sink of N2O to the atmosphere. Multivariate analysis of data describing water column biogeochemistry and its links to N-cycling activity failed to explain the observed variance in rates of N-regeneration and N-assimilation, possibly due to the limited number of process rate observations. In the surface waters of 5 further stations, Ocean Acidification (OA) bioassay experiments were conducted to investigate the response of NH4+ oxidising and regenerating organisms to simulated OA conditions, including the implications for [N2O]. Multivariate analysis was undertaken which considered the complete bioassay dataset of measured variables describing changes in N-regeneration rate, [N2O] and the biogeochemical composition of seawater. While anticipating biogeochemical differences between locations, we aimed to test the hypothesis that the underlying mechanism through which pelagic N-regeneration responded to simulated OA conditions was independent of location and that a mechanistic understanding of how NH4+ oxidation, NH4+ regeneration and N2O production responded to OA could be developed. Results indicated that N-regeneration process responses to OA treatments were location specific; no mechanistic understanding of how N-regeneration processes respond to OA in the surface ocean of the NW European shelf sea could be developed.

The influence of ocean acidification on nitrogen regeneration and nitrous oxide production in the North-West European shelf sea

The assimilation and regeneration of dissolved inorganic nitrogen, and the concentration of N2O, was investigated at stations located in the NW European shelf sea during June/July 2011. These observational measurements within the photic zone demonstrated the simultaneous regeneration and assimilation of NH4+, NO2− and NO3−. NH4+ was assimilated at 1.82–49.12 nmol N L−1 h−1 and regenerated at 3.46–14.60 nmol N L−1 h−1; NO2− was assimilated at 0–2.08 nmol N L−1 h−1 and regenerated at 0.01–1.85 nmol N L−1 h−1; NO3− was assimilated at 0.67–18.75 nmol N L−1 h−1 and regenerated at 0.05–28.97 nmol N L−1 h−1. Observations implied that these processes were closely coupled at the regional scale and nitrogen recycling played an important role in sustaining phytoplankton growth during the summer. The [N2O], measured in water column profiles, was 10.13 ± 1.11 nmol L−1 and did not strongly diverge from atmospheric equilibrium indicating that sampled marine regions where neither a strong source nor sink of N2O to the atmosphere. Multivariate analysis of data describing water column biogeochemistry and its links to N-cycling activity failed to explain the observed variance in rates of N-regeneration and N-assimilation, possibly due to the limited number of process rate observations. In the surface waters of 5 further stations, Ocean Acidification (OA) bioassay experiments were conducted to investigate the response of NH4+ oxidising and regenerating organisms to simulated OA conditions, including the implications for [N2O]. Multivariate analysis was undertaken which considered the complete bioassay dataset of measured variables describing changes in N-regeneration rate, [N2O] and the biogeochemical composition of seawater. While anticipating biogeochemical differences between locations, we aimed to test the hypothesis that the underlying mechanism through which pelagic N-regeneration responded to simulated OA conditions was independent of location and that a mechanistic understanding of how NH4+ oxidation, NH4+ regeneration and N2O production responded to OA could be developed. Results indicated that N-regeneration process responses to OA treatments were location specific; no mechanistic understanding of how N-regeneration processes respond to OA in the surface ocean of the NW European shelf sea could be developed.

The influence of ocean acidification on nitrogen regeneration and nitrous oxide production in the northwest European shelf sea

The assimilation and regeneration of dissolved inorganic nitrogen, and the concentration of N₂O, was investigated at stations located in the NW European shelf sea during June/July 2011. These observational measurements within the photic zone demonstrated the simultaneous regeneration and assimilation of NH₄⁺, NO₂⁻ and NO₃⁻. NH₄⁺ was assimilated at 1.82–49.12 nmol N L⁻¹ h⁻¹ and regenerated at 3.46–14.60 nmol N L⁻¹ h⁻¹; NO₂^- was assimilated at 0–2.08 nmol N L⁻¹ h⁻¹ and regenerated at 0.01–1.85 nmol N L⁻¹ h⁻¹; NO₃⁻ was assimilated at 0.67–18.75 nmol N L⁻¹ h⁻¹ and regenerated at 0.05–28.97 nmol N L⁻¹ h⁻¹. Observations implied that these processes were closely coupled at the regional scale and that nitrogen recycling played an important role in sustaining phytoplankton growth during the summer. The [N₂O], measured in water column profiles, was 10.13 ± 1.11 nmol L⁻¹ and did not strongly diverge from atmospheric equilibrium indicating that sampled marine regions were neither a strong source nor sink of N₂O to the atmosphere. Multivariate analysis of data describing water column biogeochemistry and its links to N-cycling activity failed to explain the observed variance in rates of N-regeneration and N-assimilation, possibly due to the limited number of process rate observations. In the surface waters of five further stations, ocean acidification (OA) bioassay experiments were conducted to investigate the response of NH₄⁺ oxidising and regenerating organisms to simulated OA conditions, including the implications for [N₂O]. Multivariate analysis was undertaken which considered the complete bioassay data set of measured variables describing changes in N-regeneration rate, [N₂O] and the biogeochemical composition of seawater. While anticipating biogeochemical differences between locations, we aimed to test the hypothesis that the underlying mechanism through which pelagic N-regeneration responded to simulated OA conditions was independent of location. Our objective was to develop a mechanistic understanding of how NH₄⁺ regeneration, NH₄⁺ oxidation and N₂O production responded to OA. Results indicated that N-regeneration process responses to OA treatments were location specific; no mechanistic understanding of how N-regeneration processes respond to OA in the surface ocean of the NW European shelf sea could be developed.

Late winter under ice pelagic microbial communities in the high Arctic Ocean and the impact of short-term exposure to elevated CO2 levels

Polar Oceans are natural CO2 sinks because of the enhanced solubility of CO2 in cold water. The Arctic Ocean is at additional risk of accelerated ocean acidification (OA) because of freshwater inputs from sea ice and rivers, which influence the carbonate system. Winter conditions in the Arctic are of interest because of both cold temperatures and limited CO2 venting to the atmosphere when sea ice is present. Earlier OA experiments on Arctic microbial communities conducted in the absence of ice cover, hinted at shifts in taxa dominance and diversity under lowered pH. The Catlin Arctic Survey provided an opportunity to conduct in situ, under-ice, OA experiments during late Arctic winter. Seawater was collected from under the sea ice off Ellef Ringnes Island, and communities were exposed to three CO2 levels for 6 days. Phylogenetic diversity was greater in the attached fraction compared to the free-living fraction in situ, in the controls and in the treatments. The dominant taxa in all cases were Gammaproteobacteria but acidification had little effect compared to the effects of containment. Phylogenetic net relatedness indices suggested that acidification may have decreased the diversity within some bacterial orders, but overall there was no clear trend. Within the experimental communities, alkalinity best explained the variance among samples and replicates, suggesting subtle changes in the carbonate system need to be considered in such experiments. We conclude that under ice communities have the capacity to respond either by selection or phenotypic plasticity to heightened CO2 levels over the short term.

Spatio-temporal Bayesian network models with latent variables for revealing trophic dynamics and functional networks in fisheries ecology

Ecosystems consist of complex dynamic interactions among species and the environment, the understanding of which has implications for predicting the environmental response to changes in climate and biodiversity. However, with the recent adoption of more explorative tools, like Bayesian networks, in predictive ecology, few assumptions can be made about the data and complex, spatially varying interactions can be recovered from collected field data. In this study, we compare Bayesian network modelling approaches accounting for latent effects to reveal species dynamics for 7 geographically and temporally varied areas within the North Sea. We also apply structure learning techniques to identify functional relationships such as prey–predator between trophic groups of species that vary across space and time. We examine if the use of a general hidden variable can reflect overall changes in the trophic dynamics of each spatial system and whether the inclusion of a specific hidden variable can model unmeasured group of species. The general hidden variable appears to capture changes in the variance of different groups of species biomass. Models that include both general and specific hidden variables resulted in identifying similarity with the underlying food web dynamics and modelling spatial unmeasured effect. We predict the biomass of the trophic groups and find that predictive accuracy varies with the models' features and across the different spatial areas thus proposing a model that allows for spatial autocorrelation and two hidden variables. Our proposed model was able to produce novel insights on this ecosystem's dynamics and ecological interactions mainly because we account for the heterogeneous nature of the driving factors within each area and their changes over time. Our findings demonstrate that accounting for additional sources of variation, by combining structure learning from data and experts' knowledge in the model architecture, has the potential for gaining deeper insights into the structure and stability of ecosystems. Finally, we were able to discover meaningful functional networks that were spatially and temporally differentiated with the particular mechanisms varying from trophic associations through interactions with climate and commercial fisheries.