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.

Free atmospheric CO2 enrichment increased above ground biomass but did not affect symbiotic N2-fixation and soil carbon dynamics in a mixed deciduous stand in Wales

Through increases in net primary production (NPP), elevated CO2 is hypothesizes to increase the amount of plant litter entering the soil. The fate of this extra carbon on the forest floor or in mineral soil is currently not clear. Moreover, increased rates of NPP can be maintained only if forests can escape nitrogen limitation. In a Free atmospheric CO2 Enrichment (FACE) experiment near Bangor, Wales, 4 ambient CO2 and 4 FACE plots were planted with patches of Betula pendula, Alnus glutinosa and Fagus sylvatica on a former arable field. Four years after establishment, only a shallow L forest floor litter layer had formed due to intensive bioturbation. Total soil C and N contents increased irrespective of treatment and species as a result of afforestation. We could not detect an additional C sink in the soil, nor were soil C stabilization processes affected by FACE. We observed a decrease of leaf N content in Betula and Alnus under FACE, while the soil C/N ratio decreased regardless of CO2 treatment. The ratio of N taken up from the soil and by N2-fixation in Alnus was not affected by FACE. We infer that increased nitrogen use efficiency is the mechanism by which increased NPP is sustained under elevated CO2 at this site.

The role of soil microbes in the global carbon cycle: tracking the below-ground microbial processing of plant-derived carbon for manipulating carbon dynamics in agricultural systems.

It is well known that atmospheric concentrations of carbon dioxide (CO2) (and other greenhouse gases) have increased markedly as a result of human activity since the industrial revolution. It is perhaps less appreciated that natural and managed soils are an important source and sink for atmospheric CO2 and that, primarily as a result of the activities of soil microorganisms, there is a soil-derived respiratory flux of CO2 to the atmosphere that overshadows by tenfold the annual CO2 flux from fossil fuel emissions. Therefore small changes in the soil carbon cycle could have large impacts on atmospheric CO2 concentrations. Here we discuss the role of soil microbes in the global carbon cycle and review the main methods that have been used to identify the microorganisms responsible for the processing of plant photosynthetic carbon inputs to soil. We discuss whether application of these techniques can provide the information required to underpin the management of agro-ecosystems for carbon sequestration and increased agricultural sustainability. We conclude that, although crucial in enabling the identification of plant-derived carbon-utilising microbes, current technologies lack the high-throughput ability to quantitatively apportion carbon use by phylogentic groups and its use efficiency and destination within the microbial metabolome. It is this information that is required to inform rational manipulation of the plant–soil system to favour organisms or physiologies most important for promoting soil carbon storage in agricultural soil.

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.

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.

Multi-season eddy covariance observations of energy, water and carbon fluxes over a suburban area in Swindon, UK

Eddy covariance measurements of the turbulent sensible heat, latent heat and carbon dioxide fluxes for 12 months (2011–2012) are reported for the first time for a suburban area in the UK. The results from Swindon are comparable to suburban studies of similar surface cover elsewhere but reveal large seasonal variability. Energy partitioning favours turbulent sensible heat during summer (midday Bowen ratio 1.4–1.6) and latent heat in winter (0.05–0.7). A significant proportion of energy is stored (and released) by the urban fabric and the estimated anthropogenic heat flux is small but non-negligible (0.5–0.9 MJ m−2 day−1). The sensible heat flux is negative at night and for much of winter daytimes, reflecting the suburban nature of the site (44% vegetation) and relatively low built fraction (16%). Latent heat fluxes appear to be water limited during a dry spring in both 2011 and 2012, when the response of the surface to moisture availability can be seen on a daily timescale. Energy and other factors are more relevant controls at other times; at night the wind speed is important. On average, surface conductance follows a smooth, asymmetrical diurnal course peaking at around 6–9 mm s−1, but values are larger and highly variable in wet conditions. The combination of natural (vegetative) and anthropogenic (emission) processes is most evident in the temporal variation of the carbon flux: significant photosynthetic uptake is seen during summer, whilst traffic and building emissions explain peak release in winter (9.5 g C m−2 day−1). The area is a net source of CO2 annually. Analysis by wind direction highlights the role of urban vegetation in promoting evapotranspiration and offsetting CO2 emissions, especially when contrasted against peak traffic emissions from sectors with more roads. Given the extent of suburban land use, these results have important implications for understanding urban energy, water and carbon dynamics.

Contributions of matric and osmotic potentials to the unfrozen water content of frozen soils

Recent reports show that biogeochemical processes continue when the soil is frozen, but are limited by water availability. However, there is little knowledge about the interactive effects of soil and environmental variables on amounts of unfrozen water in frozen soils. The aims of this study were to determine the contributions of matric and osmotic potentials to the unfrozen water content of frozen soil. We determined the effects of matric and osmotic potential on unfrozen water contents of frozen mineral soil fractions (ranging from coarse sand to fine silt) at -7 degrees C, and estimated the contributions of these potentials to liquid water contents in samples from organic surface layers of boreal soils frozen at -4 degrees C. In the mineral soil fractions the unfrozen water contents appeared to be governed solely by the osmotic potential, but in the humus layers of the sampled boreal soils both the osmotic and matric potentials control unfrozen water content, with osmotic potential contributing 20 to 69% of the total water potential. We also determined pore size equivalents, where unfrozen water resides at -4 degrees C, and found a strong correlation between these equivalents and microbial CO2 production. The larger the pores in which the unfrozen water is found the larger the microbial activity that can be sustained. The osmotic potential may therefore be a key determinant of unfrozen water and carbon dynamics in frozen soil. (C) 2008 Elsevier B.V. All rights reserved.

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.

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.

Robust dynamics of Amazon dieback to climate change with perturbed ecosystem model parameters

Climate change science is increasingly concerned with methods for managing and integrating sources of uncertainty from emission storylines, climate model projections, and ecosystem model parameterizations. In tropical ecosystems, regional climate projections and modeled ecosystem responses vary greatly, leading to a significant source of uncertainty in global biogeochemical accounting and possible future climate feedbacks. Here, we combine an ensemble of IPCC-AR4 climate change projections for the Amazon Basin (eight general circulation models) with alternative ecosystem parameter sets for the dynamic global vegetation model, LPJmL. We evaluate LPJmL simulations of carbon stocks and fluxes against flux tower and aboveground biomass datasets for individual sites and the entire basin. Variability in LPJmL model sensitivity to future climate change is primarily related to light and water limitations through biochemical and water-balance-related parameters. Temperature-dependent parameters related to plant respiration and photosynthesis appear to be less important than vegetation dynamics (and their parameters) for determining the magnitude of ecosystem response to climate change. Variance partitioning approaches reveal that relationships between uncertainty from ecosystem dynamics and climate projections are dependent on geographic location and the targeted ecosystem process. Parameter uncertainty from the LPJmL model does not affect the trajectory of ecosystem response for a given climate change scenario and the primary source of uncertainty for Amazon 'dieback' results from the uncertainty among climate projections. Our approach for describing uncertainty is applicable for informing and prioritizing policy options related to mitigation and adaptation where long-term investments are required.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

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.