Model constraints on the anthropogenic carbon budget of the Arctic Ocean
Abstract. The Arctic Ocean is projected to experience not only amplified climate change but also amplified ocean acidification. Modeling future acidification depends on our ability to simulate baseline conditions and changes over the industrial era. Such centennial-scale changes require a global model to account for exchange between the Arctic and surrounding regions. Yet the coarse resolution of typical global models may poorly resolve that exchange as well as critical features of Arctic Ocean circulation. Here we assess how simulations of Arctic Ocean storage of anthropogenic carbon (Cant), the main driver of open- ocean acidification, differ when moving from coarse to eddy admitting resolution in a global ocean circulation-biogeochemistry model (NEMO-PISCES). The Arctic's regional storage of Cant is enhanced as model resolution increases. While the coarse- resolution model configuration ORCA2 (2°) stores 2.0 Pg C in the Arctic Ocean between 1765 and 2005, the eddy-admitting versions ORCA05 and ORCA025 (1/2° and 1/4°) store 2.4 and 2.6 Pg C. That result from ORCA025 falls within the uncertainty range from a previous data-based Cant storage estimate (2.5 to 3.3 Pg C). Yet those limits may each need to be reduced by about 10 % because data-based Cant concentrations in deep waters remain at ∼ 6 μmol kg−1, while they should be almost negligible by analogy to the near-zero observed CFC-12 concentrations from which they are calculated. Across the three resolutions, there was roughly three times as much anthropogenic carbon that entered the Arctic Ocean through lateral transport than via the flux of CO2 across the air-sea interface. Wider comparison to nine earth system models that participated in the Coupled Model Intercomparison Project Phase 5 (CMIP5) reveals much larger diversity of stored anthropogenic carbon and lateral transport. Only the CMIP5 models with higher lateral transport obtain Cant inventories that are close to the data-based estimates. Increasing resolution also enhances acidification, e.g., with greater shoaling of the Arctic's average depth of the aragonite saturation horizon during 1960–2012, from 50 m in ORCA2 to 210 m in ORCA025. To assess the potential to further refine modeled estimates of the Arctic Ocean's Cant storage and acidification, sensitivity tests that adjust model parameters are needed given that century-scale global ocean biogeochemical simulations still cannot be run routinely at high resolution.