Predicted levels of future ocean acidification and temperature rise could alter community structure and biodiversity in marine benthic communities

A mesocosm experiment was conducted to quantify the effects of reduced pH and elevated temperature on an intact marine invertebrate community. Standardised faunal communities, collected from the extreme low intertidal zone using artificial substrate units, were exposed to one of eight nominal treatments (four pH levels: 8.0, 7.7, 7.3 and 6.7, crossed with two temperature levels: 12 and 16°C). After 60 days exposure communities showed significant changes in structure and lower diversity in response to reduced pH. The response to temperature was more complex. At higher pH levels (8.0 and 7.7) elevated temperature treatments contained higher species abundances and diversity than the lower temperature treatments. In contrast, at lower pH levels (7.3 and 6.7), elevated temperature treatments had lower species abundances and diversity than lower temperature treatments. The species losses responsible for these changes in community structure and diversity were not randomly distributed across the different phyla examined. Molluscs showed the greatest reduction in abundance and diversity in response to low pH and elevated temperature, whilst annelid abundance and diversity was mostly unaffected by low pH and was higher at the elevated temperature. The arthropod response was between these two extremes with moderately reduced abundance and diversity at low pH and elevated temperature. Nematode abundance increased in response to low pH and elevated temperature, probably due to the reduction of ecological constraints, such as predation and competition, caused by a decrease in macrofaunal abundance. This community-based mesocosm study supports previous suggestions, based on observations of direct physiological impacts, that ocean acidification induced changes in marine biodiversity will be driven by differential vulnerability within and between different taxonomical groups. This study also illustrates the importance of considering indirect effects that occur within multispecies assemblages when attempting to predict the consequences of ocean acidification and global warming on marine communities.

Space-borne study of seasonal, multi-year, and decadal phytoplankton dynamics in the Bay of Biscay

Seasonal and inter-annual variations in phytoplankton community abundance in the Bay of Biscay are studied. Preliminarily processed by the National Aeronautics and Space Administration (NASA) to yield normalized water-leaving radiance and the top-of-the-atmosphere solar radiance, Sea-viewing Wide Field-of-View Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Coastal Zone Color Scanner (CZCS) data are further supplied to our dedicated retrieval algorithms to infer the sought for parameters. By applying the National Oceanic and Atmospheric Administration's (NOAA's) Advanced Very High Resolution Radiometer (AVHRR) data, the surface reflection coefficient in the only band in the visible spectrum is derived and employed for analysis. Decadal bridged time series of variations of diatom-dominated phytoplankton and green dinoflagellate Lepidodinium chlorophorum within the shelf zone and the coccolithophore Emiliania huxleyi in the pelagic area of the Bay are documented and analysed in terms of impacts of some biogeochemical and geophysical forcing factors.

Effect of elevated CO2 on the dynamics of particle-attached and free-living bacterioplankton communities in an Arctic fjord

In the frame of the European Project on Ocean Acidification (EPOCA), the response of an Arctic pelagic community (<3 mm) to a gradient of seawater pCO(2) was investigated. For this purpose 9 large-scale in situ mesocosms were deployed in Kongsfjorden, Svalbard (78 degrees 56.2' N, 11 degrees 53.6' E), in 2010. The present study investigates effects on the communities of particle-attached (PA; >3 mu m) and free-living (FL; <3 mu m > 0.2 mu m) bacteria by Automated Ribosomal Intergenic Spacer Analysis (ARISA) in 6 of the mesocosms, ranging from 185 to 1050 mu atm initial pCO(2), and the surrounding fjord. ARISA was able to resolve, on average, 27 bacterial band classes per sample and allowed for a detailed investigation of the explicit richness and diversity. Both, the PA and the FL bacterioplankton community exhibited a strong temporal development, which was driven mainly by temperature and phytoplankton development. In response to the breakdown of a picophytoplankton bloom, numbers of ARISA band classes in the PA community were reduced at low and medium CO2 (similar to 185-685 mu atm) by about 25 %, while they were more or less stable at high CO2 (similar to 820-1050 mu atm). We hypothesise that enhanced viral lysis and enhanced availability of organic substrates at high CO2 resulted in a more diverse PA bacterial community in the post-bloom phase. Despite lower cell numbers and extracellular enzyme activities in the post-bloom phase, bacterial protein production was enhanced in high CO2 mesocosms, suggesting a positive effect of community richness on this function and on carbon cycling by bacteria.

Transient dynamics of an altered large marine ecosystem

Overfishing of large-bodied benthic fishes and their subsequent population collapses on the Scotian Shelf of Canada’s east coast1, 2 and elsewhere3, 4 resulted in restructuring of entire food webs now dominated by planktivorous, forage fish species and macroinvertebrates. Despite the imposition of strict management measures in force since the early 1990s, the Scotian Shelf ecosystem has not reverted back to its former structure. Here we provide evidence of the transient nature of this ecosystem and its current return path towards benthic fish species domination. The prolonged duration of the altered food web, and its current recovery, was and is being governed by the oscillatory, runaway consumption dynamics of the forage fish complex. These erupting forage species, which reached biomass levels 900% greater than those prevalent during the pre-collapse years of large benthic predators, are now in decline, having outstripped their zooplankton food supply. This dampening, and the associated reduction in the intensity of predation, was accompanied by lagged increases in species abundances at both lower and higher trophic levels, first witnessed in zooplankton and then in large-bodied predators, all consistent with a return towards the earlier ecosystem structure. We conclude that the reversibility of perturbed ecosystems can occur and that this bodes well for other collapsed fisheries.

Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link

Exploring climate and anthropogenic impacts on marine ecosystems requires an understanding of how trophic components interact. However, integrative end-to-end ecosystem studies (experimental and/or modelling) are rare. Experimental investigations often concentrate on a particular group or individual species within a trophic level, while tropho-dynamic field studies typically employ either a bottom-up approach concentrating on the phytoplankton community or a top-down approach concentrating on the fish community. Likewise the emphasis within modelling studies is usually placed upon phytoplankton-dominated biogeochemistry or on aspects of fisheries regulation. In consequence the roles of zooplankton communities (protists and metazoans) linking phytoplankton and fish communities are typically under-represented if not (especially in fisheries models) ignored. Where represented in ecosystem models, zooplankton are usually incorporated in an extremely simplistic fashion, using empirical descriptions merging various interacting physiological functions governing zooplankton growth and development, and thence ignoring physiological feedback mechanisms. Here we demonstrate, within a modelled plankton food-web system, how trophic dynamics are sensitive to small changes in parameter values describing zooplankton vital rates and thus the importance of using appropriate zooplankton descriptors. Through a comprehensive review, we reveal the mismatch between empirical understanding and modelling activities identifying important issues that warrant further experimental and modelling investigation. These include: food selectivity, kinetics of prey consumption and interactions with assimilation and growth, form of voided material, mortality rates at different age-stages relative to prior nutrient history. In particular there is a need for dynamic data series in which predator and prey of known nutrient history are studied interacting under varied pH and temperature regimes.

Macroalgal blooms alter community structure and primary productivity in marine ecosystems

Eutrophication, coupled with loss of herbivory due to habitat degradation and overharvesting, has increased the frequency and severity of macroalgal blooms worldwide. Macroalgal blooms interfere with human activities in coastal areas, and sometimes necessitate costly algal removal programs. They also have many detrimental effects on marine and estuarine ecosystems, including induction of hypoxia, release of toxic hydrogen sulfide into the sediments and atmosphere, and the loss of ecologically and economically important species. However, macroalgal blooms can also increase habitat complexity, provide organisms with food and shelter, and reduce other problems associated with eutrophication. These contrasting effects make their overall ecological impacts unclear. We conducted a systematic review and meta-analysis to estimate the overall effects of macroalgal blooms on several key measures of ecosystem structure and functioning in marine ecosystems. We also evaluated some of the ecological and methodological factors that might explain the highly variable effects observed in different studies. Averaged across all studies, macroalgal blooms had negative effects on the abundance and species richness of marine organisms, but blooms by different algal taxa had different consequences, ranging from strong negative to strong positive effects. Blooms' effects on species richness also depended on the habitat where they occurred, with the strongest negative effects seen in sandy or muddy subtidal habitats and in the rocky intertidal. Invertebrate communities also appeared to be particularly sensitive to blooms, suffering reductions in their abundance, species richness, and diversity. The total net primary productivity, gross primary productivity, and respiration of benthic ecosystems were higher during macroalgal blooms, but blooms had negative effects on the productivity and respiration of other organisms. These results suggest that, in addition to their direct social and economic costs, macroalgal blooms have ecological effects that may alter their capacity to deliver important ecosystem services.

High-resolution view of the spring bloom initiation and net community production in the Subantarctic Southern Ocean using glider data

In the Southern Ocean, there is increasing evidence that seasonal to subseasonal temporal scales, and meso- to submesoscales play an important role in understanding the sensitivity of ocean primary productivity to climate change. This drives the need for a high-resolution approach to re- solving biogeochemical processes. In this study, 5.5 months of continuous, high-resolution (3 h, 2 km horizontal resolution) glider data from spring to summer in the Atlantic Subantarctic Zone is used to investigate: (i) the mechanisms that drive bloom initiation and high growth rates in the region and (ii) the seasonal evolution of water column production and respiration. Bloom initiation dates were analysed in the context of upper ocean boundary layer physics highlighting sensitivities of different bloom detection methods to different environmental processes. Model results show that in early spring (September to mid-November) increased rates of net community production (NCP) are strongly affected by meso- to submesoscale features. In late spring/early summer (late-November to mid-December) seasonal shoaling of the mixed layer drives a more spatially homogenous bloom with maximum rates of NCP and chlorophyll biomass. A comparison of biomass accumulation rates with a study in the North Atlantic highlights the sensitivity of phytoplankton growth to fine-scale dynamics and emphasizes the need to sample the ocean at high resolution to accurately resolve phytoplankton phenology and improve our ability to estimate the sensitivity of the biological carbon pump to climate change.

ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

The ERSEM model is one of the most established ecosystem models for the lower trophic levels of the marine food-web in the scientific literature. Since its original development in the early nineties it has evolved significantly from a coastal ecosystem model for the North-Sea to a generic tool for ecosystem simulations from shelf seas to the global ocean. The current model release contains all essential elements for the pelagic and benthic part of the marine ecosystem, including the microbial food-web, the carbonate system and calcification. Its distribution is accompanied by a testing framework enabling the analysis of individual parts of the model. Here we provide a detailed mathematical description of all ERSEM components along with case-studies of mesocosm type simulations, water column implementations and a brief example of a full-scale application for the North-West European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.

Analyses of sublittoral macrobenthic community change in a marine nature reserve using similarity profiles (SIMPROF).

Sublittoral macrobenthic communities in the Skomer Marine Nature Reserve (SMNR), Pembrokeshire, Wales, were sampled at 10 stations in 1993, 1996, 1998, 2003, 2007 and 2009 using a Day grab and a 0.5 mm mesh. The time series is analysed using Similarities Profiles (SIMPROF) tests and associated methods. Q-mode analysis using clustering with Type 1 SIMPROF addresses multivariate structure among samples, showing that there is clear structure associated with differences among years. Inverse (r-mode) analysis using Type 2 SIMPROF decisively rejects a hypothesis that species are not associated with each other. Clustering of the variables (species) with Type 3 SIMPROF identifies groups of species which covary coherently through the time-series. The time-series is characterised by a dramatic decline in abundances and diversity between the 1993 and 1996 surveys. By 1998 there had been a shift in community composition from the 1993 situation, with different species dominating. Communities had recovered in terms of abundance and species richness, but different species dominated the community. No single factor could be identified which unequivocally explained the dramatic changes observed in the SMNR. Possible causes were the effects of dispersed oil and dispersants from the Sea Empress oil spill in February 1996 and the cessation of dredge-spoil disposal off St Anne’s Head in 1995, but the most likely cause was severe weather. With many species, and a demonstrable recovery from an impact, communities within the SMNR appear to be diverse and resilient. If attributable to natural storms, the changes observed here indicate that natural variability may be much more important than is generally taken into account in the design of monitoring programmes.

ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

The European Regional Seas Ecosystem Model (ERSEM) is one of the most established ecosystem models for the lower trophic levels of the marine food web in the scientific literature. Since its original development in the early nineties it has evolved significantly from a coastal ecosystem model for the North Sea to a generic tool for ecosystem simulations from shelf seas to the global ocean. The current model release contains all essential elements for the pelagic and benthic parts of the marine ecosystem, including the microbial food web, the carbonate system, and calcification. Its distribution is accompanied by a testing framework enabling the analysis of individual parts of the model. Here we provide a detailed mathematical description of all ERSEM components along with case studies of mesocosm-type simulations, water column implementations, and a brief example of a full-scale application for the north-western European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.