Parameters Optimization Method of Ecosystem Model Based on Bayesian Machine Learning
Parameter optimization is an effective means for the accurate estimation of ecosystem model parameters and the reduction of the uncertainty in model predictions. We proposed a method for parameter optimization of the ecosystem model, which is based on the Bayesian machine learning and called No-U-Turn-Sampler (NUTS). As an efficient means of parameter optimization, NUTS uses a recursive algorithm to build a set of candidate points to obtain the posterior information of the parameters. If the constraint condition of"Non-U-Turn"is met, subtrees will be built to update parameters. Otherwise,"the optimal"set of parameters from current sample will be recorded, and then the next sampling begin to run until enough samples are taken. This algorithm avoids sampling redundancy caused by random walk and thus improves the efficiency of parameter optimization. Taking the carbon flux simulations of the Qianyanzhou subtropical coniferous plantation as an example, we implemented the parameter inversion of the carbon flux (Net Ecosystem Exchange, NEE) model using the NUTS method based on the Pymc3 framework. The comparison between the inversion results of NUTS and Metropolis-Hastings (MH) shows that the sampling frequency reduces about 85%, and the optimization efficiency increases about 3 times when the parameter values of the NUTS algorithm reaches convergence. The uncertainties of the seven parameters estimated by NUTS in the two NEE models are reduced by 10%-53%compared to MH. The NEE simulation improved significantly, with the R2 between the simulated values and the observed values increased by 23%and 17%, respectively and the RMSE decreased by 3%and 4%, respectively. In sum, the NUTS parameter optimization method proposed in this paper provides an efficient approach for the parameter optimization in ecosystem modeling.