Description: The seasonal cycle is the dominant mode of variability in the air-sea CO2 flux in most regions of the global ocean, yet discrepancies between different seasonality estimates are rather large. As part of the Regional Carbon Cycle Assessment and Processes Phase 2 project (RECCAP2), we synthesize surface ocean pCO2 and air-sea CO2 flux seasonality from models and observation-based estimates, focusing on both a present-day climatology and decadal changes between the 1980s and 2010s. Four main findings emerge: First, global ocean biogeochemistry models (GOBMs) and observation-based estimates (pCO2 products) of surface pCO2 seasonality disagree in amplitude and phase, primarily due to discrepancies in the seasonal variability in surface DIC. Second, the seasonal cycle in pCO2 has increased in amplitude over the last three decades in both pCO2 products and GOBMs. Third, decadal increases in pCO2 seasonal cycle amplitudes in subtropical biomes for both pCO2 products and GOBMs are driven by increasing DIC concentrations stemming from the uptake of anthropogenic CO2 (Cant). In subpolar and Southern Ocean biomes, however, the seasonality change for GOBMs is dominated by Cant invasion, whereas for pCO2 products an indeterminate combination of Cant invasion and climate change modulates the changes. Fourth, biome-aggregated decadal changes in the amplitude of pCO2 seasonal variability are largely detectable against both mapping uncertainty (reducible) and natural variability uncertainty (irreducible), but not at the gridpoint scale over much of the northern subpolar oceans and over the Southern Ocean, underscoring the importance of sustained high-quality seasonally resolved measurements over these regions.

Description: Idealised and hindcast simulations performed with the stand-alone ocean carbon-cycle configuration of the Norwegian Earth System Model (NorESM-OC) are described and evaluated. We present simulation results of two different model versions at different grid resolutions and using two different atmospheric forcing data sets. Model version NorESM-OC1 corresponds to the version that is included in the fully coupled model NorESM-ME1, which participated in CMIP5. The main update between NorESM-OC1 and NorESM-OC1.2 is the addition of two new options for the treatment of sinking particles. We find that using a constant sinking speed, which has been the standard in NorESM's ocean carbon cycle module HAMOCC (HAMburg Ocean Carbon Cycle model) does not transport enough particulate organic carbon (POC) into the deep ocean below approximately 2000 m depth. The two newly implemented parameterisations, a particle aggregation scheme with prognostic sinking speed, and a simpler scheme prescribing a linear increase of sinking speed with depth, provide better agreement with observed POC fluxes. Additionally, reduced deep ocean biases of oxygen and remineralised phosphate indicate a better performance of the new parameterisations. For model version 1.2, a re-tuning of the ecosystem parameterisation has been performed, which (i) reduces previously too high primary production in high latitudes, (ii) consequently improves model results for surface nutrients, and (iii) reduces alkalinity and dissolved inorganic carbon biases at low latitudes. We use hindcast simulations with prescribed observed and constant (pre-industrial) atmospheric CO2 concentrations to derive the past and contemporary ocean carbon sink. For the period 1990–1999 we find an average ocean carbon uptake ranging from 2.01 to 2.58 Pg C yr-1 depending on model version, grid resolution and atmospheric forcing data set.