(Appendix) Vertebral number of topsmelt, Atherinops affinis, adults and offspring from different latitudes along the North-American Pacific coast

Organisms that are distributed across spatial climate gradients often exhibit adaptive local variations in morphological and physiological traits, but to what extent such gradients shape evolutionary responses is still unclear. Given the strong natural contrast in latitudinal temperature gradients between the North-American Pacific and Atlantic coast, we asked how increases in vertebral number (VN, known as Jordan's Rule) with latitude would differ between Pacific (Atherinops affinis) and Atlantic Silversides (Menidia menidia), two ecologically equivalent and taxonomically similar fishes with similar latitudinal distributions. VN was determined from radiographs of wild-caught adults (genetic + environmental differences) and its genetic basis confirmed by rearing offspring in common garden experiments. Compared to published data on VN variation in M. menidia (a mean increase of 7.0 vertebrae from 32 to 46°N, VN slope = 0.42/lat), the latitudinal VN increase in Pacific Silversides was approximately half as strong (a mean increase of 3.3 vertebrae from 28 to 43°N, VN slope = 0.23/lat). This mimicked the strong Atlantic (1.11°C/lat) versus weak Pacific latitudinal gradient (0.40°C/lat) in median annual sea surface temperature (SST). Importantly, the relationship of VN to SST was not significantly different between the two species (average slope = -0.39 vertebrae/°C), thus suggesting a common thermal dependency of VN in silverside fishes. Our findings provide novel support for the hypothesis that temperature gradients are the ultimate cause of Jordan's Rule, even though its exact adaptive significance remains speculative. A second investigated trait, the mode of sex determination in Atlantic versus Pacific Silversides, revealed patterns that were inconsistent with our expectation: M. menidia displays temperature-dependent sex determination (TSD) at low latitudes, where growing seasons are long or unconstrained, but also a gradual shift to genetic sex determination (GSD) with increasing latitude due to more and more curtailed growing seasons. Sex ratios in A. affinis, on the other hand, were independent of latitude and rearing temperature (indicating GSD), even though growing seasons are thermally unconstrained across most of the geographical distribution of A. affinis. This suggests that additional factors (e.g., longevity) play an important role in shaping the mode of sex determination in silverside fishes.

Database of Ice-Rich Yedoma Permafrost (IRYP)

Vast portions of Arctic and sub-Arctic Siberia, Alaska and the Yukon Territory are covered by ice-rich silty to sandy deposits that are containing large ice wedges, resulting from syngenetic sedimentation and freezing. Accompanied by wedge-ice growth in polygonal landscapes, the sedimentation process was driven by cold continental climatic and environmental conditions in unglaciated regions during the late Pleistocene, inducing the accumulation of the unique Yedoma deposits up to >50 meters thick. Because of fast incorporation of organic material into syngenetic permafrost during its formation, Yedoma deposits include well-preserved organic matter. Ice-rich deposits like Yedoma are especially prone to degradation triggered by climate changes or human activity. When Yedoma deposits degrade, large amounts of sequestered organic carbon as well as other nutrients are released and become part of active biogeochemical cycling. This could be of global significance for future climate warming as increased permafrost thaw is likely to lead to a positive feedback through enhanced greenhouse gas fluxes. Therefore, a detailed assessment of the current Yedoma deposit coverage and its volume is of importance to estimate its potential response to future climate changes. We synthesized the map of the coverage and thickness estimation, which will provide critical data needed for further research. In particular, this preliminary Yedoma map is a great step forward to understand the spatial heterogeneity of Yedoma deposits and its regional coverage. There will be further applications in the context of reconstructing paleo-environmental dynamics and past ecosystems like the mammoth-steppe-tundra, or ground ice distribution including future thermokarst vulnerability. Moreover, the map will be a crucial improvement of the data basis needed to refine the present-day Yedoma permafrost organic carbon inventory, which is assumed to be between 83±12 (Strauss et al., 2013, doi:10.1002/2013GL058088) and 129±30 (Walter Anthony et al., 2014, doi:10.1038/nature13560) gigatonnes (Gt) of organic carbon in perennially-frozen archives. Hence, here we synthesize data on the circum-Arctic and sub-Arctic distribution and thickness of Yedoma for compiling a preliminary circum-polar Yedoma map. For compiling this map, we used (1) maps of the previous Yedoma coverage estimates, (2) included the digitized areas from Grosse et al. (2013) as well as extracted areas of potential Yedoma distribution from additional surface geological and Quaternary geological maps (1.: 1:500,000: Q-51-V,G; P-51-A,B; P-52-A,B; Q-52-V,G; P-52-V,G; Q-51-A,B; R-51-V,G; R-52-V,G; R-52-A,B; 2.: 1:1,000,000: P-50-51; P-52-53; P-58-59; Q-42-43; Q-44-45; Q-50-51; Q-52-53; Q-54-55; Q-56-57; Q-58-59; Q-60-1; R-(40)-42; R-43-(45); R-(45)-47; R-48-(50); R-51; R-53-(55); R-(55)-57; R-58-(60); S-44-46; S-47-49; S-50-52; S-53-55; 3.: 1:2,500,000: Quaternary map of the territory of Russian Federation, 4.: Alaska Permafrost Map). The digitalization was done using GIS techniques (ArcGIS) and vectorization of raster Images (Adobe Photoshop and Illustrator). Data on Yedoma thickness are obtained from boreholes and exposures reported in the scientific literature. The map and database are still preliminary and will have to undergo a technical and scientific vetting and review process. In their current form, we included a range of attributes for Yedoma area polygons based on lithological and stratigraphical information from the original source maps as well as a confidence level for our classification of an area as Yedoma (3 stages: confirmed, likely, or uncertain). In its current version, our database includes more than 365 boreholes and exposures and more than 2000 digitized Yedoma areas. We expect that the database will continue to grow. In this preliminary stage, we estimate the Northern Hemisphere Yedoma deposit area to cover approximately 625,000 km². We estimate that 53% of the total Yedoma area today is located in the tundra zone, 47% in the taiga zone. Separated from west to east, 29% of the Yedoma area is found in North America and 71 % in North Asia. The latter include 9% in West Siberia, 11% in Central Siberia, 44% in East Siberia and 7% in Far East Russia. Adding the recent maximum Yedoma region (including all Yedoma uplands, thermokarst lakes and basins, and river valleys) of 1.4 million km² (Strauss et al., 2013, doi:10.1002/2013GL058088) and postulating that Yedoma occupied up to 80% of the adjacent formerly exposed and now flooded Beringia shelves (1.9 million km², down to 125 m below modern sea level, between 105°E - 128°W and >68°N), we assume that the Last Glacial Maximum Yedoma region likely covered more than 3 million km² of Beringia. Acknowledgements: This project is part of the Action Group "The Yedoma Region: A Synthesis of Circum-Arctic Distribution and Thickness" (funded by the International Permafrost Association (IPA) to J. Strauss) and is embedded into the Permafrost Carbon Network (working group Yedoma Carbon Stocks). We acknowledge the support by the European Research Council (Starting Grant #338335), the German Federal Ministry of Education and Research (Grant 01DM12011 and "CarboPerm" (03G0836A)), the Initiative and Networking Fund of the Helmholtz Association (#ERC-0013) and the German Federal Environment Agency (UBA, project UFOPLAN FKZ 3712 41 106).

Incorporation of Mg and Sr in calcite of cultured benthic foraminifera (Heterostegina depressa and Ammonia tepida) and seawater carbonate chemistry, 2010

We investigated the effect of the calcium concentration in seawater and thereby the calcite saturation state (omega) on the magnesium and strontium incorporation into benthic foraminiferal calcite under laboratory conditions. For this purpose individuals of the shallow-water species Heterostegina depressa (precipitating high-Mg calcite, symbiont-bearing) and Ammonia tepida (low-Mg calcite, symbiont-barren) were cultured in media under a range of [Ca2+], but similar Mg/Ca ratios. Trace element/Ca ratios of newly formed calcite were analysed with Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) and normalized to the seawater elemental composition using the equation DTE=(TE/Cacalcite)/(TE/Caseawater). The culturing study shows that DMg of A. tepida significantly decreases with increasing omega at a gradient of -4.3x10-5 per omega unit. The DSr value of A. tepida does not change with omega, suggesting that fossil Sr/Ca in this species may be a potential tool to reconstruct past variations in seawater Sr/Ca. Conversely, DMg of H. depressa shows only a minor decrease with increasing omega, while DSr increases considerably with omega at a gradient of 0.009 per omega unit. The different responses to seawater chemistry of the two species may be explained by a difference in the calcification pathway that is, at the same time, responsible for the variation in the total Mg incorporation between the two species. Since the Mg/Ca ratio in H. depressa is 50-100 times higher than that of A. tepida, it is suggested that the latter exhibits a mechanism that decreases the Mg/Ca ratio of the calcification fluid, while the high-Mg calcite forming species may not have this physiological tool. If the dependency of Mg incorporation on seawater [Ca2+] is also valid for deep-sea benthic foraminifera typically used for paleostudies, the higher Ca concentrations in the past may potentially bias temperature reconstructions to a considerable degree. For instance, 25 Myr ago Mg/Ca ratios in A. tepida would have been 0.2 mmol/mol lower than today, due to the 1.5 times higher [Ca2+] of seawater, which in turn would lead to a temperature underestimation of more than 2 °C.

Combined effects of ocean acidification and light or nitrogen availabilities on 13C fractionation in marine dinoflagellates

Along with increasing oceanic CO2 concentrations, enhanced stratification constrains phytoplankton to shallower upper mixed layers with altered light regimes and nutrient concentrations. Here, we investigate the effects of elevated pCO2 in combination with light or nitrogen-limitation on 13C fractionation (epsilon p) in four dinoflagellate species. We cultured Gonyaulax spinifera and Protoceratium reticulatum in dilute batches under low-light (LL) and high-light (HL) conditions, and grew Alexandrium fundyense and Scrippsiella trochoidea in nitrogen-limited continuous cultures (LN) and nitrogen-replete batches (HN). The observed CO2-dependency of epsilon p remained unaffected by the availability of light for both G. spinifera and P. reticulatum, though at HL epsilon p was consistently lower by about 2.7 per mil over the tested CO2 range for P. reticulatum. This may reflect increased uptake of (13C-enriched) bicarbonate fueled by increased ATP production under HL conditions. The observed CO2-dependency of epsilon p disappeared under LN conditions in both A. fundyense and S. trochoidea. The generally higher epsilon p under LN may be associated with lower organic carbon production rates and/or higher ATP:NADPH ratios. CO2-dependent epsilon p under non-limiting conditions has been observed in several dinoflagellate species, showing potential for a new CO2-proxy. Our results however demonstrate that light- and nitrogen-limitation also affect epsilon p, thereby illustrating the need to carefully consider prevailing environmental conditions.

Seawater carbonate chemistry, calcification and shell size of hard clam Mercenaria mercenaria during experiments, 2011

Increasing atmospheric carbon dioxide threatens to decrease pH in the world's oceans. Coastal and estuarine calcifying organisms of significant ecological and economical importance are at risk; however, several biogeochemical processes drive pH in these habitats. In particular, coastal and estuarine sediments are frequently undersaturated with respect to calcium carbonate due to high rates of organic matter remineralization, even when overlying waters are saturated. As a result, the post-larval stages of infaunal marine bivalves must be able to deposit new shell material in conditions that are corrosive to shell. We measured calcification rates on the hard clam, Mercenaria spp.,in 5 post-larval size classes (0.39, 0.56, 0.78, 0.98, and 2.90 mm shell height) using the alkalinity anomaly method. Acidity of experimental water was controlled by bubbling with air-CO2 blends to obtain pH values of 8.02, 7.64, and 7.41, corresponding to pCO2 values of 424, 1120, and 1950 µatm. These pH values are typical of those found in many near-shore terrigenous marine sediments. Our results show that calcification rate decreased with lower pH in all 5 size classes measured. We also found a significant effect of size on calcification rate, with the smaller post-larval sizes unable to overcome dissolution pressure. Increased calcification rate with size allowed the larger sizes to overcome dissolution pressure and deposit new shell material under corrosive conditions. Size dependency of pH effects on calcification is likely due to organogenesis and developmental shifts in shell mineralogy occurring through the post-larval stage. Furthermore, we found significantly different calcification rates between the 2 sources of hard clams we used for these experiments, most likely due to genotypic differences. Our findings confirm the susceptibility of the early life stages of this important bivalve to decreasing pH and reveal mechanisms behind the increased mortality in post-larval juvenile hard clams related to dissolution pressure, that has been found in previous studies.

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.

30 years of synoptic observations from Neumayer Station with links to datasets

The analysis of research data plays a key role in data-driven areas of science. Varieties of mixed research data sets exist and scientists aim to derive or validate hypotheses to find undiscovered knowledge. Many analysis techniques identify relations of an entire dataset only. This may level the characteristic behavior of different subgroups in the data. Like automatic subspace clustering, we aim at identifying interesting subgroups and attribute sets. We present a visual-interactive system that supports scientists to explore interesting relations between aggregated bins of multivariate attributes in mixed data sets. The abstraction of data to bins enables the application of statistical dependency tests as the measure of interestingness. An overview matrix view shows all attributes, ranked with respect to the interestingness of bins. Complementary, a node-link view reveals multivariate bin relations by positioning dependent bins close to each other. The system supports information drill-down based on both expert knowledge and algorithmic support. Finally, visual-interactive subset clustering assigns multivariate bin relations to groups. A list-based cluster result representation enables the scientist to communicate multivariate findings at a glance. We demonstrate the applicability of the system with two case studies from the earth observation domain and the prostate cancer research domain. In both cases, the system enabled us to identify the most interesting multivariate bin relations, to validate already published results, and, moreover, to discover unexpected relations.

(Tables 1,2) Site characteristics and phenotypic attributes of King Penguins (Aptenodytes patagonicus) on Possession Island

In colonial species, it is often assumed that locations in the center of the colony are of highest quality and provide highest breeding success. We tested this prediction, known as the "central-periphery model," in a King Penguin colony in the subantarctic Crozet Archipelago. Breeding activity and survival of 150 penguins, fitted with transponder tags, were monitored over an entire breeding season. Among these 150 birds, 50 bred on the slope at the upper periphery of the colony, where the rates of predation and parasitism by ticks were high. Fifty birds bred in the center of the colony, where rates of predation and tick parasitism were low, and 50 bred at the lower end of the colony, where the rate of tick parasitism was low but predation and flooding were important risks. We predicted that the center of the colony should provide the safest breeding place and consequently be characterized by the highest breeding success and be used by the highest-quality individuals. Yet we found that penguins breeding in the center of the colony had the same breeding success as those at both peripheral locations. In addition, penguins breeding on the upper slope had a higher survival rate than penguins breeding at the center or bottom of the slope and were likely of higher quality. Our study does not support the central-periphery model and emphasizes the complexity behind the relationships among breeding site, breeding success, and individual quality.