|
Minet, J., Laloy, E., Tychon, B., & François, L. (2015). Bayesian inversions of a dynamic vegetation model at four European grassland sites. Biogeosciences, 12(9), 2809–2829.
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.
|
|
|
Kraus, D., Weller, S., Klatt, S., Haas, E., Wassmann, R., Kiese, R., et al. (2015). A new LandscapeDNDC biogeochemical module to predict CH4 and N2O emissions from lowland rice and upland cropping systems. Plant Soil, 386(1-2), 125–149.
Abstract: Replacing paddy rice by upland systems such as maize cultivation is an on-going trend in SE Asia caused by increasing water scarcity and higher demand for meat. How such land management changes will feedback on soil C and N cycles and soil greenhouse gas emissions is not well understood at present. A new LandscapeDNDC biogeochemical module was developed that allows the effect of land management changes on soil C and N cycle to be simulated. The new module is applied in combination with further modules simulating microclimate and crop growth and evaluated against observations from field experiments. The model simulations agree well with observed dynamics of CH (4) emissions in paddy rice depending on changes in climatic conditions and agricultural management. Magnitude and peak emission periods of N (2) O from maize cultivation are simulated correctly, though there are still deficits in reproducing day-to-day dynamics. These shortcomings are most likely related to simulated soil hydrology and may only be resolved if LandscapeDNDC is coupled to more complex hydrological models. LandscapeDNDC allows for simulation of changing land management practices in SE Asia. The possibility to couple LandscapeDNDC to more complex hydrological models is a feature needed to better understand related effects on soil-atmosphere-hydrosphere interactions.
|
|
|
Höglind, M., Van Oijen, M., Cameron, D., & Persson, T. (2016). Process-based simulation of growth and overwintering of grassland using the BASGRA model. Ecol. Model., 335, 1–15.
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.
|
|
|
Ventrella, D., Stellacci, A. M., Castrignanò, A., Charfeddine, M., & Castellini, M. (2016). Effects of crop residue management on winter durum wheat productivity in a long term experiment in Southern Italy. European Journal of Agronomy, 77, 188–198.
Abstract: A long-term experiment comparing different crop residue (CR) managements was established in 1977 in Foggia (Apulia region, southern Italy). The objective of this study was to investigate the long-term effects of different types of crop residue management on main yield response parameters in a continuous cropping system of winter durum wheat. In order to correctly interpret the results, models accounting for spatial error autocorrelation were used and compared with ordinary least square models. Eight crop residue management treatments, based on burning of wheat straw and stubble or their incorporation with or without N fertilization and irrigation, were compared. The experimental design was a complete randomized block with five replicates. Results indicated that the dynamics of yield, grain protein content and hectolitric weight of winter durum wheat did not show any decline as usually expected when a monoculture is carried out for a long time. In addition, the temporal variability of productivity was more affected by meteorological factors, such as air temperature and rainfall, than CR management treatments. Higher wheat grain yields and hectolitric weights quite frequently occurred after burning of wheat straw compared with straw incorporation without nitrogen fertilization and autumn irrigation and this was attributed to temporary mineral N immobilization in the soil. The rate of 50 kg ha(-1) of N seemed to counterbalance this negative effect when good condition of soil moisture occurred in the autumn period, so yielding the same productive level of straw burning treatment. (C) 2016 Elsevier B.V. All rights reserved.
|
|
|
Carabano, M. J., Logar, B., Bormann, J., Minet, J., Vanrobays, M. L., Diaz, C., et al. (2016). Modeling heat stress under different environmental conditions. J. Dairy Sci., 99(5), 3798–3814.
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.
|
|