Mid- to late Holocene Indian Ocean Monsoon variability recorded in four speleothems from Socotra Island, Yemen

Four stalagmites covering the last 7.0 ka were sampled on Socotra, an island in the northern Indian Ocean to investigate the evolution of the northeast Indian Ocean Monsoon (IOM) since the mid Holocene. On Socotra, rain is delivered at the start of the southwest IOM in May–June and at the start of the northeast IOM from September to December. The Haggeher Mountains act as a barrier forcing precipitation brought by the northeast winds to fall preferentially on the eastern side of the island, where the studied caves are located. δ18O and δ13C and Mg/Ca and Sr/Ca signals in the stalagmites reflect precipitation amounts brought by the northeast winds. For stalagmite STM6, this amount effect is amplified by kinetic effects during calcite deposition. Combined interpretation of the stalagmites' signals suggest a weakening of the northeast precipitation between 6.0 and 3.8 ka. After 3.8 ka precipitation intensities remain constant with two superimposed drier periods, between 0 and 0.6 ka and from 2.2 to 3.8 ka. No link can be established with Greenland ice cores and with the summer IOM variability. In contrast to the stable northeast rainy season suggested by the records in this study, speleothem records from western Socotra indicate a wettening of the southwest rainy season on Socotra after 4.4 ka. The local wettening of western Socotra could relate to a more southerly path (more over the Indian Ocean) taken by the southwest winds. Stalagmite STM5, sampled at the fringe between both rain areas displays intermediate δ18O values. After 6.2 ka, similar precipitation changes are seen between eastern Socotra and northern Oman indicating that both regions are affected similarly by the monsoon. Different palaeoclimatologic records from the Arabian Peninsula currently located outside the ITCZ migration pathway display an abrupt drying around 6 ka due to their disconnection from the southwest rain influence. Records that are nowadays still receiving rain by the southwest winds, suggest a more gradual drying reflecting the weakening of the southwest monsoon.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global three-dimensional ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during the early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

Fingerprints of changes in the terrestrial carbon cycle in response to large reorganizations in ocean circulation

CO2 and carbon cycle changes in the land, ocean and atmosphere are investigated using the comprehensive carbon cycle-climate model NCAR CSM1.4-carbon. Ensemble simulations are forced with freshwater perturbations applied at the North Atlantic and Southern Ocean deep water formation sites under pre-industrial climate conditions. As a result, the Atlantic Meridional Overturning Circulation reduces in each experiment to varying degrees. The physical climate fields show changes qualitatively in agreement with results documented in the literature, but there is a clear distinction between northern and southern perturbations. Changes in the physical variables, in turn, affect the land and ocean biogeochemical cycles and cause a reduction, or an increase, in the atmospheric CO2 concentration by up to 20 ppmv, depending on the location of the perturbation. In the case of a North Atlantic perturbation, the land biosphere reacts with a strong reduction in carbon stocks in some tropical locations and in high northern latitudes. In contrast, land carbon stocks tend to increase in response to a southern perturbation. The ocean is generally a sink of carbon although large reorganizations occur throughout various basins. The response of the land biosphere is strongest in the tropical regions due to a shift of the Intertropical Convergence Zone. The carbon fingerprints of this shift, either to the south or to the north depending on where the freshwater is applied, can be found most clearly in South America. For this reason, a compilation of various paleoclimate proxy records of Younger Dryas precipitation changes are compared with our model results. The proxy records, in general, show good agreement with the model's response to a North Atlantic freshwater perturbation.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global 3-D ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

Sensitivity of pelagic calcification to ocean acidification

Ocean acidification might reduce the ability of calcifying plankton to produce and maintain their shells of calcite, or of aragonite, the more soluble form of CaCO3. In addition to possibly large biological impacts, reduced CaCO3 production corresponds to a negative feedback on atmospheric CO2. In order to explore the sensitivity of the ocean carbon cycle to increasing concentrations of atmospheric CO2, we use the new biogeochemical Bern3D/PISCES model. The model reproduces the large scale distributions of biogeochemical tracers. With a range of sensitivity studies, we explore the effect of (i) using different parameterizations of CaCO3 production fitted to available laboratory and field experiments, of (ii) letting calcite and aragonite be produced by auto- and heterotrophic plankton groups, and of (iii) using carbon emissions from the range of the most recent IPCC Representative Concentration Pathways (RCP). Under a high-emission scenario, the CaCO3 production of all the model versions decreases from ~1 Pg C yr−1 to between 0.36 and 0.82 Pg C yr−1 by the year 2100. The changes in CaCO3 production and dissolution resulting from ocean acidification provide only a small feedback on atmospheric CO2 of −1 to −11 ppm by the year 2100, despite the wide range of parameterizations, model versions and scenarios included in our study. A potential upper limit of the CO2-calcification/dissolution feedback of −30 ppm by the year 2100 is computed by setting calcification to zero after 2000 in a high 21st century emission scenario. The similarity of feedback estimates yielded by the model version with calcite produced by nanophytoplankton and the one with calcite, respectively aragonite produced by mesozooplankton suggests that expending biogeochemical models to calcifying zooplankton might not be needed to simulate biogeochemical impacts on the marine carbonate cycle. The changes in saturation state confirm previous studies indicating that future anthropogenic CO2 emissions may lead to irreversible changes in ΩA for several centuries. Furthermore, due to the long-term changes in the deep ocean, the ratio of open water CaCO3 dissolution to production stabilizes by the year 2500 at a value that is 30–50% higher than at pre-industrial times when carbon emissions are set to zero after 2100.