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.