Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Collection of data on particulate and dissolved soil phosphorus in the Jena Experiment (Main Experiment, time series since 2002)

This data set contains two time series of measurements of dissolved phosphorus (organic, inorganic and total with a biweekly resolution) and dissolved inorganic phosphorus with a seasonal resolution. In addition, data on phosphorus from soil samples measured in 2007 and fractionated by different acid-extrations (Hedley fractions) are provided. All data measured at 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. 1. Dissolved phosphorus in soil solution: Suction plates installed on the field site in 10, 20, 30 and 60 cm depth were used to sample soil pore water. Cumulatively extracted soil solution was collected every two weeks from October 2002 to May 2006. The biweekly samples from 2002, 2003 and 2004 were analyzed for dissolved organic phosphorus (DOP), dissolved inorganic phosphorus (PO4P) and dissolved total phosphorus (TDP) by Continuous Flow Analyzer (CFA SAN ++, SKALAR [Breda, The Netherlands]). 2. Seasonal values of dissolved inorganic phosphorus in soil solution were calculated as volume-weighted mean values of the biweekly measurements (spring = March to May, summer = June to August, fall = September to November, winter = December to February). 3. Phosphorus fractions in soil: Five independent soil samples per plot were taken in a depth of 0-15 cm using a soil corer with an inner diameter of 1 cm. The five samples per plot were combined to one composite sample per plot. A four-step sequential P fractionation (Hedley fractions) was applied and concentrations of P fractions in soil were measured photometrically (molybdenum blue-reactive P) with a Continuous Flow Analyzer (Bran&Luebbe, Germany).

Proxy data and model simulations of sediment core GeoB7920-2

The evolution of the northwest African hydrological balance throughout the Pleistocene epoch influenced the migration of prehistoric humans**1. The hydrological balance is also thought to be important to global teleconnection mechanisms during Dansgaard-Oeschger and Heinrich events**2. However, most high-resolution African climate records do not span the millennial-scale climate changes of the last glacial-interglacial cycle**1, 3, 4, 5, or lack an accurate chronology**6. Here, we use grain-size analyses of siliciclastic marine sediments from off the coast of Mauritania to reconstruct changes in northwest African humidity over the past 120,000 years. We compare this reconstruction to simulations of palaeo-humidity from a coupled atmosphere-ocean-vegetation model. These records are in good agreement, and indicate the reoccurrence of precession-forced humid periods during the last interglacial period similar to the Holocene African Humid Period. We suggest that millennial-scale arid events are associated with a reduction of the North Atlantic meridional overturning circulation and that millennial-scale humid events are linked to a regional increase of winter rainfall over the coastal regions of northwest Africa.

PumaPlex 1.0: Data generated during the development of the SNP markers and SNP genotypes of pumas

Pumas are one of the most studied terrestrial mammals because of their widespread distribution, substantial ecological impacts, and conflicts with humans. Extensive efforts, often employing genetic methods, are undertaken to manage this species. However, the comparison of population genetic data is difficult because few of the microsatellite loci chosen are shared across research programs. Here, we describe the development of PumaPlex, a high-throughput assay to genotype 25 single nucleotide polymorphisms in pumas. We validated PumaPlex in more than 700 North American pumas (Puma concolor couguar), and demonstrated its ability to generate reproducible genotypes and accurately identify individuals. Furthermore, we compared PumaPlex with traditional genotyping of 12 microsatellite loci in fecal DNA samples and found that PumaPlex produced significantly more genotypes with fewer false alleles. PumaPlex promotes the cross-laboratory comparison of genotypes, is easily expandable in the future, and is a valuable tool for the genetic monitoring and management of North American puma populations.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

Collection of data on particulate and dissolved soil nitrogen in the Jena Experiment (Main Experiment, time series since 2002)

This data set contains four time series of particulate and dissolved soil nitrogen measurements from 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. 1. Total nitrogen from solid phase: Stratified soil sampling was performed every two years since before sowing in April 2002 and was repeated in April 2004, 2006 and 2008 to a depth of 30 cm segmented to a depth resolution of 5 cm giving six depth subsamples per core. In 2002 five samples per plot were taken and analyzed independently. Averaged values per depth layer are reported. In later years, three samples per plot were taken, pooled in the field, and measured as a combined sample. Sampling locations were less than 30 cm apart from sampling locations in other years. All soil samples were passed through a sieve with a mesh size of 2 mm in 2002. In later years samples were further sieved to 1 mm. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany). 2. Total nitrogen from solid phase (high intensity sampling): In block 2 of the Jena Experiment, soil samples were taken to a depth of 1m (segmented to a depth resolution of 5 cm giving 20 depth subsamples per core) with three replicates per block ever 5 years starting before sowing in April 2002. Samples were processed as for the more frequent sampling but were always analyzed independently and never pooled. 3. Mineral nitrogen from KCl extractions: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m (and between 2002 and 2004 also at a depth of 0.15 to 0.3 m) of the mineral soil from each of the experimental plots at various times over the years. In addition also plots of the management experiment, that altered mowing frequency and fertilized subplots (see further details below) were sampled in some later years. Samples of the soil cores per plot (subplots in case of the management experiment) were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, 2003-2005: Skalar, Breda, Netherlands; 2006-2007: AutoAnalyzer, Seal, Burgess Hill, United Kingdom). 4. Dissolved nitrogen in soil solution: Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for nitrate (NO3-), ammonium (NH4+) and total dissolved nitrogen concentrations with a continuous flow analyzer (CFA, Skalar, Breda, The Netherlands). Nitrate was analyzed photometrically after reduction to NO2- and reaction with sulfanilamide and naphthylethylenediamine-dihydrochloride to an azo-dye. Our NO3- concentrations contained an unknown contribution of NO2- that is expected to be small. Simultaneously to the NO3- analysis, NH4+ was determined photometrically as 5-aminosalicylate after a modified Berthelot reaction. The detection limits of NO3- and NH4+ were 0.02 and 0.03 mg N L-1, respectively. Total dissolved N in soil solution was analyzed by oxidation with K2S2O8 followed by reduction to NO2- as described above for NO3-. Dissolved organic N (DON) concentrations in soil solution were calculated as the difference between TDN and the sum of mineral N (NO3- + NH4+).

Adverse effects of ocean acidification on early development of squid (Doryteuthis pealeii)

Anthropogenic carbon dioxide (CO2) is being absorbed into the ocean, altering seawater chemistry, with potentially negative impacts on a wide range of marine organisms. The early life stages of invertebrates with internal and external aragonite structures may be particularly vulnerable to this ocean acidification. Impacts to cephalopods, which form aragonite cuttlebones and statoliths, are of concern because of the central role they play in many ocean ecosystems and because of their importance to global fisheries. Atlantic longfin squid (Doryteuthis pealeii), an ecologically and economically valuable taxon, were reared from eggs to hatchlings (paralarvae) under ambient and elevated CO2 concentrations in replicated experimental trials. Animals raised under elevated pCO2 demonstrated significant developmental changes including increased time to hatching and shorter mantle lengths, although differences were small. Aragonite statoliths, critical for balance and detecting movement, had significantly reduced surface area and were abnormally shaped with increased porosity and altered crystal structure in elevated pCO2-reared paralarvae. These developmental and physiological effects could alter squid paralarvae behavior and survival in the wild, directly and indirectly impacting marine food webs and commercial fisheries.

Effects of pre-hatching pCO2 on development and survival of zoea stages of Arctic spider crab Hyas araneus

Sensitivity of marine crustaceans to anthropogenic CO2 emissions and the associated acidification of the oceans may be less than that of other, especially lower, invertebrates. However, effects on critical transition phases or carry-over effects between life stages have not comprehensively been explored. Here we report the impact of elevated seawater PCO2 values (3100 µatm) on Hyas araneus during the last 2 weeks of their embryonic development (pre-hatching phase) and during development while in the consecutive zoea I and zoea II larval stages (post-hatching phase). We measured oxygen consumption, dry weight, developmental time and mortality in zoea I to assess changes in performance. Feeding rates and survival under starvation were investigated at different temperatures to detect differences in thermal sensitivities of zoea I and zoea II larvae depending on pre-hatch history. When embryos were pre-exposed to elevated PCO2 during maternal care, mortality increased about 60% under continued CO2 exposure during the zoea I phase. The larvae that moulted into zoea II, displayed a developmental delay by about 20 days compared to larvae exposed to control PCO2 during embryonic and zoeal phases. Elevated PCO2 caused a reduction in zoea I dry weight and feeding rates, while survival of the starved larvae was not affected by the seawater CO2 concentration. In conclusion, CO2 effects on egg masses under maternal care carried over to the first larval stages of crustaceans and reduced their survival and development to levels below those previously reported in studies exclusively focussing on acute PCO2 effects on the larval stages.

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.