Sea urchins in a high CO2 world: partitioned effects of body-size, ocean warming and acidification on metabolic rate

Body-size and temperature are the major factors explaining metabolic rate, and the additional factor of pH is a major driver at the biochemical level. These three factors have frequently been found to interact, complicating the formulation of broad models predicting metabolic rates and hence ecological functioning. In this first study of the effects of warming and ocean acidification, and their potential interaction, on metabolic rate across a broad body-size range (two-to-three orders of magnitude difference in body mass) we addressed the impact of climate change on the sea urchin Heliocidaris erythrogramma in context with climate projections for east Australia, an ocean warming hotspot. Urchins were gradually introduced to two temperatures (18 and 23 °C) and two pH (7.5 and 8.0), and maintained for two months. That a new physiological steady-state had been reached, otherwise know as acclimation, was validated through identical experimental trials separated by several weeks. The relationship between body-size, temperature and acidification on the metabolic rate of H. erythrogramma was strikingly stable. Both stressors caused increases in metabolic rate; 20% for temperature and 19% for pH. Combined effects were additive; a 44% increase in metabolism. Body-size had a highly stable relationship with metabolic rate regardless of temperature or pH. None of these diverse drivers of metabolism interacted or modulated the effects of the others, highlighting the partitioned nature of how each influences metabolic rate, and the importance of achieving a full acclimation state. Despite these increases in energetic demand there was very limited capacity for compensatory modulating of feeding rate; food consumption increased only in the very smallest specimens, and only in response to temperature, and not pH. Our data show that warming, acidification and body-size all substantially affect metabolism and are highly consistent and partitioned in their effects, and for H. erythrogramma near-future climate change will incur a substantial energetic cost.

Impact of atmospheric pCO2, seawater temperature, and calcification rate on the delta 18O and delta 13C composition of echinoid calcite (Echinometra viridis)

The tropical echinoid Echinometra viridis was reared in controlled laboratory experiments at temperatures of approximately 20°C and 30°C to mimic winter and summer temperatures and at carbon dioxide (CO2) partial pressures of approximately 487 ppm-v and 805 ppm-v to simulate current and predicted-end-of-century levels. Spine material produced during the experimental period and dissolved inorganic carbon (DIC) of the corresponding culture solutions were then analyzed for stable oxygen (delta 18Oe, delta 18ODIC) and carbon (The tropical echinoid Echinometra viridis was reared in controlled laboratory experiments at temperatures of approximately 20°C and 30°C to mimic winter and summer temperatures and at carbon dioxide (CO2) partial pressures of approximately 487 ppm-v and 805 ppm-v to simulate current and predicted-end-of-century levels. Spine material produced during the experimental period and dissolved inorganic carbon (DIC) of the corresponding culture solutions were then analyzed for stable oxygen (delta18Oe, delta18ODIC) and carbon (delta13Ce, delta13CDIC) isotopic composition. Fractionation of oxygen stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta18O = delta18Oe - delta18ODIC) was significantly inversely correlated with seawater temperature but not significantly correlated with atmospheric pCO2. Fractionation of carbon stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (Delta delta13C = delta13Ce - delta13CDIC) was significantly positively correlated with pCO2 and significantly inversely correlated with temperature, with pCO2 functioning as the primary factor and temperature moderating the pCO2-delta13C relationship. Echinoid calcification rate was significantly inversely correlated with both delta18O and delta13C, both within treatments (i.e., pCO2 and temperature fixed) and across treatments (i.e., with effects of pCO2 and temperature controlled for through ANOVA). Therefore, calcification rate and potentially the rate of co-occurring dissolution appear to be important drivers of the kinetic isotope effects observed in the echinoid spines. Study results suggest that echinoid delta18O monitors seawater temperature, but not atmospheric pCO2, and that echinoid delta13C monitors atmospheric pCO2, with temperature moderating this relationship. These findings, coupled with echinoids' long and generally high-quality fossil record, supports prior assertions that fossil echinoid delta18O is a viable archive of paleo-seawater temperature throughout Phanerozoic time, and that delta13C merits further investigation as a potential proxy of paleo-atmospheric pCO2. However, the apparent impact of calcification rate on echinoid delta18O and delta13C suggests that paleoceanographic reconstructions derived from these proxies in fossil echinoids could be improved by incorporating the effects of growth rate.13Ce, delta13CDIC) isotopic composition. Fractionation of oxygen stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta18O = delta18Oe - delta18ODIC) was significantly inversely correlated with seawater temperature but not significantly correlated with atmospheric pCO2. Fractionation of carbon stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta13C = delta13Ce - delta13CDIC) was significantly positively correlated with pCO2 and significantly inversely correlated with temperature, with pCO2 functioning as the primary factor and temperature moderating the pCO2-delta13C relationship. Echinoid calcification rate was significantly inversely correlated with both delta18O and delta13C, both within treatments (i.e., pCO2 and temperature fixed) and across treatments (i.e., with effects of pCO2 and temperature controlled for through ANOVA). Therefore, calcification rate and potentially the rate of co-occurring dissolution appear to be important drivers of the kinetic isotope effects observed in the echinoid spines. Study results suggest that echinoid delta18O monitors seawater temperature, but not atmospheric pCO2, and that echinoid delta13C monitors atmospheric pCO2, with temperature moderating this relationship. These findings, coupled with echinoids' long and generally high-quality fossil record, supports prior assertions that fossil echinoid delta18O is a viable archive of paleo-seawater temperature throughout Phanerozoic time, and that delta13C merits further investigation as a potential proxy of paleo-atmospheric pCO2. However, the apparent impact of calcification rate on echinoid delta18O and delta13C suggests that paleoceanographic reconstructions derived from these proxies in fossil echinoids could be improved by incorporating the effects of growth rate.

Registry of selected samples from the Tara Oceans Expedition (2009-2013)

The present data set is a registry of samples from the Tara Oceans Expedition (2009-2013) that were selected for publication in a special issue of the SCIENCE journal (see related references below). The registry provides details about the sampling location and methodology of each sample. Uniform resource locators (URLs) offer direct links to additional contextual environmental data and to the corresponding sequence runs used for analysis in the related literature publications in the SCIENCE journal.

Biodiversity context of selected samples from the Tara Oceans Expedition (2009-2013)

The present data set provides contextual data for samples from the Tara Oceans Expedition (2009-2013) that were selected for publication in a special issue of the SCIENCE journal (see related references below). Contextual data include various diversity indexes calculated for the sampling location using satellite and model climatologies (Darwin project, Physat) and results from the sequencing of Tara Oceans samples.

Middle Eocene Fe and Ca counts, Ca/Fe ratios, revised splices and astronomical age models form different Sites of IODP Exp342

Accurate age models are a tool of utmost important in paleoclimatology. Constraining the rate and pace of past climate change are at the core of paleoclimate research, as such knowledge is crucial to our understanding of the climate system. Indeed, it allows for the disentanglement of the various drivers of climate change. The scarcity of highly resolved sedimentary records from the middle Eocene (Bartonian - Lutetian Stages; 47.8 - 37.8 Ma) has led to the existence of the "Eocene astronomical time scale gap" and hindered the establishment of a comprehensive astronomical time scale (ATS) for the entire Cenozoic. Sediments from the Newfoundland Ridge drilled during Integrated Ocean Drilling Program (IODP) Expedition 342 span the Eocene gap at an unprecedented stratigraphic resolution with carbonate bearing sediments. Moreover, these sediments exhibit cyclic lithological changes that allow for an astronomical calibration of geologic time. In this study, we use the dominant obliquity imprint in XRF-derived calcium-iron ratio series (Ca/Fe) from three sites drilled during IODP Expedition 342 (U1408, U1409, U1410) to construct a floating astrochronology. We then anchor this chronology to numerical geological time by tuning 173-kyr cycles in the amplitude modulation pattern of obliquity to an astronomical solution. This study is one of the first to use the 173-kyr obliquity amplitude cycle for astrochronologic purposes, as previous studies primarily use the 405-kyr long eccentricity cycle as a tuning target to calibrate the Paleogene geologic time scale. We demonstrate that the 173-kyr cycles in obliquity's amplitude are stable between 40 and 50 Ma, which means that one can use the 173-kyr cycle for astrochronologic calibration in the Eocene. Our tuning provides new age estimates for magnetochron reversals C18n.1n - C21r and a stratigraphic framework for key sites from Expedition 342 for the Eocene. Some disagreements emerge when we compare our tuning for the interval between C19r and C20r with previous tuning attempts from the South Atlantic. We therefore present a revision of the original astronomical interpretations for the latter records, so that the various astrochronologic age models for the middle Eocene in the North- and South-Atlantic are consistent.

Projecting the fate of fish populations using ecological-economic modeling

Four marine fish species are among the most important on the world market: cod, salmon, tuna, and sea bass. While the supply of North American and European markets for two of these species - Atlantic salmon and European sea bass - mainly comes from fish farming, Atlantic cod and tunas are mainly caught from wild stocks. We address the question what will be the status of these wild stocks in the midterm future, in the year 2048, to be specific. Whereas the effects of climate change and ecological driving forces on fish stocks have already gained much attention, our prime interest is in studying the effects of changing economic drivers, as well as the impact of variable management effectiveness. Using a process-based ecological-economic multispecies optimization model, we assess the future stock status under different scenarios of change. We simulate (i) technological progress in fishing, (ii) increasing demand for fish, and (iii) increasing supply of farmed fish, as well as the interplay of these driving forces under different sce- narios of (limited) fishery management effectiveness. We find that economic change has a substantial effect on fish populations. Increasing aquaculture production can dampen the fishing pressure on wild stocks, but this effect is likely to be overwhelmed by increasing demand and technological progress, both increasing fishing pressure. The only solution to avoid collapse of the majority of stocks is institutional change to improve management effectiveness significantly above the current state. We conclude that full recognition of economic drivers of change will be needed to successfully develop an integrated ecosystem management and to sustain the wild fish stocks until 2048 and beyond.

Abundance, biomass and secondary production of benthic megafauna on the Barents Sea shelf in 2008 and 2009

Megabenthos plays a major role in the overall energy flow on Arctic shelves, but information on megabenthic secondary production on large spatial scales is scarce. Here, we estimated for the first time megabenthic secondary production for the entire Barents Sea shelf by applying a species-based empirical model to an extensive dataset from the joint Norwegian? Russian ecosystem survey. Spatial patterns and relationships were analyzed within a GIS. The environmental drivers behind the observed production pattern were identified by applying an ordinary least squares regression model. Geographically weighted regression (GWR) was used to examine the varying relationship of secondary production and the environment on a shelfwide scale. Significantly higher megabenthic secondary production was found in the northeastern, seasonally ice-covered regions of the Barents Sea than in the permanently ice-free southwest. The environmental parameters that significantly relate to the observed pattern are bottom temperature and salinity, sea ice cover, new primary production, trawling pressure, and bottom current speed. The GWR proved to be a versatile tool for analyzing the regionally varying relationships of benthic secondary production and its environmental drivers (R² = 0.73). The observed pattern indicates tight pelagic? benthic coupling in the realm of the productive marginal ice zone. Ongoing decrease of winter sea ice extent and the associated poleward movement of the seasonal ice edge point towards a distinct decline of benthic secondary production in the northeastern Barents Sea in the future.

Collection of aboveground community and species-specific plant biomass from the Jena Experiment (time series since 2002).

This data set comprises time series of aboveground community plant 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 several experiments at the field site of a large grassland biodiversity experiment (the Jena Experiment; see further details below). Aboveground community biomass was normally harvested twice a year just prior to mowing (during peak standing biomass twice a year, generally in May and August; in 2002 only once in September) on all experimental plots in the Jena Experiment. This was done by clipping the vegetation at 3 cm above ground in up to four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned by random selection of new coordinates every year within the core area of the plots. 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 following series of datasets are contained in this collection: 1. Plant biomass form the Main Experiment: 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). 2. Plant biomass from the Dominance Experiment: In the Dominance Experiment, 206 grassland plots of 3.5 x 3.5 m were established from a pool of 9 species that can be dominant in semi-natural grassland communities of the study region. 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, 3, 4, 6, and 9 species). 3. Plant biomass from the monoculture plots: In the monoculture plots the sown plant community contains only a single species per plot and this species is a different one for each plot. Which species has been sown in which plot is stated in the plot information table for monocultures (see further details below). The monoculture plots of 3.5 x 3.5 m were established for all of the 60 plant species of the Jena Experiment species pool with two replicates per species like the other experiments in May 2002. All plots were maintained by bi-annual weeding and mowing.

Long-term chemical characterization of aerosols at CVAO from 2007 to 2011

The first long-term aerosol sampling and chemical characterization results from measurements at the Cape Verde Atmospheric Observatory (CVAO) on the island of São Vicente are presented and are discussed with respect to air mass origin and seasonal trends. In total 671 samples were collected using a high-volume PM10 sampler on quartz fiber filters from January 2007 to December 2011. The samples were analyzed for their aerosol chemical composition, including their ionic and organic constituents. Back trajectory analyses showed that the aerosol at CVAO was strongly influenced by emissions from Europe and Africa, with the latter often responsible for high mineral dust loading. Sea salt and mineral dust dominated the aerosol mass and made up in total about 80% of the aerosol mass. The 5-year PM10 mean was 47.1 ± 55.5 µg/m**2, while the mineral dust and sea salt means were 27.9 ± 48.7 and 11.1 ± 5.5 µg/m**2, respectively. Non-sea-salt (nss) sulfate made up 62% of the total sulfate and originated from both long-range transport from Africa or Europe and marine sources. Strong seasonal variation was observed for the aerosol components. While nitrate showed no clear seasonal variation with an annual mean of 1.1 ± 0.6 µg/m**3, the aerosol mass, OC (organic carbon) and EC (elemental carbon), showed strong winter maxima due to strong influence of African air mass inflow. Additionally during summer, elevated concentrations of OM were observed originating from marine emissions. A summer maximum was observed for non-sea-salt sulfate and was connected to periods when air mass inflow was predominantly of marine origin, indicating that marine biogenic emissions were a significant source. Ammonium showed a distinct maximum in spring and coincided with ocean surface water chlorophyll a concentrations. Good correlations were also observed between nss-sulfate and oxalate during the summer and winter seasons, indicating a likely photochemical in-cloud processing of the marine and anthropogenic precursors of these species. High temporal variability was observed in both chloride and bromide depletion, differing significantly within the seasons, air mass history and Saharan dust concentration. Chloride (bromide) depletion varied from 8.8 ± 8.5% (62 ± 42%) in Saharan-dust-dominated air mass to 30 ± 12% (87 ± 11%) in polluted Europe air masses. During summer, bromide depletion often reached 100% in marine as well as in polluted continental samples. In addition to the influence of the aerosol acidic components, photochemistry was one of the main drivers of halogenide depletion during the summer; while during dust events, displacement reaction with nitric acid was found to be the dominant mechanism. Positive matrix factorization (PMF) analysis identified three major aerosol sources: sea salt, aged sea salt and long-range transport. The ionic budget was dominated by the first two of these factors, while the long-range transport factor could only account for about 14% of the total observed ionic mass.