Search

Atmospheric carbon dioxide (CO₂) accounts for the largest radiative forcing among anthropogenic greenhouse gases. There is, therefore, a pressing need to understand the rate at which CO₂ accumulates in the atmosphere, including the interannual variations (IAVs) in this rate. IAV in the CO₂ growth rate is a small signal relative to the long-term trend and the mean annual cycle of atmospheric CO₂, and IAV is tied to climatic variations that may provide insights into long-term carbon–climate feedbacks. Observations from the Orbiting Carbon Observatory-2 (OCO-2) mission offer a new opportunity to refine our understanding of atmospheric CO₂ IAV since the satellite can measure over remote terrestrial regions and the open ocean, where traditional in situ CO₂ monitoring is difficult, providing better spatial coverage compared to ground-based monitoring techniques. In this study, we analyze the IAV of column-averaged dry-air CO₂ mole fraction (XCO₂) from OCO-2 between September 2014 and June 2021. The amplitude of the IAV, which is calculated as the standard deviation of the time series, is up to 1.2 ppm over the continents and around 0.4 ppm over the open ocean. Across all latitudes, the OCO-2-detected XCO₂ IAV shows a clear relationship with El Niño–Southern Oscillation (ENSO)-driven variations that originate in the tropics and are transported poleward. Similar, but smoother, zonal patterns of OCO-2 XCO₂ IAV time series compared to ground-based in situ observations and with column observations from the Total Carbon Column Observing Network (TCCON) and the Greenhouse Gases Observing Satellite (GOSAT) show that OCO-2 observations can be used reliably to estimate IAV. Furthermore, the extensive spatial coverage of the OCO-2 satellite data leads to smoother IAV time series than those from other datasets, suggesting that OCO-2 provides new capabilities for revealing small IAV signals despite sources of noise and error that are inherent to remote-sensing datasets.

GCAM model version,and processing scripts for figures.&nbsp

Marine heterotrophic Bacteria (or referred to as bacteria) play an important role in the ocean carbon cycle by utilizing, respiring, and remineralizing organic matter exported from the surface to deep ocean. Here, we investigate the responses of bacteria to climate change using a three-dimensional coupled ocean biogeochemical model with explicit bacterial dynamics as part of the Coupled Model Intercomparison Project Phase 6. First, we assess the credibility of the century-scale projections (2015–2099) of bacterial carbon stock and rates in the upper 100 m layer using skill scores and compilations of the measurements for the contemporary period (1988–2011). Second, we demonstrate that across different climate scenarios, the simulated bacterial biomass trends (2076–2099) are sensitive to the regional trends in temperature and organic carbon stocks. Bacterial carbon biomass declines by 5–10% globally, while it increases by 3–5% in the Southern Ocean where semi-labile dissolved organic carbon (DOC) stocks are relatively low and particle-attached bacteria dominate. While a full analysis of drivers underpinning the simulated changes in all bacterial stock and rates is not possible due to data constraints, we investigate the mechanisms of the changes in DOC uptake rates of free-living bacteria using the first-order Taylor decomposition. The results demonstrate that the increase in semi-labile DOC stocks drives the increase in DOC uptake rates in the Southern Ocean, while the increase in temperature drives the increase in DOC uptake rates in the northern high and low latitudes. Our study provides a systematic analysis of bacteria at global scale and a critical step toward a better understanding of how bacteria affect the functioning of the biological carbon pump and partitioning of organic carbon pools between surface and deep layers.

China's pledge to reach carbon neutrality before 2060 is an ambitious goal and could provide the world with much-needed leadership on how to limit warming to +1.5C warming above pre-industrial levels by the end of the century. But the pathways that would achieve net zero by 2060 are still unclear, including the role of negative emissions technologies. We use the Global Change Analysis Model to simulate how negative emissions technologies, in general, and direct air capture (DAC) in particular, could contribute to China's meeting this target. Our results show that negative emissions could play a large role, offsetting on the order of 3 GtCO2 per year from difficult-to-mitigate sectors such as freight transportation and heavy industry. This includes up to a 1.6 GtCO2 per year contribution from DAC, constituting up to 60% of total projected negative emissions in China. But DAC, like bioenergy with carbon capture and storage and afforestation, has not yet been demonstrated at anywhere approaching the scales required to meaningfully contribute to climate mitigation. Deploying NETs at these scales will have widespread impacts on financial systems and natural resources such as water, land, and energy in China.

We analyzed 20 years (1993–2013) of observations of dissolved inorganic macronutrients (nitrate, N; phosphate, P; and silicate, Si) and chlorophyll a (Chl) at Palmer Station, Antarctica (64.8°S, 64.1°W) to elucidate how large-scale climate and local physical forcing affect the interannual variability in the seasonal phytoplankton bloom and associated drawdown of nutrients. The leading modes of nutrients (N, P, and Si empirical orthogonal functions 1, EOF1) represent overall negative anomalies throughout growing seasons, showing a mixed signal of variability in the initial levels and drawdown thereafter (low-frequency dynamics). The second most common seasonal patterns of nitrate and phosphate (N and P EOF2) capture prolonged drawdown events during December–March, which are correlated to Chl EOF1. Si EOF2 captures a drawdown event during November–December, which is correlated to Chl EOF2. These different drawdown patterns are shaped by different sets of physical and climate forcing mechanisms. N and P drawdown events during December–March are influenced by the winter and spring Southern Annular Mode (SAM) phase, where nutrient utilization is enhanced in a stabilized upper water column as a consequence of SAM-driven winter sea ice and spring wind dynamics. Si drawdown during November–December is influenced by early sea ice retreat, where ice breakup may induce abrupt water column stratification and a subsequent diatom bloom or release of diatom cells from within the sea ice. Our findings underscore that seasonal nutrient dynamics in the coastal WAP are coupled to large-scale climate forcing and related physics, understanding of which may enable improved projections of biogeochemical responses to climate change.

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., Biogeosciences Discussions, 12, ISSN:1810-6277, ISSN:1810-6285

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). ISSN:1726-4170 ISSN:1726-4170

Here we analyze the impact of projected climate change on plankton ecology in all major ocean biomes over the 21st century, using a multidecade (1880–2090) experiment conducted with the Community Climate System Model (CCSM-3.1) coupled ocean-atmosphere-land-sea ice model. The climate response differs fundamentally in the Northern and Southern Hemispheres for diatom and small phytoplankton biomass and consequently for total biomass, primary, and export production. Increasing vertical stratification in the Northern Hemisphere oceans decreases the nutrient supply to the ocean surface. Resulting decreases in diatom and small phytoplankton biomass together with a relative shift from diatoms to small phytoplankton in the Northern Hemisphere result in decreases in the total primary and export production and export ratio, and a shift to a more oligotrophic, more efficiently recycled, lower biomass euphotic layer. By contrast, temperature and stratification increases are smaller in the Southern compared to the Northern Hemisphere. Additionally, a southward shift and increase in strength of the Southern Ocean westerlies act against increasing temperature and freshwater fluxes to destratify the water-column. The wind-driven, poleward shift in the Southern Ocean subpolar-subtropical boundary results in a poleward shift and increase in the frontal diatom bloom. This boundary shift, localized increases in iron supply, and the direct impact of warming temperatures on phytoplankton growth result in diatom increases in the Southern Hemisphere. An increase in diatoms and decrease in small phytoplankton partly compensate such that while total production and the efficiency of organic matter export to the deep ocean increase, total Southern Hemisphere biomass does not change substantially. The impact of ecological shifts on the global carbon cycle is complex and varies across ecological biomes, with Northern and Southern Hemisphere effects on the biological production and export partially compensating. The net result of climate change is a small Northern Hemisphere-driven decrease in total primary production and efficiency of organic matter export to the deep ocean.

publication Equipe CO2 qui intègre des analyses issues de OISO/CARAUS; International audience; The Indian Ocean (44° S-30° N) plays an important role in the global carbon cycle, yet 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 Pg C yr-1, is consistent with the -0.24 ± 0.12 Pg C yr-1 calculated from observations. The fluxes from the Southern Indian Ocean (18° S-44° S; -0.43 ± 0.07 Pg C 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), over estimation in the northeastern region (Bay of Bengal) and underestimation of the CO2 sink in the subtropical convergence zone. These differences were mainly driven by a 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 over estimate the magnitude. The predicted sea-air CO2 fluxes by Ocean BioGeochemical Models (OBGM) respond to seasonal variability with strong phase lags with reference to climatological CO2 flux, whereas the atmospheric inversions predict an order of magnitude higher seasonal flux than OBGMs. The simulated interannual variability by the OBGMs is weaker than atmospheric inversions. Prediction of such weak interannual variability in CO2 fluxes by atmospheric inversions was mainly caused by 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 Pg C decade-1. This is inconsistent with the observations in the southwest 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 Annual Mode (SAM) during the 1990s.

This study used an ocean general circulation model to simulate the marine iron cycle in an investigation of how simulated distributions of weak iron-binding ligands would be expected to control dissolved iron concentrations in the ocean, with a particular focus on deep ocean waters. The distribution of apparent oxygen utilization was used as a proxy for humic substances that have recently been hypothesized to account for the bulk of weak iron-binding ligands in seawater. Compared to simulations using a conventional approach with homogeneous ligand distributions, the simulations that incorporated spatially variable ligand concentrations exhibited substantial improvement in the simulation of global dissolved iron distributions as revealed by comparisons with available field data. The improved skill of the simulations resulted largely because the spatially variable ligand distributions led to a more reasonable basin-scale variation of the residence time of iron when present at high concentrations. The model results, in conjunction with evidence from recent field studies, suggest that humic substances play an important role in the iron cycle in the ocean.

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., Biogeosciences Discussions, 10, ISSN:1810-6277, ISSN:1810-6285

Detection and attribution of hydrographic and biogeochemical changes in the deep ocean are challenging due to the small magnitude of their signals and to limitations in the accuracy of available data. However, there are indications that anthropogenic and climate change signals are starting to manifest at depth. The deep ocean below 2000 m comprises about 50% of the total ocean volume, and changes in the deep ocean should be followed over time to accurately assess the partitioning of anthropogenic carbon dioxide (CO₂) between the ocean, terrestrial biosphere, and atmosphere. Here we determine the changes in the interior deep-water inorganic carbon content by a novel means that uses the partial pressure of CO₂ measured at 20 °C, pCO₂(20), along three meridional transects in the Atlantic and Pacific oceans. These changes are measured on decadal time scales using observations from the World Ocean Circulation Experiment (WOCE)/World Hydrographic Program (WHP) of the 1980s and 1990s and the CLIVAR/CO₂ Repeat Hydrography Program of the past decade. The pCO₂(20) values show a consistent increase in deep water over the time period. Changes in total dissolved inorganic carbon (DIC) content in the deep interior are not significant or consistent, as most of the signal is below the level of analytical uncertainty. Using an approximate relationship between pCO₂(20) and DIC change, we infer DIC changes that are at the margin of detectability. However, when integrated on the basin scale, the increases range from 8–40% of the total specific water column changes over the past several decades. Patterns in chlorofluorocarbons (CFCs), along with output from an ocean model, suggest that the changes in pCO₂(20) and DIC are of anthropogenic origin.

Estimates of the anthropogenic global-integrated sea-air carbon dioxide (CO2) fluxfrom 1990 to 2009, based on different models and measurements, range from−1.4to−2.6 PgC yr−1. The median values of anthropogenic CO2for each method showbetter agreement and are:−1.9 for Pg C yr−1for numerical ocean general circulation hind cast models (OGCMs) with parameterized biogeochemistry;−2.1 PgC yr−1for at-mospheric inverse models;−1.9 PgC yr−1for global atmospheric constraints based onO2/ N2ratios for 1990–2000; and−2.4 PgC yr−1for oceanic inverse models. An up-dated estimate of this anthropogenic CO2flux based on a climatology of sea-air partialpressure of CO2differences (∆pCO2) (Takahashi et al., 2009) and a bulk formulation of gas transfer with wind speed for year 2000 is−2.0 PgC yr−1. Using this∆pCO2climatology and empirical relationships ofpCO2with sea-surface temperature (SST)anomalies (Park et al., 2010a), the interannual variability of the contemporary CO2flux is estimated to be 0.20Pg C yr−1(1σ) from 1990 through 2009. This is similar tothe variability estimated by the OGCMs of 0.16 Pg Cyr−1but smaller than the interannual variability from atmospheric inverse estimates of 0.40Pg C yr−1. The variability islargely driven by large-scale climate re-organizations. The decadal trends for differentmethods range from−0.13 (Pg Cyr−1) decade−1to−0.50 (Pg Cyr−1) decade−1. TheOGCMs and the data based sea-air CO2flux estimates show smaller uptakes and ap-preciably smaller decadal trends than estimates based on changes in carbon inventory suggesting that methods capable of resolving shorter timescales are showing a slow-ing of the rate of ocean CO2uptake. It is not clear if this large difference in trend is amethodological issue or a real natural feedback, Biogeosciences Discussions, 9, ISSN:1810-6277, ISSN:1810-6285

ISSN:1810-6277 ISSN:1810-6285

he global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data-based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air–sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on forward ocean models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 ± 31 PgC (±20% uncertainty). This estimate includes a broad range of values, suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2. ISSN:1810-6277 ISSN:1810-6285

We present a summary of biomass data for 11 Plankton Functional Types (PFTs)plus phytoplankton pigment data, compiled as part of the MARine Ecosystem biomassDATa (MAREDAT) initiative. The goal of the MAREDAT initiative is to provide global gridded data products with coverage of all biological components of the global ocean ecosystem. This special issue is the first step towards achieving this. The PFTs presented here include picophytoplankton, diazotrophs, coccolithophores,Phaeocystis, diatoms, picoheterotrophs, microzooplankton, foraminifers, mesozooplankton, pteropodsand macrozooplankton. All variables have been gridded onto a World Ocean Atlas (WOA) grid (1° × 1° × 33 vertical levels × monthly climatologies). The data show that (1) the global total heterotrophic biomass (2.0–6.4Pg C) is at least as high as the totalautotrophic biomass (0.5–2.6 Pg C excluding nanophytoplankton and autotrophic dinoflagellates), (2) the biomass of zooplankton calcifiers (0.9–2.3Pg C) is substantiallyhigher than that of coccolithophores (0.01–0.14 Pg C), (3) patchiness of biomass distribution increases with organism size, and (4) although zooplankton biomass measurements below 200m are rare, the limited measurements available suggest thatBacteriaandArchaeaare not the only heterotrophs in the deep sea. More data will be neededto characterize ocean ecosystem functioning and associated biogeochemistry in theSouthern Hemisphere and below 200 m. ISSN:1866-3591

Biogeosciences Discussions, 9 (8), ISSN:1810-6277, ISSN:1810-6285

International audience; In the North Pacific Subtropical Gyre (NPSG), the regular occurrence of summer phytoplankton blooms contributes to marine ecosystem productivity and the annual carbon export. The mechanisms underlying the formation, maintenance, and decay of these blooms remain largely unknown; nitrogen fixation, episodic vertical mixing of nutrients, and meso- (

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., Biogeosciences Discussions, 9 (8), ISSN:1810-6277, ISSN:1810-6285

The performance of 36 models (22 ocean color models and 14 biogeochemical ocean circulation models (BOGCMs)) that estimate depth-integrated marine net primary productivity (NPP) was assessed by comparing their output to in situ 14C data at the Bermuda Atlantic Time series Study (BATS) and the Hawaii Ocean Time series (HOT) over nearly two decades. Specifically, skill was assessed based on the models' ability to estimate the observed mean, variability, and trends of NPP. At both sites, more than 90% of the models underestimated mean NPP, with the average bias of the BOGCMs being nearly twice that of the ocean color models. However, the difference in overall skill between the best BOGCM and the best ocean color model at each site was not significant. Between 1989 and 2007, in situ NPP at BATS and HOT increased by an average of nearly 2% per year and was positively correlated to the North Pacific Gyre Oscillation index. The majority of ocean color models produced in situ NPP trends that were closer to the observed trends when chlorophyll-a was derived from high-performance liquid chromatography (HPLC), rather than fluorometric or SeaWiFS data. However, this was a function of time such that average trend magnitude was more accurately estimated over longer time periods. Among BOGCMs, only two individual models successfully produced an increasing NPP trend (one model at each site). We caution against the use of models to assess multiannual changes in NPP over short time periods. Ocean color model estimates of NPP trends could improve if more high quality HPLC chlorophyll-a time series were available. © 2010 by the American Geophysical Union.

We present results from transient sensitivity studies with the Biogeochemical Elemental Cycling (BEC) ocean model to increasing anthropogenic atmospheric inorganic nitrogen (N) and soluble iron (Fe) deposition over the industrial era. Elevated N deposition results from fossil fuel combustion and agriculture, and elevated soluble Fe deposition results from increased atmospheric processing in the presence of anthropogenic pollutants and soluble Fe from combustion sources. Simulations with increasing Fe and increasing Fe and N inputs raised simulated marine nitrogen fixation, with the majority of the increase in the subtropical North and South Pacific, and raised primary production and export in the high-nutrient low-chlorophyll (HNLC) regions. Increasing N inputs alone elevated small phytoplankton and diatom production, resulting in increased phosphorus (P) and Fe limitation for diazotrophs, hence reducing nitrogen fixation (∼6%). Globally, the simulated primary production, sinking particulate organic carbon (POC) export. and atmospheric CO2 uptake were highest under combined increase in Fe and N inputs compared to preindustrial control. Our results suggest that increasing combustion iron sources and aerosol Fe solubility along with atmospheric anthropogenic nitrogen deposition are perturbing marine biogeochemical cycling and could partially explain the observed trend toward increased P limitation at station ALOHA in the subtropical North Pacific. Excess inorganic nitrogen ([NO3 −] + [NH4 +] − 16[PO4 3−]) distributions may offer useful insights for understanding changing ocean circulation and biogeochemistry.

We synthesize estimates of the contemporary net air-sea CO2 flux on the basis of an inversion of interior ocean carbon observations using a suite of 10 ocean general circulation models (Mikaloff Fletcher et al., 2006, 2007) and compare them to estimates based on a new climatology of the air-sea difference of the partial pressure of CO2 (pCO(2)) (Takahashi et al., 2008). These two independent flux estimates reveal a consistent description of the regional distribution of annual mean sources and sinks of atmospheric CO2 for the decade of the 1990s and the early 2000s with differences at the regional level of generally less than 0.1 Pg C a(-1). This distribution is characterized by outgassing in the tropics, uptake in midlatitudes, and comparatively small fluxes in the high latitudes. Both estimates point toward a small(similar to -0.3 Pg C a(-1)) contemporary CO2 sink in the Southern Ocean (south of 44 degrees S), a result of the near cancellation between a substantial outgassing of natural CO2 and a strong uptake of anthropogenic CO2. A notable exception in the generally good agreement between the two estimates exists within the Southern Ocean: the ocean inversion suggests a relatively uniform uptake, while the pCO(2)-based estimate suggests strong uptake in the region between 58 degrees S and 44 degrees S, and a source in the region south of 58 degrees S. Globally and for a nominal period between 1995 and 2000, the contemporary net air-sea flux of CO2 is estimated to be -1.7 +/- 0.4 Pg C a(-1) (inversion) and -1.4 +/- 0.7 Pg C a(-1) (pCO(2)-climatology), respectively, consisting of an outgassing flux of river-derived carbon of similar to+0.5 Pg C a(-1), and an uptake flux of anthropogenic carbon of -2.2 +/- 0.3 Pg C a(-1) (inversion) and -1.9 +/- 0.7 Pg C a(-1) (pCO(2)-climatology). The two flux estimates also imply a consistent description of the contemporary meridional transport of carbon with southward ocean transport throughout most of the Atlantic basin, and strong equatorward convergence in the Indo-Pacific basins. Both transport estimates suggest a small hemispheric asymmetry with a southward transport of between -0.2 and -0.3 Pg C a(-1) across the equator. While the convergence of these two independent estimates is encouraging and suggests that it is now possible to provide relatively tight constraints for the net air-sea CO2 fluxes at the regional basis, both studies are limited by their lack of consideration of long-term changes in the ocean carbon cycle, such as the recent possible stalling in the expected growth of the Southern Ocean carbon sink.

International audience; Observational studies report a rapid decline of ocean CO2 uptake in the temperate North Atlantic during the last decade. We analyze these findings using ocean physical-biological numerical simulations forced with interannually varying atmospheric conditions for the period 1979-2004. In the simulations, surface ocean water mass properties and CO2 system variables exhibit substantial multiannual variability on sub-basin scales in response to wind-driven reorganization in ocean circulation and surface warming/cooling. The simulated temporal evolution of the ocean CO2 system is broadly consistent with reported observational trends and is influenced substantially by the phase of the North Atlantic Oscillation (NAO). Many of the observational estimates cover a period after 1995 of mostly negative or weakly positive NAO conditions, which are characterized in the simulations by reduced North Atlantic Current transport of subtropical waters into the eastern basin and by a decline in CO2 uptake. We suggest therefore that air-sea CO2 uptake may rebound in the eastern temperate North Atlantic during future periods of more positive NAO, similar to the patterns found in our model for the sustained positive NAO period in the early 1990s. Thus, our analysis indicates that the recent rapid shifts in CO2 flux reflect decadal perturbations superimposed on more gradual secular trends. The simulations highlight the need for long-term ocean carbon observations and modeling to fully resolve multiannual variability, which can obscure detection of the long-term changes associated with anthropogenic CO2 uptake and climate change.

Seasonal and interannual variability in atmospheric carbon dioxide (CO2) concentrations was simulated using fluxes from fossil fuel, ocean and terrestrial biogeochemical models, and a tracer transport model with time-varying winds. The atmospheric CO2 variability resulting from these surface fluxes was compared to observations from 89 GLOBALVIEW monitoring stations. At northern hemisphere stations, the model simulations captured most of the observed seasonal cycle in atmospheric CO2, with the land tracer accounting for the majority of the signal. The ocean tracer was 3–6 months out of phase with the observed cycle at these stations and had a seasonal amplitude only ∼10% on average of observed. Model and observed interannual CO2 growth anomalies were only moderately well correlated in the northern hemisphere (R ∼ 0.4–0.8), and more poorly correlated in the southern hemisphere (R < 0.6). Land dominated the interannual variability (IAV) in the northern hemisphere, and biomass burning in particular accounted for much of the strong positive CO2 growth anomaly observed during the 1997–1998 El Niño event. The signals in atmospheric CO2 from the terrestrial biosphere extended throughout the southern hemisphere, but oceanic fluxes also exerted a strong influence there, accounting for roughly half of the IAV at many extratropical stations. However, the modeled ocean tracer was generally uncorrelated with observations in either hemisphere from 1979–2004, except during the weak El Niño/post-Pinatubo period of the early 1990s. During that time, model results suggested that the ocean may have accounted for 20–25% of the observed slowdown in the atmospheric CO2 growth rate

La principal preocupación de muchos científicos hoy en día es que el episodio de acidificación actual está ocurriendo mucho más rápidamente de los que tuvieron lugar en el pasado, por lo que no da tiempo a que las especies oceánicas se adapten a estos cambios. Aunque los efectos puedan estar ocultos a la vista de la gente, las dramáticas alteraciones del ambiente marino parecen inevitables., Resumen del siguiente artículo: Doney, Scott C. (2006). "The dangers of ocean acidification". Scientific American 294(3): 38-45., Fundación Museo La Plata

New observations from the North Sea, a NW European shelf sea, show that between 2001 and 2005 the CO2 partial pressure (pCO(2)) in surface waters rose by 22 mu atm, thus faster than atmospheric pCO(2), which in the same period rose approximately 11 matm. The surprisingly rapid decline in air-sea partial pressure difference (Delta pCO(2)) is primarily a response to an elevated water column inventory of dissolved inorganic carbon (DIC), which, in turn, reflects mostly anthropogenic CO2 input rather than natural interannual variability. The resulting decline in the buffering capacity of the inorganic carbonate system (increasing Revelle factor) sets up a theoretically predicted feedback loop whereby the invasion of anthropogenic CO2 reduces the ocean's ability to uptake additional CO2. Model simulations for the North Atlantic Ocean and thermodynamic principles reveal that this feedback should be stronger, at present, in colder midlatitude and subpolar waters because of the lower present-day buffer capacity and elevated DIC levels driven either by northward advected surface water and/or excess local air-sea CO2 uptake. This buffer capacity feedback mechanism helps to explain at least part of the observed trend of decreasing air-sea Delta pCO(2) over time as reported in several other recent North Atlantic studies.

International audience

Recent upward revisions in key sink/source terms for fixed nitrogen (N) in the oceans imply a short residence time and strong negative feedbacks involving denitrification and N fixation to prevent large swings in the ocean N inventory over timescales of a few centuries. We tested the strength of these feedbacks in a global biogeochemical elemental cycling (BEC) ocean model that includes water column denitrification and an explicit N fixing phytoplankton group. In the northern Indian Ocean and over longer timescales in the tropical Atlantic, we find strong stabilizing feedbacks that minimize changes in marine N inventory over timescales of ∼30–200 years. In these regions high atmospheric dust/iron inputs lead to phosphorus limitation of diazotrophs, and thus a tight link between N fixation and surface water N/P ratios. Maintenance of the oxygen minimum zones in these basins depends on N fixation driven export. The stabilizing feedbacks in other regions are significant but weaker owing to iron limitation of the diazotrophs. Thus Fe limitation appears to restrict the ability of N fixation to compensate for changes in denitrification in the current climate, perhaps leading the oceans to lose fixed N. We suggest that iron is the ultimate limiting nutrient leading to nitrogen being the proximate limiting nutrient over wide regions today. Iron stress was at least partially alleviated during more dusty, glacial times, leading to a higher marine N inventory, increased export production, and perhaps widespread phosphorus limitation of the phytoplankton community. The increased efficiency of the biological pump would have contributed to the glacial drawdown in atmospheric CO2.

Estimates of sources/sinks of carbon dioxide (CO2) at the Earth's surface are commonly made using atmospheric CO2 inverse modeling, terrestrial and oceanic biogeochemical modeling, and inventory-based studies. In this study, we compare sea-air CO2 fluxes from the Time-Dependent Inverse (TDI) atmosphere model and the marine Biogeochemical Elemental Cycling (BEC) model to study the processes involved in ocean carbon cycling at subbasin scales. A dust generation and transport model, based on analyzed meteorology and terrestrial vegetation cover, is also used to estimate the interannual variability in dust and iron deposition to different ocean basins. Overall, a fairly good agreement is established between the TDI and BEC model results for the net annual patterns and seasonal cycle of sea-air CO2 exchange. Sensitivity studies with the ocean biogeochemical model using increased or reduced atmospheric iron inputs indicate the relative sensitivity of air-sea CO2 exchange. The simulated responses to changes in iron inputs are not instantaneous (peak response after ∼2−3 years). The TDI model derived seasonal cycles for the Southern Ocean (South Atlantic) are better matched by the BEC model by increasing (decreasing) iron inputs through atmospheric aerosols. Our results suggest that some of the interannual variability in TDI model air-sea CO2 fluxes during the past decade may be explainable by dust variability that relaxes/increases iron limitation in high-nitrate, low-chlorophyll (HNLC) ocean regions.

The sources and sinks of important climatic trace gases such as carbon dioxide (CO2) are often deduced from spatial and temporal variations in atmospheric concentrations. Reducing uncertainties in our understanding of the contemporary carbon budget and its underlying dynamics, however, requires significantly denser observations globally than is practical with in situ measurements. Space-based measurements appear technically feasible but require innovations in data analysis approaches. We develop a variational data assimilation scheme to estimate surface CO2 fluxes at fine time/space scales from such dense atmospheric data. Global flux estimates at a daily time step and model-grid spatial resolution (4°× 5° here) are rapidly achieved after only a few dozen minimization steps. We quantify the flux errors from existing, planned and hypothetical surface and space-borne observing systems. Simulations show that the planned NASA Orbital Carbon Observatory (OCO) satellite should provide significant additional information beyond that from existing and proposed in situ observations. Improvements in data assimilation techniques and in mechanistic process models are both needed to fully exploit the emerging global carbon observing system.DOI: 10.1111/j.1600-0889.2006.00218.x

Satellite remote sensing estimates of surface chlorophyll, temperature, wind speed, and sea ice cover are examined in the region of the Southern Ocean Iron Experiment (SOFeX). Our objectives are to place SOFeX into a regional context and highlight regional mesoscale spatial and monthly temporal variability. SOFeX fertilized two patches with iron, one south of the Antarctic Polar Front (PF) and one north of the PF but south of the Subantarctic Front (SAF). Satellite-observable phytoplankton blooms developed in both patches. The spring sea ice retreat near the south patch site was delayed in the 2001–2002 season, in turn delaying the naturally occurring, modest spring bloom in this region. Ambient surface chlorophyll concentrations for the area surrounding the southern patch during January 2002 are low (mean 0.26 mg/m3) compared with climatological January values (0.42 mg/m3). Regions east and west at similar latitudes exhibited higher mean chlorophyll concentrations (0.79 and 0.74 mg/m3, respectively). These modest phytoplankton blooms were likely stimulated by melting sea ice via changes in the light-mixing regime and release of iron and were smaller in magnitude than the iron-induced bloom within the SOFeX southern patch (>3 mg/m3). Iron inputs from melting ice may drive much of the natural spatial and temporal variability within the seasonal ice zone. Mean chlorophyll concentrations surrounding the SOFeX northern patch site were similar to climatological values during the SOFeX season. The northern patch was stretched into a long, thin filament along the southern boundary of the SAF, likely increasing the mixing/dilution rate with surrounding waters.

A global three-dimensional marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary and export production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed high nitrate, low chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 55 Tg N. Diazotrophs directly account for a small fraction of primary production (0.5%) but indirectly support 10% of primary production and 8% of sinking particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in the upper ocean is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, owing to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.

We study the dynamics of the planktonic ecosystem in the coastal upwelling zone within the California Current System using a three-dimensional (3-D), eddy-resolving circulation model coupled to an ecosystem/biogeochemistry model. The physical model is based on the Regional Oceanic Modeling System (ROMS), configured at a resolution of 15 km for a domain covering the entire US West Coast, with an embedded child grid covering the central California upwelling region at a resolution of 5 km. The model is forced with monthly mean boundary conditions at the open lateral boundaries as well as at the surface. The ecological/biogeochemical model is nitrogen based, includes single classes for phytoplankton and zooplankton, and considers two detrital pools with different sinking speeds. The model also explicitly simulates a variable chlorophyll-to-carbon ratio. Comparisons of model results with either remote sensing observations (AVHRR, SeaWiFS) or in-situ measurements from the CalCOFI program indicate that our model is capable of replicating many of the large-scale, time-averaged features of the coastal upwelling system. An exception is the underestimation of the chlorophyll levels in the northern part of the domain, perhaps because of the lack of short-term variations in the atmospheric forcing. Another shortcoming is that the modeled thermocline is too diffuse, and that the upward slope of the isolines toward the coast is too small. Detailed time-series comparisons with observations from Monterey Bay reveal similar agreements and discrepancies. We attribute the good agreement between the modeled and observed ecological properties in large part to the accuracy of the physical fields. In turn, many of the discrepancies can be traced back to our use of monthly mean forcing. Analysis of the ecosystem structure and dynamics reveal that the magnitude and pattern of phytoplankton biomass in the nearshore region are determined largely by the balance of growth and zooplankton grazing, while in the offshore region, growth is balanced by mortality. The latter appears to be inconsistent with in situ observations and is a result of our consideration of only one zooplankton size class (mesozooplankton), neglecting the importance of microzooplankton grazing in the offshore region. A comparison of the allocation of nitrogen into the different pools of the ecosystem in the 3-D results with those obtained from a box model configuration of the same ecosystem model reveals that only a few components of the ecosystem reach a local steady-state, i.e. where biological sources and sinks balance each other. The balances for the majority of the components are achieved by local biological source and sink terms balancing the net physical divergence, confirming the importance of the 3-D nature of circulation and mixing in a coastal upwelling system.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – 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 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 and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, 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 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. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. 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 (2006–2015), EFF was 9.3 ± 0.5 GtC yr−1, ELUC 1.0 ± 0.5 GtC yr−1, GATM 4.5 ± 0.1 GtC yr−1, SOCEAN 2.6 ± 0.5 GtC yr−1, and SLAND 3.1 ± 0.9 GtC yr−1. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr−1 that took place during 2006–2015. Also, for 2015, ELUC was 1.3 ± 0.5 GtC yr−1, GATM was 6.3 ± 0.2 GtC yr−1, SOCEAN was 3.0 ± 0.5 GtC yr−1, and SLAND was 1.9 ± 0.9 GtC yr−1. GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565 ± 55 GtC (2075 ± 205 GtCO2) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2016).

We investigated the mechanisms of phytoplankton competition during the spring bloom,one of the most dramatic seasonal events in lower-trophic level ecosystems, in fourstate-of-the-art Plankton Functional Type (PFTs) models: PISCES, NEMURO, Plank-TOM5 and CCSM-BEC. In particular, we investigated the relative importance of dif-5ferent ecophysiological processes on the determination of the community structure,focusing both on the bottom-up and the top-down controls. The models reasonablyreproduced the observed global distribution and seasonal variation of phytoplanktonbiomass. The fraction of diatoms with respect to the total phytoplankton biomass in-creases with the magnitude of the spring bloom in all models. However, the governing10mechanisms differ between models, despite the fact that current PFT models repre-sent ecophysiological processes using the same types of parameterizations. The in-creasing trend in the percentage of diatoms with increasing bloom magnitude is mainlycaused by a stronger nutrient dependence of photosynthesis for diatoms compared tonanophytoplankton (bottom-up control). The difference in the maximum photosynthesis rate plays an important role in NEMURO and PlankTOM5 and determines the abso-lute values of the percentage of diatoms during the bloom. In CCSM-BEC, the lightdependency of photosynthesis plays an important role in the North Atlantic and theSouthern Ocean. The grazing pressure by zooplankton (top-down control), however,strongly contributes to the dominance of diatoms in PISCES and CCSM-BEC. The re-20gional differences in the percentage of diatoms in PlankTOM5 are mainly determinedby top-down control. These differences in the mechanisms suggest that the responseof marine ecosystems to climate change could significantly differ among models, evenif the present-day ecosystem is reproduced to a similar degree of confidence. For fur-ther understanding of plankton competition and for the prediction of future change in marine ecosystems, it is important to understand the relative differences in each phys-iological rate and life history rate in the bottom-up and the top-down controls betweenPFTs, Biogeosciences Discussions, 9 (12), ISSN:1810-6277, ISSN:1810-6285

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., Biogeosciences, 13 (13), ISSN:1726-4170

We present a summary of biomass data for 11 plankton functional types (PFTs) plus phytoplankton pigment data, compiled as part of the MARine Ecosystem biomass DATa (MAREDAT) initiative. The goal of the MAREDAT initiative is to provide, in due course, global gridded data products with coverage of all planktic components of the global ocean ecosystem. This special issue is the first step towards achieving this. The PFTs presented here include picophytoplankton, diazotrophs, coccolithophores, Phaeocystis, diatoms, picoheterotrophs, microzooplankton, foraminifers, mesozooplankton, pteropods and macrozooplankton. All variables have been gridded onto a World Ocean Atlas (WOA) grid (1° × 1° × 33 vertical levels × monthly climatologies). The results show that abundance is much better constrained than their carbon content/elemental composition, and coastal seas and other high productivity regions have much better coverage than the much larger volumes where biomass is relatively low. The data show that (1) the global total heterotrophic biomass (2.0–4.6 Pg C) is at least as high as the total autotrophic biomass (0.5–2.4 Pg C excluding nanophytoplankton and autotrophic dinoflagellates); (2) the biomass of zooplankton calcifiers (0.03–0.67 Pg C) is substantially higher than that of coccolithophores (0.001–0.03 Pg C); (3) patchiness of biomass distribution increases with organism size; and (4) although zooplankton biomass measurements below 200 m are rare, the limited measurements available suggest that Bacteria and Archaea are not the only important heterotrophs in the deep sea. More data will be needed to characterise ocean ecosystem functioning and associated biogeochemistry in the Southern Hemisphere and below 200 m. Future efforts to understand marine ecosystem composition and functioning will be helped both by further archiving of historical data and future sampling at new locations., Earth System Science Data, 5 (2), ISSN:1866-3516, ISSN:1866-3508

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., Biogeosciences, 12 (3), ISSN:1726-4170

The extent, duration, and seasonality of sea ice and glacial discharge strongly influence Antarctic marine ecosystems. Most organisms' life cycles in this region are attuned to ice seasonality. The annual retreat and melting of sea ice in the austral spring stratifies the upper ocean, triggering large phytoplankton blooms. The magnitude of the blooms is proportional to the winter extent of ice cover, which can act as a barrier to wind mixing. Antarctic krill, one of the most abundant metazoan populations on Earth, consume phytoplankton blooms dominated by large diatoms. Krill, in turn, support a large biomass of predators, including penguins, seals, and whales. Human activity has altered even these remote ecosystems. The western Antarctic Peninsula region has warmed by 7°C over the past 50 years, and sea ice duration has declined by almost 100 days since 1978, causing a decrease in phytoplankton productivity in the northern peninsula region. Besides climate change, Antarctic marine systems have been greatly altered by harvesting of the great whales and now krill. It is unclear to what extent the ecosystems we observe today differ from the pristine state.

In the North Pacific Subtropical Gyre (NPSG), the regular occurrence of summer phytoplankton blooms contributes to marine ecosystem productivity and the annual carbon export. The mechanisms underlying the formation, maintenance, and decay of these blooms remain largely unknown; nitrogen fixation, episodic vertical mixing of nutrients, and meso- (

Accurate assessmen of anthropogenic carbon dioxide (CO₂) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to betterunderstand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and 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 and consistency within and among components, alongside methodology and data limitations. CO₂ emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respec- tively, 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 CO₂ concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO₂ sink (SOCEAN) is based on observations from the 1990s, while the an- nual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO₂ measurements. The global residual terrestrial CO₂ 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. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. 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 (2006–2015), EFF was 9.3± 0.5 GtC yr⁻¹, ELUC 1.0± 0.5 GtC yr⁻¹,GATM 4.5± 0.1 GtC yr⁻¹, SOCEAN 2.6± 0.5 GtC yr⁻¹, and SLAND 3.1± 0.9 GtC yr⁻¹. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9± 0.5 GtC yr⁻¹, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr⁻¹ that took place during 2006–2015. Also, for 2015, ELUC was 1.3± 0.5 GtC yr⁻¹,GATM was 6.3± 0.2 GtC yr⁻¹, SOCEAN was 3.0± 0.5 GtC yr⁻¹, and SLAND was 1.9± 0.9 GtC yr⁻¹ GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO₂ concentration reached 399.4± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to + 1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO₂ concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO₂ will reach 565± 55 GtC (2075± 205 GtCO₂) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013).

The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF₆) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and ¹⁴CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO₂ as well as carbon isotope ratios of CO₂. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

We present a summary of biomass data for 11 Plankton Functional Types (PFTs)plus phytoplankton pigment data, compiled as part of the MARine Ecosystem biomassDATa (MAREDAT) initiative. The goal of the MAREDAT initiative is to provide global gridded data products with coverage of all biological components of the global ocean ecosystem. This special issue is the first step towards achieving this. The PFTs presented here include picophytoplankton, diazotrophs, coccolithophores,Phaeocystis, diatoms, picoheterotrophs, microzooplankton, foraminifers, mesozooplankton, pteropodsand macrozooplankton. All variables have been gridded onto a World Ocean Atlas (WOA) grid (1° × 1° × 33 vertical levels × monthly climatologies). The data show that (1) the global total heterotrophic biomass (2.0–6.4Pg C) is at least as high as the totalautotrophic biomass (0.5–2.6 Pg C excluding nanophytoplankton and autotrophic dinoflagellates), (2) the biomass of zooplankton calcifiers (0.9–2.3Pg C) is substantiallyhigher than that of coccolithophores (0.01–0.14 Pg C), (3) patchiness of biomass distribution increases with organism size, and (4) although zooplankton biomass measurements below 200m are rare, the limited measurements available suggest thatBacteriaandArchaeaare not the only heterotrophs in the deep sea. More data will be neededto characterize ocean ecosystem functioning and associated biogeochemistry in theSouthern Hemisphere and below 200 m., Earth System Science Data Discussions, 5 (2), ISSN:1866-3591

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., Biogeosciences, 10 (1), ISSN:1726-4170

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., Biogeosciences, 10 (1), ISSN:1726-4170

The global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data-based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air–sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on forward ocean models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 ± 31 PgC (±20% uncertainty). This estimate includes a broad range of values, suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2., Biogeosciences, 10 (4), ISSN:1726-4170

We propose to add dissolved oxygen sensors to the Argo (Global array of free-drifting profiling floats)float program in order to determine, on a global-scale, seasonal to decadal time-scale variations in dissolved oxygen concentrations throughout the upper ocean. Such observations are especially important to document the ocean’s expected loss of oxygen as a result of ocean warming, but there are many other benefits including the opportunity to estimate net community and export production, the assessment of changes in low oxygen regions, and improved estimates of theoceanic uptake of anthropogenic CO2. The proposed joint Argo-Oxygen program is made possible by the recent development of dissolved oxygen sensors that are both precise and stable over extended periods and can be easily integrated with the currently usedArgo floats. Results from the more than 200 oxygen equipped Argo float have not only demonstrated the feasibility of the program, but also produced already many insights and discoveries. Achieving the main goal of the Argo-Oxygenprogram does not require any appreciable changes in the deployment and operating strategies of the current Argo program and can therefore be implemented without a significant impact on Argo’s core mission. A two-phase implementation is proposed, consisting of a small set of regional pilot-studies, followed by a build-up toward a global implementation. The cost of adding oxygen sensors to all floats of the Argo program is estimated to be about USD 4.2 million per year. The proposed Argo-Oxygenprogram will add substantial value to Argo by expanding the number of Argo data users, as well as by creating new synergies between the physical and the biogeochemical ocean research communities. The new observations will also contribute to the activities of various international networks and partnerships for Earth Observing Systems.

he global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data-based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air–sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on forward ocean models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 ± 31 PgC (±20% uncertainty). This estimate includes a broad range of values, suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2., Biogeosciences Discussions, 9, ISSN:1810-6277, ISSN:1810-6285