Assessing the accuracy and robustness of a process-based model for coffee agroforestry systems in Central America

Coffee is often grown in production systems associated with shade trees that provide different ecosystem services. Management, weather and soil conditions are spatially variable production factors. CAF2007 is a dynamic model for coffee agroforestry systems that takes these factors as inputs and simulates the processes underlying berry production at the field scale. There remain, however, uncertainties about process rates that need to be reduced through calibration. Bayesian statistics using Markov chain Monte Carlo algorithms is increasingly used for calibration of parameter-rich models. However, very few studies have employed multi-site calibration, which aims to reduce parameter uncertainties using data from multiple sites simultaneously. The main objectives of this study were to calibrate the coffee agroforestry model using data gathered in long-term experiments in Costa Rica and Nicaragua, and to test the calibrated model against independent data from commercial coffee-growing farms. Two sub-models were improved: calculation of flowering date and the modelling of biennial production patterns. The modified model, referred to as CAF2014, can be downloaded at https://doi.org/10.5281/zenodo.3608877. Calibration improved model performance (higher R2, lower RMSE) for Turrialba (Costa Rica) and Masatepe (Nicaragua), including when all experiments were pooled together. Multi-site and single-site Bayesian calibration led to similar RMSE. Validation on new data from coffee-growing farms revealed that both calibration methods improved simulation of yield and its bienniality. The thus improved model was used to test the effect of N fertilizer and shade in different locations on coffee yield.