Abstract: Buried heating-cable method for manipulating soil temperature was designed and tested its performance in large concrete lysimeters grown with the wheat crop in Denmark. Soil temperature in heated plots was elevated by 5℃ compared with that in control by burying heating-cable at 0.1 m depth in a plough layer. Temperature sensors were placed at 0.05, 0.1 and 0.25 m depths in soil, and 0.1 m above the soil surface in all plots, which were connected to an automated data logger. Soil-warming setup was able to maintain a mean seasonal temperature difference of 5.0 ± 0.005℃ between heated and control plots at 0.1 m depth while the mean seasonal rise in soil temperature in the top 0.25 m depth (plough layer) was 3℃. Soil temperature in control plots froze (≤ 0℃) for 15 and 13 days respectively at 0.05 and 0.1 m depths while it did not in heated plots during the coldest period (Nov-Apr). This study clearly showed the efficacy of buried heating-cable technique in simulating soil temperature, and thus offers a simple, effective and alternative technique to study soil biogeochemical processes under warmer climates. This technique, however, decouples below-ground soil responses from that of above-ground vegetation response as this method heats only the soil. Therefore, using infrared heaters seems to represent natural climate warming (both air and soil) much more closely and may be used for future climate manipulation field studies.

Abstract: As the degradation of global ecosystem services (ES) continues in the last five decades, maintaining or even enhancing the ES of agro-ecosystem is one of the approaches to mitigate the global ES loss. This study provides the first estimate of an economic valuation of ES provided by organic cereal crop production systems with different management practices in relation to organic matter input (low, medium and high). Our results show that organic matter inputs significantly affect the total ES value on organic cereal crop production systems. The system with high organic matter input has the highest gross total ES value (US$ 1969 ha(-1) yr(-1)), followed by the low organic matter input system (US$ 1688 ha(-1) yr(-1)), and the lowest ES value are found in the medium organic matter input system (US$ 1492 ha(-1) yr(-1)). Organic matter inputs have strong positive relationship with non-marketable ES values, while this relationship was not found in marketable ES values. Monetizing the ES can be used by land managers and policy makers to adjust management practices in terms of organic matter input in cereal production system with a long term goal for sustainable agriculture.

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.