Incorporating the stable carbon isotope ¹³C in the ocean biogeochemical component of the Max Planck Institute Earth System Model

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.

Incorporating the stable carbon isotope ¹³C in the ocean biogeochemical component of the Max Planck Institute Earth System Model

Direct comparison between paleo oceanic δ13C records and model results facilitates assessing simulated distributions and properties of water masses in the past. To accomplish this, 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 existing oceanic carbon pools. 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 in dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to δ13C in particulate organic carbon εpPopp shows a considerably improved performance than εpLaws, because the latter results in a too strong preference for 12C. The model also well reproduces the oceanic 13C Suess effect, i.e. the intrusion of the isotopically light anthropogenic CO2 into the ocean, based on comparison to other existing 13C models and to observation-based oceanic carbon uptake estimates over the industrial period. We further apply the approach of Eide et al. (2017a), who derived the first global oceanic 13C Suess effect estimate based on observations, to our model data that has ample spatial and temporal coverage. With this we are able to analyse in detail the underestimation of 13C Suess effect by this approach as it has been noted by Eide et al. (2017a). Based on our model we find underestimations of 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 the underestimation to two assumptions in Eide et al. (2017a)'s approach: a spatially-constant preformed component of δ13CDIC in year 1940 and neglecting 13C Suess effect in CFC-12 free water.

Stable isotopes of amino acids from reef fishes uncover Suess and nitrogen enrichment effects on local ecosystems

In 1979, the Suess effect was described as decreasing δ13C in the oceans linked to anthropogenic CO2 emissions. After years of over-fertilization of farming soils and runoff, we hypothesized that δ15N in coastal environments would also decline, whereby synthetic fertilizers lead to depletion of the heavy isotope 15N. We used museum-preserved and modern samples of 3 fishes from Otago, New Zealand, to reconstruct the isotopic baselines of C and N and assess specific trophic positions through time (1955-present) based on bulk and amino acid stable isotope values. Our sample set included Odax pullus, a strictly herbivorous species, and 2 commercially important species: Nemadactylus macropterus and Parapercis colias. Muscle tissue of the fishes recorded the change in δ13CBulk through time, which matched estimated Suess effect values for New Zealand. We also resolved the effects on the C isotopic baseline from natural changes in the food web using analysis of the δ13C of essential amino acids and found that while P. colias maintained a steady diet, the food web position of N. macropterus likely changed. Analysis of δ15N of phenylalanine in O. pullus indicated a decrease of 0.023‰ yr-1 since 1955, which corroborates our coastal N-enrichment hypothesis. Furthermore, we found that isotopic changes for N. macropterus were consistent with overfishing and habitat degradation in the region. These data provide vital information for our resolution and understanding of how past environments have changed in terms of both anthropogenic influences on coastal food web structure and biogeochemical cycles of C and N in marine ecosystems.

A ~1000-year 13C Suess correction model for the study of past ecosystems

Inferences about how an ecosystem has changed through time often rely on longitudinal records of species characteristics or niche parameters, and stable isotope analysis is a common tool employed to study changes in an organism’s niche. One of the most frequently used stable isotope measures is δ13C, a ratio of 13C to 12C. However, applying δ13C to historical samples comes with some methodological hurdles. One such hurdle is correcting for the 13C Suess effect or the change in atmospheric δ13C due to increased anthropogenic CO2 emissions. The change in the amount of carbon isotopes in the atmosphere through time can confound the study of historical shifts in species characteristics. No standard way of correcting for the 13C Suess effect has been suggested despite this problem. Here, I propose a standard 13C Suess correction model for the past ~1000 years using three prehistoric/historic records of atmospheric δ13C.