How can soil monitoring networks be used to improve predictions of organic carbon pool dynamics and CO2 fluxes in agricultural soils?

Abstract As regional and continental carbon balances of terrestrial ecosystems become available, it becomes clear that the soils are the largest source of uncertainty. Repeated inventories of soil organic carbon (SOC) organized in soil monitoring networks (SMN) are being implemented in a number of countries. This paper reviews the concepts and design of SMNs in ten countries, and discusses the contribution of such networks to reducing the uncertainty of soil carbon balances. Some SMNs are designed to estimate country-specific land use or management effects on SOC stocks, while others collect soil carbon and ancillary data to provide a nationally consistent assessment of soil carbon condition across the major land-use/soil type combinations. The former use a single sampling campaign of paired sites, while for the latter both systematic (usually grid based) and stratified repeated sampling campaigns (5–10 years interval) are used with densities of one site per 10–1,040 km². For paired sites, multiple samples at each site are taken in order to allow statistical analysis, while for the single sites, composite samples are taken. In both cases, fixed depth increments together with samples for bulk density and stone content are recommended. Samples should be archived to allow for re-measurement purposes using updated techniques. Information on land management, and where possible, land use history should be systematically recorded for each site. A case study of the agricultural frontier in Brazil is presented in which land use effect factors are calculated in order to quantify the CO2 fluxes from national land use/management conversion matrices. Process-based SOC models can be run for the individual points of the SMN, provided detailed land management records are available. These studies are still rare, as most SMNs have been implemented recently or are in progress. Examples from the USA and Belgium show that uncertainties in SOC change range from 1.6–6.5 Mg C ha−1 for the prediction of SOC stock changes on individual sites to 11.72 Mg C ha−1 or 34% of the median SOC change for soil/land use/climate units. For national SOC monitoring, stratified sampling sites appears to be the most straightforward attribution of SOC values to units with similar soil/land use/climate conditions (i.e. a spatially implicit upscaling approach). Keywords Soil monitoring networks - Soil organic carbon - Modeling - Sampling design

N2O emissions from agricultural lands: a synthesis of simulation approaches

Nitrous oxide (N2O) is primarily produced by the microbially-mediated nitrification and denitrification processes in soils. It is influenced by a suite of climate (i.e. temperature and rainfall) and soil (physical and chemical) variables, interacting soil and plant nitrogen (N) transformations (either competing or supplying substrates) as well as land management practices. It is not surprising that N2O emissions are highly variable both spatially and temporally. Computer simulation models, which can integrate all of these variables, are required for the complex task of providing quantitative determinations of N2O emissions. Numerous simulation models have been developed to predict N2O production. Each model has its own philosophy in constructing simulation components as well as performance strengths. The models range from those that attempt to comprehensively simulate all soil processes to more empirical approaches requiring minimal input data. These N2O simulation models can be classified into three categories: laboratory, field and regional/global levels. Process-based field-scale N2O simulation models, which simulate whole agroecosystems and can be used to develop N2O mitigation measures, are the most widely used. The current challenge is how to scale up the relatively more robust field-scale model to catchment, regional and national scales. This paper reviews the development history, main construction components, strengths, limitations and applications of N2O emissions models, which have been published in the literature. The three scale levels are considered and the current knowledge gaps and challenges in modelling N2O emissions from soils are discussed.

Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide

Aboveground–belowground interactions exert critical controls on the composition and function of terrestrial ecosystems, yet the fundamental relationships between plant diversity and soil microbial diversity remain elusive. Theory predicts predominantly positive associations but tests within single sites have shown variable relationships, and associations between plant and microbial diversity across broad spatial scales remain largely unexplored. We compared the diversity of plant, bacterial, archaeal and fungal communities in one hundred and forty-five 1 m2 plots across 25 temperate grassland sites from four continents. Across sites, the plant alpha diversity patterns were poorly related to those observed for any soil microbial group. However, plant beta diversity (compositional dissimilarity between sites) was significantly correlated with the beta diversity of bacterial and fungal communities, even after controlling for environmental factors. Thus, across a global range of temperate grasslands, plant diversity can predict patterns in the composition of soil microbial communities, but not patterns in alpha diversity.

Potential soil carbon sequestration in overgrazed grassland ecosystems

Excessive grazing pressure is detrimental to plant productivity and may lead to declines in soil organic matter. Soil organic matter is an important source of plant nutrients and can enhance soil aggregation, limit soil erosion, and can also increase cation exchange and water holding capacities, and is, therefore, a key regulator of grassland ecosystem processes. Changes in grassland management which reverse the process of declining productivity can potentially lead to increased soil C. Thus, rehabilitation of areas degraded by overgrazing can potentially sequester atmospheric C. We compiled data from the literature to evaluate the influence of grazing intensity on soil C. Based on data contained within these studies, we ascertained a positive linear relationship between potential C sequestration and mean annual precipitation which we extrapolated to estimate global C sequestration potential with rehabilitation of overgrazed grassland. The GLASOD and IGBP DISCover data sets were integrated to generate a map of overgrazed grassland area for each of four severity classes on each continent. Our regression model predicted losses of soil C with decreased grazing intensity in drier areas (precipitation less than 333 mm yr(-1)), but substantial sequestration in wetter areas. Most (93%) C sequestration potential occurred in areas with MAP less than 1800 mm. Universal rehabilitation of overgrazed grasslands can sequester approximately 45 Tg C yr(-1), most of which can be achieved simply by cessation of overgrazing and implementation of moderate grazing intensity. Institutional level investments by governments may be required to sequester additional C.

Stabilization mechanisms of soil organic matter: Implications for C-saturation of soils

The relationship between soil structure and the ability of soil to stabilize soil organic matter (SOM) is a key element in soil C dynamics that has either been overlooked or treated in a cursory fashion when developing SOM models. The purpose of this paper is to review current knowledge of SOM dynamics within the framework of a newly proposed soil C saturation concept. Initially, we distinguish SOM that is protected against decomposition by various mechanisms from that which is not protected from decomposition. Methods of quantification and characteristics of three SOM pools defined as protected are discussed. Soil organic matter can be: (1) physically stabilized, or protected from decomposition, through microaggregation, or (2) intimate association with silt and clay particles, and (3) can be biochemically stabilized through the formation of recalcitrant SOM compounds. In addition to behavior of each SOM pool, we discuss implications of changes in land management on processes by which SOM compounds undergo protection and release. The characteristics and responses to changes in land use or land management are described for the light fraction (LF) and particulate organic matter (POM). We defined the LF and POM not occluded within microaggregates (53-250 mum sized aggregates as unprotected. Our conclusions are illustrated in a new conceptual SOM model that differs from most SOM models in that the model state variables are measurable SOM pools. We suggest that physicochemical characteristics inherent to soils define the maximum protective capacity of these pools, which limits increases in SOM (i.e. C sequestration) with increased organic residue inputs.

Effect of biochar amendment on the soil-atmosphere exchange of greenhouse gases from an intensive subtropical pasture in northern New South Wales, Australia

We assessed the effect of biochar incorporation into the soil on the soil-atmosphere exchange of the greenhouse gases (GHG) from an intensive subtropical pasture. For this, we measured N2O, CH4 and CO2 emissions with high temporal resolution from April to June 2009 in an existing factorial experiment where cattle feedlot biochar had been applied at 10 t ha-1 in November 2006. Over the whole measurement period, significant emissions of N2O and CO2 were observed, whereas a net uptake of CH4 was measured. N2O emissions were found to be highly episodic with one major emission pulse (up to 502 µg N2O-N m-2 h 1) following heavy rainfall. There was no significant difference in the net flux of GHGs from the biochar amended vs. the control plots. Our results demonstrate that intensively managed subtropical pastures on ferrosols in northern New South Wales of Australia can be a significant source of GHG. Our hypothesis that the application of biochar would lead to a reduction in emissions of GHG from soils was not supported in this field assessment. Additional studies with longer observation periods are needed to clarify the long term effect of biochar amendment on soil microbial processes and the emission of GHGs under field conditions.

Oxygen transport in soil and the vertical distribution of roots

An analytical solution for steady-state oxygen transport in soils including 2 sink terms, viz roots and microbes with the corresponding vertical distribution scaling lengths forming a ratio p, showed p governed the critical air-filled porosity, θ_c, needed by most plants. For low temperature and p, θ_c was <0.1 but at higher temperatures and p = 1, θ_c was >0.15 m³/m³. When root length density at the surface was 10⁴ m/m³ and p > 3, θ_c was 0.25 m³/m³, more than half the pore space. Few combinations of soil and climate regularly meet this condition. However, for sandy soils and seasonally warm, arid regions, the theory is consistent with observation, in that plants may have some deep roots. Critical θ_c values are used to formulate theoretical solutions in a forward mode, so different levels of oxygen uptake by roots may be compared to microbial activity. The proportion of respiration by plant roots increases rapidly with p up to p ≈2.

Calibration functions for estimating soil moisture from GPR dielectric constant measurements

Ground-penetrating radar (GPR) is widely used for assessment of soil moisture variability in field soils. Because GPR does not measure soil water content directly, it is common practice to use calibration functions that describe its relationship with the soil dielectric properties and textural parameters. However, the large variety of models complicates the selection of the appropriate function. In this article an overview is presented of the different functions available, including volumetric models, empirical functions, effective medium theories, and frequency-specific functions. Using detailed information presented in summary tables, the choice for which calibration function to use can be guided by the soil variables available to the user, the frequency of the GPR equipment, and the desired level of detail of the output. This article can thus serve as a guide for GPR practitioners to obtain soil moisture values and to estimate soil dielectric properties.

Herbivores and nutrients control grassland plant diversity via light limitation

Human alterations to nutrient cycles1, 2 and herbivore communities3, 4, 5, 6, 7 are affecting global biodiversity dramatically2. Ecological theory predicts these changes should be strongly counteractive: nutrient addition drives plant species loss through intensified competition for light, whereas herbivores prevent competitive exclusion by increasing ground-level light, particularly in productive systems8, 9. Here we use experimental data spanning a globally relevant range of conditions to test the hypothesis that herbaceous plant species losses caused by eutrophication may be offset by increased light availability due to herbivory. This experiment, replicated in 40 grasslands on 6 continents, demonstrates that nutrients and herbivores can serve as counteracting forces to control local plant diversity through light limitation, independent of site productivity, soil nitrogen, herbivore type and climate. Nutrient addition consistently reduced local diversity through light limitation, and herbivory rescued diversity at sites where it alleviated light limitation. Thus, species loss from anthropogenic eutrophication can be ameliorated in grasslands where herbivory increases ground-level light.

Using the Gini coefficient with BIOLOG substrate utilisation data to provide an alternative quantitative measure for comparing bacterial soil communities

A measure quantifying unequal use of carbon sources, the Gini coefficient (G), has been developed to allow comparisons of the observed functional diversity of bacterial soil communities. This approach was applied to the analysis of substrate utilisation data obtained from using BIOLOG microtiter plates in a study which compared decomposition processes in two contrasting plant substrates in two different soils. The relevance of applying the Gini coefficient as a measure of observed functional diversity, for soil bacterial communities is evaluated against the Shannon index (H) and average well colour development (AWCD), a measure of the total microbial activity. Correlation analysis and analysis of variance of the experimental data show that the Gini coefficient, the Shannon index and AWCD provided similar information when used in isolation. However, analyses based on the Gini coefficient and the Shannon index, when total activity on the microtiter plates was maintained constant (i.e. AWCD as a covariate), indicate that additional information about the distribution of carbon sources being utilised can be obtained. We demonstrate that the Lorenz curve and its measure of inequality, the Gini coefficient, provides not only comparable information to AWCD and the Shannon index but when used together with AWCD encompasses measures of total microbial activity and absorbance inequality across all the carbon sources. This information is especially relevant for comparing the observed functional diversity of soil microbial communities.

Diverse urban plantings managed with sufficient resource availability can increase plant productivity and arthropod diversity

Buildings structures and surfaces are explicitly being used to grow plants, and these “urban plantings” are generally designed for aesthetic value. Urban plantings also have the potential to contribute significant “ecological values” by increasing urban habitat for animals such as arthropods and by increasing plant productivity. In this study, we evaluated how the provision of these additional ecological values is affected by plant species richness; the availability of essential resources for plants, such as water, light, space; and soil characteristics. We sampled 33 plantings located on the exterior of three buildings in the urban center of Brisbane, Australia (subtropical climatic region) over 2, 6 week sampling periods characterized by different temperature and rainfall conditions. Plant cover was estimated as a surrogate for productivity as destructive sampling of biomass was not possible. We measured weekly light levels (photosynthetically active radiation), plant CO2 assimilation, soil CO2 efflux, and arthropod diversity. Differences in plant cover were best explained by a three-way interaction of plant species richness, management water regime and sampling period. As the richness of plant species increased in a planter, productivity and total arthropod richness also increased significantly—likely due to greater habitat heterogeneity and quality. Overall we found urban plantings can provide additional ecological values if essential resources are maintained within a planter such as water, light and soil temperature. Diverse urban plantings that are managed with these principles in mind can contribute to the attraction of diverse arthropod communities, and lead to increased plant productivity within a dense urban context.

Hydrological consequences of land-cover change: Quantifying the influence of plants on soil moisture with time-lapse electrical resistivity

Electrical resistivity of soils and sediments is strongly influenced by the presence of interstitial water. Taking advantage of this dependency, electrical-resistivity imaging (ERI) can be effectively utilized to estimate subsurface soil-moisture distributions. The ability to obtain spatially extensive data combined with time-lapse measurements provides further opportunities to understand links between land use and climate processes. In natural settings, spatial and temporal changes in temperature and porewater salinity influence the relationship between soil moisture and electrical resistivity. Apart from environmental factors, technical, theoretical, and methodological ambiguities may also interfere with accurate estimation of soil moisture from ERI data. We have examined several of these complicating factors using data from a two-year study at a forest-grassland ecotone, a boundary between neighboring but different plant communities.At this site, temperature variability accounts for approximately 20-45 of resistivity changes from cold winter to warm summer months. Temporal changes in groundwater conductivity (mean=650 S/cm =57.7) and a roughly 100-S/cm spatial difference between the forest and grassland had only a minor influence on the moisture estimates. Significant seasonal fluctuations in temperature and precipitation had negligible influence on the basic measurement errors in data sets. Extracting accurate temporal changes from ERI can be hindered by nonuniqueness of the inversion process and uncertainties related to time-lapse inversion schemes. The accuracy of soil moisture obtained from ERI depends on all of these factors, in addition to empirical parameters that define the petrophysical soil-moisture/resistivity relationship. Many of the complicating factors and modifying variables to accurately quantify soil moisture changes with ERI can be accounted for using field and theoretical principles.

Soil N2O emissions under N2-fixing legumes and N-fertilised canola: A reappraisal of emissions factor calculations

Introducing nitrogen (N)-fixing legumes into cereal-based crop rotations reduces synthetic fertiliser-N use and may mitigate soil emissions of nitrous oxide (N2O). Current IPCC calculations assume 100% of legume biomass N as the anthropogenic N input and use 1% of this as an emission factor (EF)—the percentage of input N emitted as N2O. However, legumes also utilise soil inorganic N, so legume-fixed N is typically less than 100% of legume biomass N. In two field experiments, we measured soil N2O emissions from a black Vertosol in sub-tropical Australia for 12 months after sowing of chickpea (Cicer arietinum L.), canola (Brassica napus L.), faba bean (Vicia faba L.), and field pea (Pisum sativum L.). Cumulative N2O emissions from N-fertilised canola (624 g N2O-N ha−1) greatly exceeded those from chickpea (127 g N2O-N ha−1) in Experiment 1. Similarly, N2O emitted from canola (385 g N2O-N ha−1) in Experiment 2 was significantly greater than chickpea (166 g N2O-N ha−1), faba bean (166 g N2O-N ha−1) or field pea (135 g N2O-N ha−1). Highest losses from canola were recorded during the growing season, whereas 75% of the annual N2O losses from the legumes occurred post-harvest. Legume N2-fixation provided 37–43% (chickpea), 54% (field pea) and 64% (faba bean) of total plant biomass N. Using only fixed-N inputs, we calculated EFs for chickpea (0.13–0.31%), field pea (0.18%) and faba bean (0.04%) that were significantly less than N-fertilised canola (0.48–0.78%) (P < 0.05), suggesting legume-fixed N is a less emissive form of N input to the soil than fertiliser N. Inputs of legume-fixed N should be more accurately quantified to properly gauge the potential for legumes to mitigate soil N2O emissions. EF’s from legume crops need to be revised and should include a factor for the proportion of the legume’s N derived from the atmosphere.

Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe

Soil microorganisms are critical to ecosystem functioning and the maintenance of soil fertility. However, despite global increases in the inputs of nitrogen (N) and phosphorus (P) to ecosystems due to human activities, we lack a predictive understanding of how microbial communities respond to elevated nutrient inputs across environmental gradients. Here we used high-throughput sequencing of marker genes to elucidate the responses of soil fungal, archaeal, and bacterial communities using an N and P addition experiment replicated at 25 globally distributed grassland sites. We also sequenced metagenomes from a subset of the sites to determine how the functional attributes of bacterial communities change in response to elevated nutrients. Despite strong compositional differences across sites, microbial communities shifted in a consistent manner with N or P additions, and the magnitude of these shifts was related to the magnitude of plant community responses to nutrient inputs. Mycorrhizal fungi and methanogenic archaea decreased in relative abundance with nutrient additions, as did the relative abundances of oligotrophic bacterial taxa. The metagenomic data provided additional evidence for this shift in bacterial life history strategies because nutrient additions decreased the average genome sizes of the bacterial community members and elicited changes in the relative abundances of representative functional genes. Our results suggest that elevated N and P inputs lead to predictable shifts in the taxonomic and functional traits of soil microbial communities, including increases in the relative abundances of faster-growing, copiotrophic bacterial taxa, with these shifts likely to impact belowground ecosystems worldwide.