An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Climate extremes and the carbon cycle

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Ice core-based isotopic constraints on past carbon cycle changes

High-precision ice core data on both atmospheric CO2 concentrations and their carbon isotopic composition (δ13Catm) provide improved constraints on the marine and terrestrial processes responsible for carbon cycle changes during the last two interglacials and the preceding glacial/interglacial transitions.

Holocene biomass burning and global dynamics of the carbon cycle

Fire regimes have changed during the Holocene due to changes in climate, vegetation, and in human practices. Here, we hypothesise that changes in fire regime may have affected the global CO2 concentration in the atmosphere through the Holocene. Our data are based on quantitative reconstructions of biomass burning deduced from stratified charcoal records from Europe, and South-, Central- and North America, and Oceania to test the fire-carbon release hypothesis. In Europe the significant increase of fire activity is dated ≈6000 cal. yr ago. In north-eastern North America burning activity was greatest before 7500 years ago, very low between 7500–3000 years, and has been increasing since 3000 years ago. In tropical America, the pattern is more complex and apparently latitudinally zonal. Maximum burning occurred in the southern Amazon basin and in Central America during the middle Holocene, and during the last 2000 years in the northern Amazon basin. In Oceania, biomass burning has decreased since a maximum 5000 years ago. Biomass burning has broadly increased in the Northern and Southern hemispheres throughout the second half of the Holocene associated with changes in climate and human practices. Global fire indices parallel the increase of atmospheric CO2 concentration recorded in Antarctic ice cores. Future issues on carbon dynamics relatively to biomass burning are discussed to improve the quantitative reconstructions.

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 Impact of Whaling on the Ocean Carbon Cycle: Why Bigger Was Better

Background: Humans have reduced the abundance of many large marine vertebrates, including whales, large fish, and sharks, to only a small percentage of their pre-exploitation levels. Industrial fishing and whaling also tended to preferentially harvest the largest species and largest individuals within a population. We consider the consequences of removing these animals on the ocean's ability to store carbon. Methodology/Principal Findings: Because body size is critical to our arguments, our analysis focuses on populations of baleen whales. Using reconstructions of pre-whaling and modern abundances, we consider the impact of whaling on the amount of carbon stored in living whales and on the amount of carbon exported to the deep sea by sinking whale carcasses. Populations of large baleen whales now store 9.1 x 10(6) tons less carbon than before whaling. Some of the lost storage has been offset by increases in smaller competitors; however, due to the relative metabolic efficiency of larger organisms, a shift toward smaller animals could decrease the total community biomass by 30% or more. Because of their large size and few predators, whales and other large marine vertebrates can efficiently export carbon from the surface waters to the deep sea. We estimate that rebuilding whale populations would remove 1.6 x 10(5) tons of carbon each year through sinking whale carcasses. Conclusions/Significance: Even though fish and whales are only a small portion of the ocean's overall biomass, fishing and whaling have altered the ocean's ability to store and sequester carbon. Although these changes are small relative to the total ocean carbon sink, rebuilding populations of fish and whales would be comparable to other carbon management schemes, including ocean iron fertilization.

North Pacific Carbon Cycle Response to Climate Variability on Seasonal to Decadal Timescales

Climate variability drives significant changes in the physical state of the North Pacific, and there may be important impacts of this variability on the upper ocean carbon balance across the basin. We address this issue by considering the response of seven biogeochemical ocean models to climate variability in the North Pacific. The models' upper ocean pCO(2) and air-sea CO(2) flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO(2) and its temperature- and non-temperature-driven components at three contrasting oceanographic sites capture the basic features found in observations (Takahashi et al., 2002, 2006; Keeling et al., 2004; Brix et al., 2004). However, particularly in the Western Subarctic Gyre, the models have difficulty representing the temporal structure of the total pCO(2) seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO(2) flux interannual variability (1 sigma) in the North Pacific is smaller ( ranges across models from 0.03 to 0.11 PgC/yr) than in the Tropical Pacific ( ranges across models from 0.08 to 0.19 PgC/yr), and the time series of the first or second EOF of the air-sea CO(2) flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO(2) flux anomalies are correlated with the PDO, their magnitudes are small ( up to +/- 0.025 PgC/yr ( 1 sigma)). Flux anomalies are damped because anomalies in the key drivers of pCO(2) ( temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO(2) anomalies.

Holocene carbon cycle dynamics, links to model files

We are investigating the late Holocene rise in CO2 by performing four experiments with the climate-carbon-cycle model CLIMBER2-LPJ. Apart from the deep sea sediments, important carbon cycle processes considered are carbon uptake or release by the vegetation, carbon uptake by peatlands, and CO 2 release due to shallow water sedimentation of CaCO3. Ice core data of atmospheric CO2 between 8 ka BP and preindustrial climate can only be reproduced if CO2 outgassing due to shallow water sedimentation of CaCO3 is considered. In this case the model displays an increase of nearly 20 ppmv CO2 between 8 ka BP and present day. Model configurations that do not contain this forcing show a slight decrease in atmospheric CO2. We can therefore explain the late Holocene rise in CO2 by invoking natural forcing factors only, and anthropogenic forcing is not required to understand preindustrial CO2 dynamics.

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.

Carbon cycle change in the Eocene

The transition from the extreme global warmth of the early Eocene 'greenhouse' climate ~55 million years ago to the present glaciated state is one of the most prominent changes in Earth's climatic evolution. It is widely accepted that large ice sheets first appeared on Antarctica ~34 million years ago, coincident with decreasing atmospheric carbon dioxide concentrations and a deepening of the calcite compensation depth in the world's oceans, and that glaciation in the Northern Hemisphere began much later, between 10 and 6 million years ago. Here we present records of sediment and foraminiferal geochemistry covering the greenhouse-icehouse climate transition. We report evidence for synchronous deepening and subsequent oscillations in the calcite compensation depth in the tropical Pacific and South Atlantic oceans from ~42 million years ago, with a permanent deepening 34 million years ago. The most prominent variations in the calcite compensation depth coincide with changes in seawater oxygen isotope ratios of up to 1.5 per mil, suggesting a lowering of global sea level through significant storage of ice in both hemispheres by at least 100 to 125 metres. Variations in benthic carbon isotope ratios of up to ~1.4 per mil occurred at the same time, indicating large changes in carbon cycling. We suggest that the greenhouse-icehouse transition was closely coupled to the evolution of atmospheric carbon dioxide, and that negative carbon cycle feedbacks may have prevented the permanent establishment of large ice sheets earlier than 34 million years ago.