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.

Carbon dioxide induced stomatal closure increases radiative forcing of climate via a rapid reduction in low cloud

We performed an ensemble of twelve five-year experiments using a coupled climate-carbon-cycle model with scenarios of prescribed atmospheric carbon dioxide concentration; CO2 was instantaneously doubled or quadrupled at the start of the experiments. Within these five years, climate feedback is not significantly influenced by the effects of climate change on the carbon system. However, rapid changes take place, within much less than a year, due to the physiological effect of CO2 on plant stomatal conductance, leading to adjustment in the shortwave cloud radiative effect over land, due to a reduction in low cloud cover. This causes a 10% enhancement to the radiative forcing due to CO2, which leads to an increase in the equilibrium warming of 0.4 and 0.7 K for doubling and quadrupling. The implications for calibration of energy-balance models are discussed.

Alternative explanations for rising dissolved organic carbon export from organic soils

Since 1988, there has been, on average, a 91% increase in dissolved organic carbon (DOC) concentrations of UK lakes and streams in the Acid Waters Monitoring Network (AWMN). Similar DOC increases have been observed in surface waters across much of Europe and North America. Much of the debate about the causes of rising DOC has, as in other studies relating to the carbon cycle, focused on factors related to climate change. Data from our peat-core experiments support an influence of climate on DOC, notably an increase in production with temperature under aerobic, and to a lesser extent anaerobic, conditions. However, we argue that climatic factors may not be the dominant drivers of DOC change. DOC solubility is suppressed by high soil water acidity and ionic strength, both of which have decreased as a result of declining sulphur deposition since the 1980s, augmented during the 1990s in the United Kingdom by a cyclical decline in sea-salt deposition. Our observational and experimental data demonstrate a clear, inverse and quantitatively important link between DOC and sulphate concentrations in soil solution. Statistical analysis of 11 AWMN lakes suggests that rising temperature, declining sulphur deposition and changing sea-salt loading can account for the majority of the observed DOC trend. This combination of evidence points to the changing chemical composition of atmospheric deposition, particularly the substantial reduction in anthropogenic sulphur emissions during the last 20 years, as a key cause of rising DOC. The implications of rising DOC export for the carbon cycle will be very different if linked primarily to decreasing acid deposition, rather than to changes in climate, suggesting that these systems may be recovering rather than destabilising.

Global climate change and soil carbon stocks: predictions from two contrasting models for the turnover of organic carbon in soil

Enhanced release of CO2 to the atmosphere from soil organic carbon as a result of increased temperatures may lead to a positive feedback between climate change and the carbon cycle, resulting in much higher CO2 levels and accelerated lobal warming. However, the magnitude of this effect is uncertain and critically dependent on how the decomposition of soil organic C (heterotrophic respiration) responds to changes in climate. Previous studies with the Hadley Centre’s coupled climate–carbon cycle general circulation model (GCM) (HadCM3LC) used a simple, single-pool soil carbon model to simulate the response. Here we present results from numerical simulations that use the more sophisticated ‘RothC’ multipool soil carbon model, driven with the same climate data. The results show strong similarities in the behaviour of the two models, although RothC tends to simulate slightly smaller changes in global soil carbon stocks for the same forcing. RothC simulates global soil carbon stocks decreasing by 54 GtC by 2100 in a climate change simulation compared with an 80 GtC decrease in HadCM3LC. The multipool carbon dynamics of RothC cause it to exhibit a slower magnitude of transient response to both increased organic carbon inputs and changes in climate. We conclude that the projection of a positive feedback between climate and carbon cycle is robust, but the magnitude of the feedback is dependent on the structure of the soil carbon model.

Biotic carbon feedbacks in a materially-closed soil-vegetation-atmosphere system

The magnitude and direction of the coupled feedbacks between the biotic and abiotic components of the terrestrial carbon cycle is a major source of uncertainty in coupled climate–carbon-cycle models1, 2, 3. Materially closed, energetically open biological systems continuously and simultaneously allow the two-way feedback loop between the biotic and abiotic components to take place4, 5, 6, 7, but so far have not been used to their full potential in ecological research, owing to the challenge of achieving sustainable model systems6, 7. We show that using materially closed soil–vegetation–atmosphere systems with pro rata carbon amounts for the main terrestrial carbon pools enables the establishment of conditions that balance plant carbon assimilation, and autotrophic and heterotrophic respiration fluxes over periods suitable to investigate short-term biotic carbon feedbacks. Using this approach, we tested an alternative way of assessing the impact of increased CO2 and temperature on biotic carbon feedbacks. The results show that without nutrient and water limitations, the short-term biotic responses could potentially buffer a temperature increase of 2.3 °C without significant positive feedbacks to atmospheric CO2. We argue that such closed-system research represents an important test-bed platform for model validation and parameterization of plant and soil biotic responses to environmental changes.

The proportionality of global warming to cumulative carbon emissions

The global temperature response to increasing atmospheric CO2 is often quantified by metrics such as equilibrium climate sensitivity and transient climate response1. These approaches, however, do not account for carbon cycle feedbacks and therefore do not fully represent the net response of the Earth system to anthropogenic CO2 emissions. Climate–carbon modelling experiments have shown that: (1) the warming per unit CO2 emitted does not depend on the background CO2 concentration2; (2) the total allowable emissions for climate stabilization do not depend on the timing of those emissions3, 4, 5; and (3) the temperature response to a pulse of CO2 is approximately constant on timescales of decades to centuries3, 6, 7, 8. Here we generalize these results and show that the carbon–climate response (CCR), defined as the ratio of temperature change to cumulative carbon emissions, is approximately independent of both the atmospheric CO2 concentration and its rate of change on these timescales. From observational constraints, we estimate CCR to be in the range 1.0–2.1 °C per trillion tonnes of carbon (Tt C) emitted (5th to 95th percentiles), consistent with twenty-first-century CCR values simulated by climate–carbon models. Uncertainty in land-use CO2 emissions and aerosol forcing, however, means that higher observationally constrained values cannot be excluded. The CCR, when evaluated from climate–carbon models under idealized conditions, represents a simple yet robust metric for comparing models, which aggregates both climate feedbacks and carbon cycle feedbacks. CCR is also likely to be a useful concept for climate change mitigation and policy; by combining the uncertainties associated with climate sensitivity, carbon sinks and climate–carbon feedbacks into a single quantity, the CCR allows CO2-induced global mean temperature change to be inferred directly from cumulative carbon emissions.

Warming caused by cumulative carbon emissions towards the trillionth tonne

Global efforts to mitigate climate change are guided by projections of future temperatures1. But the eventual equilibrium global mean temperature associated with a given stabilization level of atmospheric greenhouse gas concentrations remains uncertain1, 2, 3, complicating the setting of stabilization targets to avoid potentially dangerous levels of global warming4, 5, 6, 7, 8. Similar problems apply to the carbon cycle: observations currently provide only a weak constraint on the response to future emissions9, 10, 11. Here we use ensemble simulations of simple climate-carbon-cycle models constrained by observations and projections from more comprehensive models to simulate the temperature response to a broad range of carbon dioxide emission pathways. We find that the peak warming caused by a given cumulative carbon dioxide emission is better constrained than the warming response to a stabilization scenario. Furthermore, the relationship between cumulative emissions and peak warming is remarkably insensitive to the emission pathway (timing of emissions or peak emission rate). Hence policy targets based on limiting cumulative emissions of carbon dioxide are likely to be more robust to scientific uncertainty than emission-rate or concentration targets. Total anthropogenic emissions of one trillion tonnes of carbon (3.67 trillion tonnes of CO2), about half of which has already been emitted since industrialization began, results in a most likely peak carbon-dioxide-induced warming of 2 °C above pre-industrial temperatures, with a 5–95% confidence interval of 1.3–3.9 °C.

Whole life carbon – a building case study

There is potential to reduce both operational and embodied greenhouse gas emission from buildings. To date the focus has been on reducing the operational element, although given the urgency of carbon reductions, it may be more beneficial to consider upfront embodied carbon reductions. This paper describes a case study on the whole life carbon cycle of a warehouse building in Swindon, UK. It examines the relationship between embodied carbon (Ec) and operational carbon (Oc), the proportions of Ec from the structural and non-structural elements, carbon benchmarking of the structure, the value of ‘cradle to site’ or ‘cradle to grave’ assessments and the significance of the timing of emissions during the life of the building. The case study indicates that Ec was dominant for the building and that the structure was responsible for more than half of the Ec. Weighting of future emissions appears to be an important factor to consider. The PAS 2050 reduction factors had only a modest effect but weighting to allow for future decarbonisation of the national grid energy supply had a large effect. This suggests that future operational carbon emissions are being overestimated compared to embodied.

Climate-related changes in peatland carbon accumulation during the last millennium.

Peatlands are a major terrestrial carbon store and a persistent natural carbon sink during the Holocene, but there is considerable uncertainty over the fate of peatland carbon in a changing climate. It is generally assumed that higher temperatures will increase peat decay, causing a positive feedback to climate warming and contributing to the global positive carbon cycle feedback. Here we use a new extensive database of peat profiles across northern high latitudes to examine spatial and temporal patterns of carbon accumulation over the past millennium. Opposite to expectations, our results indicate a small negative carbon cycle feedback from past changes in the long-term accumulation rates of northern peatlands. Total carbon accumulated over the last 1000 yr is linearly related to contemporary growing season length and photosynthetically active radiation, suggesting that variability in net primary productivity is more important than decomposition in determining long-term carbon accumulation. Furthermore, northern peatland carbon sequestration rate declined over the climate transition from the Medieval Climate Anomaly (MCA) to the Little Ice Age (LIA), probably because of lower LIA temperatures combined with increased cloudiness suppressing net primary productivity. Other factors including changing moisture status, peatland distribution, fire, nitrogen deposition, permafrost thaw and methane emissions will also influence future peatland carbon cycle feedbacks, but our data suggest that the carbon sequestration rate could increase over many areas of northern peatlands in a warmer future.

Optimising the FAMOUS climate model: inclusion of global carbon cycling

FAMOUS fills an important role in the hierarchy of climate models, both explicitly resolving atmospheric and oceanic dynamics yet being sufficiently computationally efficient that either very long simulations or large ensembles are possible. An improved set of carbon cycle parameters for this model has been found using a perturbed physics ensemble technique. This is an important step towards building the "Earth System" modelling capability of FAMOUS, which is a reduced resolution, and hence faster running, version of the Hadley Centre Climate model, HadCM3. Two separate 100 member perturbed parameter ensembles were performed; one for the land surface and one for the ocean. The land surface scheme was tested against present-day and past representations of vegetation and the ocean ensemble was tested against observations of nitrate. An advantage of using a relatively fast climate model is that a large number of simulations can be run and hence the model parameter space (a large source of climate model uncertainty) can be more thoroughly sampled. This has the associated benefit of being able to assess the sensitivity of model results to changes in each parameter. The climatologies of surface and tropospheric air temperature and precipitation are improved relative to previous versions of FAMOUS. The improved representation of upper atmosphere temperatures is driven by improved ozone concentrations near the tropopause and better upper level winds.

Large inert carbon pool in the terrestrial biosphere at the Last Glacial Maximum

During each of the late Pleistocene glacial–interglacial transitions, atmospheric carbon dioxide concentrations rose by almost 100 ppm. The sources of this carbon are unclear, and efforts to identify them are hampered by uncertainties in the magnitude of carbon reservoirs and fluxes under glacial conditions. Here we use oxygen isotope measurements from air trapped in ice cores and ocean carbon-cycle modelling to estimate terrestrial and oceanic gross primary productivity during the Last Glacial Maximum. We find that the rate of gross terrestrial primary production during the Last Glacial Maximum was about 40±10 Pg C yr−1, half that of the pre-industrial Holocene. Despite the low levels of photosynthesis, we estimate that the late glacial terrestrial biosphere contained only 330 Pg less carbon than pre-industrial levels. We infer that the area covered by carbon-rich but unproductive biomes such as tundra and cold steppes was significantly larger during the Last Glacial Maximum, consistent with palaeoecological data. Our data also indicate the presence of an inert carbon pool of 2,300 Pg C, about 700 Pg larger than the inert carbon locked in permafrost today. We suggest that the disappearance of this carbon pool at the end of the Last Glacial Maximum may have contributed to the deglacial rise in atmospheric carbon dioxide concentrations.

Modelling fire and the terrestrial carbon balance

Four CO2 concentration inversions and the Global Fire Emissions Database (GFED) versions 2.1 and 3 are used to provide benchmarks for climate-driven modeling of the global land-atmosphere CO2 flux and the contribution of wildfire to this flux. The Land surface Processes and exchanges (LPX) model is introduced. LPX is based on the Lund-Potsdam-Jena Spread and Intensity of FIRE (LPJ-SPITFIRE) model with amended fire probability calculations. LPX omits human ignition sources yet simulates many aspects of global fire adequately. It captures the major features of observed geographic pattern in burnt area and its seasonal timing and the unimodal relationship of burnt area to precipitation. It simulates features of geographic variation in the sign of the interannual correlations of burnt area with antecedent dryness and precipitation. It simulates well the interannual variability of the global total land-atmosphere CO2 flux. There are differences among the global burnt area time series from GFED2.1, GFED3 and LPX, but some features are common to all. GFED3 fire CO2 fluxes account for only about 1/3 of the variation in total CO2 flux during 1997–2005. This relationship appears to be dominated by the strong climatic dependence of deforestation fires. The relationship of LPX-modeled fire CO2 fluxes to total CO2 fluxes is weak. Observed and modeled total CO2 fluxes track the El Niño–Southern Oscillation (ENSO) closely; GFED3 burnt area and global fire CO2 flux track the ENSO much less so. The GFED3 fire CO2 flux-ENSO connection is most prominent for the El Niño of 1997–1998, which produced exceptional burning conditions in several regions, especially equatorial Asia. The sign of the observed relationship between ENSO and fire varies regionally, and LPX captures the broad features of this variation. These complexities underscore the need for process-based modeling to assess the consequences of global change for fire and its implications for the carbon cycle.

Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum

We use a state-of-the-art ocean general circulation and biogeochemistry model to examine the impact of changes in ocean circulation and biogeochemistry in governing the change in ocean carbon-13 and atmospheric CO2 at the last glacial maximum (LGM). We examine 5 different realisations of the ocean's overturning circulation produced by a fully coupled atmosphere-ocean model under LGM forcing and suggested changes in the atmospheric deposition of iron and phytoplankton physiology at the LGM. Measured changes in carbon-13 and carbon-14, as well as a qualitative reconstruction of the change in ocean carbon export are used to evaluate the results. Overall, we find that while a reduction in ocean ventilation at the LGM is necessary to reproduce carbon-13 and carbon-14 observations, this circulation results in a low net sink for atmospheric CO2. In contrast, while biogeochemical processes contribute little to carbon isotopes, we propose that most of the change in atmospheric CO2 was due to such factors. However, the lesser role for circulation means that when all plausible factors are accounted for, most of the necessary CO2 change remains to be explained. This presents a serious challenge to our understanding of the mechanisms behind changes in the global carbon cycle during the geologic past.

Time valued life cycle greenhouse gas emissions from buildings

The UK has adopted legally binding carbon reduction targets of 34% by 2020 and 80% by 2050 (measured against the 1990 baseline). Buildings are estimated to be responsible for more than 50% of greenhouse gas (GHG) emissions in the UK. These consist of both operational, produced during use, and embodied, produced during manufacture of materials and components, and during construction, refurbishments and demolition. A brief assessment suggests that it is unlikely that UK emission reduction targets can be met without substantial reductions in both Oc and Ec. Oc occurs over the lifetime of a building whereas the bulk of Ec occurs at the start of a building’s life. A time value for emissions could influence the decision making process when it comes to comparing mitigation measures which have benefits that occur at different times. An example might be the choice between building construction using low Ec construction materials versus building construction using high Ec construction materials but with lower Oc, although the use of high Ec materials does not necessarily imply a lower Oc. Particular time related issues examined here are: the urgency of the need to achieve large emissions reductions during the next 10 to 20 years; the earlier effective action is taken, the less costly it will be; future reduction in carbon intensity of energy supply; the carbon cycle and relationship between the release of GHG’s and their subsequent concentrations in the atmosphere. An equation is proposed, which weights emissions according to when they occur during the building life cycle, and which effectively increases Ec as a proportion of the total, suggesting that reducing Ec is likely to be more beneficial, in terms of climate change, for most new buildings. Thus, giving higher priority to Ec reductions is likely to result in a bigger positive impact on climate change and mitigation costs.

Direct soil moisture controls of future global soil carbon changes: An important source of uncertainty

The nature of the climate–carbon cycle feedback depends critically on the response of soil carbon to climate, including changes in moisture. However, soil moisture–carbon feedback responses have not been investigated thoroughly. Uncertainty in the response of soil carbon to soil moisture changes could arise from uncertainty in the relationship between soil moisture and heterotrophic respiration. We used twelve soil moisture–respiration functions (SMRFs) with a soil carbon model (RothC) and data from a coupled climate–carbon cycle general circulation model to investigate the impact of direct heterotrophic respiration dependence on soil moisture on the climate carbon cycle feedback. Global changes in soil moisture acted to oppose temperature‐driven decreases in soil carbon and hence tended to increase soil carbon storage. We found considerable uncertainty in soil carbon changes due to the response of soil respiration to soil moisture. The use of different SMRFs resulted in both large losses and small gains in future global soil carbon stocks, whether considering all climate forcings or only moisture changes. Regionally, the greatest range in soil carbon changes across SMRFs was found where the largest soil carbon changes occurred. Further research is needed to constrain the soil moisture–respiration relationship and thus reduce uncertainty in climate–carbon cycle feedbacks. There may also be considerable uncertainty in the regional responses of soil carbon to soil moisture changes since climate model predictions of regional soil moisture changes are less coherent than temperature changes.