• European Marine Science
• Other research productsclose
• Open Access close
• European Commission close
• EC|FP7 close
• CH close
• JP close

Changes in marine net primary productivity (PP) and export of particulate organic carbon (EP) are projected over the 21st century with four global coupled carbon cycle-climate models. These include representations of marine ecosystems and the carbon cycle of different structure and complexity. All four models show a decrease in global mean PP and EP between 2 and 20% by 2100 relative to preindustrial conditions, for the SRES A2 emission scenario. Two different regimes for productivity changes are consistently identified in all models. The first chain of mechanisms is dominant in the low- and mid-latitude ocean and in the North Atlantic: reduced input of macro-nutrients into the euphotic zone related to enhanced stratification, reduced mixed layer depth, and slowed circulation causes a decrease in macro-nutrient concentrations and in PP and EP. The second regime is projected for parts of the Southern Ocean: an alleviation of light and/or temperature limitation leads to an increase in PP and EP as productivity is fueled by a sustained nutrient input. A region of disagreement among the models is the Arctic, where three models project an increase in PP while one model projects a decrease. Projected changes in seasonal and interannual variability are modest in most regions. Regional model skill metrics are proposed to generate multi-model mean fields that show an improved skill in representing observation-based estimates compared to a simple multi-model average. Model results are compared to recent productivity projections with three different algorithms, usually applied to infer net primary production from satellite observations.

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

This study presents simulations of Greenland surface melt for the Eemian interglacial period (∼130 000 to 115 000 years ago) derived from regional climate simulations with a coupled surface energy balance model. Surface melt is of high relevance due to its potential effect on ice core observations, e.g., lowering the preserved total air content (TAC) used to infer past surface elevation. An investigation of surface melt is particularly interesting for warm periods with high surface melt, such as the Eemian interglacial period. Furthermore, Eemian ice is the deepest and most compressed ice preserved on Greenland, resulting in our inability to identify melt layers visually. Therefore, simulating Eemian melt rates and associated melt layers is beneficial to improve the reconstruction of past surface elevation. Estimated TAC, based on simulated melt during the Eemian, could explain the lower TAC observations. The simulations show Eemian surface melt at all deep Greenland ice core locations and an average of up to ∼30 melt days per year at Dye-3, corresponding to more than 600 mm water equivalent (w.e.) of annual melt. For higher ice sheet locations, between 60 and 150 mmw.e.yr-1 on average are simulated. At the summit of Greenland, this yields a refreezing ratio of more than 25 % of the annual accumulation. As a consequence, high melt rates during warm periods should be considered when interpreting Greenland TAC fluctuations as surface elevation changes. In addition to estimating the influence of melt on past TAC in ice cores, the simulated surface melt could potentially be used to identify coring locations where Greenland ice is best preserved.

The stable carbon isotope ratio of atmospheric CO2 (d13Catm) is a key parameter in deciphering past carbon cycle changes. Here we present d13Catm data for the past 24,000 years derived from three independent records from two Antarctic ice cores. We conclude that a pronounced 0.3 per mil decrease in d13Catm during the early deglaciation can be best explained by upwelling of old, carbon-enriched waters in the Southern Ocean. Later in the deglaciation, regrowth of the terrestrial biosphere, changes in sea surface temperature, and ocean circulation governed the d13Catm evolution. During the Last Glacial Maximum, d13Catm and atmospheric CO2 concentration were essentially constant, which suggests that the carbon cycle was in dynamic equilibrium and that the net transfer of carbon to the deep ocean had occurred before then.

We quantify the CO2 source/sink nature of the California Current System (CalCS) and determine the drivers and processes behind the mean and spatiotemporal variability of the partial pressure of CO2 (pCO2) in the surface ocean. To this end, we analyze eddy-resolving, climatological simulations of a coupled physical–biogeochemical oceanic model on the basis of the Regional Oceanic Modeling System (ROMS). In the annual mean, the entire CalCS within 800 km of the coast and from ∼33° N to 46° N is essentially neutral with regard to atmospheric CO2: the model simulates an integrated uptake flux of −0.9 ± 3.6 Tg C yr−1, corresponding to an average flux density of −0.05 ± 0.20 mol C m−2 yr−1. This near zero flux is a consequence of an almost complete regional compensation between (i) strong outgassing in the nearshore region (first 100 km) that brings waters with high concentrations of dissolved inorganic carbon (DIC) to the surface and (ii) and a weaker, but more widespread uptake flux in the offshore region due to an intense biological reduction of this DIC, driven by the nutrients that are upwelled together with the DIC. The air–sea CO2 fluxes vary substantially in time, both on seasonal and sub-seasonal timescales, largely driven by variations in surface ocean pCO2. Most of the variability in pCO2 is associated with the seasonal cycle, with the exception of the nearshore region, where sub-seasonal variations driven by mesoscale processes dominate. In the regions offshore of 100 km, changes in surface temperature are the main driver, while in the nearshore region, changes in surface temperature, as well as anomalies in DIC and alkalinity (Alk) owing to changes in circulation, biological productivity and air–sea CO2 fluxes dominate. The prevalence of eddy-driven variability in the nearshore 100 km leads to a complex spatiotemporal mosaic of surface ocean pCO2 and air–sea CO2 fluxes that require a substantial observational effort to determine the source/sink nature of this region reliably.

Understanding the role of atmospheric CO2 during past climate changes requires clear knowledge of how it varies in time relative to temperature. Antarctic ice cores preserve highly resolved records of atmospheric CO2 and Antarctic temperature for the past 800,000 years. Here we propose a revised relative age scale for the concentration of atmospheric CO2 and Antarctic temperature for the last deglacial warming, using data from five Antarctic ice cores. We infer the phasing between CO2 concentration and Antarctic temperature at four times when their trends change abruptly. We find no significant asynchrony between them, indicating that Antarctic temperature did not begin to rise hundreds of years before the concentration of atmospheric CO2, as has been suggested by earlier studies.

In this paper we provide an overview of new knowledge on oxygen depletion (hypoxia) and related phenomena in aquatic systems resulting from the EU-FP7 project HYPOX ("In situ monitoring of oxygen depletion in hypoxic ecosystems of coastal and open seas, and landlocked water bodies", http://www.hypox.net). In view of the anticipated oxygen loss in aquatic systems due to eutrophication and climate change, HYPOX was set up to improve capacities to monitor hypoxia as well as to understand its causes and consequences. Temporal dynamics and spatial patterns of hypoxia were analyzed in field studies in various aquatic environments, including the Baltic Sea, the Black Sea, Scottish and Scandinavian fjords, Ionian Sea lagoons and embayments, and Swiss lakes. Examples of episodic and rapid (hours) occurrences of hypoxia, as well as seasonal changes in bottom-water oxygenation in stratified systems, are discussed. Geologically driven hypoxia caused by gas seepage is demonstrated. Using novel technologies, temporal and spatial patterns of water-column oxygenation, from basin-scale seasonal patterns to meter-scale sub-micromolar oxygen distributions, were resolved. Existing multidecadal monitoring data were used to demonstrate the imprint of climate change and eutrophication on long-term oxygen distributions. Organic and inorganic proxies were used to extend investigations on past oxygen conditions to centennial and even longer timescales that cannot be resolved by monitoring. The effects of hypoxia on faunal communities and biogeochemical processes were also addressed in the project. An investigation of benthic fauna is presented as an example of hypoxia-devastated benthic communities that slowly recover upon a reduction in eutrophication in a system where naturally occurring hypoxia overlaps with anthropogenic hypoxia. Biogeochemical investigations reveal that oxygen intrusions have a strong effect on the microbially mediated redox cycling of elements. Observations and modeling studies of the sediments demonstrate the effect of seasonally changing oxygen conditions on benthic mineralization pathways and fluxes. Data quality and access are crucial in hypoxia research. Technical issues are therefore also addressed, including the availability of suitable sensor technology to resolve the gradual changes in bottom-water oxygen in marine systems that can be expected as a result of climate change. Using cabled observatories as examples, we show how the benefit of continuous oxygen monitoring can be maximized by adopting proper quality control. Finally, we discuss strategies for state-of-the-art data archiving and dissemination in compliance with global standards, and how ocean observations can contribute to global earth observation attempts.

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Computer models are necessary for understanding and predicting marine ice sheet behaviour. However, there is uncertainty over implementation of physical processes at the ice base, both for grounded and floating glacial ice. Here we implement several sliding relations in a marine ice sheet flow-line model accounting for all stress components and demonstrate that model resolution requirements are strongly dependent on both the choice of basal sliding relation and the spatial distribution of ice shelf basal melting.Sliding relations that reduce the magnitude of the step change in basal drag from grounded ice to floating ice (where basal drag is set to zero) show reduced dependence on resolution compared to a commonly used relation, in which basal drag is purely a power law function of basal ice velocity. Sliding relations in which basal drag goes smoothly to zero as the grounding line is approached from inland (due to a physically motivated incorporation of effective pressure at the bed) provide further reduction in resolution dependence.A similar issue is found with the imposition of basal melt under the floating part of the ice shelf: melt parameterisations that reduce the abruptness of change in basal melting from grounded ice (where basal melt is set to zero) to floating ice provide improved convergence with resolution compared to parameterisations in which high melt occurs adjacent to the grounding line.Thus physical processes, such as sub-glacial outflow (which could cause high melt near the grounding line), impact on capability to simulate marine ice sheets. If there exists an abrupt change across the grounding line in either basal drag or basal melting, then high resolution will be required to solve the problem. However, the plausible combination of a physical dependency of basal drag on effective pressure, and the possibility of low ice shelf basal melt rates next to the grounding line, may mean that some marine ice sheet systems can be reliably simulated at a coarser resolution than currently thought necessary.

The impact of viral lysis and grazing by flagellates on bacterioplankton production was assessed during a mesocosm experiment in the Eastern Mediterranean Sea, in response to Saharan dust (SD) vs. mixed aerosols (A) addition. The results highlight a positive effect on bacterial abundance, production and growth rate (~1.2, ~2.4, and ~1.9-fold higher than the controls) in both SD and A, which was also confirmed by the increased portion of high DNA content bacteria (up to 48% of the bacterial community). Lytic viral production and the portion of bacterial production lost due to viral lysis were lower in SD and A after dust addition than in the controls (0.33 ± 0.17 × 106 virus-like particles mL-1 h-1 and 6 ± 4%, respectively). Potential ingestion rate of bacteria by flagellates increased upon dust enrichment, but did not differ between mesocosms. Larger predators possibly down regulated flagellate abundance, and the calculated portion of bacterial production lost due to flagellate grazing was probably an artifact. Higher frequency of lysogenic cells in A compared to SD and the controls four days after dust addition may reflect faster phosphorus limitation in A, due to receiving less dissolved inorganic phosphorus and more dissolved inorganic nitrogen than SD. Science Citation Index Expanded WOS: 000457690600057