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: Process-based models (PBM) for simulation of weather dependent grass growth can assist farmers and plant breeders in addressing the challenges of climate change by simulating alternative roads of adaptation. They can also provide management decision support under current conditions. A drawback of existing grass models is that they do not take into account the effect of winter stresses, limiting their use for full-year simulations in areas where winter survival is a key factor for yield security. Here, we present a novel full-year PBM for grassland named BASGRA. It was developed by combining the LINGRA grassland model (Van Oijen et al., 2005a) with models for cold hardening and soil physical winter processes. We present the model and show how it was parameterized for timothy (Phleum pratense L.), the most important forage grass in Scandinavia and parts of North America and Asia. Uniquely, BASGRA simulates the processes taking place in the sward during the transition from summer to winter, including growth cessation and gradual cold hardening, and functions for simulating plant injury due to low temperatures, snow and ice affecting regrowth in spring. For the calibration, we used detailed data from five different locations in Norway, covering a wide range of agroclimatic regions, day lengths (latitudes from 59 degrees to 70 degrees N) and soil conditions. The total dataset included 11 variables, notably above-ground dry matter, leaf area index, tiller density, content of C reserves, and frost tolerance. All data were used in the calibration. When BASGRA was run with the maximum a-posteriori (MAP) parameter vector from the single, Bayesian calibration, nearly all measured variables were simulated to an overall normalized root mean squared error (NRMSE) <0.5. For many site x experiment combinations, NRMSE was <0.3. The temporal dynamics were captured well for most variables, as evaluated by comparing simulated time courses versus data for the individual sites. The results may suggest that BASGRA is a reasonably robust model, allowing for simulation of growth and several important underlying processes with acceptable accuracy for a range of agroclimatic conditions. However, the robustness of the model needs to be tested further using independent data from a wide range of growing conditions. Finally we show an example of application of the model, comparing overwintering risks in two climatically different sites, and discuss future model applications. Further development work should include improved simulation of the dynamics of C reserves, and validation of winter tiller dynamics against independent data. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Renewed interest in heat stress effects on livestock productivity derives from climate change, which is expected to increase temperatures and the frequency of extreme weather events. This study aimed at evaluating the effect of temperature and humidity on milk production in highly selected dairy cattle populations across 3 European regions differing in climate and production systems to detect differences and similarities that can be used to optimize heat stress (HS) effect modeling. Milk, fat, and protein test day data from official milk recording for 1999 to 2010 in 4 Holstein populations located in the Walloon Region of Belgium (BEL), Luxembourg (LUX), Slovenia (SLO), and southern Spain (SPA) were merged with temperature and humidity data provided by the state meteorological agencies. After merging, the number of test day records/cows per trait ranged from 686,726/49,655 in SLO to 1,982,047/136,746 in BEL. Values for the daily average and maximum temperature-humidity index (THIavg and THImax) ranges for THIavg/THImax were largest in SLO (22-74/28-84) and shortest in SPA (39-76/46-83). Change point techniques were used to determine comfort thresholds, which differed across traits and climatic regions. Milk yield showed an inverted U-shaped pattern of response across the THI scale with a HS threshold around 73 THImax units. For fat and protein, thresholds were lower than for milk yield and were shifted around 6 THI units toward larger values in SPA compared with the other countries. Fat showed lower HS thresholds than protein traits in all countries. The traditional broken line model was compared with quadratic and cubic fits of the pattern of response in production to increasing heat loads. A cubic polynomial model allowing for individual variation in patterns of response and THIavg as heat load measure showed the best statistical features. Higher/lower producing animals showed less/more persistent production (quantity and quality) across the THI scale. The estimated correlations between comfort and THIavg values of 70 (which represents the upper end of the THIavg scale in BEL-LUX) were lower for BEL-LUX (0.70-0.80) than for SPA (0.83-0.85). Overall, animals producing in the more temperate climates and semi-extensive grazing systems of BEL and LUX showed HS at lower heat loads and more re-ranking across the THI scale than animals producing in the warmer climate and intensive indoor system of SPA.