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.

Seawater carbonate chemistry and calcification at Anse des Cuivres, NW Mediterranean, 2004

The community metabolism of a shallow infralittoral ecosystem dominated by the calcareous macroalgae Corallina elongata was investigated in Marseilles (NW Mediterranean), by monitoring hourly changes of seawater pH and total alkalinity over 6 d in February 2000. Fair weather conditions prevailed over the study period as indicated by oceanographic (temperature, salinity, and current velocity and direction) and meteorological variables, which validated the standing water hypothesis. This temperate ecosystem exhibited high community gross primary production (GPP = 519 ± 106 mmol C m-2 d-1, n = 6) and also supported high rates of community respiration (R). As a result, the system was slightly autotrophic (net community production, NCP = 20 mmol C m-2 d-1), with a GPP/R ratio of 1.06. NCP exhibited circadian variations with 2- to 3-fold changes in community respiration, both in the light and in the dark. Rates of net community calcification also exhibited circadian variations, with positive rates (up to 24 mmol CaCO3 m-2 h-1) for irradiance values >300 W m-2 (about 1380 µmol photon m-2 s-1). Below this irradiance threshold, net community dissolution prevailed. Daily net calcification (G) was on average 8 mmol CaCO3 m-2 d-1. CO2 fluxes generated by primary production, respiration, and calcification suggest that the study site was a potential atmospheric CO2 sink of 15 mmol CO2 m-2 d-1 at the time of measurement.

Seawater carbonate chemistry and calcification during a tropical lagoon study at SW lagoon of New Caledonia, 1994

Carbon dioxide and oxygen fluxes were measured in 0.2 m2 enclosures placed at the water sediment interface in the SW lagoon of New Caledonia. Experiments, performed at several stations in a wide range of environments, were carried out both in darkness to estimate respiration and at ambient light, to assess the effects of primary production. The community respiratory quotient (CRQ = CO2 production rate/02 consumption rate) and the community photosynthetic quotient (CPQ= gross O2 production rate/gross CO2 consumption rate) were calculated by functional regressions. The CRQ value, calculated from 61 incubations, was 1.14 (S.E. 0.05) and the CPQ value, obtained from 18 incubations, was 1.03 (S.E. 0.08). The linearity of the relationship between the O2 and the CO2 fluxes suggests that these values are representative for the whole lagoon

Changes in the storage and sink of carbon dioxide in subsurface atmospheres controlled by climate-driven processes: the case of the Ojo Guareña karst system

A comprehensive environmental monitoring program was conducted in the Ojo Guareña cave system (Spain), one of the longest cave systems in Europe, to assess the magnitude of the spatiotemporal changes in carbon dioxide gas (CO2) in the cave–soil–atmosphere profile. The key climate-driven processes involved in gas exchange, primarily gas diffusion and cave ventilation due to advective forces, were characterized. The spatial distributions of both processes were described through measurements of CO2 and its carbon isotopic signal (δ13C[CO2]) from exterior, soil and cave air samples analyzed by cavity ring-down spectroscopy (CRDS). The trigger mechanisms of air advection (temperature or air density differences or barometric imbalances) were controlled by continuous logging systems. Radon monitoring was also used to characterize the changing airflow that results in a predictable seasonal or daily pattern of CO2 concentrations and its carbon isotopic signal. Large daily oscillations of CO2 levels, ranging from 680 to 1900 ppm day−1 on average, were registered during the daily oscillations of the exterior air temperature around the cave air temperature. These daily variations in CO2 concentration were unobservable once the outside air temperature was continuously below the cave temperature and a prevailing advective-renewal of cave air was established, such that the daily-averaged concentrations of CO2 reached minimum values close to atmospheric background. The daily pulses of CO2 and other tracer gases such as radon (222Rn) were smoothed in the inner cave locations, where fluctuation of both gases was primarily correlated with medium-term changes in air pressure. A pooled analysis of these data provided evidence that atmospheric air that is inhaled into dynamically ventilated caves can then return to the lower troposphere as CO2-rich cave air.

CO2, CH4 and N2O fluxes in an Ultisol treated with sewage sludge and cultivated with castor bean

Organic residue application into soil alter the emission of gases to atmosphere and CO2, CH4, N2O may contribute to increase the greenhouse effect. This experiment was carried out in a restoration area on a dystrophic Ultisol (PVAd) to quantify greenhouse gas (GHG) emissions from soil under castor bean cultivation, treated with sewage sludge (SS) or mineral fertilizer. The following treatments were tested: control without N; FertMin = mineral fertilizer; SS5 = 5 t ha-1 SS (37.5 kg ha-1 N); SS10 = 10 t ha-1 SS (75 kg ha-1 N); and SS20 = 20 t ha-1 SS (150 kg ha-1 N). The amount of sludge was based on the recommended rate of N for castor bean (75 kg ha-1), the N level of SS and the mineralization fraction of N from SS. Soil gas emission was measured for 21 days. Sewage sludge and mineral fertilizers altered the CO2, CH4 and N2O fluxes. Soil moisture had no effect on GHG emissions and the gas fluxes was statistically equivalent after the application of FertMin and of 5 t ha-1 SS. The application of the entire crop N requirement in the form of SS practically doubled the Global Warming Potential (GWP) and the C equivalent emissions in comparison with FertMin treatments.

Improved understanding of drought controls on seasonal variation in Mediterranean forest canopy CO2 and water fluxes through combined in situ measurements and ecosystem modelling

Water stress is a defining characteristic of Mediterranean ecosystems, and is likely to become more severe in the coming decades. Simulation models are key tools for making predictions, but our current understanding of how soil moisture controls ecosystem functioning is not sufficient to adequately constrain parameterisations. Canopy-scale flux data from four forest ecosystems with Mediterranean-type climates were used in order to analyse the physiological controls on carbon and water flues through the year. Significant non-stomatal limitations on photosynthesis were detected, along with lesser changes in the conductance-assimilation relationship. New model parameterisations were derived and implemented in two contrasting modelling approaches. The effectiveness of two models, one a dynamic global vegetation model ('ORCHIDEE'), and the other a forest growth model particularly developed for Mediterranean simulations ('GOTILWA+'), was assessed and modelled canopy responses to seasonal changes in soil moisture were analysed in comparison with in situ flux measurements. In contrast to commonly held assumptions, we find that changing the ratio of conductance to assimilation under natural, seasonally-developing, soil moisture stress is not sufficient to reproduce forest canopy CO2 and water fluxes. However, accurate predictions of both CO2 and water fluxes under all soil moisture levels encountered in the field are obtained if photosynthetic capacity is assumed to vary with soil moisture. This new parameterisation has important consequences for simulated responses of carbon and water fluxes to seasonal soil moisture stress, and should greatly improve our ability to anticipate future impacts of climate changes on the functioning of ecosystems in Mediterranean-type climates.

Simulated and observed fluxes of sensible and latent heat and CO2 at the WLEF-TV tower using SiB2.5

Three years of meteorological data collected at the WLEF-TV tower were used to drive a revised version of the Simple Biosphere (SiB 2.5) Model. Physiological properties and vegetation phenology were specified from satellite imagery. Simulated fluxes of heat, moisture, and carbon were compared to eddy covariance measurements taken onsite as a means of evaluating model performance on diurnal, synoptic, seasonal, and interannual time scales. The model was very successful in simulating variations of latent heat flux when compared to observations, slightly less so in the simulation of sensible heat flux. The model overestimated peak values of sensible heat flux on both monthly and diurnal scales. There was evidence that the differences between observed and simulated fluxes might be linked to wetlands near the WLEF tower, which were not present in the SiB simulation. The model overestimated the magnitude of the net ecosystem exchange of CO2 in both summer and winter. Mid-day maximum assimilation was well represented by the model, but late afternoon simulations showed excessive carbon uptake due to misrepresentation of within-canopy shading in the model. Interannual variability was not well simulated because only a single year of satellite imagery was used to parameterize the model.

Identification of micro-scale anthropogenic CO2, heat and moisture sources – processing eddy covariance fluxes for a dense urban environment

Anthropogenic emissions of heat and exhaust gases play an important role in the atmospheric boundary layer, altering air quality, greenhouse gas concentrations and the transport of heat and moisture at various scales. This is particularly evident in urban areas where emission sources are integrated in the highly heterogeneous urban canopy layer and directly linked to human activities which exhibit significant temporal variability. It is common practice to use eddy covariance observations to estimate turbulent surface fluxes of latent heat, sensible heat and carbon dioxide, which can be attributed to a local scale source area. This study provides a method to assess the influence of micro-scale anthropogenic emissions on heat, moisture and carbon dioxide exchange in a highly urbanized environment for two sites in central London, UK. A new algorithm for the Identification of Micro-scale Anthropogenic Sources (IMAS) is presented, with two aims. Firstly, IMAS filters out the influence of micro-scale emissions and allows for the analysis of the turbulent fluxes representative of the local scale source area. Secondly, it is used to give a first order estimate of anthropogenic heat flux and carbon dioxide flux representative of the building scale. The algorithm is evaluated using directional and temporal analysis. The algorithm is then used at a second site which was not incorporated in its development. The spatial and temporal local scale patterns, as well as micro-scale fluxes, appear physically reasonable and can be incorporated in the analysis of long-term eddy covariance measurements at the sites in central London. In addition to the new IMAS-technique, further steps in quality control and quality assurance used for the flux processing are presented. The methods and results have implications for urban flux measurements in dense urbanised settings with significant sources of heat and greenhouse gases.

Integrated Data Sets of the EU FP5 Research Project ORFOIS: Origin and fate of biogenic particle fluxes in the ocean and their interactions with atmospheric CO2 concentrations as well as the marine sediment

For a reliable simulation of the time and space dependent CO2 redistribution between ocean and atmosphere an appropriate time dependent simulation of particle dynamics processes is essential but has not been carried out so far. The major difficulties were the lack of suitable modules for particle dynamics and early diagenesis (in order to close the carbon and nutrient budget) in ocean general circulation models, and the lack of an understanding of biogeochemical processes, such as the partial dissolution of calcareous particles in oversaturated water. The main target of ORFOIS was to fill in this gap in our knowledge and prediction capability infrastructure. This goal has been achieved step by step. At first comprehensive data bases (already existing data) of observations of relevance for the three major types of biogenic particles, organic carbon (POC), calcium carbonate (CaCO3), and biogenic silica (BSi or opal), as well as for refractory particles of terrestrial origin were collated and made publicly available.

Impact of elevated levels of CO2 on animal mediated ecosystem function: The modification of sediment nutrient fluxes by burrowing urchins

A mesocosm experiment was conducted to quantify the relationships between the presence and body size of two burrowing heart urchins (Brissopsis lyrifera and Echinocardium cordatum) and rates of sediment nutrient flux. Furthermore, the impact of seawater acidification on these relationships was determined during this 40-day exposure experiment. Using carbon dioxide (CO2) gas, seawater was acidified to pHNBS 7.6, 7.2 or 6.8. Control treatments were maintained in natural seawater (pH = 8.0). Under normocapnic conditions, burrowing urchins were seen to reduce the sediment uptake of nitrite or nitrate whilst enhancing the release of silicate and phosphate. In acidified (hypercapnic) treatments, the biological control of biogeochemical cycles by urchins was significantly affected, probably through the combined impacts of high CO2 on nitrifying bacteria, benthic algae and urchin behaviour. This study highlights the importance of considering biological interactions when predicting the consequences of seawater acidification on ecosystem function.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Nitrous oxide and methane fluxes in south Brazilian gleysol as affected by nitrogen fertilizers

Nitrogen fertilizers increase the nitrous oxide (N2O) emission and can reduce the methane (CH4) oxidation from agricultural soils. However, the magnitude of this effect is unknown in Southern Brazilian edaphoclimatic conditions, as well as the potential of different sources of mineral N fertilizers in such an effect. The aim of this study was to investigate the effects of different mineral N sources (urea, ammonium sulphate, calcium nitrate, ammonium nitrate, Uran, controlled- release N fertilizer, and urea with urease inhibitor) on N2O and CH4 fluxes from Gleysol in the South of Brazil (Porto Alegre, RS), in comparison to a control treatment without a N application. The experiment was arranged in a randomized block with three replications, and the N fertilizer was applied to corn at the V5 growth stage. Air samples were collected from a static chambers for 15 days after the N application and the N2O and CH4 concentration were determined by gas chromatography. The topmost emissions occurred three days after the N fertilizer application and ranged from 187.8 to 8587.4 µg m-2 h-1 N. The greatest emissions were observed for N-nitric based fertilizers, while N sources with a urease inhibitor and controlled release N presented the smallest values and the N-ammonium and amidic were intermediate. This peak of N2O emissions was related to soil NO3--N (R² = 0.56, p < 0.08) when the soil water-filled pore space was up to 70 % and it indicated that N2O was predominantly produced by a denitrification process in the soil. Soil CH4 fluxes ranged from -30.1 µg m-2 h-1 C (absorption) to +32.5 µg m-2 h-1 C (emission), and the accumulated emission in the period was related to the soil NH4+-N concentration (R² = 0.82, p < 0.001), probably due to enzymatic competition between nitrification and metanotrophy processes. Despite both of the gas fluxes being affected by N fertilizers, in the average of the treatments, the impact on CH4 emission (0.2 kg ha-1 equivalent CO2-C ) was a hundredfold minor than for N2O (132.8 kg ha-1 equivalent CO2-C). Accounting for the N2O and CH4 emissions plus energetic costs of N fertilizers of 1.3 kg CO2-C kg-1 N regarding the manufacture, transport and application, we estimated an environmental impact of N sources ranging from 220.4 to 664.5 kg ha-1 CO2 -C , which can only be partially offset by C sequestration in the soil, as no study in South Brazil reported an annual net soil C accumulation rate larger than 160 kg ha-1 C due to N fertilization. The N2O mitigation can be obtained by the replacement of N-nitric sources by ammonium and amidic fertilizers. Controlled release N fertilizers and urea with urease inhibitor are also potential alternatives to N2O emission mitigation to atmospheric and systematic studies are necessary to quantify their potential in Brazilian agroecosystems.

Measurements of glial metabolic fluxes with C-11-acetate using positron emission and an adapted NMR-based metabolic modeling approach

Objectives: Acetate brain metabolism has the particularity to occur specifically in glial cells. Labeling studies, using acetate labeled either with 13C (NMR) or 11C (PET), are governed by the same biochemical reactions and thus follow the same mathematical principles. In this study, the objective was to adapt an NMR acetate brain metabolism model to analyse [1-11C]acetate infusion in rats. Methods: Brain acetate infusion experiments were modeled using a two-compartment model approach used in NMR.1-3 The [1-11C]acetate labeling study was done using a beta scintillator.4 The measured radioactive signal represents the time evolution of the sum of all labeled metabolites in the brain. Using a coincidence counter in parallel, an arterial input curve was measured. The 11C at position C-1 of acetate is metabolized in the first turn of the TCA cycle to the position 5 of glutamate (Figure 1A). Through the neurotransmission process, it is further transported to the position 5 of glutamine and the position 5 of neuronal glutamate. After the second turn of the TCA cycle, tracer from [1-11C]acetate (and also a part from glial [5-11C]glutamate) is transferred to glial [1-11C]glutamate and further to [1-11C]glutamine and neuronal glutamate through the neurotransmission cycle. Brain poster session: oxidative mechanisms S460 Journal of Cerebral Blood Flow & Metabolism (2009) 29, S455-S466 Results: The standard acetate two-pool PET model describes the system by a plasma pool and a tissue pool linked by rate constants. Experimental data are not fully described with only one tissue compartment (Figure 1B). The modified NMR model was fitted successfully to tissue time-activity curves from 6 single animals, by varying the glial mitochondrial fluxes and the neurotransmission flux Vnt. A glial composite rate constant Kgtg=Vgtg/[Ace]plasma was extracted. Considering an average acetate concentration in plasma of 1 mmol/g5 and the negligible additional amount injected, we found an average Vgtg = 0.08±0.02 (n = 6), in agreement with previous NMR measurements.1 The tissue time-activity curve is dominated by glial glutamate and later by glutamine (Figure 1B). Labeling of neuronal pools has a low influence, at least for the 20 mins of beta-probe acquisition. Based on the high diffusivity of CO2 across the blood-brain barrier; 11CO2 is not predominant in the total tissue curve, even if the brain CO2 pool is big compared with other metabolites, due to its strong dilution through unlabeled CO2 from neuronal metabolism and diffusion from plasma. Conclusion: The two-compartment model presented here is also able to fit data of positron emission experiments and to extract specific glial metabolic fluxes. 11C-labeled acetate presents an alternative for faster measurements of glial oxidative metabolism compared to NMR, potentially applicable to human PET imaging. However, to quantify the relative value of the TCA cycle flux compared to the transmitochondrial flux, the chemical sensitivity of NMR is required. PET and NMR are thus complementary.