A First Intercomparison of the Simulated LGM Carbon Results Within PMIP-Carbon: Role of the Ocean Boundary Conditions

F. Lhardy*, N. Bouttes, D. M. Roche, A. Abe-Ouchi, Z. Chase, K. A. Crichton, T. Ilyina, R. Ivanovic, M. Jochum, M. Kageyama, H. Kobayashi, B. Liu, L. Menviel, J. Muglia, R. Nuterman, A. Oka, G. Vettoretti, A. Yamamoto

*Corresponding author for this work

Research output: Contribution to JournalArticleAcademicpeer-review

72 Downloads (Pure)

Abstract

Model intercomparison studies of coupled carbon-climate simulations have the potential to improve our understanding of the processes explaining the (Formula presented.) drawdown at the Last Glacial Maximum (LGM) and to identify related model biases. Models participating in the Paleoclimate Modeling Intercomparison Project (PMIP) now frequently include the carbon cycle. The ongoing PMIP-carbon project provides the first opportunity to conduct multimodel comparisons of simulated carbon content for the LGM time window. However, such a study remains challenging due to differing implementation of ocean boundary conditions (e.g., bathymetry and coastlines reflecting the low sea level) and to various associated adjustments of biogeochemical variables (i.e., alkalinity, nutrients, dissolved inorganic carbon). After assessing the ocean volume of PMIP models at the pre-industrial and LGM, we investigate the impact of these modeling choices on the simulated carbon at the global scale, using both PMIP-carbon model outputs and sensitivity tests with the iLOVECLIM model. We show that the carbon distribution in reservoirs is significantly affected by the choice of ocean boundary conditions in iLOVECLIM. In particular, our simulations demonstrate a (Formula presented.) GtC effect of an alkalinity adjustment on carbon sequestration in the ocean. Finally, we observe that PMIP-carbon models with a freely evolving (Formula presented.) and no additional glacial mechanisms do not simulate the (Formula presented.) drawdown at the LGM (with concentrations as high as 313, 331, and 315 ppm), especially if they use a low ocean volume. Our findings suggest that great care should be taken on accounting for large bathymetry changes in models including the carbon cycle.

Original languageEnglish
Article numbere2021PA004302
Pages (from-to)1-15
Number of pages15
JournalPaleoceanography and Paleoclimatology
Volume36
Issue number10
Early online date17 Sept 2021
DOIs
Publication statusPublished - Oct 2021

Bibliographical note

Funding Information:
F. Lhardy, N. Bouttes, and D. M. Roche designed the research. N. Bouttes coordinated the PMIP‐carbon project and obtained funding. Participating modeling groups all performed a PI and a LGM simulation, provided their model outputs and the relevant metadata and computed the equilibrated carbon content in reservoirs. These modeling groups included AA‐O, HK, and AO (MIROC4m‐COCO); KC (CLIMBER‐2); MJ, RN, GV, and ZC (CESM); BL and TI (MPI‐ESM); MK (IPSL‐CM5A2); AY (MIROC‐ES2L); LM (LOVECLIM); JM and AS (UVic). F. Lhardy, D. M. Roche, and N. Bouttes generated new boundary conditions in the iLOVECLIM model. N. Bouttes and D. M. Roche developed the automated adjustments to allow for a change of ocean boundary conditions. F. Lhardy ran the iLOVECLIM simulations and analyzed both the iLOVECLIM and the PMIP‐carbon outputs under supervision of N. Bouttes and D. M. Roche. F. Lhardy wrote the manuscript with the inputs from all co‐authors. This study was supported by the French National program LEFE (Les Enveloppes Fluides et l'Environnement). F. Lhardy is supported by the Université Versailles Saint‐Quentin‐en‐Yvelines (UVSQ). N. Bouttes and D. M. Roche are supported by the Centre national de la recherche scientifique (CNRS). In addition, D. M. Roche is supported by the Vrije Universiteit Amsterdam. B. Liu and T. Ilyina are supported by the German Federal Ministry of Education and Research (PalMod initiative, FKZ: 420 grant no. 01LP1919B). L. Menviel acknowledges funding from the Australian Research Council grant FT180100606. A. Yamamoto acknowledges funding from the Integrated Research Program for Advancing Climate Models (TOUGOU) Grant Number JPMXD0717935715 from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. The authors acknowledge the use of the LSCE storage and computing facilities. The authors thank Théo Mandonnet for his preliminary work on the PMIP‐carbon project. Last but not least, the authors thank the two anonymous reviewers for their help with this manuscript.

Funding Information:
F. Lhardy, N. Bouttes, and D. M. Roche designed the research. N. Bouttes coordinated the PMIP-carbon project and obtained funding. Participating modeling groups all performed a PI and a LGM simulation, provided their model outputs and the relevant metadata and computed the equilibrated carbon content in reservoirs. These modeling groups included AA-O, HK, and AO (MIROC4m-COCO); KC (CLIMBER-2); MJ, RN, GV, and ZC (CESM); BL and TI (MPI-ESM); MK (IPSL-CM5A2); AY (MIROC-ES2L); LM (LOVECLIM); JM and AS (UVic). F. Lhardy, D. M. Roche, and N. Bouttes generated new boundary conditions in the iLOVECLIM model. N. Bouttes and D. M. Roche developed the automated adjustments to allow for a change of ocean boundary conditions. F. Lhardy ran the iLOVECLIM simulations and analyzed both the iLOVECLIM and the PMIP-carbon outputs under supervision of N. Bouttes and D. M. Roche. F. Lhardy wrote the manuscript with the inputs from all co-authors. This study was supported by the French National program LEFE (Les Enveloppes Fluides et l'Environnement). F. Lhardy is supported by the Universit? Versailles Saint-Quentin-en-Yvelines (UVSQ). N. Bouttes and D. M. Roche are supported by the Centre national de la recherche scientifique (CNRS). In addition, D. M. Roche is supported by the Vrije Universiteit Amsterdam. B. Liu and T. Ilyina are supported by the German Federal Ministry of Education and Research (PalMod initiative, FKZ: 420 grant no. 01LP1919B). L. Menviel acknowledges funding from the Australian Research Council grant FT180100606. A. Yamamoto acknowledges funding from the Integrated Research Program for Advancing Climate Models (TOUGOU) Grant Number JPMXD0717935715 from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. The authors acknowledge the use of the LSCE storage and computing facilities. The authors thank Th?o Mandonnet for his preliminary work on the PMIP-carbon project. Last but not least, the authors thank the two anonymous reviewers for their help with this manuscript.

Publisher Copyright:
© 2021. American Geophysical Union. All Rights Reserved.

Funding

F. Lhardy, N. Bouttes, and D. M. Roche designed the research. N. Bouttes coordinated the PMIP‐carbon project and obtained funding. Participating modeling groups all performed a PI and a LGM simulation, provided their model outputs and the relevant metadata and computed the equilibrated carbon content in reservoirs. These modeling groups included AA‐O, HK, and AO (MIROC4m‐COCO); KC (CLIMBER‐2); MJ, RN, GV, and ZC (CESM); BL and TI (MPI‐ESM); MK (IPSL‐CM5A2); AY (MIROC‐ES2L); LM (LOVECLIM); JM and AS (UVic). F. Lhardy, D. M. Roche, and N. Bouttes generated new boundary conditions in the iLOVECLIM model. N. Bouttes and D. M. Roche developed the automated adjustments to allow for a change of ocean boundary conditions. F. Lhardy ran the iLOVECLIM simulations and analyzed both the iLOVECLIM and the PMIP‐carbon outputs under supervision of N. Bouttes and D. M. Roche. F. Lhardy wrote the manuscript with the inputs from all co‐authors. This study was supported by the French National program LEFE (Les Enveloppes Fluides et l'Environnement). F. Lhardy is supported by the Université Versailles Saint‐Quentin‐en‐Yvelines (UVSQ). N. Bouttes and D. M. Roche are supported by the Centre national de la recherche scientifique (CNRS). In addition, D. M. Roche is supported by the Vrije Universiteit Amsterdam. B. Liu and T. Ilyina are supported by the German Federal Ministry of Education and Research (PalMod initiative, FKZ: 420 grant no. 01LP1919B). L. Menviel acknowledges funding from the Australian Research Council grant FT180100606. A. Yamamoto acknowledges funding from the Integrated Research Program for Advancing Climate Models (TOUGOU) Grant Number JPMXD0717935715 from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. The authors acknowledge the use of the LSCE storage and computing facilities. The authors thank Théo Mandonnet for his preliminary work on the PMIP‐carbon project. Last but not least, the authors thank the two anonymous reviewers for their help with this manuscript. F. Lhardy, N. Bouttes, and D. M. Roche designed the research. N. Bouttes coordinated the PMIP-carbon project and obtained funding. Participating modeling groups all performed a PI and a LGM simulation, provided their model outputs and the relevant metadata and computed the equilibrated carbon content in reservoirs. These modeling groups included AA-O, HK, and AO (MIROC4m-COCO); KC (CLIMBER-2); MJ, RN, GV, and ZC (CESM); BL and TI (MPI-ESM); MK (IPSL-CM5A2); AY (MIROC-ES2L); LM (LOVECLIM); JM and AS (UVic). F. Lhardy, D. M. Roche, and N. Bouttes generated new boundary conditions in the iLOVECLIM model. N. Bouttes and D. M. Roche developed the automated adjustments to allow for a change of ocean boundary conditions. F. Lhardy ran the iLOVECLIM simulations and analyzed both the iLOVECLIM and the PMIP-carbon outputs under supervision of N. Bouttes and D. M. Roche. F. Lhardy wrote the manuscript with the inputs from all co-authors. This study was supported by the French National program LEFE (Les Enveloppes Fluides et l'Environnement). F. Lhardy is supported by the Universit? Versailles Saint-Quentin-en-Yvelines (UVSQ). N. Bouttes and D. M. Roche are supported by the Centre national de la recherche scientifique (CNRS). In addition, D. M. Roche is supported by the Vrije Universiteit Amsterdam. B. Liu and T. Ilyina are supported by the German Federal Ministry of Education and Research (PalMod initiative, FKZ: 420 grant no. 01LP1919B). L. Menviel acknowledges funding from the Australian Research Council grant FT180100606. A. Yamamoto acknowledges funding from the Integrated Research Program for Advancing Climate Models (TOUGOU) Grant Number JPMXD0717935715 from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. The authors acknowledge the use of the LSCE storage and computing facilities. The authors thank Th?o Mandonnet for his preliminary work on the PMIP-carbon project. Last but not least, the authors thank the two anonymous reviewers for their help with this manuscript.

FundersFunder number
Les Enveloppes Fluides et l'Environnement
MPI-ESMIPSL-CM5A2
Australian Research CouncilFT180100606, JPMXD0717935715
Ministry of Education, Culture, Sports, Science and Technology
Bundesministerium für Bildung und Forschung01LP1919B
Centre National de la Recherche Scientifique
Université de Versailles Saint-Quentin-en-Yvelines

    Keywords

    • alkalinity
    • atmospheric CO
    • carbon cycle
    • climate models
    • glacial-interglacial cycles
    • ocean volume

    Fingerprint

    Dive into the research topics of 'A First Intercomparison of the Simulated LGM Carbon Results Within PMIP-Carbon: Role of the Ocean Boundary Conditions'. Together they form a unique fingerprint.

    Cite this