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.

Inverse modeling of the 14C bomb pulse in stalagmites to constrain the dynamics of soil carbon cycling at selected European cave sites

The decomposition of soil organic matter (SOM) is temperature dependent, but its response to a future warmer climate remains equivocal. Enhanced rates of decomposition of SOM under increased global temperatures might cause higher CO2 emissions to the atmosphere, and could therefore constitute a strong positive feedback. The magnitude of this feedback however remains poorly understood, primarily because of the difficulty in quantifying the temperature sensitivity of stored, recalcitrant carbon that comprises the bulk (>90%) of SOM in most soils. In this study we investigated the effects of climatic conditions on soil carbon dynamics using the attenuation of the 14C ‘bomb’ pulse as recorded in selected modern European speleothems. These new data were combined with published results to further examine soil carbon dynamics, and to explore the sensitivity of labile and recalcitrant organic matter decomposition to different climatic conditions. Temporal changes in 14C activity inferred from each speleothem was modelled using a three pool soil carbon inverse model (applying a Monte Carlo method) to constrain soil carbon turnover rates at each site. Speleothems from sites that are characterised by semi-arid conditions, sparse vegetation, thin soil cover and high mean annual air temperatures (MAATs), exhibit weak attenuation of atmospheric 14C ‘bomb’ peak (a low damping effect, D in the range: 55–77%) and low modelled mean respired carbon ages (MRCA), indicating that decomposition is dominated by young, recently fixed soil carbon. By contrast, humid and high MAAT sites that are characterised by a thick soil cover and dense, well developed vegetation, display the highest damping effect (D = c. 90%), and the highest MRCA values (in the range from 350 ± 126 years to 571 ± 128 years). This suggests that carbon incorporated into these stalagmites originates predominantly from decomposition of old, recalcitrant organic matter. SOM turnover rates cannot be ascribed to a single climate variable, e.g. (MAAT) but instead reflect a complex interplay of climate (e.g. MAAT and moisture budget) and vegetation development.

OzFACE: the Australian savanna free air CO2 enrichment facility and its relevance to carbon-cycling issues in a tropical savanna

A free air CO2 enrichment (FACE) facility has recently been constructed in a tropical savanna in north-eastern Queensland, Australia. The system has a novel and cost-effective design and uses an industrial source of pure CO2 piped directly to the site. We describe the design details of this facility and assess the likely contribution it will make towards advancing our understanding of the direct impacts of rising atmospheric CO2 on savannas. These include addressing uncertainties about future shifts in the tree–grass balance and associated changes in carbon stocks, responses of C4 grasses in dry tropical environments, potential sequestration of soil carbon, and the modifications of CO2 responses by moisture and nutrient interactions. Tropical regions have been poorly represented in climate change research, and the work at the OzFACE facility will complement existing and ongoing FACE studies at temperate latitudes.

Holocene vegetation dynamics in the northeastern Rub' al-Khali desert, Arabian Peninsula: a phytolith, pollen and carbon isotope study

The Holocene vegetation history of the Arabian Peninsula is poorly understood, with few palaeobotanical studies to date. At Awafi, Ras al-Khaimah, UAE, a 3.3 m lake sediment sequence records the vegetation development for the period 8500 cal. yr BP to similar to3000 cal. yr BP. delta(13)C isotope, pollen and phytolith analyses indicate that C3 Pooid grassland with a strong woody element existed during the early Holocene (between 8500 and 6000 cal. yr BP) and became replaced by mixed C3 and C4 grasses with a strong C4 Panicoid tall grass element between 5900 and 5400 cal. yr BP. An intense, arid event Occurred at 4100 cal. yr BP when the lake desiccated and was infilled by Aeolian sand. From 4100 cal. yr BP the vegetation was dominated by C4 Chloridoid types and Cyperaceae, suggesting an incomplete vegetation cover and Aeolian dune reactivation owing to increased regional aridity. These data outline the ecosystem dynamics and carbon cycling in response to palaeomon-soon and north-westerly variability during the Holocene. Copyright (C) 2004 John Wiley Sons, Ltd.

Seasonal pattern of regional carbon balance in the central Rocky Mountains from surface and airborne measurements

[1] High-elevation forests represent a large fraction of potential carbon uptake in North America, but this uptake is not well constrained by observations. Additionally, forests in the Rocky Mountains have recently been severely damaged by drought, fire, and insect outbreaks, which have been quantified at local scales but not assessed in terms of carbon uptake at regional scales. The Airborne Carbon in the Mountains Experiment was carried out in 2007 partly to assess carbon uptake in western U.S. mountain ecosystems. The magnitude and seasonal change of carbon uptake were quantified by (1) paired upwind-downwind airborne CO2 observations applied in a boundary layer budget, (2) a spatially explicit ecosystem model constrained using remote sensing and flux tower observations, and (3) a downscaled global tracer transport inversion. Top-down approaches had mean carbon uptake equivalent to flux tower observations at a subalpine forest, while the ecosystem model showed less. The techniques disagreed on temporal evolution. Regional carbon uptake was greatest in the early summer immediately following snowmelt and tended to lessen as the region experienced dry summer conditions. This reduction was more pronounced in the airborne budget and inversion than in flux tower or upscaling, possibly related to lower snow water availability in forests sampled by the aircraft, which were lower in elevation than the tower site. Changes in vegetative greenness associated with insect outbreaks were detected using satellite reflectance observations, but impacts on regional carbon cycling were unclear, highlighting the need to better quantify this emerging disturbance effect on montane forest carbon cycling.

Simulated last glacial maximum D14Catm and the deep glacial ocean carbon reservoir

∆14Catm has been estimated as 420 ± 80‰ (IntCal09) during the Last Glacial Maximum (LGM) compared to preindustrial times (0‰), but mechanisms explaining this difference are not yet resolved. ∆14Catm is a function of both cosmogenic production in the high atmosphere and of carbon cycling and partitioning in the Earth system. 10Be-based reconstructions show a contribution of the cosmogenic production term of only 200 ± 200‰ in the LGM. The remaining 220‰ have thus to be explained by changes in the carbon cycle. Recently, Bouttes et al. (2010, 2011) proposed to explain most of the difference in pCO2atm and δ13C between glacial and interglacial times as a result of brine-induced ocean stratification in the Southern Ocean. This mechanism involves the formation of very saline water masses that contribute to high carbon storage in the deep ocean. During glacial times, the sinking of brines is enhanced and more carbon is stored in the deep ocean, lowering pCO2atm. Moreover, the sinking of brines induces increased stratification in the Southern Ocean, which keeps the deep ocean well isolated from the surface. Such an isolated ocean reservoir would be characterized by a low ∆14C signature. Evidence of such 14C-depleted deep waters during the LGM has recently been found in the Southern Ocean (Skinner et al. 2010). The degassing of this carbon with low ∆14C would then reduce ∆14Catm throughout the deglaciation. We have further developed the CLIMBER-2 model to include a cosmogenic production of 14C as well as an interactive atmospheric 14C reservoir. We investigate the role of both the sinking of brine and cosmogenic production, alongside iron fertilization mechanisms, to explain changes in ∆14Catm during the last deglaciation. In our simulations, not only is the sinking of brine mechanism consistent with past ∆14C data, but it also explains most of the differences in pCO2atm and ∆14Catm between the LGM and preindustrial times. Finally, this study represents the first time to our knowledge that a model experiment explains glacial-interglacial differences in pCO2atm, δ13C, and ∆14C together with a coherent LGM climate.

Enhanced Australian carbon sink despite increased wildfire during the 21st century

Climate projections show Australia becoming significantly warmer during the 21st century, and precipitation decreasing over much of the continent. Such changes are conventionally considered to increase wildfire risk. Nevertheless, we show that burnt area increases in southern Australia, but decreases in northern Australia. Overall the projected increase in fire is small (0.72–1.31% of land area, depending on the climate scenario used), and does not cause a decrease in carbon storage. In fact, carbon storage increases by 3.7–5.6 Pg C (depending on the climate scenario used). Using a process-based model of vegetation dynamics, vegetation–fire interactions and carbon cycling, we show increased fire promotes a shift to more fire-adapted trees in wooded areas and their encroachment into grasslands, with an overall increase in forested area of 3.9–11.9%. Both changes increase carbon uptake and storage. The increase in woody vegetation increases the amount of coarse litter, which decays more slowly than fine litter hence leading to a relative reduction in overall heterotrophic respiration, further reducing carbon losses. Direct CO2 effects increase woody cover, water-use efficiency and productivity, such that carbon storage is increased by 8.5–14.8 Pg C compared to simulations in which CO2 is held constant at modern values. CO2 effects tend to increase burnt area, fire fluxes and therefore carbon losses in arid areas, but increase vegetation density and reduce burnt area in wooded areas.

Earthworms produce granules of intricately zoned calcite

Earthworms of the family Lumbricidae, which includes many common species, produce and secrete up to millimeter-sized calcite granules, and the intricate fine-scale zoning of their constituent crystals is unique for a biomineral. Granule calcite is produced by crystallization of amorphous calcium carbonate (ACC) that initially precipitates within the earthworm calciferous glands, then forms protogranules by accretion on quartz grain cores. Crystallization of ACC is mediated by migrating fluid films and is largely complete within 24 11 of ACC production and before granules leave the earthworm. Variations in the density of defects formed as a byproduct of trace element incorporation during calcite crystall growth have generated zoning that can be resolved by cathodoluminescence imaging at ultraviolet to blue wavelengths and using the novel technique of scanning electron microscope charge contrast imaging. Mapping of calcite crystal orientations by electron backscatter diffraction reveals an approximate radial fabric to the granules that reflects crystal growth from internal nucleation sites toward their margins. The survival within granules of ACC inclusions for months after they enter soils indicates that they crystallize only within the earthworm and in the presence of fluids containing biochemical catalysts. The earthworm probably promotes crystallization of ACC in order to prevent remobilization of the calcium carbonate by dissolution. Calcite granules vividly illustrate the role of transient precursors in biomineralization, but the underlying question of why earth-worms produce granules in volumes sufficient to have a measurable impact on soil carbon cycling remains to be answered.

Fine-root turnover rates of European forests revisited: an analysis of data from sequential coring and ingrowth cores

Background and Aims Forest trees directly contribute to carbon cycling in forest soils through the turnover of their fine roots. In this study we aimed to calculate root turnover rates of common European forest tree species and to compare them with most frequently published values. Methods We compiled available European data and applied various turnover rate calculation methods to the resulting database. We used Decision Matrix and Maximum-Minimum formula as suggested in the literature. Results Mean turnover rates obtained by the combination of sequential coring and Decision Matrix were 0.86 yr−1 for Fagus sylvatica and 0.88 yr−1 for Picea abies when maximum biomass data were used for the calculation, and 1.11 yr−1 for both species when mean biomass data were used. Using mean biomass rather than maximum resulted in about 30 % higher values of root turnover. Using the Decision Matrix to calculate turnover rate doubled the rates when compared to the Maximum-Minimum formula. The Decision Matrix, however, makes use of more input information than the Maximum-Minimum formula. Conclusions We propose that calculations using the Decision Matrix with mean biomass give the most reliable estimates of root turnover rates in European forests and should preferentially be used in models and C reporting.

Interacting controls on innate sources of CO2 efflux from a calcareous arid zone soil under experimental acidification and wetting

More than half of global soil carbon is stored as carbonates, primarily in arid and semi-arid zones. Climate change models predict more frequent and severe rainfall events in some parts of the globe, many of which are dominated by calcareous soils. Such events trigger substantial increases in soil CO2 efflux. We hypothesised that the primary source of CO2 emissions from calcareous, arid zone soil during a single wetting event is abiotic and that soil acidification and wetting have a positive, potentially interacting, effect. We manipulated soil pH, soil moisture, and controlled soil respiration by gamma irradiating half of an 11 day incubation experiment. All manipulated experimental treatments had a rapid and enormous effect on CO2 emission. Respiration contributed ca. 5% of total CO2 efflux; the major source (carbonate buffering) varied depending on the extent of acidification and wetting. Maximum CO2 efflux occurred when pH was lowest and at intermediate matric potential. CO2 efflux was lowest at native pH when soil was air dry. Our data suggest that there may be an underestimate of soil-atmosphere carbon fluxes in arid ecosystems with calcareous soils. There is also a clear potential that these soils may become net carbon sources depending on changes in rainfall patterns, rainfall acidity, and future land management. Our findings have major implications for carbon cycling in arid zone soil and further study of carbon dynamics in these terrestrial systems at a landscape level will be required if we are to improve global climate and carbon cycling models.

Linking carbon and sulphur cycling during simulated drought cycles in peat from six sites across the UK

Water table draw-down is thought to increase peat decomposition and, therefore, DOC release. However, several studies have shown lower DOC concentrations during droughts relative to ‘normal’ periods with high water table. We carried out controlled incubation experiments at 10°C on 10x10 cm peat soil cores collected from six UK sites across a sulphur deposition gradient. Our aim was to quantify the balance between microbial consumption and chemical precipitation of DOC due to episodic acidification driven by sulphur redox reactions by comparing changes in soil water chemistry to microbial activity (i.e. soil respiration and trace gas fluxes). During dry periods, all sites showed a concurrent increase in SO4 and soil respiration and a decline in DOC. However, the magnitude of change in both DOC and SO4 varied considerably between sites according to historical sulphur deposition loads and the variation in acid/base chemistry.

Modeling the mechanisms that control in-stream dissolved organic carbon dynamics in upland and forested catchments

[1] We present a new, process-based model of soil and stream water dissolved organic carbon (DOC): the Integrated Catchments Model for Carbon (INCA-C). INCA-C is the first model of DOC cycling to explicitly include effects of different land cover types, hydrological flow paths, in-soil carbon biogeochemistry, and surface water processes on in-stream DOC concentrations. It can be calibrated using only routinely available monitoring data. INCA-C simulates daily DOC concentrations over a period of years to decades. Sources, sinks, and transformation of solid and dissolved organic carbon in peat and forest soils, wetlands, and streams as well as organic carbon mineralization in stream waters are modeled. INCA-C is designed to be applied to natural and seminatural forested and peat-dominated catchments in boreal and temperate regions. Simulations at two forested catchments showed that seasonal and interannual patterns of DOC concentration could be modeled using climate-related parameters alone. A sensitivity analysis showed that model predictions were dependent on the mass of organic carbon in the soil and that in-soil process rates were dependent on soil moisture status. Sensitive rate coefficients in the model included those for organic carbon sorption and desorption and DOC mineralization in the soil. The model was also sensitive to the amount of litter fall. Our results show the importance of climate variability in controlling surface water DOC concentrations and suggest the need for further research on the mechanisms controlling production and consumption of DOC in soils.

Utilizing chromophoric dissolved organic matter measurements to derive export and reactivity of dissolved organic carbon exported to the Arctic Ocean: A case study of the Yukon River, Alaska

The quality and quantity of dissolved organic matter (DOM) exported by Arctic rivers is known to vary with hydrology and this exported material plays a fundamental role in the biogeochemical cycling of carbon at high latitudes. We highlight the potential of optical measurements to examine DOM quality across the hydrograph in Arctic rivers. Furthermore, we establish chromophoric DOM (CDOM) relationships to dissolved organic carbon (DOC) and lignin phenols in the Yukon River and model DOC and lignin loads from CDOM measurements, the former in excellent agreement with long-term DOC monitoring data. Intensive sampling across the historically under-sampled spring flush period highlights the importance of this time for total export of DOC and particularly lignin. Calculated riverine DOC loads to the Arctic Ocean show an increase from previous estimates, especially when new higher discharge data are incorporated. Increased DOC loads indicate decreased residence times for terrigenous DOM in the Arctic Ocean with important implications for the reactivity and export of this material to the Atlantic Ocean.

Amino-acid cycling drives nitrogen fixation in the legume - Rhizobium symbiosis

The biological reduction of atmospheric N-2 to ammonium (nitrogen fixation) provides about 65% of the biosphere's available nitrogen. Most of this ammonium is contributed by legume rhizobia symbioses(1), which are initiated by the infection of legume hosts by bacteria (rhizobia), resulting in formation of root nodules. Within the nodules, rhizobia are found as bacteroids, which perform the nitrogen fixation: to do this, they obtain sources of carbon and energy from the plant, in the form of dicarboxylic acids(2,3). It has been thought that, in return, bacteroids simply provide the plant with ammonium. But here we show that a more complex amino-acid cycle is essential for symbiotic nitrogen fixation by Rhizobium in pea nodules. The plant provides amino acids to the bacteroids, enabling them to shut down their ammonium assimilation. In return, bacteroids act like plant organelles to cycle amino acids back to the plant for asparagine synthesis. The mutual dependence of this exchange prevents the symbiosis being dominated by the plant, and provides a selective pressure for the evolution of mutualism.