Description: During the Last Glacial Maximum (LGM, ∼21,000 years ago) the cold climate was strongly tied to low atmospheric CO〈inf〉2〈/inf〉 concentration (∼190 ppm). Although it is generally assumed that this low CO〈inf〉2〈/inf〉 was due to an expansion of the oceanic carbon reservoir, simulating the glacial level has remained a challenge especially with the additional δ13C constraint. Indeed the LGM carbon cycle was also characterized by a modern-like δ13C in the atmosphere and a higher surface to deep Atlantic δ13C gradient indicating probable changes in the thermohaline circulation. Here we show with a model of intermediate complexity, that adding three oceanic mechanisms: brine induced stratification, stratification-dependant diffusion and iron fertilization to the standard glacial simulation (which includes sea level drop, temperature change, carbonate compensation and terrestrial carbon release) decreases CO〈inf〉2〈/inf〉 down to the glacial value of ∼190 ppm and simultaneously matches glacial atmospheric and oceanic δ13C inferred from proxy data. LGM CO〈inf〉2〈/inf〉 and δ13C can at last be successfully reconciled.

Description: During the cold period of the Last Glacial Maximum (LGM, about 21 000 years ago) atmospheric CO2 was around 190 ppm, much lower than the pre-industrial concentration of 280 ppm. The causes of this substantial drop remain partially unresolved, despite intense research. Understanding the origin of reduced atmospheric CO2 during glacial times is crucial to comprehend the evolution of the different carbon reservoirs within the Earth system (atmosphere, terrestrial biosphere and ocean). In this context, the ocean is believed to play a major role as it can store large amounts of carbon, especially in the abyss, which is a carbon reservoir that is thought to have expanded during glacial times. To create this larger reservoir, one possible mechanism is to produce very dense glacial waters, thereby stratifying the deep ocean and reducing the carbon exchange between the deep and upper ocean. The existence of such very dense waters has been inferred in the LGM deep Atlantic from sediment pore water salinity and δ18O inferred temperature. Based on these observations, we study the impact of a brine mechanism on the glacial carbon cycle. This mechanism relies on the formation and rapid sinking of brines, very salty water released during sea ice formation, which brings salty dense water down to the bottom of the ocean. It provides two major features: a direct link from the surface to the deep ocean along with an efficient way of setting a strong stratification. We show with the CLIMBER-2 carbon-climate model that such a brine mechanism can account for a significant decrease in atmospheric CO2 and contribute to the glacial-interglacial change. This mechanism can be amplified by low vertical diffusion resulting from the brine-induced stratification. The modeled glacial distribution of oceanic δ13C as well as the deep ocean salinity are substantially improved and better agree with reconstructions from sediment cores, suggesting that such a mechanism could have played an important role during glacial times.

Description: Model intercomparison studies of coupled carbon-climate simulations have the potential to improve our understanding of the processes explaininEg the pCO2 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 demonstErate a ~ 250 GtC effect of an alkalinity adjustment on carbon sequestration in the ocean. Finally, we observe that PMIP-carbon models with a freely evolving CO2 and no additional glacial mechanisms do not simulaEte the pCO2 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.