The calcium looping cycle study for capturing carbon dioxide applied to the energy generation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Toward explaining the Holocene carbon dioxide and carbon isotope records: Results from transient ocean carbon cycle-climate simulations

[1] The Bern3D model was applied to quantify the mechanisms of carbon cycle changes during the Holocene (last 11,000 years). We rely on scenarios from the literature to prescribe the evolution of shallow water carbonate deposition and of land carbon inventory changes over the glacial termination (18,000 to 11,000 years ago) and the Holocene and modify these scenarios within uncertainties. Model results are consistent with Holocene records of atmospheric CO2 and δ13C as well as the spatiotemporal evolution of δ13C and carbonate ion concentration in the deep sea. Deposition of shallow water carbonate, carbonate compensation of land uptake during the glacial termination, land carbon uptake and release during the Holocene, and the response of the ocean-sediment system to marine changes during the termination contribute roughly equally to the reconstructed late Holocene pCO2 rise of 20 ppmv. The 5 ppmv early Holocene pCO2 decrease reflects terrestrial uptake largely compensated by carbonate deposition and ocean sediment responses. Additional small contributions arise from Holocene changes in sea surface temperature, ocean circulation, and export productivity. The Holocene pCO2 variations result from the subtle balance of forcings and processes acting on different timescales and partly in opposite direction as well as from memory effects associated with changes occurring during the termination. Different interglacial periods with different forcing histories are thus expected to yield different pCO2 evolutions as documented by ice cores.

Modeling the marine aragonite cycle: changes under rising carbon dioxide and its role in shallow water CaCO3 dissolution

The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by 29% and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.

Imminent ocean acidification in the Arctic projected with the NCAR global coupled carbon cycle-climate model

Ocean acidification from the uptake of anthropogenic carbon is simulated for the industrial period and IPCC SRES emission scenarios A2 and B1 with a global coupled carbon cycle-climate model. Earlier studies identified seawater saturation state with respect to aragonite, a mineral phase of calcium carbonate, as a key variable governing impacts on corals and other shell-forming organisms. Globally in the A2 scenario, water saturated by more than 300%, considered suitable for coral growth, vanishes by 2070 AD (CO2≈630 ppm), and the ocean volume fraction occupied by saturated water decreases from 42% to 25% over this century. The largest simulated pH changes worldwide occur in Arctic surface waters, where hydrogen ion concentration increases by up to 185% (ΔpH=−0.45). Projected climate change amplifies the decrease in Arctic surface mean saturation and pH by more than 20%, mainly due to freshening and increased carbon uptake in response to sea ice retreat. Modeled saturation compares well with observation-based estimates along an Arctic transect and simulated changes have been corrected for remaining model-data differences in this region. Aragonite undersaturation in Arctic surface waters is projected to occur locally within a decade and to become more widespread as atmospheric CO2 continues to grow. The results imply that surface waters in the Arctic Ocean will become corrosive to aragonite, with potentially large implications for the marine ecosystem, if anthropogenic carbon emissions are not reduced and atmospheric CO2 not kept below 450 ppm.

Comparative carbon cycle dynamics of the present and last interglacial

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO₂ remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO₂ fertilization, land use, wild fire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO₂ dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO₂ and δ¹³ CO₂ changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO₂ dynamics from 8ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO₂ changes after 122 ka BP. This failure to simulate late-Eemian CO₂ dynamics could be a result of the imposed forcings such as prescribed CaCO₃ accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO₂ dynamics eshallow water CaCO₃ accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO₂ dynamics.

The unstable CO2 feedback cycle on ocean planets

Ocean planets are volatile-rich planets, not present in our Solar system, which are thought to be dominated by deep, global oceans. This results in the formation of high-pressure water ice, separating the planetary crust from the liquid ocean and, thus, also from the atmosphere. Therefore, instead of a carbonate-silicate cycle like on the Earth, the atmospheric carbon dioxide concentration is governed by the capability of the ocean to dissolve carbon dioxide (CO2). In our study, we focus on the CO2 cycle between the atmosphere and the ocean which determines the atmospheric CO2 content. The atmospheric amount of CO2 is a fundamental quantity for assessing the potential habitability of the planet's surface because of its strong greenhouse effect, which determines the planetary surface temperature to a large degree. In contrast to the stabilizing carbonate-silicate cycle regulating the long-term CO2 inventory of the Earth atmosphere, we find that the CO2 cycle feedback on ocean planets is negative and has strong destabilizing effects on the planetary climate. By using a chemistry model for oceanic CO2 dissolution and an atmospheric model for exoplanets, we show that the CO2 feedback cycle can severely limit the extension of the habitable zone for ocean planets.

Late Neogene carbonate sedimentation and stable isotope record of ODP Site 121-758

The evolution of oceanic and climatic conditions the northeast Indian Ocean during the last 7 m.y. is revealed in the sediments from Site 758. We present detailed and continuous records of d18O and d13C from planktonic foraminifers, weight percent calcium carbonate, weight percent coarse fraction, magnetic susceptibility, and geomagnetic reversals. Sample spacing of the records ranges from 3 to 10 cm and is equivalent to an average time interval of 2000 to 6000 yr. Despite the fact that core recovery ranged between 100% and 105%, recovery gaps as large as 2.7 m occurred at nearly every break between advanced hydraulic piston cores. Approximately 12% of the late Neogene sequence was not recovered in each of the two holes drilled at Site 758. To circumvent the discontinuity introduced by the gaps, a composite depth section was constructed from multiple cores taken from offset holes at Site 758. The resulting composite depth section extends continuously from 0 to 116 mbsf, from the Holocene to the upper Miocene. A detailed chronostratigraphy is based on geomagnetic reversals which extend from the Brunhes Chron to Chron 6, and on d18O stages 1 through 105, which span from 0 to 2.5 Ma. The d18O record is dominated by a ~40-k.y. cycle in the late Pliocene and early Pleistocene, and is followed by a change to a ~100-k.y. cycle in the late Pleistocene. The mid-Pleistocene transition between these two modes of variability occurs between d18O stages 25 and 22 (between 860 and 800 Ka). Thirteen major volcanic ash horizons from the Indonesian arc are observed throughout the sedimentary section and are dated by their relative position within the geomagnetic reversals and the d18O chronostratigraphy. Since 5 Ma, there has been a long-term decline in weight percent CaCO3 and CaCO3 mass accumulation rates, and an associated rise in non-CaCO3 mass accumulation rates. We attribute these changes to a decrease in CaCO3 productivity and an increase in terrigenous sedimentation through enhanced riverine input. Such input may be linked to rapid tectonic uplift of the Himalayas and the Tibetan Plateau via mechanisms such as the intensification of the monsoonal rains, increased fluvial erosion, and regional glaciation. The long-term increase in percent coarse fraction since 5 Ma suggests a gradual increase in CaCO3 preservation. Higher frequency fluctuations in CaCO3 preservation are superimposed on the long-term trend and are related to climate fluctuations. The abrupt drop (-50%) in CaCO3 accumulation at 3.4 Ma signals a dramatic decrease in CaCO3 production that occurred over much of the Indian Ocean.

Eocene sedimentary calcium carbonate contents and stable isotope composition of benthic foraminifera

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

An early Eocene carbon cycle perturbation at in the Southern Alps

At least two transient events of extreme global warming occurred superimposed on the long-term latest Paleocene and early Eocene warming trend in the Paleocene-Eocene thermal maximum (PETM) (or ETM1 ~55.5 Ma) and the Elmo (or ETM2 ?53.6 Ma). Other than warmth, the best known PETM is characterized by (1) significant injection of 13C-depleted carbon into the ocean-atmosphere system, (2) deep-sea carbonate dissolution, (3) strong biotic responses, and (4) perturbations of the hydrological cycle. Documentation of the other documented and suspected "hyperthermals" is, as yet, insufficient to assess whether they are similar in nature to the PETM. Here we present and discuss biomagnetostratigraphic data and geochemical records across two lower Eocene successions deposited on a continental margin of the western Tethys: the Farra and Possagno sections in the Venetian pre-Alps. We recognize four negative carbon isotope excursions within chron C24. Three of these shifts correlate to known or suspected hyperthermals: the PETM, the Eocene thermal maximum 2 (~53.6 Ma), and the informally named "X event" (~52.5 Ma). The fourth excursion lies within a reverse subchron and occurred between the latter two. In the Farra section, the X event is marked by a ~0.6 per mil negative carbon isotope excursion and carbonate dissolution. Furthermore, the event exhibits responses among calcareous nannofossils, planktic foraminifera, and dinoflagellates that are similar to, though less intense than, those observed across the PETM. Sedimentological and quantitative micropaleontological data from the Farra section also suggest increased weathering and runoff as well as sea surface eutrophication during this event.

Calcium carbonate, isotopes and Rock-Eval pyrolysis at DSDP Holes 77-535 and 77-540

At DSDP Sites 534 (Central Atlantic) and 535 and 540 (Gulf of Mexico), and in the Vocontian Basin (France), Lower Cretaceous deposits show a very pronounced alternation of limestone and marl. This rhythm characterizes the pelagic background sedimentation and is independent of detritic intercalations related to contour and turbidity currents. Bed-scale cycles, estimated to be 6000-26,000 yr. long, comprise major and minor units. Their biological and mineralogic components, burrowing, heavy isotopes C and O, and some geochemical indicators, vary in close correlation with CaCO3 content. Vertical changes of frequency and asymmetry of the cycles are connected with fluctuations of the sedimentation rate. Plots of cycle thickness ("cyclograms") permit detailed correlations of the three areas and improve the stratigraphic subdivision of Neocomian deposits at the DSDP sites. Small-scale alternations, only observed in DSDP cores, comprise centimetric to millimetric banding and millimetric to micrometric lamination, here interpreted as varvelike alternations between laminae that are rich in calcareous plankton and others rich in clay. The laminations are estimated to correspond to cycles approximately 1,3, and 13 yr. in duration. The cyclic patterns appear to be governed by an interplay of continental and oceanic processes. Oceanic controls express themselves in variations of the biogenic carbonate flux, which depends on variations of such elements as temperature, oxygenation, salinity, and nutrient content. Continental controls modulate the influxes of terrigenous material, organic matter, and nutrients derived from cyclic erosion on land. Among the possible causes of cyclic sedimentation, episodic carbonate dissolution has been ruled out in favor of climatic fluctuations with a large range of periods. Such fluctuations are consistent with the great geographic extension shown by alternation controls and with the continuous spectrum of scales that characterizes limestone-marl cycles. The climatic variations induced by the Earth's orbital parameters (Milankovitch cycles) could be connected to bed-interbed alternations.

Changes in bottom water corrosiveness and carbonate ion concentrations in the sub-Antarctic Atlantic during the last glacial period

The glacial climate system transitioned rapidly between cold (stadial) and warm (interstadial) conditions in the Northern Hemisphere. This variability, referred to as Dansgaard-Oeschger variability, is widely believed to arise from perturbations of the Atlantic Meridional Overturning Circulation. Evidence for such changes during the longer Heinrich stadials has been identified, but direct evidence for overturning circulation changes during Dansgaard-Oeschger events has proven elusive. Here we reconstruct bottom water [CO3]2- variability from B/Ca ratios of benthic foraminifera and indicators of sedimentary dissolution, and use these reconstructions to infer the flow of northern-sourced deep water to the deep central sub-Antarctic Atlantic Ocean. We find that nearly every Dansgaard-Oeschger interstadial is accompanied by a rapid incursion of North Atlantic Deep Water into the deep South Atlantic. Based on these results and transient climate model simulations, we conclude that North Atlantic stadial-interstadial climate variability was associated with significant Atlantic overturning circulation changes that were rapidly transmitted across the Atlantic. However, by demonstrating the persistent role of Atlantic overturning circulation changes in past abrupt climate variability, our reconstructions of carbonate chemistry further indicate that the carbon cycle response to abrupt climate change was not a simple function of North Atlantic overturning.