Abstract: Apart from measurements, numerical models are the most convenient instruments to analyze the carbon and water balance of terrestrial ecosystems and their interactions with changing environmental conditions. The process-based Biome-BGC model is widely used to simulate the storage and flux of water, carbon, and nitrogen within the vegetation, litter, and soil of unmanaged terrestrial ecosystems. Considering herbaceous vegetation related simulations with Biome-BGC, soil moisture and growing season control on ecosystem functioning is inaccurate due to the simple soil hydrology and plant phenology representation within the model. Consequently, Biome-BGC has limited applicability in herbaceous ecosystems because (1) they are usually managed; (2) they are sensitive to soil processes, most of all hydrology; and (3) their carbon balance is closely connected with the growing season length. Our aim was to improve the applicability of Biome-BGC for managed herbaceous ecosystems by implementing several new modules, including management. A new index (heatsum growing season index) was defined to accurately estimate the first and the final days of the growing season. Instead of a simple bucket soil sub-model, a multilayer soil sub-model was implemented, which can handle the processes of runoff, diffusion and percolation. A new module was implemented to simulate the ecophysiological effect of drought stress on plant mortality. Mowing and grazing modules were integrated in order to quantify the functioning of managed ecosystems. After modifications, the Biome-BGC model was calibrated and validated using eddy covariance-based measurement data collected in Hungarian managed grassland ecosystems. Model calibration was performed based on the Bayes theorem. As a result of these developments and calibration, the performance of the model was substantially improved. Comparison with measurement-based estimate showed that the start and the end of the growing season are now predicted with an average accuracy of 5 and 4 days instead of 46 and 85 days as in the original model. Regarding the different sites and modeled fluxes (gross primary production, total ecosystem respiration, evapotranspiration), relative errors were between 18-60% using the original model and 10-18% using the developed model; squares of the correlation coefficients were between 0.02-0.49 using the original model and 0.50-0.81 using the developed model. (c) 2011 Elsevier B.V. All rights reserved.

Abstract: Predicting regional and global carbon (C) and water dynamics on grasslands has become of major interest, as grasslands are one of the most widespread vegetation types worldwide, providing a number of ecosystem services (such as forage production and C storage). The present study is a contribution to a regional-scale analysis of the C and water cycles on managed grasslands. The mechanistic biogeochemical model PaSim (Pasture Simulation model) was evaluated at 12 grassland sites in Europe. A new parameterization was obtained on a common set of eco-physiological parameters, which represented an improvement of previous parameterization schemes (essentially obtained via calibration at specific sites). We found that C and water fluxes estimated with the parameter set are in good agreement with observations. The model with the new parameters estimated that European grassland are a sink of C with 213 g C m(-2) yr(-1), which is close to the observed net ecosystem exchange (NEE) flux of the studied sites (185 g C m(-2) yr(-1) on average). The estimated yearly average gross primary productivity (GPP) and ecosystem respiration (RECO) for all of the study sites are 1220 and 1006 g C m(-2) yr(-1), respectively, in agreement with observed average GPP (1230 g C m(-2) yr(-1)) and RECO (1046 g C m(-2) yr(-1)). For both variables aggregated on a weekly basis, the root mean square error (RMSE) was similar to 5-16 g C week(-1) across the study sites, while the goodness of fit (R-2) was similar to 0.4-0.9. For evapotranspiration (ET), the average value of simulated ET (415 mmyr(-1)) for all sites and years is close to the average value of the observed ET (451 mm yr(-1)) by flux towers (on a weekly basis, RMSE similar to 2-8 mm week(-1); R-2 = 0.3-0.9). However, further model development is needed to better represent soil water dynamics under dry conditions and soil temperature in winter. A quantification of the uncertainties introduced by spatially generalized parameter values in C and water exchange estimates is also necessary. In addition, some uncertainties in the input management data call for the need to improve the quality of the observational system.

Abstract: Croplands play an important role in the carbon budget of many regions. However, the estimation of their carbon balance remains difficult due to diversity and complexity of the processes involved. We report the coupling of a one-dimensional soil water, heat, and CO2 flux model (SOILCO2), a pool concept of soil carbon turnover (RothC), and a crop growth module (SUCROS) to predict the net ecosystem exchange (NEE) of carbon. The coupled model, further referred to as AgroC, was extended with routines for managed grassland as well as for root exudation and root decay. In a first step, the coupled model was applied to two winter wheat sites and one upland grassland site in Germany. The model was calibrated based on soil water content, soil temperature, biometric, and soil respiration measurements for each site, and validated in terms of hourly NEE measured with the eddy covariance technique. The overall model performance of AgroC was sufficient with a model efficiency above 0.78 and a correlation coefficient above 0.91 for NEE. In a second step, AgroC was optimized with eddy covariance NEE measurements to examine the effect of different objective functions, constraints, and data-transformations on estimated NEE. It was found that NEE showed a distinct sensitivity to the choice of objective function and the inclusion of soil respiration data in the optimization process. In particular, both positive and negative day- and nighttime fluxes were found to be sensitive to the selected optimization strategy. Additional consideration of soil respiration measurements improved the simulation of small positive fluxes remarkably. Even though the model performance of the selected optimization strategies did not diverge substantially, the resulting cumulative NEE over simulation time period differed substantially. Therefore, it is concluded that data transformations, definitions of objective functions, and data sources have to be considered cautiously when a terrestrial ecosystem model is used to determine NEE by means of eddy covariance measurements. (C) 2017 Elsevier B.V. All rights reserved.