Experimental data on photophysiology and growth rates of Pararotalia calcariformata

The eastern Mediterranean is a hotspot of biological invasions. Numerous species of Indo-pacific origin have colonized the Mediterranean in recent times, including tropical symbiont-bearing foraminifera. Among these is the species Pararotalia calcariformata. Unlike other invasive foraminifera, this species has been discovered only two decades ago and is restricted to the eastern Mediterranean coast. Combining ecological, genetic and physiological observations, we attempt to explain the recent invasion of this species in the Mediterranean Sea. Using morphological and genetic data, we confirm the species attribution to P. calcariformata McCulloch 1977 and identify its symbionts as a consortium of diatom species dominated by Minutocellus polymorphus. We document photosynthetic activity of its endosymbionts using Pulse Amplitude Modulated Fluorometry and test the effects of elevated temperatures on growth rates of asexual offspring. The culturing of asexual offspring for 120 days shows a 30-day period of rapid growth followed by a period of slower growth. A subsequent 48-day temperature sensitivity experiment indicates a similar developmental pathway and high growth rate at 28°C, whereas an almost complete inhibition of growth was observed at 20°C and 35°C. This indicates that the offspring of this species may have lower tolerance to cold temperatures than what would be expected for species native to the Mediterranean. We expand this hypothesis by applying a Species Distribution Model (SDM) based on modern occurrences in the Mediterranean using three environmental variables: irradiance, turbidity and yearly minimum temperature. The model reproduces the observed restricted distribution and indicates that the range of the species will drastically expand westwards under future global change scenarios. We conclude that P. calcariformata established a population in the Levant because of the recent warming in the region. In line with observations from other groups of organisms, our results indicate that continued warming of the eastern Mediterranean will facilitate the invasion of more tropical marine taxa into the Mediterranean, disturbing local biodiversity and ecosystem structure.

Simulated Methyl iodide emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by OMIP data

Production pathways of the prominent volatile organic halogen compound methyl iodide (CH3I) are not fully understood. Based on observations, production of CH3I via photochemical degradation of organic material or via phytoplankton production has been proposed. Additional insights could not be gained from correlations between observed biological and environmental variables or from biogeochemical modeling to identify unambiguously the source of methyl iodide. In this study, we aim to address this question of source mechanisms with a three-dimensional global ocean general circulation model including biogeochemistry (MPIOM-HAMOCC (MPIOM - Max Planck Institute Ocean Model HAMOCC - HAMburg Ocean Carbon Cycle model)) by carrying out a series of sensitivity experiments. The simulated fields are compared with a newly available global data set. Simulated distribution patterns and emissions of CH3I differ largely for the two different production pathways. The evaluation of our model results with observations shows that, on the global scale, observed surface concentrations of CH3I can be best explained by the photochemical production pathway. Our results further emphasize that correlations between CH3I and abiotic or biotic factors do not necessarily provide meaningful insights concerning the source of origin. Overall, we find a net global annual CH3I air-sea flux that ranges between 70 and 260 Gg/yr. On the global scale, the ocean acts as a net source of methyl iodide for the atmosphere, though in some regions in boreal winter, fluxes are of the opposite direction (from the atmosphere to the ocean).

Geographical variation in shell morphology of juvenile snails (Concholepas concholepas) along the physical-chemical gradient of the Chilean coast

Changes in phenotypic traits, such as mollusc shells, are indicative of variations in selective pressure along environmental gradients. Recently, increased sea surface temperature (SST) and ocean acidification (OA) due to increased levels of carbon dioxide in the seawater have been described as selective agents that may affect the biological processes underlying shell formation in calcifying marine organisms. The benthic snail Concholepas concholepas (Muricidae) is widely distributed along the Chilean coast, and so is naturally exposed to a strong physical-chemical latitudinal gradient. In this study, based on elliptical Fourier analysis, we assess changes in shell morphology (outlines analysis) in juvenile C. concholepas collected at northern (23°S), central (33°S) and southern (39°S) locations off the Chilean coast. Shell morphology of individuals collected in northern and central regions correspond to extreme morphotypes, which is in agreement with both the observed regional differences in the shell apex outlines, the high reclassification success of individuals (discriminant function analysis) collected in these regions, and the scaling relationship in shell weight variability among regions. However, these extreme morphotypes showed similar patterns of mineralization of calcium carbonate forms (calcite and aragonite). Geographical variability in shell shape of C. concholepas described by discriminant functions was partially explained by environmental variables (pCO2, SST). This suggests the influence of corrosive waters, such as upwelling and freshwaters penetrating into the coastal ocean, upon spatial variation in shell morphology. Changes in the proportion of calcium carbonate forms precipitated by C. concholepas across their shells and its susceptibility to corrosive coastal waters are discussed.

(Table S1) Diatom abundance in coastal upwelling systems

Coastal upwelling systems account for approximately half of global ocean primary production and contribute disproportionately to biologically driven carbon sequestration. Diatoms, silica-precipitating microalgae, constitute the dominant phytoplankton in these productive regions, and their abundance and assemblage composition in the sedimentary record is considered one of the best proxies for primary production. The study of the sedimentary diatom abundance (SDA) and total organic carbon content (TOC) in the five most important coastal upwelling systems of the modern ocean (Iberia-Canary, Benguela, Peru-Humboldt, California and Somalia-Oman) reveals a global-scale positive relationship between diatom production and organic carbon burial. The analysis of SDA in conjunction with environmental variables of coastal upwelling systems such as upwelling strength, satellite-derived net primary production and surface water nutrient concentrations shows different relations between SDA and primary production on the regional scale. At the global-scale, SDA appears modulated by the capacity of diatoms to take up silicic acid, which ultimately sets an upper limit to global export production in these ocean regions.

Temperature is the evil twin: Effects of increased temperature and ocean acidification on reproduction in a reef fish

Reproduction in many organisms can be disrupted by changes to the physical environment, such as those predicted to occur during climate change. Marine organisms face the dual climate change threats of increasing temperature and ocean acidification, yet no studies have examined the potential interactive effects of these stressors on reproduction in marine fishes. We used a long-term experiment to test the interactive effects of increased temperature and CO2 on the reproductive performance of the anemonefish, Amphiprion melanopus. Adult breeding pairs were kept for 10 months at three temperatures, 28.5°C (+0.0°C), 30.0°C (+1.5°C) and 31.5°C (+3.0°C), cross-factored with 3 CO2 levels, a current day control (417 µatm) and moderate (644 µatm) and high (1134 µatm) treatments consistent with the range of CO2 projections for the year 2100 under RCP8.5. We recorded each egg clutch produced during the breeding season, the number of eggs laid per clutch, average egg size, fertilization success, survival to hatching, hatchling length and yolk provisioning. Adult body condition, hepatosomatic index, gonadosomatic index, and plasma 17beta-estradiol concentrations were measured at the end of the breeding season to determine the effect of prolonged exposure to increased temperature and elevated CO2 on adults, and to examine potential physiological mechanisms for changes in reproduction. Temperature had by far the stronger influence on reproduction, with clear declines in reproduction occurring in the +1.5°C treatment and ceasing altogether in the +3.0°C treatment. In contrast, CO2 had a minimal effect on the majority of reproductive traits measured, but caused a decline in offspring quality in combination with elevated temperature. We detected no significant effect of temperature or CO2 on adult body condition or hepatosomatic index. Elevated temperature had a significant negative effect on plasma 17beta-estradiol concentrations, suggesting that declines in reproduction with increasing temperature were due to the thermal sensitivity of reproductive hormones rather than a reduction in energy available for reproduction. Our results show that elevated temperature exerts a stronger influence than high CO2 on reproduction in A. melanopus. Understanding how these two environmental variables interact to affect the reproductive performance of marine organisms will be important for predicting the future impacts of climate change.

Dynamics of halocarbons in coastal surface waters during short term mesocosm experiments

The aim of the present study was to evaluate the influence of different light quality, especially ultraviolet radiation (UVR), on the dynamics of volatile halogenated organic compounds (VHOCs) at the sea surface. Short term experiments were conducted with floating gas-tight mesocosms of different optical qualities. Six halocarbons (CH3I, CHCl3, CH2Br2, CH2ClI, CHBr3 and CH2I2), known to be produced by phytoplankton, together with a variety of biological and environmental variables were measured in the coastal southern Baltic Sea and in the Raunefjord (North Sea). These experiments showed that ambient levels of UVR have no significant influence on VHOC dynamics in the natural systems. We attribute it to the low radiation doses that phytoplankton cells receive in a normal turbulent surface mixed layer. The VHOC concentrations were influenced by their production and removal processes, but they were not correlated with biological or environmental parameters investigated. Diatoms were most likely the dominant biogenic source of VHOCs in the Baltic Sea experiment, whereas in the Raunefjord experiment macroalgae probably contributed strongly to the production of VHOCs. The variable stable carbon isotope signatures (d13C values) of bromoform (CHBr3) also indicate that different autotrophic organisms were responsible for CHBr3 production in the two coastal environments. In the Raunefjord, despite strong daily variations in CHBr3 concentration, the carbon isotopic ratio was fairly stable with a mean value of -26 per mil. During the declining spring phytoplankton bloom in the Baltic Sea, the d13C values of CHBr3 were enriched in 13C and showed noticeable diurnal changes (-12 per mil ± 4). These results show that isotope signature analysis is a useful tool to study both the origin and dynamics of VHOCs in natural systems.

(Table B2) Relative abundances of diatom species used in the Canoco analysis at sediment surface samples from the Chinese inshore waters

Canonical correspondence analysis has been used to analyze and to visualize the relationships between the main species and selected environmental variables in a study of diatoms from surface sediment samples in Chinese inshore waters. The result shows that the diatom distribution in Chinese inshore waters is closely correlated with the environmental variables and that the measured environmental variables account for the major changes of the diatom composition. Winter sea-surface temperature (WST), winter sea-surface salinity (WSS), water depth and summer sea-surface salinity (SSS) play an important role for the diatom distribution. Among the environmental factors, winter sea-surface temperature is the most important, controlling the distribution of diatoms in the surface sediments in Chinese inshore waters, and therefore, it may be potentially reconstructed in palaeoceanographic studies. Three diatom assemblages are distinguished, representing environments with different hydrological characteristics. The temperate-water diatom assemblage may be used as an indicator of the coastal circulation system of Bohai Sea and Yellow Sea. While the warm-temperate water diatom assemblage is closely related to Shanghai-Zhejiang-Fujian coastal currents and Northern Bay coastal currents of South China Sea. The deep water diatom assemblage is a response to that the waters are less controlled by coastal currents, but are more influenced by open sea currents, such as Kuroshio.

Spatial distribution patterns of ascidians (Ascidiacea: Tunicata) in combination with information on bathymetry, oceanography, chlorophyll-a, and sea-ice on the continental shelves off the northern Antarctic Peninsula during POLARSTERN cruise ANT-XXIX/3

Ascidians (Ascidiacea: Tunicata) are sessile suspension feeders that represent dominant epifaunal components of the Southern Ocean shelf benthos and play a significant role in the pelagic-benthic coupling. Here, we report the results of a first study on the relationship between the distribution patterns of eight common and/or abundant (putative) ascidian species, and environmental drivers in the waters off the northern Antarctic Peninsula. During RV Polarstern cruise XXIX/3 (PS81) in January-March 2013, we used seabed imaging surveys along 28 photographic transects of 2 km length each at water depths from 70 to 770 m in three regions (northwestern Weddell Sea, southern Bransfield Strait and southern Drake Passage), differing in their general environmental setting, primarily oceanographic characteristics and sea-ice dynamics, to comparatively analyze the spatial patterns in the abundance of the selected ascidians, reliably to be identified in the photographs, at three nested spatial scales. At a regional (100-km) scale, the ascidian assemblages of the Weddell Sea differed significantly from those of the other two regions, whereas at an intermediate 10-km scale no such differences were detected among habitat types (bank, upper slope, slope, deep/canyon) on the shelf and at the shelf break within each region. These spatial patterns were superimposed by a marked small-scale (10-m) patchiness of ascidian distribution within the 2-km-long transects. Among the environmental variables considered in our study, a combination of water-mass characteristics, sea-ice dynamics (approximated by 5-year averages in sea-ice cover in the region of or surrounding the photographic stations), as well as the seabed ruggedness, was identified as explaining best the distribution patterns of the ascidians.

Growth and carbon metabolism of the seaweed Fucus vesiculosus in the western Baltic Sea ("Kiel Outdoor Benthocosms")

Warming and acidification of the oceans as a consequence of increasing CO2-concentrations occur at large scales. Numerous studies have shown the impact of single stressors on individual species. However, studies on the combined effect of multiple stressors on a multi-species assemblage, which is ecologically much more realistic and relevant, are still scarce. Therefore, we orthogonally crossed the two factors warming and acidification in mesocosm experiments and studied their single and combined impact on the brown alga Fucus vesiculosus associated with its natural community (epiphytes and mesograzers) in the Baltic Sea in all seasons (from April 2013 to April 2014). We superimposed our treatment factors onto the natural fluctuations of all environmental variables present in the Benthocosms in so-called delta-treatments. Thereby we compared the physiological responses of F. vesiculosus (growth and metabolites) to the single and combined effects of natural Kiel Fjord temperatures and pCO2 conditions with a 5 °C temperature increase and/or pCO2 increase treatment (1100 ppm in the headspace above the mesocosms). Responses were also related to the factor photoperiod which changes over the course of the year. Our results demonstrate complex seasonal pattern. Elevated pCO2 positively affected growth of F. vesiculosus alone and/or interactively with warming. The response direction (additive, synergistic or antagonistic), however, depended on season and daylength. The effects were most obvious when plants were actively growing during spring and early summer. Our study revealed for the first time that it is crucial to always consider the impact of variable environmental conditions throughout all seasons. In summary, our study indicates that in future F. vesiculosus will be more affected by detrimental summer heat-waves than by ocean acidification although the latter consequently enhances growth throughout the year. The mainly negative influence of rising temperatures on the physiology of this keystone macroalga may alter and/or hamper its ecological functions in the shallow coastal ecosystem of the Baltic Sea.

Drivers of multi-scale plant community structure in agricultural field margins of South Korea

This data set describes the distribution of a total of 90 plant species growing on field margins of an agricultural landscape in the Haean-myun catchment in South Korea. We conducted our survey between July and August 2011 in 100 sampling plots, covering the whole catchment. In each plot we measured three environmental variables: slope, width of the field margin, and management type (i.e. "managed" for field margins that had signs of management activities from the ongoing season such as cutting or spraying herbicides and "unmanaged" for field margins that had been left untouched in the season). For the botanical survey each plot was sampled using three subplots of one square meter per subplot; subplots were 4 m apart from each other. In each subplot, we estimated three different vegetation characteristics: vegetation cover (i.e. the percentage of ground covered by vegetation), species richness (i.e. the number of observed species) and species abundance (i.e. the number of observed individuals / species). We calculated the percentage of the non-farmed habitats by creating buffer zones of 100, 200, 300, 400 and 500 m radii around each plot using data provided by (Seo et al. 2014). Non-farmed habitats included field margins, fallows, forest, riparian areas, pasture and grassland.

Oxygen content and paleoproductivity from surface sediment samples Y69-86 to PLDS-004G

Quantitative estimation of surface ocean productivity and bottom water oxygen concentration with benthic foraminifera was attempted using 70 samples from equatorial and North Pacific surface sediments. These samples come from a well defined depth range in the ocean, between 2200 and 3200 m, so that depth related factors do not interfere with the estimation. Samples were selected so that foraminifera were well preserved in the sediments and temperature and salinity were nearly uniform (T = 1.5° C; S = 34.6 per mil). The sample set was also assembled so as to minimize the correlation often seen between surface ocean productivity and bottom water oxygen values (r**2 = 0.23 for prediction purposes in this case). This procedure reduced the chances of spurious results due to correlations between the environmental variables. The samples encompass a range of productivities from about 25 to >300 gC m**-2 yr**-1, and a bottom water oxygen range from 1.8 to 3.5 ml/L. Benthic foraminiferal assemblages were quantified using the >62 µm fraction of the sediments and 46 taxon categories. MANOVA multivariate regression was used to project the faunal matrix onto the two environmental dimensions using published values for productivity and bottom water oxygen to calibrate this operation. The success of this regression was measured with the multivariate r? which was 0.98 for the productivity dimension and 0.96 for the oxygen dimension. These high coefficients indicate that both environmental variables are strongly imbedded in the faunal data matrix. Analysis of the beta regression coefficients shows that the environmental signals are carried by groups of taxa which are consistent with previous work characterizing benthic foraminiferal responses to productivity and bottom water oxygen. The results of this study suggest that benthic foraminiferal assemblages can be used for quantitative reconstruction of surface ocean productivity and bottom water oxygen concentrations if suitable surface sediment calibration data sets are developed and appropriate means for detecting no-analog samples are found.

Diatom abundance in surface sediments of the northern North Pacific

In order to map the modern distribution of diatoms and to establish a reliable reference data set for paleoenvironmental reconstruction in the northern North Pacific, a new data set including the relative abundance of diatom species preserved in a total of 422 surface sediments was generated, which covers a broad range of environmental variables characteristic of the subarctic North Pacific, the Sea of Okhotsk and the Bering Sea between 30° and 70°N. The biogeographic distribution patterns as well as the preferences in sea surface temperature of 38 diatom species and species groups are documented. A Q-mode factor analysis yields a three-factor model representing assemblages associated with the Arctic, Subarctic and Subtropical water mass, indicating a close relationship between the diatom composition and the sea surface temperatures. The relative abundance pattern of 38 diatom species and species groups was statistically compared with nine environmental variables, i.e. the summer sea surface temperature and salinity, annual surface nutrient concentration (nitrate, phosphate, silicate), summer and winter mixed layer depth and summer and winter sea ice concentrations. Canonical Correspondence Analysis (CCA) indicates 32 species and species groups have strong correspondence with the pattern of summer sea surface temperature. In addition, the total diatom flux data compiled from ten sediment traps reveal that the seasonal signals preserved in the surface sediments are mostly from spring through autumn. This close relationship between diatom composition and the summer sea surface temperature will be useful in deriving a transfer function in the subarctic North Pacific for the quantitative paleoceanographic and paleoenvironmental studies. The relative abundance of the sea-ice indicator diatoms Fragilariopsis cylindrus and F. oceanica of >20% in the diatom composition is used to represent the winter sea ice edge in the Bering Sea. The northern boundary of the distribution of F. doliolus in the open ocean is suggested to be an indicator of the Subarctic Front, while the abundance of Chaetoceros resting spores may indicate iron input from nearby continents and shelves and induced productivity events in the study area.