Proglacial river dataset from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008 - 2010, version 1.0

Pressing scientific questions concerning the Greenland ice sheet's climatic sensitivity, hydrology, and contributions to current and future sea level rise require hydrological datasets to resolve. While direct observations of ice sheet meltwater losses can be obtained in terrestrial rivers draining the ice sheet and from lake levels, few such datasets exist. We present a new dataset of meltwater river discharge for the vicinity of Kangerlussuaq, Southwest Greenland. The dataset contains measurements of river stage and discharge for three sites along the Akuliarusiarsuup Kuua (Watson) River's northern tributary, with 30 minute temporal resolution between June 2008 and August 2010. Additional data of water temperature, air pressure, and lake water depth and temperature are also provided. Discharge data were measured at sites with near-ideal properties for such data collection. Regardless, high water bedload and turbulent flow introduce considerable uncertainty. These were constrained and quantified using statistical techniques, thereby providing a high quality dataset from this important site. The greatest data uncertainties are associated with streambed elevation change and measurements. Large portions of stream channels deepened according to statistical tests, but poor precision of streambed depth measurements also added uncertainty. Quality checked data are freely available for scientific use as supplementary online material.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Effects of various pollutant mixtures on immune responses of the blue mussel (Mytilus edulis) collected at a salinity gradient in Danish coastal waters, links to supplementary material

The Baltic Sea is a semi-enclosed sea with a steady salinity gradient (3 per mil-30 per mil). Organisms have adapted to such low salinities, but are suspected to be more susceptible to stress. Within the frame of the integrated environmental monitoring BONUS + project "BEAST" the applicability of immune responses of the blue mussel was investigated in Danish coastal waters. The sampling sites were characterised by a salinity range (11-19 per mil) and different mixtures of contaminants (metals, PAHs and POPs), according to chemical analysis of mussel tissues. Variation partitioning (redundancy analysis) was applied to decompose salinity and contamination effects. The results indicated that cellular immune responses (total and differential haemocyte count, phagocytic activity and apoptosis) were mainly influenced by contaminants, whereas humoral factors (haemolytic activity) were mainly impacted by salinity. Hence, cellular immune functions may be suitable as biomarkers in monitoring programmes for the Baltic Sea and other geographic regions with salinity variances of the studied range.

Experimental evaluation of seaweeds as a vector for microplastics into marine food webs, supplementary material

The ingestion of microplastics has been shown for a great variety of marine organisms. However, benthic marine mesoherbivores such as the common periwinkle Littorina littorea have been largely disregarded in studies about the effects of microplastics on the marine biota, probably because the pathway for microplastics to this functional group of organisms was not obvious. In laboratory experiments we showed that the seaweed Fucus vesiculosus retains suspended microplastics on its surface. The numbers of microplastics that adhered to the algae correlated with the concentrations of suspended particles in the water. In choice feeding assays L. littorea did not distinguish between algae with adherent microplastics and clean algae without microplastics, indicating that the snails do not recognize solid nonfood particles in the submillimeter size range as deleterious. In periwinkles that were feeding on contaminated algae, microplastics were found in the stomach and in the gut. However, no microplastics were found in the midgut gland, which is the principle digestive organ of gastropods. Microplastics in the fecal pellets of the periwinkles indicate that the particles do not accumulate rapidly inside the animals but are mostly released with the feces. Our results provide the first evidence that seaweeds may represent an efficient pathway for microplastics from the water to marine benthic herbivores.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Intraspecific facilitation by allelochemical mediated grazing protection within a toxigenic dinoflagellate population, link to supplementary material

Dinoflagellates are a major cause of harmful algal blooms, with consequences for coastal marine ecosystem functioning and services. Alexandrium tamarense is one of the most abundant and widespread toxigenic species in the temperate northern and southern hemisphere, and produces paralytic shellfish poisoning toxins as well as lytic allelochemical substances. These bioactive compounds may support the success of A. tamarense and its ability to form blooms. Here we investigate the impact of grazing on monoclonal and mixed set-ups of highly (Alex2) and moderately (Alex4) allelochemically active A. tamarense strains and on a non-allelochemically active conspecific (Alex5) by the heterotrophic dinoflagellate Polykrikos kofoidii. While Alex4 and particularly Alex5 were strongly grazed by P. kofoidii when offered alone, both strains grew well in the mixed assemblages (Alex4+Alex5 and Alex2+Alex5). Hence, the allelochemical active strains facilitated growth of the non-active strain by protecting the population as a whole against grazing. Based on our results, we argue that facilitation among clonal lineages within a species may partly explain the high genotypic and phenotypic diversity of Alexandrium populations. Populations of Alexandrium may comprise multiple cooperative traits that act in concert with intraspecific facilitation, and hence promote the success of this notorious harmful algal bloom species.

Spatial compartmentalization of free radical formation and mitochondrial heterogeneity in bivalve gills revealed by live-imaging techniques, supplementary material

Live-imaging techniques (LIT) utilize target-specific fluorescent dyes to visualize biochemical processes using confocal and multiphoton scanning microscopy, which are increasingly employed as non-invasive approach to physiological in-vivo and ex-vivo studies. Here we report application of LIT to bivalve gills for ex-vivo analysis of gill physiology and mapping of reactive oxygen (ROS) and nitrogen (RNS) species formation in the living tissue. Our results indicate that H2O2, HOO. and ONOO- radicals (assessed through C-H2DFFDA staining) are mainly formed within the blood sinus of the filaments and are likely to be produced by hemocytes as defense against invading pathogens. The oxidative damage in these areas is controlled by enhanced CAT (catalase) activities recorded within the filaments. The outermost areas of the ciliated epithelial cells composing the filaments, concentrated the highest mitochondrial densities (MTK Deep Red 633 staining) and the most acidic pH values (as observed with ageladine-a). These mitochondria have low (depolarized) membrane potentials (D psi m) (JC-1 staining), suggesting that the high amounts of ATP required for ciliary beating may be in part produced by non-mitochondrial mechanisms, such as the enzymatic activity of an ATP-regenerating kinase. Nitric oxide (NO, DAF-2DA staining) produced in the region of the peripheral mitochondria may have an effect on mitochondrial electron transport and possibly cause the low membrane potential. High DAF-2DA staining was moreover observed in the muscle cells composing the wall of the blood vessels where NO may be involved in regulating blood vessel diameter. On the ventral bend of the gills, subepithelial mucus glands (SMG) contain large mucous vacuoles showing higher fluorescence intensities for O2.- (DHE staining) than the rest of the tissue. Given the antimicrobial properties of superoxide, release of O2.- into the mucus may help to avoid the development of microbial biofilms on the gill surface. However, cells of the ventral bends are paying a price for this antimicrobial protection, since they show significantly higher oxidative damage, according to the antioxidant enzyme activities and the carbonyl levels, than the rest of the gill tissue. This study provides the first evidence that one single epithelial cell may contain mitochondria with significantly different membrane potentials. Furthermore, we provide new insight into ROS and RNS formation in ex-vivo gill tissues which opens new perspectives for unraveling the different ecophysiological roles of ROS and RNS in multifunctional organs such as gills.

Impact of thermal stress on evolutionary trajectories of pathogen resistance in three-spined stickleback (Gasterosteus aculeatus), supplementary material

Background: Pathogens are a major regulatory force for host populations, especially under stressful conditions. Elevated temperatures may enhance the development of pathogens, increase the number of transmission stages, and can negatively influence host susceptibility depending on host thermal tolerance. As a net result, this can lead to a higher prevalence of epidemics during summer months. These conditions also apply to marine ecosystems, where possible ecological impacts and the population-specific potential for evolutionary responses to changing environments and increasing disease prevalence are, however, less known. Therefore, we investigated the influence of thermal stress on the evolutionary trajectories of disease resistance in three marine populations of three-spined sticklebacks Gasterosteus aculeatus by combining the effects of elevated temperature and infection with a bacterial strain of Vibrio sp. using a common garden experiment. Results: We found that thermal stress had an impact on fish weight and especially on survival after infection after only short periods of thermal acclimation. Environmental stress reduced genetic differentiation (QST) between populations by releasing cryptic within-population variation. While life history traits displayed positive genetic correlations across environments with relatively weak genotype by environment interactions (GxE), environmental stress led to negative genetic correlations across environments in pathogen resistance. This reversal of genetic effects governing resistance is probably attributable to changing environment-dependent virulence mechanisms of the pathogen interacting differently with host genotypes, i.e. GPathogenxGHostxE or (GPathogenxE)x(GHostxE) interactions, rather than to pure host genetic effects, i.e. GHostxE interactions. Conclusion: To cope with climatic changes and the associated increase in pathogen virulence, host species require wide thermal tolerances and pathogen-resistant genotypes. The higher resistance we found for some families at elevated temperatures showed that there is evolutionary potential for resistance to Vibrio sp. in both thermal environments. The negative genetic correlation of pathogen resistance between thermal environments, on the other hand, indicates that adaptation to current conditions can be a weak predictor for performance in changing environments. The observed feedback on selective gradients exerted on life history traits may exacerbate this effect, as it can also modify the response to selection for other vital components of fitness.

Evidence of cryptic and pseudocryptic speciation in the Paracalanus parvus species complex (Crustacea, Copepoda, Calanoida), supplementary material

Introduction Many marine planktonic crustaceans such as copepods have been considered as widespread organisms. However, the growing evidence for cryptic and pseudo-cryptic speciation has emphasized the need of re-evaluating the status of copepod species complexes in molecular and morphological studies to get a clearer picture about pelagic marine species as evolutionary units and their distributions. This study analyses the molecular diversity of the ecologically important Paracalanus parvus species complex. Its seven currently recognized species are abundant and also often dominant in marine coastal regions worldwide from temperate to tropical oceans. Results COI and Cytochrome b sequences of 160 specimens of the Paracalanus parvus complex from all oceans were obtained. Furthermore, 42 COI sequences from GenBank were added for the genetic analyses. Thirteen distinct molecular operational taxonomic units (MOTU) and two single sequences were revealed with cladistic analyses (Maximum Likelihood, Bayesian Inference), of which seven were identical with results from species delimitation methods (barcode gaps, ABDG, GMYC, Rosenberg's P(AB)). In total, 10 to 12 putative species were detected and could be placed in three categories: (1) temperate geographically isolated, (2) warm-temperate to tropical wider spread and (3) circumglobal warm-water species. Conclusions The present study provides evidence of cryptic or pseudocryptic speciation in the Paracalanus parvus complex. One major insight is that the species Paracalanus parvus s.s. is not panmictic, but may be restricted in its distribution to the northeastern Atlantic.

Working on a baseline for the Kongsfjorden food web: production and properties of macroalgal particulate organic matter (POM), supplementary material

Macroalgae, in particular kelps, produce a large amount of biomass in Kongsfjorden, which is to a great extent released into the water in an annual cycle. As an example, the brown alga Alaria esculenta loses its blade gradually, 3 ± 0.8 % of the blade area per day (August 2012), thereby adding to the pool of particulate organic matter (POM) in the fjord. Upon release small thallus pieces are "aging" in that they are prone to leaching and serving as substrate for microorganisms, thus turning into palatable food for suspension and bottom feeders. In order to define a macroalgal baseline for the Kongsfjorden food web, stable isotopes d14C and d15N were measured in individuals of A. esculenta, Saccharina latissima and Laminaria digitata directly sampled after collection and in artificially produced POM (aPOM) of A. esculenta that was allowed to age under experimental conditions. In aPOM from this species sampled in August 2012 the C/N ratios decreased between d1 and d8 of a 14-day culture period in parallel to the fading photosynthetic activity of the algal fragments as demonstrated by use of an Imaging-PAM. Microscopic observations of the aPOM in August 2012 and 2013 revealed the frequent occurrence of small brown algal endo- and epiphytes. First feeding experiments with Mysis oculata (Mysids) and Hiatella arctica (Bivalves) showed that these species can ingest macroalgal POM. The importance of kelp-derived POM for the food web is subject of the current research.

A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Supplementary material: XRF core scanning data, CaCO3 values, foraminiferal carbon and oxygen isotope data, Mg/Ca data, clay mineralogy and helium isotope results

The Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) is associated with abrupt climate change, carbon cycle perturbation, ocean acidification, as well as biogeographic shifts in marine and terrestrial biota that were largely reversed as the climatic transient waned. We report a clear exception to the behavior of the PETM as a reversing climatic transient in the eastern North Atlantic (Deep-Sea Drilling Project Site 401, Bay of Biscay) where the PETM initiates a greatly prolonged environmental change compared to other places on Earth where records exist. The observed environmental perturbation extended well past the d13C recovery phase and up to 650 kyr after the PETM onset according to our extraterrestrial 3He-based age-model. We observe a strong decoupling of planktic foraminiferal d18O and Mg/Ca values during the PETM d13C recovery phase, which in combination with results from helium isotopes and clay mineralogy, suggests that the PETM triggered a hydrologic change in western Europe that increased freshwater flux and the delivery of weathering products to the eastern North Atlantic. This state change persisted long after the carbon-cycle perturbation had stopped. We hypothesize that either long-lived continental drainage patterns were altered by enhanced hydrological cycling induced by the PETM, or alternatively that the climate system in the hinterland area of Site 401 was forced into a new climate state that was not easily reversed in the aftermath of the PETM.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.

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.