ALBERT

Description: The global sulphur cycle has been simulated using a general circulation model with a focus on the source and oxidation of atmospheric dimethylsulphide (DMS). The sensitivity of atmospheric DMS to the oceanic DMS climatology, the parameterisation of the sea-air transfer and to the oxidant fields have been studied. The importance of additional oxidation pathways (by O3 in the gas- and aqueous-phases and by BrO in the gas phase) not incorporated in global models has also been evaluated. While three different climatologies of the oceanic DMS concentration produce rather similar global DMS fluxes to the atmosphere at 24-27 Tg S yr -1, there are large differences in the spatial and seasonal distribution. The relative contributions of OH and NO3 radicals to DMS oxidation depends critically on which oxidant fields are prescribed in the model. Oxidation by O3 appears to be significant at high latitudes in both hemispheres. Oxidation by BrO could be significant even for BrO concentrations at sub-pptv levels in the marine boundary layer. The impact of such refinements on the DMS chemistry onto the indirect radiative forcing by anthropogenic sulphate aerosols is also discussed.

Description: A watermass-based framework is presented for a quantitative understanding of the processes controlling the cycling of carbon in the Southern Ocean. The approach is developed using a model simulation of the global carbon transports within the ocean and with the atmosphere. It is shown how the watermass framework sheds light on the interplay between biology, air-sea gas exchange, and internal ocean transport including diapycnal processes, and the way in which this interplay controls the large-scale ocean-atmosphere carbon exchange. The simulated pre-industrial regional patterns of DIC distribution and the global distribution of the pre-industrial air-sea CO2 fluxes compare well with other model results and with results from an ocean inversion method. The main differences are found in the Southern Ocean where the model presents a stronger CO2 outgassing south of the polar front, a result of the upwelling of DIC-rich deep waters into the surface layer. North of the subantarctic front the typical temperature-driven solubility effect produces a net ingassing of CO2. The biological controls on surface CO2 fluxes through primary production is generally smaller than the temperature effect on solubility. Novel to this study is also a Lagrangian trajectory analysis of the meridional transport of DIC. The analysis allows to evaluate the contribution of separate branches of the global thermohaline circulation (identified by watermasses) to the vertical distribution of DIC throughout the Southern Ocean and towards the global ocean. The most important new result is that the overturning associated with Subantarctic Mode Waters sustains a northward net transport of DIC (15.7×107 mol/s across 30° S). This new finding, which has also relevant implications on the prediction of anthropogenic carbon redistribution, results from the specific mechanism of SAMW formation and its source waters whose consequences on tracer transports are analyzed for the first time in this study.

Description: The El Niño Southern Oscillation (ENSO) drives important changes in the marine productivity of the Equatorial Pacific, in particular during major El Niño/La Niña transitions. Changes in environmental conditions associated with these climatic events also likely impact phytoplankton composition. In this work, the distribution of four major phytoplankton groups (nanoeucaryotes, Prochlorococcus, Synechococcus, and diatoms) was examined between 1996 and 2007 by applying the PHYSAT algorithm to the ocean color data archive from the Ocean Color and Temperature Sensor (OCTS) and Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Coincident with the decrease in chlorophyll concentrations, a large-scale shift in the phytoplankton composition of the Equatorial Pacific, that was characterized by a decrease in Synechococcus and an increase in nanoeucaryotes dominance, was observed during the early stages of both the strong El Niño of 1997 and the moderate El Niño of 2006. A significant increase in diatoms dominance was observed in the Equatorial Pacific during the 1998 La Niña and was associated with elevated marine productivity. An analysis of the environmental variables using a coupled physical-biogeochemical model (NEMO-PISCES) suggests that the Synechococcus dominance decrease during the two El Niño events was associated with an abrupt decline in nutrient availability (−0.9 to −2.5 μM NO3 month−1). Alternatively, increased nutrient availability (3 μM NO3 month−1) during the 1998 La Niña resulted in Equatorial Pacific dominance diatom increase. Despite these phytoplankton community shifts, the mean composition is restored after a few months, which suggests resilience in community structure. Such rapid changes to the composition of phytoplankton groups should be considered in future modeling approaches to represent variability of the marine productivity in the Equatorial Pacific and to quantify its potential implications on food-web and on global carbon cycle.

Description: The land and ocean absorb on average over half of the anthropogenic emissions of carbon dioxide (CO2) every year. These CO2 "sinks" are modulated by climate change and variability. Here we use a suite of nine Dynamic Global Vegetation Models (DGVMs) and four Ocean Biogeochemical General Circulation Models (OBGCMs) to quantify the global and regional climate and atmospheric CO2 – driven trends in land and oceanic CO2 exchanges with the atmosphere over the period 1990–2009, attribute these trends to underlying processes, and quantify the uncertainty and level of model agreement. The models were forced with reconstructed climate fields and observed global atmospheric CO2; Land Use and Land Cover Changes are not included for the DGVMs. Over the period 1990–2009, the DGVMs simulate a mean global land carbon sink of −2.4 ± 0.7 Pg C yr−1 with a small significant trend of −0.06 ± 0.03 Pg C yr−2 (increasing sink). Over the more limited period 1990–2004, the ocean models simulate a mean ocean sink of –2.2 ± 0.2 Pg C yr–1 with a trend in the net C uptake that is indistinguishable from zero (−0.01 ± 0.02 Pg C yr−2). The two ocean models that extended the simulations until 2009 suggest a slightly stronger, but still small trend of −0.02 ± 0.01 Pg C yr−2. Trends from land and ocean models compare favourably to the land greenness trends from remote sensing, atmospheric inversion results, and the residual land sink required to close the global carbon budget. Trends in the land sink are driven by increasing net primary production (NPP) whose statistically significant trend of 0.22 ± 0.08 Pg C yr−2 exceeds a significant trend in heterotrophic respiration of 0.16 ± 0.05 Pg C yr−2 – primarily as a consequence of wide-spread CO2 fertilisation of plant production. Most of the land-based trend in simulated net carbon uptake originates from natural ecosystems in the tropics (−0.04 ± 0.01 Pg C yr−2), with almost no trend over the northern land region, where recent warming and reduced rainfall offsets the positive impact of elevated atmospheric CO2 on carbon storage. The small uptake trend in the ocean models emerges because climate variability and change, and in particular increasing sea surface temperatures, tend to counteract the trend in ocean uptake driven by the increase in atmospheric CO2. Large uncertainty remains in the magnitude and sign of modelled carbon trends in several regions, and on the influence of land use and land cover changes on regional trends.

Description: Planktonic foraminifera are a major contributor to the deep carbonate-flux and the planktonic biomass of the global ocean. Their microfossil deposits form one of the richest databases for reconstructing paleoenvironments, particularly through changes in their taxonomic and shell composition. Using an empirically-based foraminifer model that incorporates three known major physiological drivers of foraminifer biogeography – temperature, food and light – we investigate (i) the global redistribution of planktonic foraminifera under anthropogenic climate change, and (ii) the alteration of the carbonate chemistry of foraminifer habitat with ocean acidification. The present-day and future (2090–2100) 3-D distributions of foraminifera are simulated using temperature, plankton biomass, and light from an Earth system model forced with historical and a future (IPCC A2) high CO2 emission scenario. The broadscale patterns of present day foraminifer biogeography are well reproduced. Foraminifer abundance and diversity are projected to decrease in the tropics and subpolar regions and increase in the subtropics and around the poles. In the tropics, the geographical shifts are driven by temperature, while the vertical shifts are driven by both temperature and food availability. In the high-latitudes, vertical shifts are driven by food availability, while geographical shifts are driven by both food availability and temperature. Changes in the marine carbon cycle would be expected in response to (i) the large-scale rearrangements in foraminifer abundance, and (ii) the reduction of the carbonate concentration in the habitat range of planktonic foraminifers: from 10–30 μmol kg−1 in the polar/subpolar regions to 30–70 μmol kg−1 in the subtropical/tropical regions. High-latitude species are most vulnerable to anthropogenic change: their abundance and available habitat decrease and up to 10% of their habitat drops below the calcite saturation horizon.

Description: The observed seasonal cycles in atmospheric potential oxygen (APO) at a range of mid to high latitude surface monitoring sites are compared to those inferred from the output of 6 Earth System Models participating in the fifth phase of the Coupled Model Intercomparison Project (CMIP5). The simulated air–sea O2 fluxes are translated into APO seasonal cycles using a matrix method that takes into account atmospheric transport model (ATM) uncertainty among 13 different ATMs. Half of the ocean biogeochemistry models tested are able to reproduce the observed APO cycles at most sites, to within the current large ATM uncertainty, while the other half generally are not. Net Primary Production (NPP) and net community production (NCP), as estimated from satellite ocean color data, provide additional constraints, albeit more with respect to the seasonal phasing of ocean model productivity than the overall magnitude. The present analysis suggests that, of the tested ocean biogeochemistry models, CESM and GFDL ESM2M are best able to capture the observed APO seasonal cycle at both Northern and Southern Hemisphere sites. In the northern oceans, the comparison to observed APO suggests that most models tend to underestimate NPP or deep ventilation or both.

Description: The ocean is a substantial source of nitrous oxide (N2O) to the atmosphere, but little is known on how this flux might change in the future. Here, we investigate the potential evolution of marine N2O emissions in the 21st century in response to anthropogenic climate change using the global ocean biogeochemical model NEMO-PISCES. We implemented two different parameterizations of N2O production, which differ primarily at low oxygen (O2) conditions. When forced with output from a climate model simulation run under the business-as-usual high CO2 concentration scenario (RCP8.5), our simulations suggest a decrease of 4 to 12% in N2O emissions from 2005 to 2100, i.e., a reduction from 4.03/3.71 to 3.54/3.56 Tg N yr−1 depending on the parameterization. The emissions decrease strongly in the western basins of the Pacific and Atlantic oceans, while they tend to increase above the Oxygen Minimum Zones (OMZs), i.e., in the Eastern Tropical Pacific and in the northern Indian Ocean. The reduction in N2O emissions is caused on the one hand by weakened nitrification as a consequence of reduced primary and export production, and on the other hand by stronger vertical stratification, which reduces the transport of N2O from the ocean interior to the ocean surface. The higher emissions over the OMZ are linked to an expansion of these zones under global warming, which leads to increased N2O production associated primarily with denitrification. From the perspective of a global climate system, the averaged feedback strength associated with the projected decrease in oceanic N2O emissions amounts to around −0.009 W m−2 K−1, which is comparable to the potential increase from terrestrial N2O sources. However, the assesment for a compensation between the terrestrial and marine feedbacks calls for an improved representation of N2O production terms in fully coupled next generation of Earth System Models.

Description: Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO2 in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth System Models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration parthways (RCP) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in 2090s (compared to 1990s) for sea surface temperature, sea surface pH, global O2 content and integrated primary productivity amount to +2.73 °C, −0.33 pH unit, −3.45% and −8.6%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 °C, −0.07 pH unit, −1.81% and −2.0% respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns. Large decreases in O2 and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of sub-surface O2 concentrations in the tropics and global and regional changes in net primary productivity.

Description: This study focuses on an improved representation of the biological soft tissue pump in the global three-dimensional biogeochemical ocean model PISCES. We compare three parameterizations of particle dynamics: (1) the model standard version including two particle size classes, aggregation-disaggregation and prescribed sinking speed; (2) an aggregation-disaggregation model with a particle size spectrum and prognostic sinking speed; (3) a mineral ballast parameterization with no size classes, but prognostic sinking speed. In addition, the model includes a description of surface sediments and organic carbon early diagenesis. Model output is compared to data or data based estimates of ocean productivity, pe-ratios, particle fluxes, surface sediment bulk composition and benthic O2 fluxes. Model results suggest that different processes control POC fluxes at different depths. In the wind mixed layer turbulent particle coagulation appears as key process in controlling pe-ratios. Parameterization (2) yields simulated pe-ratios that compare well to observations. Below the wind mixed layer, POC fluxes are most sensitive to the intensity of zooplankton flux feeding, indicating the importance of zooplankton community composition. All model parameters being kept constant, the capability of the model to reproduce yearly mean POC fluxes below 2000 m and benthic oxygen demand does at first order not dependent on the resolution of the particle size spectrum. Aggregate formation appears essential to initiate an intense biological pump. At great depth the reported close to constant particle fluxes are most likely the result of the combined effect of aggregate formation and mineral ballasting.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr−1, ELUC 0.9 ± 0.5 GtC yr−1, GATM 4.3 ± 0.1 GtC yr−1, SOCEAN 2.5 ± 0.5 GtC yr−1, and SLAND 2.8 ± 0.8 GtC yr−1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr−1, 2.2% above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and assuming an ELUC of 1.0 ± 0.5 GtC yr−1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr−1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70% from EFF (390 ± 20 GtC) and 30% from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3).