Abstract: The DNDC, LandscapeDNDC and IAP-N-GAS models have been designed to simulate the carbon and nitrogen processes of terrestrial ecosystems. Until now, a comparison of these models using simultaneous observations has not been reported, although such a comparison is essential for further model development and application. This study aimed to evaluate the performance of the models, delineate the strengths and limitations of each model for simulating soil nitrous oxide (N2O) and nitric oxide (NO) emissions, and explore short-comings of these models that may require reconsideration. We conducted comparisons among the models using simultaneous observations of both gases and relevant variables from the winter wheat-summer maize rotation system at three field sites with calcareous soils. Simulations of N2O and NO emissions by the three models agreed well with annual observations, but not with daily observations. All models failed to correctly simulate soil moisture, which could explain some of the incorrect daily fluxes of N2O and NO, especially for intensive fluxes during the growing season. Multi-model ensembles are promising approaches to better simulate daily gas emissions. IAP-N-GAS underestimated the priming effect of straw incorporation on N2O and NO emissions, but better results were obtained with DNDC95 and LandscapeDNDC. LandscapeDNDC and IAP-N-GAS need to improve the simulation of irrigation water allocation and residue decomposition processes, respectively, and together to distinguish different irrigation methods as DNDC95 does. All three models overestimated the emissions of the nitrogenous gases for high nitrogen fertilizer (>430 kg N ha(-1) yr(-1)) addition treatments, and therefore, future research should focus more on the simulation of the limitation of soil dissolvable organic carbon on denitrification in calcareous soils.

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.

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.