ALBERT

Description: Ocean warming and increased stratification of the upper ocean caused by global climate change will likely lead to declines in dissolved O2 in the ocean interior (ocean deoxygenation) with implications for ocean productivity, nutrient cycling, carbon cycling, and marine habitat. Ocean models predict declines of 1 to 7% in the global ocean O2 inventory over the next century, with declines continuing for a thousand years or more into the future. An important consequence may be an expansion in the area and volume of so-called oxygen minimum zones, where O2 levels are too low to support many macrofauna and profound changes in biogeochemical cycling occur. Significant deoxygenation has occurred over the past 50 years in the North Pacific and tropical oceans, suggesting larger changes are looming. The potential for larger O2 declines in the future suggests the need for an improved observing system for tracking ocean O2 changes.

Description: © The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 10 (2013): 607-627, doi:10.5194/bg-10-607-2013.

Description: The Atlantic and Arctic Oceans are critical components of the global carbon cycle. Here we quantify the net sea–air CO2 flux, for the first time, across different methodologies for consistent time and space scales for the Atlantic and Arctic basins. We present the long-term mean, seasonal cycle, interannual variability and trends in sea–air CO2 flux for the period 1990 to 2009, and assign an uncertainty to each. We use regional cuts from global observations and modeling products, specifically a pCO2-based CO2 flux climatology, flux estimates from the inversion of oceanic and atmospheric data, and results from six ocean biogeochemical models. Additionally, we use basin-wide flux estimates from surface ocean pCO2 observations based on two distinct methodologies. Our estimate of the contemporary sea–air flux of CO2 (sum of anthropogenic and natural components) by the Atlantic between 40° S and 79° N is −0.49 ± 0.05 Pg C yr−1, and by the Arctic it is −0.12 ± 0.06 Pg C yr−1, leading to a combined sea–air flux of −0.61 ± 0.06 Pg C yr−1 for the two decades (negative reflects ocean uptake). We do find broad agreement amongst methodologies with respect to the seasonal cycle in the subtropics of both hemispheres, but not elsewhere. Agreement with respect to detailed signals of interannual variability is poor, and correlations to the North Atlantic Oscillation are weaker in the North Atlantic and Arctic than in the equatorial region and southern subtropics. Linear trends for 1995 to 2009 indicate increased uptake and generally correspond between methodologies in the North Atlantic, but there is disagreement amongst methodologies in the equatorial region and southern subtropics.

Description: U. Schuster has been supported by EU grants IP 511176-2 (CARBOOCEAN), 212196 (COCOS), and 264879 (CARBOCHANGE), and UK NERC grant NE/H017046/1 (UKOARP). G. A. McKinley and A. Fay thank NASA for support (NNX08AR68G, NNX11AF53G). P. Landsch¨utzer has been supported by EU grant 238366 (GREENCYCLESII). N. Metzl acknowledges the French national funding program LEFE/INSU. Support for N. Gruber has been provided by EU grants 264879 (CARBOCHANGE) and 283080 (GEO-CARBON) S. Doney acknowledges support from NOAA (NOAA-NA07OAR4310098). T. Takahashi is supported by NOAA (NAO80AR4320754).

Description: © The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 10 (2013): 193-216, doi:10.5194/bg-10-193-2013.

Description: Due to seasonal upwelling, the upper ocean waters of the California Current System (CCS) have a naturally low pH and aragonite saturation state (Ωarag), making this region particularly prone to the effects of ocean acidification. Here, we use the Regional Oceanic Modeling System (ROMS) to conduct preindustrial and transient (1995–2050) simulations of ocean biogeochemistry in the CCS. The transient simulations were forced with increasing atmospheric pCO2 and increasing oceanic dissolved inorganic carbon concentrations at the lateral boundaries, as projected by the NCAR CSM 1.4 model for the IPCC SRES A2 scenario. Our results show a large seasonal variability in pH (range of ~ 0.14) and Ωarag (~ 0.2) for the nearshore areas (50 km from shore). This variability is created by the interplay of physical and biogeochemical processes. Despite this large variability, we find that present-day pH and Ωarag have already moved outside of their simulated preindustrial variability envelopes (defined by ±1 temporal standard deviation) due to the rapidly increasing concentrations of atmospheric CO2. The nearshore surface pH of the northern and central CCS are simulated to move outside of their present-day variability envelopes by the mid-2040s and late 2030s, respectively. This transition may occur even earlier for nearshore surface Ωarag, which is projected to depart from its present-day variability envelope by the early- to mid-2030s. The aragonite saturation horizon of the central CCS is projected to shoal into the upper 75 m within the next 25 yr, causing near-permanent undersaturation in subsurface waters. Due to the model's overestimation of Ωarag, this transition may occur even earlier than simulated by the model. Overall, our study shows that the CCS joins the Arctic and Southern oceans as one of only a few known ocean regions presently approaching the dual threshold of widespread and near-permanent undersaturation with respect to aragonite and a departure from its variability envelope. In these regions, organisms may be forced to rapidly adjust to conditions that are both inherently chemically challenging and also substantially different from past conditions.

Description: C. H. was supported by the European Project of Ocean Acidification (EPOCA), which received funding from the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreement no. 211384. EPOCA is endorsed by the international programs Integrated Marine Biogeochemistry and Ecosystem Research (IMBER), Land-Ocean Interactions in the Coastal Zone (LOICZ), and Surface Ocean Lower Atmosphere Study (SOLAS). C. H., M. V., Z. L., A. M. P. M. and N. G. also acknowledge support by ETH Zurich. S. D. acknowledges support from NASA-NNX11AF55G.

Description: © The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 10 (2013): 1983-2000, doi:10.5194/bg-10-1983-2013.

Description: The globally integrated sea–air anthropogenic carbon dioxide (CO2) flux from 1990 to 2009 is determined from models and data-based approaches as part of the Regional Carbon Cycle Assessment and Processes (RECCAP) project. Numerical methods include ocean inverse models, atmospheric inverse models, and ocean general circulation models with parameterized biogeochemistry (OBGCMs). The median value of different approaches shows good agreement in average uptake. The best estimate of anthropogenic CO2 uptake for the time period based on a compilation of approaches is −2.0 Pg C yr−1. The interannual variability in the sea–air flux is largely driven by large-scale climate re-organizations and is estimated at 0.2 Pg C yr−1 for the two decades with some systematic differences between approaches. The largest differences between approaches are seen in the decadal trends. The trends range from −0.13 (Pg C yr−1) decade−1 to −0.50 (Pg C yr−1) decade−1 for the two decades under investigation. The OBGCMs and the data-based sea–air CO2 flux estimates show appreciably smaller decadal trends than estimates based on changes in carbon inventory suggesting that methods capable of resolving shorter timescales are showing a slowing of the rate of ocean CO2 uptake. RECCAP model outputs for five decades show similar differences in trends between approaches.

Description: RW, G-HP., RAF were supported in part through the Global Carbon Data Management and Synthesis Project of the NOAA Climate Program Office. NG and HG were supported by funds from ETH Zurich and through the FP7 projects CarboChange (Project reference 264879) and GeoCarbon. CS was supported by grants, NSF/OPP 0944761 and NOAA NA12OAR4310058. SCD acknowledges support through the NOAA Climate Process Team activity, NOAA grant NA07OAR4310098. CH and JS were supported through EU FP7 project COMBINE (grant agreement no. 226520), the Research Council of Norway funded project CarboSeason (185105/S30), the Norwegian Metacenter for Computational Science and Storage Infrastructure (NOTUR and Norstore, “Biogeochemical Earth system modeling” projects nn2980k and ns2980k) and the core project BIOFEEDBACK of the Centre for Climate Dynamics (SKD) within the Bjerknes Centre for Climate Research.

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 12 (2015): 653-679, doi:10.5194/bg-12-653-2015.

Description: The land and ocean absorb on average just 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 estimate trends driven by global and regional climate and atmospheric CO2 in land and oceanic CO2 exchanges with the atmosphere over the period 1990–2009, to attribute these trends to underlying processes in the models, and to quantify the uncertainty and level of inter-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 widespread 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 and changes in growing season length 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 counter\-act 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, as well as regarding the influence of land use and land cover changes on regional trends.

Description: S. Sitch acknowledges financial support by RCUK through NERC (grant no. NE/J010154/). N. Gruber and C. Heinze acknowledge financial support by the European Commission through the EU FP7 projects CARBOCHANGE (grant no. 264879) and GEOCARBON (grant no. 283080). N. Gruber was additionally supported through ETH Zurich. S. C. Doney acknowledges support from the US National Science Foundation (NSF AGS-1048827). P. Friedlingstein, A. Arneth, and S. Zaehle acknowledge support by the European Commission through the EU FP7 project EMBRACE (grant no. 282672). A. Arneth and S. Sitch acknowledge the support of the European Commission-funded project LUC4C (grant no. 603542). The research leading to these results received funding from the European Community’s Seventh Framework Programme (FP7 2007–2013) under grant agreement no. 238366. A. Ahlström and B. Smith acknowledge funding through the Mistra Swedish Research Programme on Climate, Impacts and Adaptation (SWECIA). C. Heinze acknowledges support from NOTUR/NorStore projects NN2980K and NS2980K.

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 12 (2015): 5793-5809, doi:10.5194/bg-12-5793-2015.

Description: We use an eddy-resolving, regional ocean biogeochemical model to investigate the main variables and processes responsible for the climatological spatio-temporal variability of pCO2 and the air-sea CO2 fluxes in the southwestern Atlantic Ocean. Overall, the region acts as a sink of atmospheric CO2 south of 30° S, and is close to equilibrium with the atmospheric CO2 to the north. On the shelves, the ocean acts as a weak source of CO2, except for the mid/outer shelves of Patagonia, which act as sinks. In contrast, the inner shelves and the low latitude open ocean of the southwestern Atlantic represent source regions. Observed nearshore-to-offshore and meridional pCO2 gradients are well represented by our simulation. A sensitivity analysis shows the importance of the counteracting effects of temperature and dissolved inorganic carbon (DIC) in controlling the seasonal variability of pCO2. Biological production and solubility are the main processes regulating pCO2, with biological production being particularly important on the shelves. The role of mixing/stratification in modulating DIC, and therefore surface pCO2, is shown in a vertical profile at the location of the Ocean Observatories Initiative (OOI) site in the Argentine Basin (42° S, 42° W).

Description: P. H. R. Calil acknowledges support from the Brazilian agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), grants 483112/2012-7 and 307385/2013-2, and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES Process 23038.004299/2014-53). R. Arruda acknowledges support from a CAPES scholarship. S. C. Doney and I. Lima acknowledge support from the National Science Foundation (NSF AGS-1048827). N. Gruber and G. Turi received support from ETH Zurich and from the EU FP7 project CarboChange (264879).

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 12 (2015): 6955-6984, doi:10.5194/bg-12-6955-2015.

Description: Past model studies have projected a global decrease in marine net primary production (NPP) over the 21st century, but these studies focused on the multi-model mean rather than on the large inter-model differences. Here, we analyze model-simulated changes in NPP for the 21st century under IPCC's high-emission scenario RCP8.5. We use a suite of nine coupled carbon–climate Earth system models with embedded marine ecosystem models and focus on the spread between the different models and the underlying reasons. Globally, NPP decreases in five out of the nine models over the course of the 21st century, while three show no significant trend and one even simulates an increase. The largest model spread occurs in the low latitudes (between 30° S and 30° N), with individual models simulating relative changes between −25 and +40 %. Of the seven models diagnosing a net decrease in NPP in the low latitudes, only three simulate this to be a consequence of the classical interpretation, i.e., a stronger nutrient limitation due to increased stratification leading to reduced phytoplankton growth. In the other four, warming-induced increases in phytoplankton growth outbalance the stronger nutrient limitation. However, temperature-driven increases in grazing and other loss processes cause a net decrease in phytoplankton biomass and reduce NPP despite higher growth rates. One model projects a strong increase in NPP in the low latitudes, caused by an intensification of the microbial loop, while NPP in the remaining model changes by less than 0.5 %. While models consistently project increases NPP in the Southern Ocean, the regional inter-model range is also very substantial. In most models, this increase in NPP is driven by temperature, but it is also modulated by changes in light, macronutrients and iron as well as grazing. Overall, current projections of future changes in global marine NPP are subject to large uncertainties and necessitate a dedicated and sustained effort to improve the models and the concepts and data that guide their development.

Description: C. Laufkötter and the research leading to these results have received funding from the European Community’s Seventh Framework Programme (FP7 2007–2013) under grant agreements no. 238366 (Greencycles II) and 264879 (CarboChange). M. Vogt and N. Gruber acknowledge funding by ETH Zürich. S. C. Doney and I. D. Lima acknowledge support from NSF (AGS-1048827).

Description: © The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 13 (2016): 4023-4047, doi:10.5194/bg-13-4023-2016.

Description: Accurate projections of marine particle export production (EP) are crucial for predicting the response of the marine carbon cycle to climate change, yet models show a wide range in both global EP and their responses to climate change. This is, in part, due to EP being the net result of a series of processes, starting with net primary production (NPP) in the sunlit upper ocean, followed by the formation of particulate organic matter and the subsequent sinking and remineralisation of these particles, with each of these processes responding differently to changes in environmental conditions. Here, we compare future projections in EP over the 21st century, generated by four marine ecosystem models under the high emission scenario Representative Concentration Pathways (RCP) 8.5 of the Intergovernmental Panel on Climate Change (IPCC), and determine the processes driving these changes. The models simulate small to modest decreases in global EP between −1 and −12 %. Models differ greatly with regard to the drivers causing these changes. Among them, the formation of particles is the most uncertain process with models not agreeing on either magnitude or the direction of change. The removal of the sinking particles by remineralisation is simulated to increase in the low and intermediate latitudes in three models, driven by either warming-induced increases in remineralisation or slower particle sinking, and show insignificant changes in the remaining model. Changes in ecosystem structure, particularly the relative role of diatoms matters as well, as diatoms produce larger and denser particles that sink faster and are partly protected from remineralisation. Also this controlling factor is afflicted with high uncertainties, particularly since the models differ already substantially with regard to both the initial (present-day) distribution of diatoms (between 11–94 % in the Southern Ocean) and the diatom contribution to particle formation (0.6–3.8 times higher than their contribution to biomass). As a consequence, changes in diatom concentration are a strong driver for EP changes in some models but of low significance in others. Observational and experimental constraints on ecosystem structure and how the fixed carbon is routed through the ecosystem to produce export production are urgently needed in order to improve current generation ecosystem models and their ability to project future changes.

Description: The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7 2007-2013) under grant agreement no. 238366. Meike Vogt and Nicolas Gruber acknowledge funding by ETH Zürich. Judith Hauck was funded by the Helmholtz Post- Doc Programme (Initiative and Networking Fund of the Helmholtz Association). Scott C. Doney and Ivan D. Lima acknowledge the support of the National Science Foundation through the Center for Microbial Oceanography Research and Education (C-MORE), an NSF Science and Technology Center (EF-0424599).

Description: © The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 19 (2013): 4037-4054, doi:10.5194/bg-10-4037-2013.

Description: The Southern Ocean (44–75° S) plays a critical role in the global carbon cycle, yet remains one of the most poorly sampled ocean regions. Different approaches have been used to estimate sea–air CO2 fluxes in this region: synthesis of surface ocean observations, ocean biogeochemical models, and atmospheric and ocean inversions. As part of the RECCAP (REgional Carbon Cycle Assessment and Processes) project, we combine these different approaches to quantify and assess the magnitude and variability in Southern Ocean sea–air CO2 fluxes between 1990–2009. Using all models and inversions (26), the integrated median annual sea–air CO2 flux of −0.42 ± 0.07 Pg C yr−1 for the 44–75° S region, is consistent with the −0.27 ± 0.13 Pg C yr−1 calculated using surface observations. The circumpolar region south of 58° S has a small net annual flux (model and inversion median: −0.04 ± 0.07 Pg C yr−1 and observations: +0.04 ± 0.02 Pg C yr−1), with most of the net annual flux located in the 44 to 58° S circumpolar band (model and inversion median: −0.36 ± 0.09 Pg C yr−1 and observations: −0.35 ± 0.09 Pg C yr−1). Seasonally, in the 44–58° S region, the median of 5 ocean biogeochemical models captures the observed sea–air CO2 flux seasonal cycle, while the median of 11 atmospheric inversions shows little seasonal change in the net flux. South of 58° S, neither atmospheric inversions nor ocean biogeochemical models reproduce the phase and amplitude of the observed seasonal sea–air CO2 flux, particularly in the Austral Winter. Importantly, no individual atmospheric inversion or ocean biogeochemical model is capable of reproducing both the observed annual mean uptake and the observed seasonal cycle. This raises concerns about projecting future changes in Southern Ocean CO2 fluxes. The median interannual variability from atmospheric inversions and ocean biogeochemical models is substantial in the Southern Ocean; up to 25% of the annual mean flux, with 25% of this interannual variability attributed to the region south of 58° S. Resolving long-term trends is difficult due to the large interannual variability and short time frame (1990–2009) of this study; this is particularly evident from the large spread in trends from inversions and ocean biogeochemical models. Nevertheless, in the period 1990–2009 ocean biogeochemical models do show increasing oceanic uptake consistent with the expected increase of −0.05 Pg C yr−1 decade−1. In contrast, atmospheric inversions suggest little change in the strength of the CO2 sink broadly consistent with the results of Le Quéré et al. (2007).

Description: A. Lenton, B. Tilbrook, R. J. Matear and R. M. Law were funded by the Australian Climate Change Science Program and theWealth from Oceans National Research Flagship. S. C. Doney acknowledges support from the National Science Foundation (OPP-0823101), T. Takahashi is supported by grants from United States NOAA (NA08OAR4320754) and National Science Foundation (ANT 06-36879). D. Baker, N. Gruber, M. Hoppema, N. Metzl acknowledge the support of EU FP7 project CARBOCHANGE (264879). S. C. Doney acknowledges support from the National Science Foundation (OPP-0823101). N. S. Lovenduski is grateful for support from NSF (OCE-1155240) and NOAA (NA12OAR4310058). This study is also a contribution to the international IMBER/SOLAS Projects. C. Sweeney acknowledges support from the United States NOAA (NA12OAR4310058) and National Science Foundation (0944761).

Description: © Author(s) 2008. This work is distributed under the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 5 (2008): 385-406, doi:10.5194/bg-5-385-2008

Description: Using numerical simulations, we quantify the impact of changes in the ocean's biological pump on the air-sea balance of CO2 by fertilizing a small surface patch in the high-nutrient, low-chlorophyll region of the eastern tropical Pacific with iron. Decade-long fertilization experiments are conducted in a basin-scale, eddy-permitting coupled physical/biogeochemical/ecological model. In contrast to previous studies, we find that most of the dissolved inorganic carbon (DIC) removed from the euphotic zone by the enhanced biological export is replaced by uptake of CO2 from the atmosphere. Atmospheric uptake efficiencies, the ratio of the perturbation in air-sea CO2 flux to the perturbation in export flux across 100 m, integrated over 10 years, are 0.75 to 0.93 in our patch size-scale experiments. The atmospheric uptake efficiency is insensitive to the duration of the experiment. The primary factor controlling the atmospheric uptake efficiency is the vertical distribution of the enhanced biological production and export. Iron fertilization at the surface tends to induce production anomalies primarily near the surface, leading to high efficiencies. In contrast, mechanisms that induce deep production anomalies (e.g. altered light availability) tend to have a low uptake efficiency, since most of the removed DIC is replaced by lateral and vertical transport and mixing. Despite high atmospheric uptake efficiencies, patch-scale iron fertilization of the ocean's biological pump tends to remove little CO2 from the atmosphere over the decadal timescale considered here.

Description: The majority of this work was funded by the Office of Science (BER) of the US Department of Energy through Grant No. DE-FG03-00ER63010. Additional funding was provided by the Information and Technology Research section of the US National Science Foundation (NG, HF, and SD) and ETH Zurich (NG).