Using ocean surface paleo-density to evaluate PMIP3 and PMIP4 Last Glacial Maximum climate simulations
- 1Université de Bordeaux (Bordeaux, France)
- 2Centre National de la Recherche Scientifique (Paris, France)
- 3Environnements et Paléoenvironnements Océaniques et Continentaux (Pessac, France)
Abstract
Quantitative reconstruction of ocean surface density during the Last Glacial Maximum (LGM) offers valuable insights into the ability of climate models to simulate past climate conditions, when global temperatures were about 4.5 to 6°C colder than today. We assess the performance of the LGM climate simulations, as part of the 3rd and 4th phase of the Paleoclimate Modeling Intercomparisons Project, using a recent ocean surface density reconstruction based on the δ ¹⁸ O of foraminiferal calcite ( δ ¹⁸ O c ). We consider the differences between the LGM and the preindustrial climates and each period separately, at both global and regional scales. Because surface density reflects the combined effects of temperature and salinity, we also examined sea surface temperature (SST) to better identify the processes underlying model–data differences. On a global scale, surface density reconstructions generally exhibit greater spatial variability than simulated surface density anomalies (LGM − PI), although part of this difference is reduced when reconstruction uncertainties are taken into account. Model simulations tend to underestimate the magnitude of reconstructed density anomalies and substantial differences between cumulative distribution persist. Part of the mismatch may arise from the uneven spatial distribution of reconstructions, which are mostly located near coastal areas. Density anomaly (LGM − PI) differences between data and models are largely controlled by sea surface salinity (SSS), with SST contributing to a lesser extent. This influence of SSS is directly linked to the reduction in tropical precipitation during the LGM: model simulations that best match the large-scale density anomalies also simulate the strongest reductions in reconstructed low-latitude precipitation during the LGM, highlighting the key role of hydrological cycle changes in shaping surface density. All simulations capture key features of the reconstructed surface densities when the LGM and PI periods are analysed separately at the global scale, with Taylor diagrams and complementary performance metrics indicating moderate to good overall agreement despite differences among simulations. Regional analyses show that most simulations reproduce reconstructed Indian Ocean surface density reasonably well, although model simulations performance is systematically lower in the North Indian Ocean than in the South Indian Ocean. Focusing on the Indo-Pacific Warm Pool, proxy reconstructions indicate a weakened West-East tropical Indian surface density gradient during the LGM, but only 7 out of 14 model simulations (50 %) reproduce this feature. These results highlight the need to improve and better constrain regional hydrological cycle changes in models, as improving their representation is crucial to reduce uncertainties in both paleoclimate simulations and future climate projections.
Versions
Related
- Sea-Level Change in Asia since the Last Glacial Maximum
- Strong latitudinal contrasts in global wet/dry patterns since the Last Glacial Maximum and under future warming
- Warming since the Last Glacial Maximum amplifies Northern Hemisphere drought expansion driven by the Pacific Decadal Oscillation
- Redefining Last Glacial Maximum environmental conditions in northeastern Brazil: pollen-based evidence of open and moist-adapted vegetation cover
- A framework to evaluate and identify development requirements for land‐surface models at km‐scale resolution: Application to a semi‐arid and mountainous region
BibTeX
@article{barathieu2026using,
title = {{Using ocean surface paleo-density to evaluate PMIP3 and PMIP4 Last Glacial Maximum climate simulations}},
author = {Héloïse Barathieu and Thibaut Caley and Masa Kageyama and Didier Swingedouw and Pascale Braconnot},
journal = {Climate of the Past},
year = {2026},
month = oct,
doi = {10.5194/cp-22-1803-2026},
url = {https://cp.copernicus.org/articles/22/1803/2026/},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Using ocean surface paleo-density to evaluate PMIP3 and PMIP4 Last Glacial Maximum climate simulations
AU - Héloïse Barathieu
AU - Thibaut Caley
AU - Masa Kageyama
AU - Didier Swingedouw
AU - Pascale Braconnot
PY - 2026
DA - 2026/10/06
JO - Climate of the Past
DO - 10.5194/cp-22-1803-2026
UR - https://cp.copernicus.org/articles/22/1803/2026/
AB - Quantitative reconstruction of ocean surface density during the Last Glacial Maximum (LGM) offers valuable insights into the ability of climate models to simulate past climate conditions, when global temperatures were about 4.5 to 6°C colder than today. We assess the performance of the LGM climate simulations, as part of the 3rd and 4th phase of the Paleoclimate Modeling Intercomparisons Project, using a recent ocean surface density reconstruction based on the δ ¹⁸ O of foraminiferal calcite ( δ ¹⁸ O c ). We consider the differences between the LGM and the preindustrial climates and each period separately, at both global and regional scales. Because surface density reflects the combined effects of temperature and salinity, we also examined sea surface temperature (SST) to better identify the processes underlying model–data differences. On a global scale, surface density reconstructions generally exhibit greater spatial variability than simulated surface density anomalies (LGM − PI), although part of this difference is reduced when reconstruction uncertainties are taken into account. Model simulations tend to underestimate the magnitude of reconstructed density anomalies and substantial differences between cumulative distribution persist. Part of the mismatch may arise from the uneven spatial distribution of reconstructions, which are mostly located near coastal areas. Density anomaly (LGM − PI) differences between data and models are largely controlled by sea surface salinity (SSS), with SST contributing to a lesser extent. This influence of SSS is directly linked to the reduction in tropical precipitation during the LGM: model simulations that best match the large-scale density anomalies also simulate the strongest reductions in reconstructed low-latitude precipitation during the LGM, highlighting the key role of hydrological cycle changes in shaping surface density. All simulations capture key features of the reconstructed surface densities when the LGM and PI periods are analysed separately at the global scale, with Taylor diagrams and complementary performance metrics indicating moderate to good overall agreement despite differences among simulations. Regional analyses show that most simulations reproduce reconstructed Indian Ocean surface density reasonably well, although model simulations performance is systematically lower in the North Indian Ocean than in the South Indian Ocean. Focusing on the Indo-Pacific Warm Pool, proxy reconstructions indicate a weakened West-East tropical Indian surface density gradient during the LGM, but only 7 out of 14 model simulations (50 %) reproduce this feature. These results highlight the need to improve and better constrain regional hydrological cycle changes in models, as improving their representation is crucial to reduce uncertainties in both paleoclimate simulations and future climate projections.
SN - 1814-9332
ER -