Climate sensitivity across marine domains of life: limits to evolutionary adaptation shape species interactions, supplementary material

Organisms in all domains, Archaea, Bacteria, and Eukarya will respond to climate change with differential vulnerabilities resulting in shifts in species distribution, coexistence, and interactions. The identification of unifying principles of organism functioning across all domains would facilitate a cause and effect understanding of such changes and their implications for ecosystem shifts. For example, the functional specialization of all organisms in limited temperature ranges leads us to ask for unifying functional reasons. Organisms also specialize in either anoxic or various oxygen ranges, with animals and plants depending on high oxygen levels. Here, we identify thermal ranges, heat limits of growth, and critically low (hypoxic) oxygen concentrations as proxies of tolerance in a meta-analysis of data available for marine organisms, with special reference to domain-specific limits. For an explanation of the patterns and differences observed, we define and quantify a proxy for organismic complexity across species from all domains. Rising complexity causes heat (and hypoxia) tolerances to decrease from Archaea to Bacteria to uni- and then multicellular Eukarya. Within and across domains, taxon-specific tolerance limits likely reflect ultimate evolutionary limits of its species to acclimatization and adaptation. We hypothesize that rising taxon-specific complexities in structure and function constrain organisms to narrower environmental ranges. Low complexity as in Archaea and some Bacteria provide life options in extreme environments. In the warmest oceans, temperature maxima reach and will surpass the permanent limits to the existence of multicellular animals, plants and unicellular phytoplankter. Smaller, less complex unicellular Eukarya, Bacteria, and Archaea will thus benefit and predominate even more in a future, warmer, and hypoxic ocean.

(Table A1) Species list of the high Antarctic Weddell Sea food web

Human-induced habitat destruction, overexploitation, introduction of alien species and climate change are causing species to go extinct at unprecedented rates, from local to global scales. There are growing concerns that these kinds of disturbances alter important functions of ecosystems. Our current understanding is that key parameters of a community (e.g. its functional diversity, species composition, and presence/absence of vulnerable species) reflect an ecological network's ability to resist or rebound from change in response to pressures and disturbances, such as species loss. If the food web structure is relatively simple, we can analyse the roles of different species interactions in determining how environmental impacts translate into species loss. However, when ecosystems harbour species-rich communities, as is the case in most natural systems, then the complex network of ecological interactions makes it a far more challenging task to perceive how species' functional roles influence the consequences of species loss. One approach to deal with such complexity is to focus on the functional traits of species in order to identify their respective roles: for instance, large species seem to be more susceptible to extinction than smaller species. Here, we introduce and analyse the marine food web from the high Antarctic Weddell Sea Shelf to illustrate the role of species traits in relation to network robustness of this complex food web. Our approach was threefold: firstly, we applied a new classification system to all species, grouping them by traits other than body size; secondly, we tested the relationship between body size and food web parameters within and across these groups and finally, we calculated food web robustness. We addressed questions regarding (i) patterns of species functional/trophic roles, (ii) relationships between species functional roles and body size and (iii) the role of species body size in terms of network robustness. Our results show that when analyzing relationships between trophic structure, body size and network structure, the diversity of predatory species types needs to be considered in future studies.

Seawater carbonate chemistry and growth, carbon acquisition, and species interaction of Antarctic phytoplankton species in a laboratory experiment

Despite the fact that ocean acidification is considered to be especially pronounced in the Southern Ocean, little is known about CO2-dependent physiological processes and the interactions of Antarctic phytoplankton key species. We therefore studied the effects of CO2 partial pressure (PCO2) (16.2, 39.5, and 101.3 Pa) on growth and photosynthetic carbon acquisition in the bloom-forming species Chaetoceros debilis, Pseudo-nitzschia subcurvata, Fragilariopsis kerguelensis, and Phaeocystis antarctica. Using membrane-inlet mass spectrometry, photosynthetic O2 evolution and inorganic carbon (Ci) fluxes were determined as a function of CO2 concentration. Only the growth of C. debilis was enhanced under high PCO2. Analysis of the carbon concentrating mechanism (CCM) revealed the operation of very efficient CCMs (i.e., high Ci affinities) in all species, but there were species-specific differences in CO2-dependent regulation of individual CCM components (i.e., CO2 and uptake kinetics, carbonic anhydrase activities). Gross CO2 uptake rates appear to increase with the cell surface area to volume ratios. Species competition experiments with C. debilis and P. subcurvata under different PCO2 levels confirmed the CO2-stimulated growth of C. debilis observed in monospecific incubations, also in the presence of P. subcurvata. Independent of PCO2, high initial cell abundances of P. subcurvata led to reduced growth rates of C. debilis. For a better understanding of future changes in phytoplankton communities, CO2-sensitive physiological processes need to be identified, but also species interactions must be taken into account because their interplay determines the success of a species.

Predicting changes in abundance and juvenile proportions of earthworms in the field using glmmADMB in R

The study aimed at investigating effects of three differently acting biocides; the insecticide esfenvalerate, the fungicide picoxystrobin and the bactericide triclosan, applied individually and as a mixture, on an earthworm community in the field. A concentration-response design was chosen and results were analyzed using univariate and multivariate approaches. Effects on juvenile proportions were less pronounced and more variable than effects on abundance, but effects in general were species- and chemical-specific, and temporal variations distinct. Esfenvalerate and picoxystrobin appeared to elicit stronger effects than triclosan at laboratory-based ECx values, which is in accordance with our previous laboratory study on Eisenia fetida. The mixture affected abundance and juvenile proportions, but the latter only at high mixture concentrations. Esfenvalerate and picoxystrobin appeared to be the main drivers for the mixture's toxicity. Species-specific toxicity patterns question the reliability of mixture toxicity predictions derived on E. fetida for field earthworms. Biocide concentrations equaling EC50s (reproduction) for E. fetida provoked effects on the field earthworms mainly exceeding 50%, indicating effect intensification from the laboratory to field as well as the influence of indirect effects produced by species interactions. The differing results of the present field study and the previous laboratory study imply that lower- and higher-tier studies may not be mutually exclusive, but to be used in complementary.

Ocean acidification affects growth but not nutritional quality of the seaweed Fucus vesiculosus (Phaeophyceae, Fucales)

Understanding the ecological implications of global climate change requires investigations of not only the direct effects of environmental change on species performance but also indirect effects that arise from altered species interactions. We performed CO2 perturbation experiments to investigate the effects of ocean acidification on the trophic interaction between the brown seaweed Fucus vesiculosus and the herbivorous isopod Idotea baltica. We predicted faster growth of F. vesiculosus at elevated CO2-concentrations and higher carbon content of the algal tissue. We expected that I. baltica has different consumption rates on algae that have been grown at different CO2 levels and that the isopods remove surplus carbon metabolically by enhanced respiration. Surprisingly, growth of F. vesiculosus as well as the C:N-ratio of the algal tissue were reduced at high CO2-levels. The changes in the elemental composition had no effect on the consumption rates and the respiration of the herbivores. An additional experiment showed that consumption of F. vesiculosus by the isopod Idotea emarginata was independent of ocean acidification and temperature. Our results could not reveal any effects of ocean acidification on the per capita strength of the trophic interaction between F. vesiculosus and its consumers. However, reduced growth of the algae at high CO2-concentrations might reduce the capability of the seaweed to compensate losses due to intense herbivory.

Experiment: Food availability outweighs ocean acidification effects in juvenile Mytilus edulis

Ocean acidification is expected to decrease calcification rates of bivalves. Nevertheless in many coastal areas high pCO2 variability is encountered already today. Kiel Fjord (Western Baltic Sea) is a brackish (12-20 g kg-1) and CO2 enriched habitat, but the blue mussel Mytilus edulis dominates the benthic community. In a coupled field and laboratory study we examined the annual pCO2 variability in this habitat and the combined effects of elevated pCO2 and food availability on juvenile M. edulis growth and calcification. In the laboratory experiment, mussel growth and calcification were found to chiefly depend on food supply, with only minor impacts of pCO2 up to 3350 µatm. Kiel Fjord was characterized by strong seasonal pCO2 variability. During summer, maximal pCO2 values of 2500 µatm were observed at the surface and >3000 µatm at the bottom. However, the field growth experiment revealed seven times higher growth and calcification rates of M. edulis at a high pCO2 inner fjord field station (mean pCO2 ca. 1000 µatm) in comparison to a low pCO2 outer fjord station (ca. 600 µatm). In addition, mussels were able to outcompete the barnacle Amphibalanus improvisus at the high pCO2 site. High mussel productivity at the inner fjord site was enabled by higher particulate organic carbon concentrations. Kiel Fjord is highly impacted by eutrophication, which causes bottom water hypoxia and consequently high seawater pCO2. At the same time, elevated nutrient concentrations increase the energy availability for filter feeding organisms such as mussels. Thus M. edulis can dominate over a seemingly more acidification resistant species such as A. improvisus. We conclude that benthic stages of M. edulis tolerate high ambient pCO2 when food supply is abundant and that important habitat characteristics such as species interactions and energy availability need to be considered to predict species vulnerability to ocean acidification.

Sea urchin response to rising pCO2 shows ocean acidification may fundamentally alter the chemistry of marine skeletons

Ocean acidification caused by an increase in pCO2 is expected to drastically affect marine ecosystem composition, yet there is much uncertainty about the mechanisms through which ecosystems may be affected. Here we studied sea urchins that are common and important grazers in the Mediterranean (Paracentrotus lividus and Arbacia lixula). Our study included a natural CO2 seep plus reference sites in the Aegean Sea off Greece. The distribution of A. lixula was unaffected by the low pH environment, whereas densities of P. lividus were much reduced. There was skeletal degradation in both species living in acidified waters compared to reference sites and remarkable increases in skeletal manganese levels (P. lividus had a 541% increase, A. lixula a 243% increase), presumably due to changes in mineral crystalline structure. Levels of strontium and zinc were also altered. It is not yet known whether such dramatic changes in skeletal chemistry will affect coastal systems but our study reveals a mechanism that may alter inter-species interactions.

Functional consequences of prey acclimation to ocean acidification for the prey and its predator

Ocean acidification is the suite of chemical changes to the carbonate system of seawater as a consequence of anthropogenic carbon dioxide (CO2) emissions. Despite a growing body of evidences demonstrating the negative effects of ocean acidification on marine species, the consequences at the ecosystem level are still unclear. One factor limiting our ability to upscale from species to ecosystem is the poor mechanistic understanding of the functional consequences of the observed effects on organisms. This is particularly true in the context of species interactions. The aim of this work was to investigate the functional consequence of the exposure of a prey (the mussel Brachidontes pharaonis) to ocean acidification for both the prey and its predator (the crab Eriphia verrucosa). Mussels exposed to pH 7.5 for >4 weeks showed significant decreases in condition index and in mechanical properties (65% decrease in maximum breaking load) as compared with mussels acclimated to pH 8.0. This translated into negative consequences for the mussel in presence of the predator crab. The crab feeding efficiency increased through a significant 27% decrease in prey handling time when offered mussels acclimated to the lowest pH. The predator was also negatively impacted by the acclimation of the prey, probably as a consequence of a decreased food quality. When fed with prey acclimated under decreased pH for 3 months, crab assimilation efficiency significantly decreased by 30% and its growth rate was 5 times slower as compared with crab fed with mussels acclimated under high pH. Our results highlight the important to consider physiological endpoints in the context of species interactions.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2004)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2004 just prior to mowing (during peak standing biomass in late May and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Mean cover and frequencies of vascular plant species during three decades around Abisko Scientific Research Station

Plant species distributions are expected to shift and diversity is expected to decline as a result of global climate change, particularly in the Arctic where climate warming is amplified. We have recorded the changes in richness and abundance of vascular plants at Abisko, sub-Arctic Sweden, by re-sampling five studies consisting of seven datasets; one in the mountain birch forest and six at open sites. The oldest study was initiated in 1977-1979 and the latest in 1992. Total species number increased at all sites except for the birch forest site where richness decreased. We found no general pattern in how composition of vascular plants has changed over time. Three species, Calamagrostis lapponica, Carex vaginata and Salix reticulata, showed an overall increase in cover/frequency, while two Equisetum taxa decreased. Instead, we showed that the magnitude and direction of changes in species richness and composition differ among sites.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2007)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2007 just prior to mowing (during peak standing biomass in early June and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four (May) or three (August) rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2006)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2006 just prior to mowing (during peak standing biomass in early June and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2003)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2003 just prior to mowing (during peak standing biomass in late May and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

The California current predator diet database: synthesis of common forage species

Characterization of the diets of upper-trophic predators is a key ingredient in management including the development of ecosystem-based fishery management plans, conservation efforts for top predators, and ecological and economic modeling of predator prey interactions. The California Current Predator Diet Database (CCPDD) synthesizes data from published records of predator food habits over the past century. The database includes diet information for 100+ upper-trophic level predator species, based on over 200 published citations from the California Current region of the Pacific Ocean, ranging from Baja, Mexico to Vancouver Island, Canada. We include diet data for all predators that consume forage species: seabirds, cetaceans, pinnipeds, bony and cartilaginous fishes, and a predatory invertebrate; data represent seven discrete geographic regions within the CCS (Canada, WA, OR, CA-n, CA-c, CA-s, Mexico). The database is organized around predator-prey links that represent an occurrence of a predator eating a prey or group of prey items. Here we present synthesized data for the occurrence of 32 forage species (see Table 2 in the affiliated paper) in the diet of pelagic predators (currently submitted to Ecological Informatics). Future versions of the shared-data will include diet information for all prey items consumed, not just the forage species of interest.