Hamburg Ocean Carbon Cycle Research Articles

The ICON Earth System Model Version 1.0

AbstractThis work documents the ICON‐Earth System Model (ICON‐ESM V1.0), the first coupled model based on the ICON (ICOsahedral Non‐hydrostatic) framework with its unstructured, icosahedral grid concept. The ICON‐A atmosphere uses a nonhydrostatic dynamical core and the ocean model ICON‐O builds on the same ICON infrastructure, but applies the Boussinesq and hydrostatic approximation and includes a sea‐ice model. The ICON‐Land module provides a new framework for the modeling of land processes and the terrestrial carbon cycle. The oceanic carbon cycle and biogeochemistry are represented by the Hamburg Ocean Carbon Cycle module. We describe the tuning and spin‐up of a base‐line version at a resolution typical for models participating in the Coupled Model Intercomparison Project (CMIP). The performance of ICON‐ESM is assessed by means of a set of standard CMIP6 simulations. Achievements are well‐balanced top‐of‐atmosphere radiation, stable key climate quantities in the control simulation, and a good representation of the historical surface temperature evolution. The model has overall biases, which are comparable to those of other CMIP models, but ICON‐ESM performs less well than its predecessor, the Max Planck Institute Earth System Model. Problematic biases are diagnosed in ICON‐ESM in the vertical cloud distribution and the mean zonal wind field. In the ocean, sub‐surface temperature and salinity biases are of concern as is a too strong seasonal cycle of the sea‐ice cover in both hemispheres. ICON‐ESM V1.0 serves as a basis for further developments that will take advantage of ICON‐specific properties such as spatially varying resolution, and configurations at very high resolution.

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.

Reconstructing the Preindustrial Coastal Carbon Cycle Through a Global Ocean Circulation Model: Was the Global Continental Shelf Already Both Autotrophic and a CO2 Sink?

AbstractThe contribution of continental shelves to the marine carbon cycle is still poorly understood. Their preindustrial state is, for one, essentially unknown, which strongly limits the quantitative assessment of their anthropogenic perturbation. To date, approaches developed to investigate and quantify carbon fluxes on continental shelves have strongly simplified their physical and biogeochemical features. In this study, we enhance the global ocean biogeochemistry model HAMburg Ocean Carbon Cycle by explicitly representing riverine loads of carbon and nutrients, as well as improving the representation of organic matter dynamics in the coastal ocean. Our simulations, at a resolution of ∼0.4°, reveal a globally averaged shelf water residence time (RT) of 12–17 months, which is much shorter than the global RTs previously assumed in benchmark studies (>4 years). This shorter global RT, induced primarily through outer shelf regions with large oceanic inflows, promotes an efficient offshore transport of terrestrial and marine organic carbon (0.44 PgCyr−1) and a dissolved inorganic carbon sink from the organic cycling of carbon on the global shelf (net ecosystem productivity [NEP] equal to +0.20 PgCyr−1). In turn, this autotrophic state of continental shelves contributes to a weak global preindustrial sink of atmospheric CO2 (0.04 PgCyr−1), dominated by extensive regions with large oceanic inflows and positive NEPs, such as the Patagonian shelf, the East China Sea and the outer North Sea. The contemporary global shelf CO2 uptake of 0.15 PgCyr−1 furthermoresuggests that the anthropogenic CO2 uptake (0.11 PgCyr−1) on the global continental shelf is less efficient with respect to the open ocean.

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.

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.

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.

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.

Ocean biogeochemistry in the warm climate of the late Paleocene

Abstract. The late Paleocene is characterized by warm and stable climatic conditions that served as the background climate for the Paleocene–Eocene Thermal Maximum (PETM, ~55 million years ago). With respect to feedback processes in the carbon cycle, the ocean biogeochemical background state is of major importance for projecting the climatic response to a carbon perturbation related to the PETM. Therefore, we use the Hamburg Ocean Carbon Cycle model (HAMOCC), embedded in the ocean general circulation model of the Max Planck Institute for Meteorology, MPIOM, to constrain the ocean biogeochemistry of the late Paleocene. We focus on the evaluation of modeled spatial and vertical distributions of the ocean carbon cycle parameters in a long-term warm steady-state ocean, based on a 560 ppm CO2 atmosphere. Model results are discussed in the context of available proxy data and simulations of pre-industrial conditions. Our results illustrate that ocean biogeochemistry is shaped by the warm and sluggish ocean state of the late Paleocene. Primary production is slightly reduced in comparison to the present day; it is intensified along the Equator, especially in the Atlantic. This enhances remineralization of organic matter, resulting in strong oxygen minimum zones and CaCO3 dissolution in intermediate waters. We show that an equilibrium CO2 exchange without increasing total alkalinity concentrations above today's values is achieved. However, consistent with the higher atmospheric CO2, the surface ocean pH and the saturation state with respect to CaCO3 are lower than today. Our results indicate that, under such conditions, the surface ocean carbonate chemistry is expected to be more sensitive to a carbon perturbation (i.e., the PETM) due to lower CO32− concentration, whereas the deep ocean calcite sediments would be less vulnerable to dissolution due to the vertically stratified ocean.

Methyl iodide production in the open ocean

Abstract. Production pathways of the prominent volatile organic halogen compound methyl iodide (CH3I) are not fully understood. Based on observations, production of CH3I via photochemical degradation of organic material or via phytoplankton production has been proposed. Additional insights could not be gained from correlations between observed biological and environmental variables or from biogeochemical modeling to identify unambiguously the source of methyl iodide. In this study, we aim to address this question of source mechanisms with a three-dimensional global ocean general circulation model including biogeochemistry (MPIOM–HAMOCC (MPIOM – Max Planck Institute Ocean Model HAMOCC – HAMburg Ocean Carbon Cycle model)) by carrying out a series of sensitivity experiments. The simulated fields are compared with a newly available global data set. Simulated distribution patterns and emissions of CH3I differ largely for the two different production pathways. The evaluation of our model results with observations shows that, on the global scale, observed surface concentrations of CH3I can be best explained by the photochemical production pathway. Our results further emphasize that correlations between CH3I and abiotic or biotic factors do not necessarily provide meaningful insights concerning the source of origin. Overall, we find a net global annual CH3I air–sea flux that ranges between 70 and 260 Gg yr−1. On the global scale, the ocean acts as a net source of methyl iodide for the atmosphere, though in some regions in boreal winter, fluxes are of the opposite direction (from the atmosphere to the ocean).

Global ocean biogeochemistry model HAMOCC: Model architecture and performance as component of the MPI‐Earth system model in different CMIP5 experimental realizations

Ocean biogeochemistry is a novel standard component of fifth phase of the Coupled Model Intercomparison Project (CMIP5) experiments which project future climate change caused by anthropogenic emissions of greenhouse gases. Of particular interest here is the evolution of the oceanic sink of carbon and the oceanic contribution to the climate‐carbon cycle feedback loop. The Hamburg ocean carbon cycle model (HAMOCC), a component of the Max Planck Institute for Meteorology Earth system model (MPI‐ESM), is employed to address these challenges. In this paper we describe the version of HAMOCC used in the CMIP5 experiments (HAMOCC 5.2) and its implementation in the MPI‐ESM to provide a documentation and basis for future CMIP5‐related studies. Modeled present day distributions of biogeochemical variables calculated in two different horizontal resolutions compare fairly well with observations. Statistical metrics indicate that the model performs better at the ocean surface and worse in the ocean interior. There is a tendency for improvements in the higher resolution model configuration in representing deeper ocean variables; however, there is little to no improvement at the ocean surface. An experiment with interactive carbon cycle driven by emissions of CO2 produces a 25% higher variability in the oceanic carbon uptake over the historical period than the same model forced by prescribed atmospheric CO2 concentrations. Furthermore, a climate warming of 3.5 K projected at atmospheric CO2 concentration of four times the preindustrial value, reduced the atmosphere‐ocean CO2 flux by 1 GtC yr−1. Overall, the model shows consistent results in different configurations, being suitable for the type of simulations required within the CMIP5 experimental design.

Miocene ocean circulation inferred from marine carbon cycle modeling combined with benthic isotope records

In a modeling sensitivity study we investigate the evolution of the ocean circulation and of marine carbon isotope (δ13C) records during the Miocene (about 23–5 million years ago). For this purpose we ran an ocean‐circulation carbon cycle model of intermediate complexity (Large Scale Geostrophic– Hamburg Ocean Carbon Cycle Model, version 2s) exploring various seaway configurations. Our investigations confirm that the Central American Seaway played a decisive role in the history of the Miocene ocean circulation. In simulations with a deep Central American Seaway (depth range 1–3 km), typical for the early to middle Miocene, deep water production in the North Atlantic is absent or weak, while the meridional overturning circulation is dominated by water mass formation in the Southern Ocean. Deep water formation in the North Atlantic begins when the Central American Seaway shoals to a few hundreds of meters, which is typical for the late Miocene. Our results do not support ideas that the mid‐Miocene closing of the Eastern Tethys contributed to Antarctic glaciation. On the other hand, we find some water exchange between the Indian Ocean and the Atlantic via the Eastern Tethys during the early Miocene. Our model results for the Atlantic meridional overturning circulation and for Atlantic δ13C during the late Miocene are largely independent from depth variations of the Greenland‐Scotland Ridge. To a large extent, the evolution of Miocene deep‐sea δ13C records can be explained with large‐scale ocean circulation changes. Our model‐data comparison for the middle and early Miocene suggests that during the early Neogene the seaway effect on benthic δ13C may have been superimposed by further factors such as climate regime shifts and/or terrestrial carbon cycle changes.

Sensitivity of ocean carbon tracer distributions to particulate organic flux parameterizations

Vertical fluxes of particulate organic carbon (POC) from the euphotic zone to the deep sea are an important part of the carbon cycle in the oceans. In this study, oceanic fluxes of POC below the euphotic zone were simulated with the Hamburg ocean carbon cycle model (HAMOCC5.1) by using different POC‐flux parameterizations and compared with sediment trap data. Overall, the geochemical distributions in the deep sea showed a high sensitivity to the selection of POC flux parameterization. Below 2000 m in the oceans, differences between simulated and observed carbon tracers (PO4, Alk*) and model‐data differences of POC fluxes are lowest when a regionally variable POC flux parameterization is used. Specifically, model‐data differences are lowest in the subtropics when simulated vertical POC flux considers mineral ballasting, while they are lowest in the Arabian Sea and west coast of Africa when a mineral dust parameterization is used for the vertical POC flux.

Terrestrial mechanisms of interannual CO2variability

The interannual variability of atmospheric CO2growth rate shows remarkable correlation with the El Niño Southern Oscillation (ENSO). Here we present results from mechanistically based terrestrial carbon cycle model VEgetation‐Global‐Atmosphere‐Soil (VEGAS), forced by observed climate fields such as precipitation and temperature. Land is found to explain most of the interannual CO2variability with a magnitude of about 5 PgC yr−1. The simulated land‐atmosphere flux has a detrended correlation of 0.53 (0.6 at the 2–7 year ENSO band) with the CO2growth rate observed at Mauna Loa from 1965 to 2000. We also present the total ocean flux from the Hamburg Ocean Carbon Cycle Model (HAMOCC) which shows ocean‐atmosphere flux variation of about 1 PgC yr−1, and it is largely out of phase with land flux. On land, much of the change comes from the tropical regions such as the Amazon and Indonesia where ENSO related climate anomalies are in the same direction across much of the tropics. The subcontinental variations over North America and Eurasia are comparable to the tropics but the total interannual variability is about 1 PgC yr−1due to the cancellation from the subregions. This has implication for flux measurement network distribution. The tropical dominance also results from a “conspiracy” between climate and plant/soil physiology, as precipitation and temperature changes drive opposite changes in net primary production (NPP) and heterotrophic respiration (Rh), both contributing to land‐atmosphere flux changes in the same direction. However, NPP contributes to about three fourths of the total tropical interannual variation and the rest is from heterotrophic respiration; thus precipitation appears to be a more important factor than temperature on the interannual timescales as tropical wet and dry regimes control vegetation growth. Fire, largely driven by drought, also contributes significantly to the interannual CO2variability at a rate of about 1 PgC yr−1, and it is not totally in phase with NPP or Rh. The robust variability in tropical fluxes agree well with atmospheric inverse modeling results. Even over North America and Eurasia, where ENSO teleconnection is less robust, the fluxes show general agreement with inversion results, an encouraging sign for fruitful carbon data assimilation.

Nutrient trapping in the equatorial Pacific: The ocean circulation solution

Nutrient trapping is a chronic problem found in global carbon cycle models with particle‐only remineralization schemes. It is defined as the excess of subsurface nutrient concentrations relative to observations and occurs principally in the eastern equatorial Pacific. Previous studies reduced excess simulated nutrients by increasing the complexity of modeled biogeochemistry, i.e., by adding pools for nutrients (and carbon) either in dissolved organic form or as plankton. Conversely, our study suggests that deficiencies in modeled circulation fields from global coarse‐resolution ocean models are mostly responsible. This new interpretation stems from our use of an ocean general circulation model with higher resolution, which offers a more realistic equatorial circulation. We used the same biogeochemical model Hamburg ocean carbon cycle model, version 3, as in some of the previous studies. Our model‐predicted distribution of PO43− in the equatorial Pacific agrees reasonably well with the observations both at the surface and in the subsurface. Subsurface PO43− concentrations in our model's eastern equatorial Pacific exceed observations by, at most, 15%, unlike coarser‐resolution models. Improvement is due to enhanced meridional resolution (0.5°) near the equator, which allows the model to simulate a vigorous equatorial undercurrent that brings in low‐nutrient water from the western basin. Furthermore, the model upwells no nutrient‐rich abyssal water into the surface equatorial Pacific. Our results suggest that dissolved organic carbon plays a minor role in the carbon budget of the equatorial Pacific.

Closing off the Southern Ocean surface

Equilibrium integrations with the three‐dimensional Hamburg ocean carbon cycle circulation model show the effect of a Southern Ocean surface outcrop that is closed off for gas exchange with the atmosphere and export production of particulate organic carbon and CaCO3 south of 32.5 °S. For such a close‐off, the northern Atlantic outcrop assumes a greater degree of control of atmospheric CO2 partial pressure. Contrary to expectations, the atmospheric pCO2 increases compared to a control run without close‐off. Maximizing the nutrient utilization in the North Atlantic can compensate for this increase. However, it does not draw down the pCO2 to glacial levels.

Glacial pCO2 Reduction by the World Ocean: Experiments With the Hamburg Carbon Cycle Model

Possible mechanisms for the 80 ppm reduction of atmospheric CO2 partial pressure during the last glaciation were investigated using the Hamburg Ocean Carbon Cycle Model. The three‐dimensional carbon cycle model is based on the velocity field of the Hamburg Large‐Scale Geostrophic Ocean General Circulation Model and uses the same grid as that model. The horizontal resolution (3.5° × 3.5°) is lower than the length scale of narrow upwelling belts which could not be represented adequately in this study, but the large‐scale features of the ocean carbon cycle are reproduced rather well. Sensitivity experiments were carried out to investigate the role of chemical and biological parameters (nutrient cycling, composition of biogenic particulate matter, CO2 solubility) and different circulation regimes for the atmospheric CO2 content. The model responses were compared to deep‐sea sediment core records and ice core data from the last glaciation. Each experiment was compared with observed average tracer patterns during 18–65 kyr B.P. It was found that none of the sensitivity experiments alone could explain all observed tracer changes (atmospheric pCO2, Δδ13Cplanktonic‐benthic, δ13Cbenthic differences, CaCO3 corrosivity indices) simultaneously, even in a qualitative sense. Thus according to the model none of the scenarios tested proves to be completely acceptable. The residual discrepancies between the observed and modeled tracer records can probably be attributed to the as yet inadequate reconstruction of the glacial ocean circulation. It is therefore suggested that more effort should be devoted to realistically reproducing the ice age ocean circulation field making use of the forthcoming glacial radiocarbon data base. The residuals between the realistically modeled and observed carbon cycle tracers (δ13C, CaCO3 saturation) should then reveal more clearly the real cause for the observed pCO2 decrease in the glacial atmosphere.