Abstract: Mitigation of greenhouse gas emissions from agriculture should be assessed across cropping systems and agroclimatic regions. In this study, we investigate the ability of the FASSET model to analyze differences in the magnitude of N2O emissions due to soil, climate and management factors in cereal-based cropping systems. Forage maize was grown in a conventional dairy system at Mabegondo (NW Spain) and wheat and barley in organic and conventional crop rotations at Foulum (NW Denmark). These two European sites represent agricultural areas with high and low to moderate emission levels, respectively. Field trials included plots with and without catch crops that were fertilized with either mineral N fertilizer, cattle slurry, pig slurry or digested manure. Non-fertilized treatments were also included. Measurements of N2O fluxes during the growing cycle of all the crops at both sites were performed with the static chamber method with more frequent measurements post-fertilization and biweekly measurements when high fluxes were not expected. All cropping systems were simulated with the FASSET version 2.5 simulation model. Cumulative soil seasonal N2O emissions were about ten-fold higher at Mabegondo than at Foulum when averaged across systems and treatments (8.99 and 0.71 kg N2O-N ha(-1), respectively). The average simulated cumulative soil N2O emissions were 9.03 and 1.71 kg N2O-N ha(-1) at Mabegondo and at Foulum, respectively. Fertilization, catch crops and cropping systems had lower influence on the seasonal soil N2O fluxes than the environmental factors. Overall, in its current version FASSET reproduced the effects of the different factors investigated on the cumulative seasonal soil N2O emissions but temporally it overestimated emissions from nitrification and denitrification on particular days when soil operations, ploughing or fertilization, took place. The errors associated with simulated daily soil N2O fluxes increased with the magnitude of the emissions. For resolving causes of differences in simulated and measured fluxes more intensive and temporally detailed measurements of N2O fluxes and soil C and N dynamics would be needed. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: The differential impact on soil organic carbon (SOC) of applying no tillage (NT) compared to conventional tillage (CT, i.e. mouldboard ploughing), along with three rates of nitrogen (N) fertilizer application (0,90 and 180 kg ha(-1) y(-1)), was studied under rain-fed Mediterranean conditions in a long-term experiment based on a durum wheat-maize rotation, in which crop residues were left on the soil (NT) or incorporated (CT). Observed SOC content following 8 and 12 years of continuous treatment application was significantly higher in the top 10 cm of the soil under NT than CT, but it was similar in the 10-40 cm layer. NT grain yields for both maize and durum wheat were below those attained under CT (on average 32% and 14% lower respectively) at a given rate of N fertilizer application. Soil, climate and crop data over 5 years were used to calibrate DSSAT model in order to simulate the impact of the different management practices over a 50-year period. Good agreement was obtained between observed and simulated values for crops grain yield, above-ground biomass and observed SOC values. Results from the simulations showed that under NT the weeds growing during the intercrop fallow period made a significant contribution to the observed SOC increase. When the contribution of the weed fallow was considered, NT significantly increased SOC in the top 40 cm of the soil at an average rate of 0.43, 0.31 and 0.03 t ha(-1) per year, respectively for 180,90 and 0 kg N ha(-1) year(-1), within the simulated 50 years. Under CT, a significant SOC increase was simulated under N180 and a significant decrease when no fertilizer was supplied.

Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.