DISCOVERY 2010: Spatial and temporal variability in a dynamic polar ecosystem

The Scotia Sea has been a focus of biological- and physical oceanographic study since the Discovery expeditions in the early 1900s. It is a physically energetic region with some of the highest levels of productivity in the Southern Ocean. It is also a region within which there have been greater than average levels of change in upper water column temperature. We describe the results of three cruises transecting the central Scotia Sea from south to north in consecutive years and covering spring, summer and autumn periods. We also report on some community level syntheses using both current-day and historical data from this region. A wide range of parameters were measured during the field campaigns, covering the physical oceanography of the region, air–sea CO2 fluxes, macro- and micronutrient concentrations, the composition and biomass of the nano-, micro- and mesoplankton communities, and the distribution and biomass of Antarctic krill and mesopelagic fish. Process studies examined the effect of iron-stress on the physiology of primary producers, reproduction and egestion in Antarctic krill and the transfer of stable isotopes between trophic layers, from primary consumers up to birds and seals. Community level syntheses included an examination of the biomass-spectra, food-web modelling, spatial analysis of multiple trophic layers and historical species distributions. The spatial analyses in particular identified two distinct community types: a northern warmer water community and a southern cold community, their boundary being broadly consistent with the position of the Southern Antarctic Circumpolar Current Front (SACCF). Temperature and ice cover appeared to be the dominant, over-riding factors in driving this pattern. Extensive phytoplankton blooms were a major feature of the surveys, and were persistent in areas such as South Georgia. In situ and bioassay measurements emphasised the important role of iron inputs as facilitators of these blooms. Based on seasonal DIC deficits, the South Georgia bloom was found to contain the strongest seasonal carbon uptake in the ice-free zone of the Southern Ocean. The surveys also encountered low-production, iron-limited regions, a situation more typical of the wider Southern Ocean. The response of primary and secondary consumers to spatial and temporal heterogeneity in production was complex. Many of the life-cycles of small pelagic organisms showed a close coupling to the seasonal cycle of food availability. For instance, Antarctic krill showed a dependence on early, non-ice-associated blooms to facilitate early reproduction. Strategies to buffer against environmental variability were also examined, such as the prevalence of multiyear life-cycles and variability in energy storage levels. Such traits were seen to influence the way in which Scotia Sea communities were structured, with biomass levels in the larger size classes being higher than in other ocean regions. Seasonal development also altered trophic function, with the trophic level of higher predators increasing through the course of the year as additional predator-prey interactions emerged in the lower trophic levels. Finally, our studies re-emphasised the role that the simple phytoplankton-krill-higher predator food chain plays in this Southern Ocean region, particularly south of the SACCF. To the north, alternative food chains, such as those involving copepods, macrozooplankton and mesopelagic fish, were increasingly important. Continued ocean warming in this region is likely to increase the prevalence of such alternative such food chains with Antarctic krill predicted to move southwards.

Temporal changes in plankton of the North Sea: community shifts and environmental drivers

This paper analyses long-term and seasonal changes in the North Sea plankton community during the period 1970 to 2008. Based on Continuous Plankton Recorder (CPR) data covering 38 yr, major changes in both phytoplankton and zooplankton abundance and community structure were identified. Regime changes were detected around 1978, 1989 and 1998. The first 2 changes have been discussed in the literature and are defined as a cold episodic event (1978) and a regime shift towards a warm dynamic regime (1989). The effect of these 2 regime changes on plankton indicators was assessed and checked against previous studies. The 1998 change represents a shift in the abundance and seasonal patterns of dinoflagellates and the dominant zooplankton group, the neritic copepods. Furthermore, environmental factors such as air temperature, wind speed and the North Atlantic water inflow were identified as potential drivers of change in seasonal patterns, and the most-likely environmental causes for detected changes were assessed. We suggest that a change in the balance of dissolved nutrients driven by these environmental factors was the cause of the latest change in plankton community structure, which in turn could have affected the North Sea fish community.

Acquisition of continuous plankton recorder data

The Continuous Plankton Recorder has been deployed on a seasonal basis in the north Pacific since 2000, accumulating a database of abundance measurements for over 290 planktonic taxa in over 3,500 processed samples. There is an additional archive of over 10,000 samples available for further analyses. Exxon Valdez Oil Spill Trustee Council financial support has contributed to about half of this tally, through four projects funded since 2002. Time series of zooplankton variables for sub-regions of the survey area are presented together with abstracts of eight papers published using data from these projects. The time series covers a period when the dominant climate signal in the north Pacific, the Pacific Decadal Oscillation (PDO), switched with unusual frequency between warm/positive states (pre-1999 and 2003-2006) and cool/negative states (1999-2002 and 2007). The CPR data suggest that cool negative years show higher biomass on the shelf and lower biomass in the open ocean, while the reverse is true in warm (PDO positive) years with lower shelf biomass (except 2005) and higher oceanic biomass. In addition, there was a delay in plankton increase on the Alaskan shelf in the colder spring of 2007, compared to the warmer springs of the preceding years. In warm years, smaller species of copepods which lack lipid reserves are also more common. Availability of the zooplankton prey to higher trophic levels (including those that society values highly) is therefore dependent on the timing of increase and peak abundance, ease of capture and nutritional value. Previously published studies using these data highlight the wide-ranging applicability of CPR data and include collaborative studies on; phenology in the key copepod species Neocalanus plumchrus, descriptions of distributions of decapod larvae and euphausiid species, the effects of hydrographic features such as mesoscale eddies and the North Pacific Current on plankton populations and a molecularbased investigation of macro-scale population structure in N. cristatus. The future funding situation is uncertain but the value of the data and studies so far accumulated is considerable and sets a strong foundation for further studies on plankton dynamics and interactions with higher trophic levels in the northern Gulf of Alaska.

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.

Decadal changes in climate and ecosystems in the North Atlantic Ocean and adjacent seas

Climate change is unambiguous and its effects are clearly detected in all functional units of the Earth system. This study presents new analyses of sea-surface temperature changes and show that climate change is affecting ecosystems of the North Atlantic. Changes are seen from phytoplankton to zooplankton to fish and are modifying the dominance of species and the structure, the diversity and the functioning of marine ecosystems. Changes also range from phenological to biogeographical shifts and have involved in some regions of the Atlantic abrupt ecosystem shifts. These alterations reflect a response of pelagic ecosystems to a warmer temperature regime. Mechanisms are complex because they are nonlinear exhibiting tipping points and varying in space and time. Sensitivity of organisms to temperature changes is high, implicating that a small temperature modification can have sustained ecosystem effects. Implications of these changes for biogeochemical cycles are discussed. Two observed changes detected in the North Sea that could have opposite effects on carbon cycle are discussed. Increase in phytoplankton, as inferred from the phytoplankton colour index derived from the Continuous Plankton Recorder (CPR) survey, has been detected in the North Sea. This pattern has been accompanied by a reduction in the abundance of the herbivorous species Calanus finmarchicus. This might have reduced the grazing pressure and increase diatomaceous ‘fluff’, therefore carbon export in the North Sea. Therefore, it could be argued that the biological carbon pump might increase in this region with sea warming. In the meantime, however, the mean size of organisms (calanoid copepods) has dropped. Such changes have implications for the turnover time of biogenic carbon in plankton organisms and the mean residence time of particulate carbon they produce. The system characterising the warmer period is more based on recycling and less on export. The increase in the minimum turnover time indicates an increase in the ecosystem metabolism, which can be considered as a response of the pelagic ecosystems to climate warming. This phenomenon could reduce carbon export. These two opposite patterns of change are examples of the diversity of mechanisms and pathways the ecosystems may exhibit with climate change. Oversimplification of current biogeochemical models, often due to lack of data and biological understanding, could lead to wrong projection on the direction ecosystems and therefore some biogeochemical cycles might take in a warmer world.

Marine biodiversity, ecosystem functioning and carbon cycles

Although recent studies suggest that climate change may substantially accelerate the rate of species loss in the biosphere, only a few studies have focused on the potential consequences of a spatial reorganization of biodiversity with global warming. Here, we show a pronounced latitudinal increase in phytoplanktonic and zooplanktonic biodiversity in the extratropical North Atlantic Ocean in recent decades. We also show that this rise in biodiversity paralleled a decrease in the mean size of zooplanktonic copepods and that the reorganization of the planktonic ecosystem toward dominance by smaller organisms may influence the networks in which carbon flows, with negative effects on the downward biological carbon pump and demersal Atlantic cod (Gadus morhua). Our study suggests that, contrary to the usual interpretation of increasing biodiversity being a positive emergent property promoting the stability/resilience of ecosystems, the parallel decrease in sizes of planktonic organisms could be viewed in the North Atlantic as reducing some of the services provided by marine ecosystems to humans.

Recent trends in plankton and upwelling intensity off Galicia (NW Spain)

The mean intensity of the NE Atlantic upwelling system at its northern limit (Galicia, NW Spain) decreased during the last 40 years. At the same time, warming of surface waters was detected. Plankton biomass and composition are expected to reflect such changes when integrated over large time and space scales. In this study, biomass, abundance and species composition of phyto- and zooplankton were analysed to search for significant patterns of annual change and relations with upwelling intensity. Regionally integrated, mostly offshore, data were obtained from the Continuous Plankton Recorder (since 1958) whereas coastal data from Vigo and A Coruña came from the Radiales program (since 1987). No significant trends were found in phytoplankton biomass at either regional or local scales. However, there was a significant decrease in diatom abundance at regional scales and also of large species at local scales. Zooplankton abundance (mainly copepods) significantly decreased offshore but increased near the coast. Biomass of zooplankton also increased near the coast, with the fastest rates in the south. Warm-water species, like Temora stylifera, were increasingly abundant at both regional and local scales. Significant correlations between upwelling intensity and plankton suggest that climatic effects were delayed for several years. Our results indicate that the effects of large scale climatic trends on plankton communities are being effectively modulated within the pelagic ecosystem in this upwelling region.

Comparing copepod time-series in the north of Spain: Spatial autocorrelation of community composition

Four time-series of copepod species biomass in the north of Spain were contrasted to demonstrate spatial autocorrelation of local communities and their responses to short-term local and regional variability in oceanographic conditions. The series represented coastal and oceanic environments along a marked gradient of influence of seasonal upwelling from Galicia to the Mar Cantábrico (S Bay of Biscay), and each one included at least 10 years of continuous data collected at monthly frequency. Community composition (i.e. species number and diversity) was very consistent through the region, but local variations in the presence of new species and the relative proportions of common species allowed for the characterisation of the response to the environment at each site. Small-sized species were more frequent near the coast. A few species, however, captured the main patterns of variability in all series. Calanus helgolandicus and Acartia (mainly Acartia clausi) were generally the main contributors to total biomass, while other species as Paracalanus parvus and Clausocalanus spp. were important only at some locations. Most copepod indices were positively correlated with upwelling, either considering the whole community (biomass, species richness and diversity) or individual species, but only in the coastal series analysed since 1991. Copepods in the nearby ocean, however, showed negative correlations with upwelling in the period 1960–1986. The effects of upwelling may have been modulated by local factors, as showed by the increases in biomass, number of species and diversity in associations with increases in sea surface temperature in Galicia, while in the Mar Cantábrico only the warming-tolerant species increased and those typical of upwelling decreased. Density stratification of the water column was associated with decreases in total copepod biomass in Galicia, while it favoured the increase in species richness in the Mar Cantábrico. Nearly all significant responses of copepods to environmental variability were delayed by up to 5 months, showing the importance of considering time-lags in the analysis of temporal responses of zooplankton.

Diet of auklet chicks in the Aleutian Islands, Alaska: similarity among islands, interspecies overlap, and relationships to ocean climate

Seabirds are effective samplers of the marine environment, and can be used to measure resource partitioning among species and sites via food loads destined for chicks. We examined the composition, overlap, and relationships to changing climate and oceanography of 3,216 food loads from Least, Crested, and Whiskered Auklets (Aethia pusilla, A. cristatella, A. pygmaea) breeding in Alaska during 1994–2006. Meals comprised calanoid copepods (Neocalanus spp.) and euphausiids (Thysanoessa spp.) that reflect secondary marine productivity, with no difference among Buldir, Kiska, and Kasatochi islands across 585 km of the Aleutian Islands. Meals were very similar among species (mean Least–Crested Auklet overlap C = 0.68; Least–Whiskered Auklet overlap C = 0.96) and among sites, indicating limited partitioning of prey resources for auklets feeding chicks. The biomass of copepods and euphausiids in Least and Crested Auklet food loads was related negatively to the summer (June–July–August) North Pacific Gyre Oscillation, while in Whiskered Auklet food loads, this was negatively related to the winter (December–January–February) Pacific Decadal Oscillation, both of which track basin-wide sea-surface temperature (SST) anomalies. We found a significant quadratic relationship between the biomass of calanoid copepods in Least Auklet food loads at all three study sites and summer (June–July) SST, with maximal copepod biomass between 3–6°C (r 2 = 0.71). Outside this temperature range, zooplankton becomes less available to auklets through delayed development. Overall, our results suggest that auklets are able to buffer climate-mediated bottom-up forcing of demographic parameters like productivity, as the composition of chick meals has remained constant over the course of our study.

Feeding of Calanus finmarchicus and Oithona similis on the microplankton assemblage in the Irminger Sea, North Atlantic

The present investigation reviews published data on the feeding rates and prey selection of Oithona similis females, Calanus finmarchicus nauplii and females in the Irminger Sea in April/May and July/August 2002. Our aim was to examine how the feeding rates and prey selection of these three copepod stages respond to concomitant changes in microplankton community composition and prey abundance. Copepods typically ingested prey overall according to its ambient concentration although significant species and stage-specific differences in prey-type ingestion and selection were apparent. Despite being of comparable weight, the ingestion rates of C. finmarchicus nauplii were always higher than those of the O. similis females. Moreover, C. finmarchicus nauplii and O. similis females fed preferentially on diatoms and ciliates respectively, whereas adult female C. finmarchicus showed limited prey selectivity. Copepod grazing impact on total and on ciliates/dinoflagellates standing stock was <0.5 and <2%, respectively. We attribute this result to a combination of low grazing rates, low copepod abundance and low microplankton biomass, all of which are indicative of the non-bloom conditions under which these experiments were conducted. The differences in copepod feeding rates and prey selection we report reflect species and stage-specific eco-physiological adaptations, which may act as important driving forces for marine ecosystem structuring and functioning.

Mechanisms of long-term decline in size of lesser sandeels in the North Sea explored using a growth and phenology model

The lesser sandeel Ammodytes marinus is a key species in the North Sea ecosystem, transferring energy from planktonic producers to top predators. Previous studies have shown a long-term decline in the size of 0-group sandeels in the western North Sea, but they were unable to pinpoint the mechanism (later hatching, slower growth or changes in size-dependent mortality) or cause. To investigate the first 2 possibilities we combined 2 independent time series of sandeel size, namely data from chick-feeding Atlantic puffins Fratercula arctica and from the Continuous Plankton Recorder (CPR), in a novel statistical model implemented using Markov Chain Monte Carlo (MCMC). The model estimated annual mean length on 1 July, as well as hatching date and growth rate for sandeels from 1973 to 2006. Mean length-at-date declined by 22% over this period, corresponding to a 60% decrease in energy content, with a sharper decline since 2002. Up to the mid-1990s, the decline was associated with a trend towards later hatching. Subsequently, hatching became earlier again, and the continued trend towards smaller size appears to have been driven by lower growth rates, particularly in the most recent years, although we could not rule out changes in size-dependent mortality. Our findings point to major changes in key aspects of sandeel life history, which we consider are most likely due to direct and indirect temperature-related changes over a range of biotic factors, including the seasonal distribution of copepods and intra- and inter-specific competition with planktivorous fish. The results have implications both for the many predators of sandeels and for age and size of maturation in this aggregation of North Sea sandeels.

Regional and annual variation in black-legged kittiwake breeding productivity is related to sea surface temperature

Spatiotemporal variation in seabird demographic parameters is often pronounced and may be an important source of information on the state of marine ecosystems. Black-legged kittiwakes Rissa tridactyla in Britain and Ireland show strong regional structure in breeding productivity, and both temporal and spatial variation are probably related to abundance of the principal prey of breeding kittiwakes, the lesser sandeel Ammodytes marinus. Annual regional estimates of sandeel abundance do not exist, prohibiting direct tests of this hypothesis. We examined relationships between kittiwake breeding productivity and 2 potential proxies of sandeel abundance, winter sea surface temperature (SST) and abundance of Calanus copepods, within and among 6 regions in Britain and Ireland from 1986 to 2004. Means and trends in winter SST differed among regions, with higher means and less pronounced increasing trends in western (Atlantic) regions than in eastern (North Sea) regions. A negative relationship between breeding productivity and winter SST in the previous year was found within 2 regions (East Scotland and Orkney), as well as in a cross-regional analysis. Results were inconclusive for Calanus abundance, with a positive relationship in East Scotland and negative in Orkney. These results demonstrate that although a single environmental driver (SST) is related to both within- and between-region variation in a key demographic parameter, regional heterogeneity in SST trends as well as the importance of other factors may lead to highly variable responses. Understanding this heterogeneity is critical for predicting long-term effects of climate change or other anthropogenic drivers on marine ecosystems.

A new procedure to optimize the selection of groups in a classification tree

Agglomerative cluster analyses encompass many techniques, which have been widely used in various fields of science. In biology, and specifically ecology, datasets are generally highly variable and may contain outliers, which increase the difficulty to identify the number of clusters. Here we present a new criterion to determine statistically the optimal level of partition in a classification tree. The criterion robustness is tested against perturbated data (outliers) using an observation or variable with values randomly generated. The technique, called Random Simulation Test (RST), is tested on (1) the well-known Iris dataset [Fisher, R.A., 1936. The use of multiple measurements in taxonomic problems. Ann. Eugenic. 7, 179–188], (2) simulated data with predetermined numbers of clusters following Milligan and Cooper [Milligan, G.W., Cooper, M.C., 1985. An examination of procedures for determining the number of clusters in a data set. Psychometrika 50, 159–179] and finally (3) is applied on real copepod communities data previously analyzed in Beaugrand et al. [Beaugrand, G., Ibanez, F., Lindley, J.A., Reid, P.C., 2002. Diversity of calanoid copepods in the North Atlantic and adjacent seas: species associations and biogeography. Mar. Ecol. Prog. Ser. 232, 179–195]. The technique is compared to several standard techniques. RST performed generally better than existing algorithms on simulated data and proved to be especially efficient with highly variable datasets.

Macrophysiology of Calanus finmarchicus in the North Atlantic Ocean

Copepods represent the major part of the dry weight of the mesozooplankton in pelagic ecosystems and therefore have a central role in the secondary production of the North Atlantic Ocean. The calanoid copepod species Calanus finmarchicus is the main large copepod in subarctic waters of the North Atlantic, dominating the dry weight of the mesozooplankton in regions such as the northern North Sea and the Norwegian Sea. The objective of this work was to investigate the relationships between both the fundamental and realised niches of C. finmarchicus in order to better understand the future influence of global climate change on the abundance, the spatial distribution and the phenology of this key-structural species. Based on standardised Principal Component Analyses (PCAs), a macroecological approach was applied to determine factors affecting the spatial distribution of C. finmarchicus and to characterise its realised niche. Second, an ecophysiological model was used to calculate the Potential Egg Production Rate (PEPR) of C. finmarchicus and the centre of its fundamental niche. Relationships between the two niches were then investigated by correlation analysis. We found a close relationship between the fundamental and realised niches of C. finmarchicus at spatial, monthly and decadal scales. While the species is at the centre of its niche in the subarctic gyre, our joint macroecological and macrophysiological analyses show that it is at the edge of its niche in the North Sea, making the species in this region more vulnerable to temperature changes.

A comparison of two zooplankton time series data collected in the Gulf of Maine

Interannual and seasonal trends of zooplankton abundance and species composition were compared between the Bongo net and Continuous Plankton Recorder (CPR) time series in the Gulf of Maine. Data from 5799 Bongo and 3118 CPR samples were compared from the years 1978–2006. The two programs use different sampling methods, with the Bongo time series composed of bimonthly vertically integrated samples from locations throughout the region, while the CPR was towed monthly at 10 m depth on a transect that bisects the region. It was found that there was a significant correlation between the interannual (r = 0.67, P < 0.01) and seasonal (r = 0.95, P < 0.01) variability of total zooplankton counts. Abundance rankings of individual taxa were highly correlated and temporal trends of dominant copepods were similar between samplers. Multivariate analysis also showed that both time series equally detected major shifts in community structure through time. However, absolute abundance levels were higher in the Bongo and temporal patterns for many of the less abundant taxa groups were not similar between the two devices. The different mesh sizes of the samplers probably caused some of the discrepancies; but diel migration patterns, damage to soft bodied animals and avoidance of the small CPR aperture by some taxa likely contributed to the catch differences between the two devices. Nonetheless, Bongo data presented here confirm the previously published patterns found in the CPR data set, and both show that the abundance increase of the 1990s has been followed by average to below average levels from 2002 to 06.