Factors influencing anthropogenic carbon dioxide uptake in the North Atlantic in models of the ocean carbon cycle

The uptake and storage of anthropogenic carbon in the North Atlantic is investigated using different configurations of ocean general circulation/carbon cycle models. We investigate how different representations of the ocean physics in the models, which represent the range of models currently in use, affect the evolution of CO2 uptake in the North Atlantic. The buffer effect of the ocean carbon system would be expected to reduce ocean CO2 uptake as the ocean absorbs increasing amounts of CO2. We find that the strength of the buffer effect is very dependent on the model ocean state, as it affects both the magnitude and timing of the changes in uptake. The timescale over which uptake of CO2 in the North Atlantic drops to below preindustrial levels is particularly sensitive to the ocean state which sets the degree of buffering; it is less sensitive to the choice of atmospheric CO2 forcing scenario. Neglecting physical climate change effects, North Atlantic CO2 uptake drops below preindustrial levels between 50 and 300 years after stabilisation of atmospheric CO2 in different model configurations. Storage of anthropogenic carbon in the North Atlantic varies much less among the different model configurations, as differences in ocean transport of dissolved inorganic carbon and uptake of CO2 compensate each other. This supports the idea that measured inventories of anthropogenic carbon in the real ocean cannot be used to constrain the surface uptake. Including physical climate change effects reduces anthropogenic CO2 uptake and storage in the North Atlantic further, due to the combined effects of surface warming, increased freshwater input, and a slowdown of the meridional overturning circulation. The timescale over which North Atlantic CO2 uptake drops to below preindustrial levels is reduced by about one-third, leading to an estimate of this timescale for the real world of about 50 years after the stabilisation of atmospheric CO2. In the climate change experiment, a shallowing of the mixed layer depths in the North Atlantic results in a significant reduction in primary production, reducing the potential role for biology in drawing down anthropogenic CO2.

Including an ocean carbon cycle model into iLOVECLIM (v1.0)

The atmospheric carbon dioxide concentration plays a crucial role in the radiative balance and as such has a strong influence on the evolution of climate. Because of the numerous interactions between climate and the carbon cycle, it is necessary to include a model of the carbon cycle within a climate model to understand and simulate past and future changes of the carbon cycle. In particular, natural variations of atmospheric CO2 have happened in the past, while anthropogenic carbon emissions are likely to continue in the future. To study changes of the carbon cycle and climate on timescales of a few hundred to a few thousand years, we have included a simple carbon cycle model into the iLOVECLIM Earth System Model. In this study, we describe the ocean and terrestrial biosphere carbon cycle models and their performance relative to observational data. We focus on the main carbon cycle variables including the carbon isotope ratios δ13C and the Δ14C. We show that the model results are in good agreement with modern observations both at the surface and in the deep ocean for the main variables, in particular phosphates, dissolved inorganic carbon and the carbon isotopes.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.

Ecosystem effects of CO2 concentration: evidence from past climates

Atmospheric CO2 concentration has varied from minima of 170-200 ppm in glacials to maxima of 280-300 ppm in the recent interglacials. Photosynthesis by C-3 plants is highly sensitive to CO2 concentration variations in this range. Physiological consequences of the CO2 changes should therefore be discernible in palaeodata. Several lines of evidence support this expectation. Reduced terrestrial carbon storage during glacials, indicated by the shift in stable isotope composition of dissolved inorganic carbon in the ocean, cannot be explained by climate or sea-level changes. It is however consistent with predictions of current process-based models that propagate known physiological CO2 effects into net primary production at the ecosystem scale. Restricted forest cover during glacial periods, indicated by pollen assemblages dominated by non-arboreal taxa, cannot be reproduced accurately by palaeoclimate models unless CO2 effects on C-3-C-4 plant competition are also modelled. It follows that methods to reconstruct climate from palaeodata should account for CO2 concentration changes. When they do so, they yield results more consistent with palaeoclimate models. In conclusion, the palaeorecord of the Late Quaternary, interpreted with the help of climate and ecosystem models, provides evidence that CO2 effects at the ecosystem scale are neither trivial nor transient.

Late Barremian-Early Aptian palaeoenvironmental change: the Cassis-La Bedoule section (SE France)

The well-dated section of Cassis-La Bédoule in the South Provencal Basin (southern France) allows for a detailed reconstruction of palaeoenvironmental change during the latest Barremian and Early Aptian. For this study, phosphorus (P) and clay-mineral contents, stable-isotope ratios on carbonate (δ13Ccarb) and organic matter (δ13Corg), and redox-sensitive trace elements (RSTE: V, U, As, Co, and Mo) have been measured in this historical stratotype. The base of the section consists of rudist limestone, which is attributed to the Urgonian platform. The presence of low P and RSTE content, and content of up to 30% kaolinite indicate deposition under oligotrophic and oxic conditions, and the presence of warm, humid climatic conditions on the adjacent continent. The top of the Urgonian succession is marked by a hardground with encrusted brachiopods and bivalves, which is interpreted as a drowning surface. The section continues with a succession of limestone and marl containing the first occurrence of planktonic foraminifera. This interval includes several laminated, organic-rich layers recording RSTE enrichments and high Corg:Ptot ratios. The deposition of these organic-rich layers was associated with oxygen-depleted conditions and a large positive excursion in δ13Corg. During this interval, a negative peak in the δ13Ccarb record is observed, which dates as latest Barremian. This excursion is coeval with negative excursions elsewhere in Tethyan platform and basin settings and is explained by the increased input of light dissolved inorganic carbon by rivers and/or volcanic activity. In this interval, an increase in P content, owing to reworking of nearshore sediments during the transgression, is coupled with a decrease in kaolinite content, which tends to be deposited in more proximal areas. The overlying hemipelagic sediments of the Early Aptian Deshayesites oglanlensis and D. weissi zones indicate rather stable palaeoenvironmental conditions with low P content and stable δ13C records. A change towards marl-dominated beds occurs close to the end of the D. weissi zone. These beds display a long decrease in their δ13Ccarb and δ13Corg records, which lasted until the end of the Deshayesites deshayesi subzone (corresponding to C3 in Menegatti et al., 1998). This is followed by a positive shift during the Roloboceras hambrovi and Deshayesites grandis subzones, which corresponds in time to oceanic anoxic event (OAE) 1a interval. This positive shift is coeval with two increases in the P content. The marly interval equivalent to OAE 1a lacks organic-rich deposits and RSTE enrichments indicating that oxic conditions prevailed in this particular part of the Tethys ocean. The clay mineralogy is dominated by smectite, which is interpreted to reflect trapping of kaolinite on the surrounding platforms rather than indicating a drier climate.

Modeling the photo-oxidation of dissolved organic matter by ultraviolet radiation in freshwater lakes: implications for mercury bioavailability

Uncertainties in projected ultraviolet (UV) radiation may lead to future increases in UV irradiation of freshwater lakes. Because dissolved organic carbon (DOC) is the main binding phase for mercury (Hg) in freshwater lakes, an increase in DOC photo-oxidation may affect Hg speciation and bioavailability. We quantified the effect of DOC concentration on the rate of abiotic DOC photo-oxidation for five lakes (DOC = 3.27–12.3 mg L−1) in Kejimkujik National Park, Canada. Samples were irradiated with UV-A or UV-B radiation over a 72-h period. UV-B radiation was found to be 2.36 times more efficient at photo-oxidizing DOC than UV-A, with energy-normalized rates of dissolved inorganic carbon (DIC) production ranging from 3.8 × 10−5 to 1.1 × 10−4 mg L−1 J−1 for UV-A, and from 6.0 × 10−5 to 3.1 × 10−4 mg L−1 J−1 for UV-B. Energy normalized rates of DIC production were positively correlated with DOC concentrations. Diffuse integrated attenuation coefficients were quantified in situ (UV-A Kd = 0.056–0.180 J cm−1; UV-B Kd = 0.015–0.165 J cm−1) and a quantitative depth-integrated model for yearly DIC photo-production in each lake was developed. The model predicts that, UV-A produces between 3.2 and 100 times more DIC (1521–2851 mg m−2 year−1) than UV-B radiation (29.17–746.7 mg m−2 year−1). Future increases in UV radiation may increase DIC production and increase Hg bioavailability in low DOC lakes to a greater extent than in high DOC lakes.