Secondary production of the central rangeland region of the United States.

Rangelands are the dominant land use across a broad swath of central North America where they span a wide gradient, from <350 to >900 mm, in mean annual precipitation. Substantial efforts have examined temporal and spatial variation in aboveground net primary production (ANPP) to precipitation (PPT) across this gradient. In contrast, net secondary productivity (NSP, e.g., primary consumer production) has not been evaluated analogously. However, livestock production, which is a form of NSP or primary consumer production supported by primary production, is the dominant non-cultivated land use and an integral economic driver in these regions. Here, we used long-term (mean length = 19 years) ANPP and NSP data from six research sites across the Central Great Plains with a history of a conservative stocking to determine resource (i.e., PPT)-productivity relationships, NSP sensitivities to dry-year precipitation, and regional trophic efficiencies (e.g., NSP:ANPP ratio). PPT-ANPP relationships were linear for both temporal (site-based) and spatial (among site) gradients. The spatial PPT-NSP model revealed that PPT mediated a saturating relationship for NSP as sites became more mesic, a finding that contrasts with many plant-based PPT-ANPP relationships. A saturating response to high growing-season precipitation suggests biogeochemical rather than vegetation growth constraints may govern NSP (i.e., large herbivore production). Differential sensitivity in NSP to dry years demonstrated that the primary consumer production response heightened as sites became more xeric. Although sensitivity generally decreased with increasing precipitation as predicted from known PPT-ANPP relationships, evidence suggests that the dominant species' identity and traits influenced secondary production efficiency. Non-native northern mixed-grass prairie was outperformed by native Central Great Plains rangeland in sensitivity to dry years and efficiency in converting ANPP to NSP. A more comprehensive understanding of the mechanisms leading to differences in producer and consumer responses will require multisite experiments to assess biotic and abiotic determinants of multi-trophic level efficiency and sensitivity.

How Modelers Model: the Overlooked Social and Human Dimensions in Model Intercomparison Studies.

There is a growing realization that the complexity of model ensemble studies depends not only on the models used but also on the experience and approach used by modelers to calibrate and validate results, which remain a source of uncertainty. Here, we applied a multi-criteria decision-making method to investigate the rationale applied by modelers in a model ensemble study where 12 process-based different biogeochemical model types were compared across five successive calibration stages. The modelers shared a common level of agreement about the importance of the variables used to initialize their models for calibration. However, we found inconsistency among modelers when judging the importance of input variables across different calibration stages. The level of subjective weighting attributed by modelers to calibration data decreased sequentially as the extent and number of variables provided increased. In this context, the perceived importance attributed to variables such as the fertilization rate, irrigation regime, soil texture, pH, and initial levels of soil organic carbon and nitrogen stocks was statistically different when classified according to model types. The importance attributed to input variables such as experimental duration, gross primary production, and net ecosystem exchange varied significantly according to the length of the modeler's experience. We argue that the gradual access to input data across the five calibration stages negatively influenced the consistency of the interpretations made by the modelers, with cognitive bias in "trial-and-error" calibration routines. Our study highlights that overlooking human and social attributes is critical in the outcomes of modeling and model intercomparison studies. While complexity of the processes captured in the model algorithms and parameterization is important, we contend that (1) the modeler's assumptions on the extent to which parameters should be altered and (2) modeler perceptions of the importance of model parameters are just as critical in obtaining a quality model calibration as numerical or analytical details.

Ensemble modelling, uncertainty and robust predictions of organic carbon in long-term bare-fallow soils.

Simulation models represent soil organic carbon (SOC) dynamics in global carbon (C) cycle scenarios to support climate-change studies. It is imperative to increase confidence in long-term predictions of SOC dynamics by reducing the uncertainty in model estimates. We evaluated SOC simulated from an ensemble of 26 process-based C models by comparing simulations to experimental data from seven long-term bare-fallow (vegetation-free) plots at six sites: Denmark (two sites), France, Russia, Sweden and the United Kingdom. The decay of SOC in these plots has been monitored for decades since the last inputs of plant material, providing the opportunity to test decomposition without the continuous input of new organic material. The models were run independently over multi-year simulation periods (from 28 to 80 years) in a blind test with no calibration (Bln) and with the following three calibration scenarios, each providing different levels of information and/or allowing different levels of model fitting: (a) calibrating decomposition parameters separately at each experimental site (Spe); (b) using a generic, knowledge-based, parameterization applicable in the Central European region (Gen); and (c) using a combination of both (a) and (b) strategies (Mix). We addressed uncertainties from different modelling approaches with or without spin-up initialization of SOC. Changes in the multi-model median (MMM) of SOC were used as descriptors of the ensemble performance. On average across sites, Gen proved adequate in describing changes in SOC, with MMM equal to average SOC (and standard deviation) of 39.2 (±15.5) Mg C/ha compared to the observed mean of 36.0 (±19.7) Mg C/ha (last observed year), indicating sufficiently reliable SOC estimates. Moving to Mix (37.5 ± 16.7 Mg C/ha) and Spe (36.8 ± 19.8 Mg C/ha) provided only marginal gains in accuracy, but modellers would need to apply more knowledge and a greater calibration effort than in Gen, thereby limiting the wider applicability of models.

Assessing uncertainties in crop and pasture ensemble model simulations of productivity and N2 O emissions.

Simulation models are extensively used to predict agricultural productivity and greenhouse gas emissions. However, the uncertainties of (reduced) model ensemble simulations have not been assessed systematically for variables affecting food security and climate change mitigation, within multi-species agricultural contexts. We report an international model comparison and benchmarking exercise, showing the potential of multi-model ensembles to predict productivity and nitrous oxide (N2 O) emissions for wheat, maize, rice and temperate grasslands. Using a multi-stage modelling protocol, from blind simulations (stage 1) to partial (stages 2-4) and full calibration (stage 5), 24 process-based biogeochemical models were assessed individually or as an ensemble against long-term experimental data from four temperate grassland and five arable crop rotation sites spanning four continents. Comparisons were performed by reference to the experimental uncertainties of observed yields and N2 O emissions. Results showed that across sites and crop/grassland types, 23%-40% of the uncalibrated individual models were within two standard deviations (SD) of observed yields, while 42 (rice) to 96% (grasslands) of the models were within 1 SD of observed N2 O emissions. At stage 1, ensembles formed by the three lowest prediction model errors predicted both yields and N2 O emissions within experimental uncertainties for 44% and 33% of the crop and grassland growth cycles, respectively. Partial model calibration (stages 2-4) markedly reduced prediction errors of the full model ensemble E-median for crop grain yields (from 36% at stage 1 down to 4% on average) and grassland productivity (from 44% to 27%) and to a lesser and more variable extent for N2 O emissions. Yield-scaled N2 O emissions (N2 O emissions divided by crop yields) were ranked accurately by three-model ensembles across crop species and field sites. The potential of using process-based model ensembles to predict jointly productivity and N2 O emissions at field scale is discussed.

Feeding incremental amounts of ground flaxseed: effects on diversity and relative abundance of ruminal microbiota and enteric methane emissions in lactating dairy cows.

We evaluated the effects of incremental amounts of ground flaxseed (GFX) on diversity and relative abundance of ruminal microbiota taxa, enteric methane (CH4) emissions, and urinary excretion of purine derivatives (PD) in lactating dairy cows in a replicated 4 × 4 Latin square design. Twenty mid-lactation Jersey cows were used in the study. Of these 20 cows, 12 were used for ruminal sampling, 16 for enteric CH4 measurements, and all for spot urine collection. Each period lasted 21 d with 14 d for diet adaptation and 7 d for data and sample collection. Diets were formulated by replacing corn meal and soybean meal with 0%, 5%, 10%, and 15% of GFX in the diet's dry matter. Ruminal fluid samples obtained via stomach tubing were used for DNA extraction. Enteric CH4 production was measured using the sulfur hexafluoride tracer technique. Diets had no effect on ruminal microbiota diversity. Similarly, the relative abundance of ruminal archaea genera was not affected by diets. In contrast, GFX decreased or increased linearly the relative abundance of Firmicutes (P < 0.01) and Bacteroidetes (P < 0.01), respectively. The relative abundance of the ruminal bacteria Ruminococcus (P < 0.01) and Clostridium (P < 0.01) decreased linearly, and that of Prevotella (P < 0.01) and Pseudobutyrivibrio (P < 0.01) increased linearly with feeding GFX. A tendency for a linear reduction (P = 0.055) in enteric CH4 production (from 304 to 256 g/d) was observed in cows fed increasing amounts of GFX. However, neither CH4 yield nor CH4 intensity was affected by treatments. Diets had no effect on the urinary excretion of uric acid, allantoin, and total PD. Overall, feeding GFX decreased linearly the relative abundance of the ruminal bacterial genera Ruminococcus and Clostridium and enteric CH4 production, but no change was seen for CH4 yield and CH4 intensity, or urinary excretion of total PD, suggesting no detrimental effect of GFX on microbial protein synthesis in the rumen.

Quantification of global and national nitrogen budgets for crop production.

Input-output estimates of nitrogen on cropland are essential for improving nitrogen management and better understanding the global nitrogen cycle. Here, we compare 13 nitrogen budget datasets covering 115 countries and regions over 1961-2015. Although most datasets showed similar spatiotemporal patterns, some annual estimates varied widely among them, resulting in large ranges and uncertainty. In 2010, global medians (in TgN yr-1) and associated minimum-maximum ranges were 73 (64-84) for global harvested crop nitrogen; 161 (139-192) for total nitrogen inputs; 86 (68-97) for nitrogen surplus; and 46% (40-53%) for nitrogen use efficiency. Some of the most uncertain nitrogen budget terms by country showed ranges as large as their medians, revealing areas for improvement. A benchmark nitrogen budget dataset, derived from central tendencies of the original datasets, can be used in model comparisons and inform sustainable nitrogen management in food systems.

Global Research Alliance N2 O chamber methodology guidelines: Guidelines for gap-filling missing measurements.

Nitrous oxide (N2 O) is a potent greenhouse gas that is primarily emitted from agriculture. Sampling limitations have generally resulted in discontinuous N2 O observations over the course of any given year. The status quo for interpolating between sampling points has been to use a simple linear interpolation. This can be problematic with N2 O emissions, since they are highly variable and sampling bias around these peak emission periods can have dramatic impacts on cumulative emissions. Here, we outline five gap-filling practices: linear interpolation, generalized additive models (GAMs), autoregressive integrated moving average (ARIMA), random forest (RF), and neural networks (NNs) that have been used for gap-filling soil N2 O emissions. To facilitate the use of improved gap-filling methods, we describe the five methods and then provide strengths and challenges or weaknesses of each method so that model selection can be improved. We then outline a protocol that details data organization and selection, splitting of data into training and testing datasets, building and testing models, and reporting results. Use of advanced gap-filling methods within a standardized protocol is likely to increase transparency, improve emission estimates, reduce uncertainty, and increase capacity to quantify the impact of mitigation practices.

Global Research Alliance N2 O chamber methodology guidelines: Summary of modeling approaches.

Measurements of nitrous oxide (N2 O) emissions from agriculture are essential for understanding the complex soil-crop-climate processes, but there are practical and economic limits to the spatial and temporal extent over which measurements can be made. Therefore, N2 O models have an important role to play. As models are comparatively cheap to run, they can be used to extrapolate field measurements to regional or national scales, to simulate emissions over long time periods, or to run scenarios to compare mitigation practices. Process-based models can also be used as an aid to understanding the underlying processes, as they can simulate feedbacks and interactions that can be difficult to distinguish in the field. However, when applying models, it is important to understand the conceptual process differences in models, how conceptual understanding changed over time in various models, and the model requirements and limitations to ensure that the model is well suited to the purpose of the investigation and the type of system being simulated. The aim of this paper is to give the reader a high-level overview of some of the important issues that should be considered when modeling. This includes conceptual understanding of widely used models, common modeling techniques such as calibration and validation, assessing model fit, sensitivity analysis, and uncertainty assessment. We also review examples of N2 O modeling for different purposes and describe three commonly used process-based N2 O models (APSIM, DayCent, and DNDC).

The use of biogeochemical models to evaluate mitigation of greenhouse gas emissions from managed grasslands.

Simulation models quantify the impacts on carbon (C) and nitrogen (N) cycling in grassland systems caused by changes in management practices. To support agricultural policies, it is however important to contrast the responses of alternative models, which can differ greatly in their treatment of key processes and in their response to management. We applied eight biogeochemical models at five grassland sites (in France, New Zealand, Switzerland, United Kingdom and United States) to compare the sensitivity of modelled C and N fluxes to changes in the density of grazing animals (from 100% to 50% of the original livestock densities), also in combination with decreasing N fertilization levels (reduced to zero from the initial levels). Simulated multi-model median values indicated that input reduction would lead to an increase in the C sink strength (negative net ecosystem C exchange) in intensive grazing systems: -64 ±â€¯74 g C m-2 yr-1 (animal density reduction) and -81 ±â€¯74 g C m-2 yr-1 (N and animal density reduction), against the baseline of -30.5 ±â€¯69.5 g C m-2 yr-1 (LSU [livestock units] ≥ 0.76 ha-1 yr-1). Simulations also indicated a strong effect of N fertilizer reduction on N fluxes, e.g. N2O-N emissions decreased from 0.34 ±â€¯0.22 (baseline) to 0.1 ±â€¯0.05 g N m-2 yr-1 (no N fertilization). Simulated decline in grazing intensity had only limited impact on the N balance. The simulated pattern of enteric methane emissions was dominated by high model-to-model variability. The reduction in simulated offtake (animal intake + cut biomass) led to a doubling in net primary production per animal (increased by 11.6 ±â€¯8.1 t C LSU-1 yr-1 across sites). The highest N2O-N intensities (N2O-N/offtake) were simulated at mown and extensively grazed arid sites. We show the possibility of using grassland models to determine sound mitigation practices while quantifying the uncertainties associated with the simulated outputs.