Ocean Carbon Cycle Model Research Articles

Periodic Solutions and Bifurcations in a Long-Timescale Oceanic Carbon Cycle Model

We study the dynamics of a recent model of Rothman for long timescale carbon cycle. We reproduce and extend variousresults of the Rothman model. We present numerical results showing that the model exhibits both stable and unstable limitcycles via Hopf bifurcations as the parameters are varied. We numerically find normal forms of Bautin bifurcations to confirmtheir criticality. We also extend the analysis of the normal form coefficients to identify where the fold limit cycle bifurcationoccurs.

Pemodelan Sistem Karbonat di Laut Jawa

<p class="Papertext"><strong>Modeling Carbonate System in the Java Sea</strong>. Besides the global fossil fuel burning activities, forest fires in Kalimantan could potentially increase atmospheric CO<sub>2</sub> concentrations, impacting air-sea CO<sub>2</sub> gas exchange in the Java Sea and changing the balance of the marine carbonate system. This study uses a marine carbonate model to examine the processes that control CO<sub>2</sub> flux in the Java Sea and their relationship to CO<sub>2</sub> increase in the atmosphere. OCMIP-2 (<em>Ocean Carbon-Cycle Model Intercomparison Model Project, Phase-2</em>) is performed in this marine carbonate model coupled with the marine ecosystem model. The model results show that the quantity of carbon air flux differs during February and October 2000. More considerable flux is produced during February 2000, where the wind speeds are higher than in October 2000. However, the wind speeds have less impact when the CO<sub>2</sub> level in the atmosphere rises significantly. Due to the influence of a relatively high surface temperature in the tropical Java sea, the Java Sea functions as a carbon source to the atmosphere in general. In this case, the role of the <em>solubility pump</em> is more significant than that of biological processes in carbon absorption. Moreover, increased CO<sub>2</sub> in the atmosphere could alter the partial pressure equilibrium. In the case of 2002 forest fires (atmospheric CO<sub>2</sub> = 460 ppm), the carbon source of the Java Sea was less than before forest fires and even became carbon sink when atmospheric CO<sub>2</sub> rose to 1135.2 ppm based on the highest SSP scenario in 2100. This modeling also reveals marine acidification issues and could rapidly assess the future changes in marine ecosystems due to CO<sub>2</sub> levels rising in the atmosphere.</p>

Incorporating the stable carbon isotope <sup>13</sup>C in the ocean biogeochemical component of the Max Planck Institute Earth System Model

Abstract. The stable carbon isotopic composition (δ13C) is an important variable to study the ocean carbon cycle across different timescales. We include a new representation of the stable carbon isotope 13C into the HAMburg Ocean Carbon Cycle model (HAMOCC), the ocean biogeochemical component of the Max Planck Institute Earth System Model (MPI-ESM). 13C is explicitly resolved for all oceanic carbon pools considered. We account for fractionation during air–sea gas exchange and for biological fractionation ϵp associated with photosynthetic carbon fixation during phytoplankton growth. We examine two ϵp parameterisations of different complexity: ϵpPopp varies with surface dissolved CO2 concentration (Popp et al., 1989), while ϵpLaws additionally depends on local phytoplankton growth rates (Laws et al., 1995). When compared to observations of δ13C of dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to δ13C in particulate organic carbon (POC) ϵpPopp shows a considerably improved performance compared to ϵpLaws. This is because ϵpLaws produces too strong a preference for 12C, resulting in δ13CPOC that is too low in our model. The model also well reproduces the global oceanic anthropogenic CO2 sink and the oceanic 13C Suess effect, i.e. the intrusion and distribution of the isotopically light anthropogenic CO2 in the ocean. The satisfactory model performance of the present-day oceanic δ13C distribution using ϵpPopp and of the anthropogenic CO2 uptake allows us to further investigate the potential sources of uncertainty of the Eide et al. (2017a) approach for estimating the oceanic 13C Suess effect. Eide et al. (2017a) derived the first global oceanic 13C Suess effect estimate based on observations. They have noted a potential underestimation, but their approach does not provide any insight about the cause. By applying the Eide et al. (2017a) approach to the model data we are able to investigate in detail potential sources of underestimation of the 13C Suess effect. Based on our model we find underestimations of the 13C Suess effect at 200 m by 0.24 ‰ in the Indian Ocean, 0.21 ‰ in the North Pacific, 0.26 ‰ in the South Pacific, 0.1 ‰ in the North Atlantic and 0.14 ‰ in the South Atlantic. We attribute the major sources of underestimation to two assumptions in the Eide et al. (2017a) approach: the spatially uniform preformed component of δ13CDIC in year 1940 and the neglect of processes that are not directly linked to the oceanic uptake and transport of chlorofluorocarbon-12 (CFC-12) such as the decrease in δ13CPOC over the industrial period. The new 13C module in the ocean biogeochemical component of MPI-ESM shows satisfying performance. It is a useful tool to study the ocean carbon sink under the anthropogenic influences, and it will be applied to investigating variations of ocean carbon cycle in the past.

The Sensitivity of the Marine Carbonate System to Regional Ocean Alkalinity Enhancement

Ocean Alkalinity Enhancement (OAE) simultaneously mitigates atmospheric concentrations of CO2 and ocean acidification; however, no previous studies have investigated the response of the non-linear marine carbonate system sensitivity to alkalinity enhancement on regional scales. We hypothesise that regional implementations of OAE can sequester more atmospheric CO2 than a global implementation. To address this, we investigate physical regimes and alkalinity sensitivity as drivers of the carbon-uptake potential response to global and different regional simulations of OAE. In this idealised ocean-only set-up, total alkalinity is enhanced at a rate of 0.25 Pmol a-1 in 75-year simulations using the Max Planck Institute Ocean Model coupled to the HAMburg Ocean Carbon Cycle model with pre-industrial atmospheric forcing. Alkalinity is enhanced globally and in eight regions: the Subpolar and Subtropical Atlantic and Pacific gyres, the Indian Ocean and the Southern Ocean. This study reveals that regional alkalinity enhancement has the capacity to exceed carbon uptake by global OAE. We find that 82–175 Pg more carbon is sequestered into the ocean when alkalinity is enhanced regionally and 156 PgC when enhanced globally, compared with the background-state. The Southern Ocean application is most efficient, sequestering 12% more carbon than the Global experiment despite OAE being applied across a surface area 40 times smaller. For the first time, we find that different carbon-uptake potentials are driven by the surface pattern of total alkalinity redistributed by physical regimes across areas of different carbon-uptake efficiencies. We also show that, while the marine carbonate system becomes less sensitive to alkalinity enhancement in all experiments globally, regional responses to enhanced alkalinity vary depending upon the background concentrations of dissolved inorganic carbon and total alkalinity. Furthermore, the Subpolar North Atlantic displays a previously unexpected alkalinity sensitivity increase in response to high total alkalinity concentrations.

Ocean carbon uptake under aggressive emission mitigation

Abstract. Nearly every nation has signed the UNFCC Paris Agreement, committing to mitigate anthropogenic carbon emissions so as to limit the global mean temperature increase above pre-industrial levels to well below 2 ∘C, and ideally to no more than 1.5 ∘C. A consequence of emission mitigation that has received limited attention is a reduced efficiency of the ocean carbon sink. Historically, the roughly exponential increase in atmospheric CO2 has resulted in a proportional increase in anthropogenic carbon uptake by the ocean. We define growth of the ocean carbon sink exactly proportional to the atmospheric growth rate to be 100 % efficient. Using a model hierarchy consisting of a common reduced-form ocean carbon cycle model and the Community Earth System Model (CESM), we assess the mechanisms of future change in the efficiency of the ocean carbon sink under three emission scenarios: aggressive mitigation (1.5 ∘C), intermediate mitigation (RCP4.5), and high emissions (RCP8.5). The reduced-form ocean carbon cycle model is tuned to emulate the global-mean behavior of the CESM and then allows for mechanistic decomposition. With intermediate or no mitigation (RCP4.5, RCP8.5), changes in efficiency through 2080 are almost entirely the result of future reductions in the carbonate buffer capacity of the ocean. Under the 1.5 ∘C scenario, the dominant driver of efficiency decline is the ocean's reduced ability to transport anthropogenic carbon from surface to depth. As the global-mean upper-ocean gradient of anthropogenic carbon reverses sign, carbon can be re-entrained in surface waters where it slows further removal from the atmosphere. Reducing uncertainty in ocean circulation is critical to better understanding the transport of anthropogenic carbon from surface to depth and to improving quantification of its role in the future ocean carbon sink.

Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models

Abstract. Results from the fully and biogeochemically coupled simulations in which CO2 increases at a rate of 1 % yr−1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm−1) and ocean (0.79 ± 0.07 PgC ppm−1), while the carbon–climate feedback over land (−45.1 ± 50.6 PgC ∘C−1) is about 3 times larger than over ocean (−17.2 ± 5.0 PgC ∘C−1). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 ∘C EgC−1 and is similar to that found in CMIP5 models (1.63 ± 0.48 ∘C EgC−1) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

Carbon Export Buffering and CO2 Drawdown by Flexible Phytoplankton C:N:P Under Glacial Conditions

AbstractModern observations indicate that variations in marine phytoplankton stoichiometry correlate with the boundaries of major surface waters. For example, phytoplankton in the oligotrophic subtropical gyres typically have much higher C:N:P ratios (i.e., higher C:P and higher N:P ratios) than those in eutrophic upwelling regions and polar regions. Such a spatial pattern points to nutrient availability as a key environmental driver of stochiometric flexibility. Environmental dependence of phytoplankton C:N:P opens unexplored possibilities for modifying the strength of the biological pump under different climate conditions. Here we present a power law formulation of C:N:P flexibility that is driven by nutrients, temperature, and light. We embed the formulation in a global ocean carbon cycle model with multiple phytoplankton types and explore biogeochemical implications under glacial conditions. We find three key controls on export C:N:P ratio: phytoplankton physiology and community structure as well as the balance in regional production at the global level. Glacial inputs of iron and sea ice expansion are important modifiers of these three controls. We also find that global export C:N:P increases substantially under glacial conditions, and this strongly buffers global carbon export against decrease and draws down approximately 20 μatm of atmospheric CO2. These results point to the importance of including phytoplankton C:N:P flexibility in a mix of mechanisms that drive atmospheric CO2 over glacial‐interglacial time scale. Finally, our simulations indicate decoupling of nutrients, which may provide a resolution to the longstanding disagreement regarding nutrient utilization in the glacial Southern Ocean derived from different nutrient proxies.

Ocean carbon pump decomposition and its application to CMIP5 earth system model simulations

The ocean stores 60 times as much carbon as the atmosphere, and the ocean carbon cycle has a critical role in controlling atmospheric CO2 concentration. The concept of ocean carbon pump is widely used for describing the ocean carbon cycle, but the term “ocean carbon pump” is not necessarily strictly defined and has been differently referred in previous studies. Here, using three dimensional distributions of dissolved inorganic carbon concentration, total alkalinity, phosphate, and salinity, four types of ocean carbon pump (organic matter, calcium carbonate, gas exchange, and freshwater flux pumps) were formulated. Although previously proposed decomposition emphasizes the enrichment in the deep ocean, my decomposition focuses on surface depletion which directly affects air-sea CO2 exchanges. Based on this decomposition, vector diagram for quantifying the individual roles of the pumps in controlling the ocean surface pCO2, which is in balance with atmospheric CO2 concentration, was demonstrated in this study. The method was applied to the observational climatology, and the contributions of the four carbon pump components to atmospheric CO2 were visualized in a single figure (the vector diagram); each carbon pump component was represented by one vector, and its contribution to CO2 concentration was measured from the difference in the contour values between the beginning and end of the vector. The same analysis was also applied to the CMIP5 earth system model simulations. All the models reproduced the same level of atmospheric CO2 concentration as the observation; however, the contributions from the four carbon pumps varied. The vector diagram was shown to quantify the differences in the contributions from the pumps between the models and against the observation. This study demonstrated that the proposed vector diagram analysis is a useful tool for quantifying the individual contributions of the ocean carbon pumps to atmospheric CO2 concentration and is helpful for evaluating the reproducibility of ocean carbon cycle models.

Microstructure and composition of marine aggregates as co-determinants for vertical particulate organic carbon transfer in the global ocean

Abstract. Marine aggregates are the vector for biogenically bound carbon and nutrients from the euphotic zone to the interior of the oceans. To improve the representation of this biological carbon pump in the global biogeochemical HAMburg Ocean Carbon Cycle (HAMOCC) model, we implemented a novel Microstructure, Multiscale, Mechanistic, Marine Aggregates in the Global Ocean (M4AGO) sinking scheme. M4AGO explicitly represents the size, microstructure, heterogeneous composition, density and porosity of aggregates and ties ballasting mineral and particulate organic carbon (POC) fluxes together. Additionally, we incorporated temperature-dependent remineralization of POC. We compare M4AGO with the standard HAMOCC version, where POC fluxes follow a Martin curve approach with (i) linearly increasing sinking velocity with depth and (ii) temperature-independent remineralization. Minerals descend separately with a constant speed. In contrast to the standard HAMOCC, M4AGO reproduces the latitudinal pattern of POC transfer efficiency, as recently constrained by Weber et al. (2016). High latitudes show transfer efficiencies of ≈0.25±0.04, and the subtropical gyres show lower values of about 0.10±0.03. In addition to temperature as a driving factor for remineralization, diatom frustule size co-determines POC fluxes in silicifier-dominated ocean regions, while calcium carbonate enhances the aggregate excess density and thus sinking velocity in subtropical gyres. Prescribing rising carbon dioxide (CO2) concentrations in stand-alone runs (without climate feedback), M4AGO alters the regional ocean atmosphere CO2 fluxes compared to the standard model. M4AGO exhibits higher CO2 uptake in the Southern Ocean compared to the standard run, while in subtropical gyres, less CO2 is taken up. Overall, the global oceanic CO2 uptake remains the same. With the explicit representation of measurable aggregate properties, M4AGO can serve as a test bed for evaluating the impact of aggregate-associated processes on global biogeochemical cycles and, in particular, on the biological carbon pump.

Enhanced Organic Carbon Burial in Sediments of Oxygen Minimum Zones Upon Ocean Deoxygenation

Oxygen minimum zones (OMZs) in the ocean are expanding. This expansion is attributed to global warming and may continue over the next 10 to 100 kyrs due to multiple climate CO$_\text{2}$-driven factors. The expansion of oxygen-deficient waters has the potential to enhance organic carbon burial in marine sediments, thereby providing a negative feedback on global warming. Here, we study the response of dissolved oxygen in the ocean to increased phosphorus and iron inputs due to CO$_\text{2}$-driven enhanced weathering and increased dust emissions, respectively. We use an ocean biogeochemical model coupled to a general ocean circulation model (the Hamburg Oceanic Carbon Cycle model, HAMOCC 2.0) to assess the impact of such regional deoxygenation on organic carbon burial in the modern ocean on time scales of up to 200~kyrs. We find that an increase in input of phosphorus and iron leads to an expansion of the area of the OMZ impinging on continental margin sediments and a significant decline in bottom water oxygen in the open ocean relative to pre-industrial conditions. The associated increase in organic carbon burial could contribute to the drawdown of $\sim$1600~Gt of carbon, which is equivalent to the total amount of CO$_\text{2}$ in the atmosphere predicted for the year 2100 in a business as usual scenario, on time scales of up to 50 kyrs. The corresponding areal extent of sediments overlain by bottom waters with little or no oxygen as estimated by the model is not very different from the minimum area estimated for two major oceanic anoxic events in Earth's past. Such events were associated with major perturbations of the oceanic carbon cycle, including high rates of organic carbon burial. We conclude that organic carbon burial in low oxygen areas in the ocean could contribute to removal of anthropogenic CO$_\text{2}$ from the atmosphere on long time scales.

The Canadian Earth System Model version 5 (CanESM5.0.3)

Abstract. The Canadian Earth System Model version 5 (CanESM5) is a global model developed to simulate historical climate change and variability, to make centennial-scale projections of future climate, and to produce initialized seasonal and decadal predictions. This paper describes the model components and their coupling, as well as various aspects of model development, including tuning, optimization, and a reproducibility strategy. We also document the stability of the model using a long control simulation, quantify the model's ability to reproduce large-scale features of the historical climate, and evaluate the response of the model to external forcing. CanESM5 is comprised of three-dimensional atmosphere (T63 spectral resolution equivalent roughly to 2.8∘) and ocean (nominally 1∘) general circulation models, a sea-ice model, a land surface scheme, and explicit land and ocean carbon cycle models. The model features relatively coarse resolution and high throughput, which facilitates the production of large ensembles. CanESM5 has a notably higher equilibrium climate sensitivity (5.6 K) than its predecessor, CanESM2 (3.7 K), which we briefly discuss, along with simulated changes over the historical period. CanESM5 simulations contribute to the Coupled Model Intercomparison Project phase 6 (CMIP6) and will be employed for climate science and service applications in Canada.

CO<sub>2</sub> drawdown due to particle ballasting by glacial aeolian dust: an estimate based on the ocean carbon cycle model MPIOM/HAMOCC version 1.6.2p3

Abstract. Despite intense efforts, the mechanisms that drive glacial–interglacial changes in atmospheric pCO2 are not fully understood. Here, we aim at quantifying the potential contribution of aeolian dust deposition changes to the atmospheric pCO2 drawdown during the Last Glacial Maximum (LGM). To this end, we use the Max Planck Institute Ocean Model (MPIOM) and the embedded Hamburg Ocean Carbon Cycle model (HAMOCC), including a new parameterization of particle ballasting that accounts for the acceleration of sinking organic soft tissue in the ocean by higher-density biogenic calcite and opal particles, as well as mineral dust. Sensitivity experiments with reconstructed LGM dust deposition rates indicate that the acceleration of detritus by mineral dust played a small role in atmospheric pCO2 variations during glacial–interglacial cycles – on the order of 5 ppmv, compared to the reconstructed ∼80 ppmv rise in atmospheric pCO2 during the last deglaciation. The additional effect of the LGM dust deposition, namely the enhanced fertilization by the iron that is associated with the glacial dust, likely played a more important role; although the full iron fertilization effect can not be estimated in the particular model version used here due to underestimated present-day non-diazotroph iron limitation, fertilization of diazotrophs in the tropical Pacific already leads to an atmospheric pCO2 drawdown of around 10 ppmv.

Changes to the Air‐Sea Flux and Distribution of Radiocarbon in the Ocean Over the 21st Century

AbstractWe investigate the spatiotemporal evolution of radiocarbon (Δ14C) in the ocean over the 21st century under different scenarios for anthropogenic CO2 emissions and atmospheric CO2 and radiocarbon changes using a 3‐D ocean carbon cycle model. Strong decreases in atmospheric Δ14C in the high‐emission scenario result in strong outgassing of 14C over 2050–2100, causing Δ14C spatial gradients in the surface ocean and vertical gradients between the surface and intermediate waters to reverse sign. Surface Δ14C in the subtropical gyres is lower than Δ14C in Pacific Deep Water and Southern Ocean surface water in 2100. In the low‐emission scenario, ocean Δ14C remains slightly higher than in 1950 and relatively constant over 2050–2100. Over the next 20 years we find decadal changes in Δ14C of −30‰ to +5‰ in the upper 2 km of the ocean, which should be detectable with continued hydrographic surveys. Our simulations can help in planning future observations, and they provide a baseline for investigating natural or anthropogenic changes in ocean circulation using ocean Δ14C observations and models.

Biological production in the Indian Ocean upwelling zones –Part 1: refined estimation via the use of a variable compensation depth in ocean carbon models

Abstract. Biological modelling approach adopted by the Ocean Carbon-Cycle Model Intercomparison Project (OCMIP-II) provided amazingly simple but surprisingly accurate rendition of the annual mean carbon cycle for the global ocean. Nonetheless, OCMIP models are known to have seasonal biases which are typically attributed to their bulk parameterisation of compensation depth. Utilising the criteria of surface Chl a-based attenuation of solar radiation and the minimum solar radiation required for production, we have proposed a new parameterisation for a spatially and temporally varying compensation depth which captures the seasonality in the production zone reasonably well. This new parameterisation is shown to improve the seasonality of CO2 fluxes, surface ocean pCO2, biological export and new production in the major upwelling zones of the Indian Ocean. The seasonally varying compensation depth enriches the nutrient concentration in the upper ocean yielding more faithful biological exports which in turn leads to accurate seasonality in the carbon cycle. The export production strengthens by ∼ 70 % over the western Arabian Sea during the monsoon period and achieves a good balance between export and new production in the model. This underscores the importance of having a seasonal balance in the model export and new productions for a better representation of the seasonality of the carbon cycle over upwelling regions. The study also implies that both the biological and solubility pumps play an important role in the Indian Ocean upwelling zones.

Large Reemergence of Anthropogenic Carbon into the Ocean’s Surface Mixed Layer Sustained by the Ocean’s Overturning Circulation

We evaluate the output from a widely used ocean carbon cycle model to identify the subduction and obduction (reemergence) rates of anthropogenic carbon (Cant) for climatological conditions during the World Ocean Circulation Experiment (WOCE) era in 1995 using a new set of Lagrangian diagnostic tools. The principal scientific value of the Lagrangian diagnostics is in providing a new means to connect Cant reemergence pathways to the relatively rapid renewal time scales of mode waters through the overturning circulation. Our main finding is that for this model with 2.04 PgC yr−1 of uptake of Cant via gas exchange, the subduction and obduction rates across the base of the mixed layer (MLbase) are 4.96 and 4.50 PgC yr−1, respectively, which are twice as large as the gas exchange at the surface. Given that there is net accumulation of 0.17 PgC yr−1 in the mixed layer itself, this implies the residual downward Cant transport of 1.40 PgC yr−1 across the MLbase is associated with diffusion. Importantly, the net patterns for subduction and obduction transports of Cant mirror the large-scale patterns for transport of water volume, thereby illustrating the processes controlling Cant uptake. Although the net transfer across the MLbase by compensating subduction and obduction is relatively smaller than the diffusion, the localized pattern of Cant subduction and obduction implies significant regional impacts. The median time scale for reemergence of obducting particles is short (<10 yr), indicating that reemergence should contribute to limiting future carbon uptake through its contribution to perturbing the Revelle factor for surface waters.

Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

Abstract. 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, SF6) 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 14CTabio, 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 CO2 as well as carbon isotope ratios of CO2. 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.

Response of export production and dissolved oxygen concentrations in oxygen minimum zones to <i>p</i>CO<sub>2</sub> and temperature stabilization scenarios in the biogeochemical model HAMOCC 2.0

Abstract. Dissolved oxygen (DO) concentration in the ocean is an important component of marine biogeochemical cycles and will be greatly altered as climate change persists. In this study a global oceanic carbon cycle model (HAMOCC 2.0) is used to address how mechanisms of oxygen minimum zone (OMZ) expansion respond to changes in CO2 radiative forcing. Atmospheric pCO2 is increased at a rate of 1 % annually and the model is stabilized at 2 ×, 4 ×, 6 ×, and 8 × preindustrial pCO2 levels. With an increase in CO2 radiative forcing, the OMZ in the Pacific Ocean is controlled largely by changes in particulate organic carbon (POC) export, resulting in increased remineralization and thus expanding the OMZs within the tropical Pacific Ocean. A potential decline in primary producers in the future as a result of environmental stress due to ocean warming and acidification could lead to a substantial reduction in POC export production, vertical POC flux, and thus increased DO concentration particularly in the Pacific Ocean at a depth of 600–800 m. In contrast, the vertical expansion of the OMZs within the Atlantic is linked to increases POC flux as well as changes in oxygen solubility with increasing seawater temperature. Changes in total organic carbon and increase sea surface temperature (SST) also lead to the formation of a new OMZ in the western subtropical Pacific Ocean. The development of the new OMZ results in dissolved oxygen concentration of ≤ 50 µmol kg−1 throughout the equatorial Pacific Ocean at 4 times preindustrial pCO2. Total ocean volume with dissolved oxygen concentrations of ≤ 50 µmol kg−1 increases by 2.4, 5.0, and 10.5 % for the 2 ×, 4 ×, and 8 × CO2 simulations, respectively.

Incorporating a prognostic representation of marine nitrogen fixers into the global ocean biogeochemical model HAMOCC

AbstractNitrogen (N2) fixation is a major source of bioavailable nitrogen to the euphotic zone, thereby exerting an important control on ocean biogeochemical cycling. This paper presents the incorporation of prognostic N2 fixers into the HAMburg Ocean Carbon Cycle model (HAMOCC), a component of the Max Planck Institute Earth System Model (MPI‐ESM). Growth dynamics of N2 fixers in the model are based on physiological characteristics of the cyanobacterium Trichodesmium. The applied temperature dependency confines diazotrophic growth and N2 fixation to the tropical and subtropical ocean roughly between 40°S and 40°N. Simulated large‐scale spatial patterns compare well with observations, and the global N2 fixation rate of 135.6 Tg N yr−1 is within the range of current estimates. The vertical distribution of N2 fixation also matches well the observations, with a major fraction of about 85% occurring in the upper 20 m. The observed seasonal variability at the stations BATS and ALOHA is reasonably reproduced, with highest fixation rates in northern summer/fall. Iron limitation was found to be an important factor in controlling the simulated distribution of N2 fixation, especially in the Pacific Ocean. The new model component considerably improves the representation of present‐day N2 fixation in HAMOCC. It provides the basis for further studies on the role of diazotrophs in global biogeochemical cycles, as well as on the response of N2 fixation to changing environmental conditions.

Fluks CO2 di Perairan Pesisir Pulau Lombok, Nusa Tenggara Barat

<strong>CO<sub>2</sub> Flux in the Coastal Waters of Lombok, West Nusa Tenggara.</strong> Coastal and marine areas play an important role in climate regulation by absorbing CO<sub>2</sub> and exchanging carbon with various carbon storing compartments as the atmosphere, land and biota. This study aims to determine the exchange of CO2 between the sea surface to the atmosphere in the coastal waters of the island of Lombok and its components sink and source, and to know the parameters that influence the process. Field observation was carried out during April, 2012 in Sekotong Bay (10 Stations) and Kodek Bay (6 Stations).The partial pressure of CO<sub>2</sub> (pCO<sub>2</sub>) in sea surface level, was calculated using ABIOTIC model of the ocean carbon cycle model intercomparison project phase-2 (OCMIP-2). The results showed that generally Kodek Bay waters and outside of Sekotong Bay act as a sink of CO2 from the atmosphere, while inside of Sekotong Bay act as source of CO<sub>2</sub> to the atmosphere.CO<sub>2</sub> flux in Sekotong bay waters and Kodek bay waters each ranging between -0.61 to 0.52 and -0.80 to -1.84 mmol.m<sup>-2</sup>.d<sup>-1</sup> with an average of 0.05 ± 0.40 and -1.29 ± 0.40 mmol.m<sup>-2</sup>.d<sup>-1</sup>. Kodek Bay act as CO<sub>2</sub> sink because influence by water mass from Lombok Strait that have high concentrations of nutrient and chlorophyll. The existence of seagrass and coral reefs that support small islands in the Kodek bay helped reduce elevated pCO<sub>2</sub> water column and strengthen the absorption of CO<sub>2</sub> in these waters. Sekotong bay inside, acts as a release of CO<sub>2</sub>, because the waters of the Sekotong bay tend to be closed, so that the residence time of water masses becomes longer resulting in enrichment of DIC and pCO<sub>2</sub>.

Evaluation of NorESM-OC (versions 1 and 1.2), the ocean carbon-cycle stand-alone configuration of the Norwegian Earth System Model (NorESM1)

Abstract. 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 three different model configurations (two different model versions at different grid resolutions) using two different atmospheric forcing data sets. Model version NorESM-OC1 corresponds to the version that is included in the NorESM-ME1 fully coupled model, 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 that uses a linear increase in the 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 at 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.