• European Marine Science
• Other research productsclose
• 2013-2022close
• European Commission close
• EC|FP7 close
• EC|FP7|SP1 close
• CH close
• RO close
• UA close

As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated pre-study experiment with interactive stratospheric aerosol models simulating the volcanic aerosol cloud from an eruption resembling the 1815 Mt. Tambora eruption (VolMIP-Tambora ISA ensemble). The pre-study provided the ancillary ability to assess intermodel diversity in the radiative forcing for a large stratospheric-injecting equatorial eruption when the volcanic aerosol cloud is simulated interactively. An initial analysis of the VolMIP-Tambora ISA ensemble showed large disparities between models in the stratospheric global mean aerosol optical depth (AOD). In this study, we now show that stratospheric global mean AOD differences among the participating models are primarily due to differences in aerosol size, which we track here by effective radius. We identify specific physical and chemical processes that are missing in some models and/or parameterized differently between models, which are together causing the differences in effective radius. In particular, our analysis indicates that interactively tracking hydroxyl radical (OH) chemistry following a large volcanic injection of sulfur dioxide (SO2) is an important factor in allowing for the timescale for sulfate formation to be properly simulated. In addition, depending on the timescale of sulfate formation, there can be a large difference in effective radius and subsequently AOD that results from whether the SO2 is injected in a single model grid cell near the location of the volcanic eruption, or whether it is injected as a longitudinally averaged band around the Earth.

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.

Concentrations of dissolved 230Th in the ocean water column increase with depth due to scavenging and downward particle flux. Due to the 230Th scavenging process, any change in the calcium carbonate (CaCO3) fraction of the marine particle flux due to changes in biological CaCO3 hard-shell production as a consequence of progressing ocean acidification would be reflected in the dissolved 230Th activity. Our prognostic simulations with a biogeochemical ocean general circulation model using different scenarios for the reduction of CaCO3 production under ocean acidification and different greenhouse gas emission scenarios – the Representative Concentration Pathways (RCPs) 8.5 to 2.6 – reveal the potential for deep 230Th measurements to detect reduced CaCO3 production at the sea surface. The time of emergence of an acidification-induced signal on dissolved 230Th is of the same order of magnitude as for alkalinity measurements. Interannual and decadal variability in factors other than a reduction in CaCO3 hard-shell production may mask the ocean-acidification-induced signal in dissolved 230Th and make detection of the pure CaCO3-induced signal more difficult so that only really strong changes in marine CaCO3 export would be unambiguously identifiable soon. Nevertheless, the impacts of changes in CaCO3 export production on marine 230Th are stronger than those for changes in POC (particulate organic carbon) or clay fluxes.

Version 2 of the Global Ocean Data Analysis Project (GLODAPv2) data product is composed of data from 724 scientific cruises covering the global ocean. It includes data assembled during the previous efforts GLODAPv1.1 (Global Ocean Data Analysis Project version 1.1) in 2004, CARINA (CARbon IN the Atlantic) in 2009/2010, and PACIFICA (PACIFic ocean Interior CArbon) in 2013, as well as data from an additional 168 cruises. Data for 12 core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have been subjected to extensive quality control, including systematic evaluation of bias. The data are available in two formats: (i) as submitted but updated to WOCE exchange format and (ii) as a merged and internally consistent data product. In the latter, adjustments have been applied to remove significant biases, respecting occurrences of any known or likely time trends or variations. Adjustments applied by previous efforts were re-evaluated. Hence, GLODAPv2 is not a simple merging of previous products with some new data added but a unique, internally consistent data product. This compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 6 µmol kg−1 in total alkalinity, 0.005 in pH, and 5 % for the halogenated transient tracers.The original data and their documentation and doi codes are available at the Carbon Dioxide Information Analysis Center (http://cdiac.ornl.gov/oceans/GLODAPv2/). This site also provides access to the calibrated data product, which is provided as a single global file or four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under the doi:10.3334/CDIAC/OTG.NDP093_GLODAPv2. The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2 methods and products and includes a broad overview of the secondary quality control results. The magnitude of and reasoning behind each adjustment is available on a per-cruise and per-variable basis in the online Adjustment Table.

We present a mapped climatology (GLODAPv2.2016b) of ocean biogeochemical variables based on the new GLODAP version 2 data product (Olsen et al., 2016; Key et al., 2015), which covers all ocean basins over the years 1972 to 2013. The quality-controlled and internally consistent GLODAPv2 was used to create global 1° × 1° mapped climatologies of salinity, temperature, oxygen, nitrate, phosphate, silicate, total dissolved inorganic carbon (TCO2), total alkalinity (TAlk), pH, and CaCO3 saturation states using the Data-Interpolating Variational Analysis (DIVA) mapping method. Improving on maps based on an earlier but similar dataset, GLODAPv1.1, this climatology also covers the Arctic Ocean. Climatologies were created for 33 standard depth surfaces. The conceivably confounding temporal trends in TCO2 and pH due to anthropogenic influence were removed prior to mapping by normalizing these data to the year 2002 using first-order calculations of anthropogenic carbon accumulation rates. We additionally provide maps of accumulated anthropogenic carbon in the year 2002 and of preindustrial TCO2. For all parameters, all data from the full 1972–2013 period were used, including data that did not receive full secondary quality control. The GLODAPv2.2016b global 1° × 1° mapped climatologies, including error fields and ancillary information, are available at the GLODAPv2 web page at the Carbon Dioxide Information Analysis Center (CDIAC; doi:10.3334/CDIAC/OTG.NDP093_GLODAPv2).

The last interglacial period (LIG, ∼ 129–116 thousand years ago) provides the most recent case study of multimillennial polar warming above the preindustrial level and a response of the Greenland and Antarctic ice sheets to this warming, as well as a test bed for climate and ice sheet models. Past changes in Greenland ice sheet thickness and surface temperature during this period were recently derived from the North Greenland Eemian Ice Drilling (NEEM) ice core records, northwest Greenland. The NEEM paradox has emerged from an estimated large local warming above the preindustrial level (7.5 ± 1.8 °C at the deposition site 126 kyr ago without correction for any overall ice sheet altitude changes between the LIG and the preindustrial period) based on water isotopes, together with limited local ice thinning, suggesting more resilience of the real Greenland ice sheet than shown in some ice sheet models. Here, we provide an independent assessment of the average LIG Greenland surface warming using ice core air isotopic composition (δ15N) and relationships between accumulation rate and temperature. The LIG surface temperature at the upstream NEEM deposition site without ice sheet altitude correction is estimated to be warmer by +8.5 ± 2.5 °C compared to the preindustrial period. This temperature estimate is consistent with the 7.5 ± 1.8 °C warming initially determined from NEEM water isotopes but at the upper end of the preindustrial period to LIG temperature difference of +5.2 ± 2.3 °C obtained at the NGRIP (North Greenland Ice Core Project) site by the same method. Climate simulations performed with present-day ice sheet topography lead in general to a warming smaller than reconstructed, but sensitivity tests show that larger amplitudes (up to 5 °C) are produced in response to prescribed changes in sea ice extent and ice sheet topography.

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.

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.

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.

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global 3-D ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.