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: © 2009 The Authors. This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 6 (2009): 515-533, doi:10.5194/bg-6-515-2009

Description: Ocean acidification from the uptake of anthropogenic carbon is simulated for the industrial period and IPCC SRES emission scenarios A2 and B1 with a global coupled carbon cycle-climate model. Earlier studies identified seawater saturation state with respect to aragonite, a mineral phase of calcium carbonate, as a key variable governing impacts on corals and other shell-forming organisms. Globally in the A2 scenario, water saturated by more than 300%, considered suitable for coral growth, vanishes by 2070 AD (CO2≈630 ppm), and the ocean volume fraction occupied by saturated water decreases from 42% to 25% over this century. The largest simulated pH changes worldwide occur in Arctic surface waters, where hydrogen ion concentration increases by up to 185% (ΔpH=−0.45). Projected climate change amplifies the decrease in Arctic surface mean saturation and pH by more than 20%, mainly due to freshening and increased carbon uptake in response to sea ice retreat. Modeled saturation compares well with observation-based estimates along an Arctic transect and simulated changes have been corrected for remaining model-data differences in this region. Aragonite undersaturation in Arctic surface waters is projected to occur locally within a decade and to become more widespread as atmospheric CO2 continues to grow. The results imply that surface waters in the Arctic Ocean will become corrosive to aragonite, with potentially large implications for the marine ecosystem, if anthropogenic carbon emissions are not reduced and atmospheric CO2 not kept below 450 ppm.

Description: This work was funded by the European Union projects CARBOOCEAN (511176-2) and EUROCEANS (511106-2) and is a contribution to the “European Project on Ocean Acidification” (EPOCA) which received funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 211384. Additional support was received from the Swiss National Science Foundation and SCD acknowledges support from the US National Science Foundation (NSF) grant ATM-0628582.

Description: © 2009 The Authors. This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 6 (2009): 2099-2120, doi:10.5194/bg-6-2099-2009

Description: Inclusion of fundamental ecological interactions between carbon and nitrogen cycles in the land component of an atmosphere-ocean general circulation model (AOGCM) leads to decreased carbon uptake associated with CO2 fertilization, and increased carbon uptake associated with warming of the climate system. The balance of these two opposing effects is to reduce the fraction of anthropogenic CO2 predicted to be sequestered in land ecosystems. The primary mechanism responsible for increased land carbon storage under radiatively forced climate change is shown to be fertilization of plant growth by increased mineralization of nitrogen directly associated with increased decomposition of soil organic matter under a warming climate, which in this particular model results in a negative gain for the climate-carbon feedback. Estimates for the land and ocean sink fractions of recent anthropogenic emissions are individually within the range of observational estimates, but the combined land plus ocean sink fractions produce an airborne fraction which is too high compared to observations. This bias is likely due in part to an underestimation of the ocean sink fraction. Our results show a significant growth in the airborne fraction of anthropogenic CO2 emissions over the coming century, attributable in part to a steady decline in the ocean sink fraction. Comparison to experimental studies on the fate of radio-labeled nitrogen tracers in temperate forests indicates that the model representation of competition between plants and microbes for new mineral nitrogen resources is reasonable. Our results suggest a weaker dependence of net land carbon flux on soil moisture changes in tropical regions, and a stronger positive growth response to warming in those regions, than predicted by a similar AOGCM implemented without land carbon-nitrogen interactions. We expect that the between-model uncertainty in predictions of future atmospheric CO2 concentration and associated anthropogenic climate change will be reduced as additional climate models introduce carbon-nitrogen cycle interactions in their land components.

Description: This work was supported in part by NASA Earth Science Enterprise, Terrestrial Ecology Program, grant #W19,953 to P. E. Thornton. Support was provided by the National Center for Atmospheric Research (NCAR) through the NCAR Community Climate System Modeling program, and through the NCAR Biogeosciences program. Additional support was provided by the US Department of Energy, Office of Science, Office of Biological and Environmental Research. I. Fung, S. Doney, N. Mahowald, and J. Randerson acknowledge support from National Science Foundation, Atmospheric Sciences Division, through the Carbon and Water Initiative.

Description: © The Authors, 2010. This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Atmospheric Chemistry and Physics 10 (2010): 10875-10893, doi:10.5194/acp-10-10875-2010.

Description: Desert dust perturbs climate by directly and indirectly interacting with incoming solar and outgoing long wave radiation, thereby changing precipitation and temperature, in addition to modifying ocean and land biogeochemistry. While we know that desert dust is sensitive to perturbations in climate and human land use, previous studies have been unable to determine whether humans were increasing or decreasing desert dust in the global average. Here we present observational estimates of desert dust based on paleodata proxies showing a doubling of desert dust during the 20th century over much, but not all the globe. Large uncertainties remain in estimates of desert dust variability over 20th century due to limited data. Using these observational estimates of desert dust change in combination with ocean, atmosphere and land models, we calculate the net radiative effect of these observed changes (top of atmosphere) over the 20th century to be −0.14 ± 0.11 W/m2 (1990–1999 vs. 1905–1914). The estimated radiative change due to dust is especially strong between the heavily loaded 1980–1989 and the less heavily loaded 1955–1964 time periods (−0.57 ± 0.46 W/m2), which model simulations suggest may have reduced the rate of temperature increase between these time periods by 0.11 °C. Model simulations also indicate strong regional shifts in precipitation and temperature from desert dust changes, causing 6 ppm (12 PgC) reduction in model carbon uptake by the terrestrial biosphere over the 20th century. Desert dust carries iron, an important micronutrient for ocean biogeochemistry that can modulate ocean carbon storage; here we show that dust deposition trends increase ocean productivity by an estimated 6% over the 20th century, drawing down an additional 4 ppm (8 PgC) of carbon dioxide into the oceans. Thus, perturbations to desert dust over the 20th century inferred from observations are potentially important for climate and biogeochemistry, and our understanding of these changes and their impacts should continue to be refined.

Description: We would like to acknowledge NASA grants NNG06G127G and NNX07AL80G, NSF grants NSF-0832782, 0932946, 0745961 and OPP-0538427, and the UK Natural Environment Research Council.

Description: © The Authors, 2011. This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 8 (2011): 387-414, doi:10.5194/bg-8-387-2011.

Description: Coupled-carbon-climate simulations are an essential tool for predicting the impact of human activity onto the climate and biogeochemistry. Here we incorporate prognostic desert dust and anthropogenic aerosols into the CCSM3.1 coupled carbon-climate model and explore the resulting interactions with climate and biogeochemical dynamics through a series of transient anthropogenic simulations (20th and 21st centuries) and sensitivity studies. The inclusion of prognostic aerosols into this model has a small net global cooling effect on climate but does not significantly impact the globally averaged carbon cycle; we argue that this is likely to be because the CCSM3.1 model has a small climate feedback onto the carbon cycle. We propose a mechanism for including desert dust and anthropogenic aerosols into a simple carbon-climate feedback analysis to explain the results of our and previous studies. Inclusion of aerosols has statistically significant impacts on regional climate and biogeochemistry, in particular through the effects on the ocean nitrogen cycle and primary productivity of altered iron inputs from desert dust deposition.

Description: This work was done under the auspices of NASA NNG06G127G, NSF grants 0748369, 0932946, 0745961 and 0832782. The work of C. J. was supported by the Joint DECC/Defra Met Office Hadley Centre Climate Programme (GA01101).

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), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 10 (2013): 6225-6245, doi:10.5194/bg-10-6225-2013.

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 pathways (RCPs) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in the 2090s (compared to the 1990s) for sea surface temperature, sea surface pH, global O2 content and integrated primary productivity amount to +2.73 (±0.72) °C, −0.33 (±0.003) pH unit, −3.45 (±0.44)% and −8.6 (±7.9)%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 (±0.45) °C, −0.07 (±0.001) pH unit, −1.81 (±0.31)% and −2.0 (±4.1)%, respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns and thus do not change coincidentally. 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 subsurface O2 concentrations in the tropics and global and regional changes in net primary productivity. These high uncertainties in projections of primary productivity and subsurface oxygen prompt us to continue inter-model comparisons to understand these model differences, while calling for caution when using the CMIP5 models to force regional impact models.

Description: This work was supported by EU FP7 project CARBOCHANGE (under grant agreement No. 264879), EU FP7 project MEECE (under grant agreement No. 212085), EU FP7 project SOCCLI (under grant agreement No. 317699), and ANR project MACROES. S. C. Doney acknowledges the US National Science Foundation (AGS-1048827). This work has been supported by the Research Council of Norway through the EarthClim (207711/E10) and NOTUR/NorStore projects. M. Vichi acknowledges the support of the Italian Ministry of Education, University and Research and the Ministry for Environment, Land and Sea through the project GEMINA.

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): 7035-7052, doi:10.5194/bg-10-7035-2013.

Description: The Indian Ocean (44° S–30° N) plays an important role in the global carbon cycle, yet it remains one of the most poorly sampled ocean regions. Several approaches have been used to estimate net sea–air CO2 fluxes in this region: interpolated observations, ocean biogeochemical models, 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 Indian Ocean sea–air CO2 fluxes between 1990 and 2009. Using all of the models and inversions, the median annual mean sea–air CO2 uptake of −0.37 ± 0.06 PgC yr−1 is consistent with the −0.24 ± 0.12 PgC yr−1 calculated from observations. The fluxes from the southern Indian Ocean (18–44° S; −0.43 ± 0.07 PgC yr−1 are similar in magnitude to the annual uptake for the entire Indian Ocean. All models capture the observed pattern of fluxes in the Indian Ocean with the following exceptions: underestimation of upwelling fluxes in the northwestern region (off Oman and Somalia), overestimation in the northeastern region (Bay of Bengal) and underestimation of the CO2 sink in the subtropical convergence zone. These differences were mainly driven by lack of atmospheric CO2 data in atmospheric inversions, and poor simulation of monsoonal currents and freshwater discharge in ocean biogeochemical models. Overall, the models and inversions do capture the phase of the observed seasonality for the entire Indian Ocean but overestimate the magnitude. The predicted sea–air CO2 fluxes by ocean biogeochemical models (OBGMs) respond to seasonal variability with strong phase lags with reference to climatological CO2 flux, whereas the atmospheric inversions predicted an order of magnitude higher seasonal flux than OBGMs. The simulated interannual variability by the OBGMs is weaker than that found by atmospheric inversions. Prediction of such weak interannual variability in CO2 fluxes by atmospheric inversions was mainly caused by a lack of atmospheric data in the Indian Ocean. The OBGM models suggest a small strengthening of the sink over the period 1990–2009 of −0.01 PgC decade−1. This is inconsistent with the observations in the southwestern Indian Ocean that shows the growth rate of oceanic pCO2 was faster than the observed atmospheric CO2 growth, a finding attributed to the trend of the Southern Annular Mode (SAM) during the 1990s.

Description: V. V. S. S. Sarma acknowledges support and encouragement from S. W. A. Naqvi, Director, CSIR-National Institute of Oceanography. A. Lenton and R. M. Law acknowledge support from the Australian Climate Change Science Program, funded by the Australian Government Department of Industry, Innovation, Climate Change, Science, Research and Tertiary Education, and by the Bureau of Meteorology and by CSIRO. S. C. Doney and I. D. Lima acknowledge support from the National Science Foundation (NSF AGS-1048827). N. Metzl acknowledges support of the EU grant 264879 CARBOCHANGE.

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.