Use of advanced modelling methods to estimate radiata pine productivity indices

•Productivity indices have been used to describe age normalised height and volume.•We predict productivity indices for radiata pine using a range of modelling methods.•Regression kriging using Random Forest, XGBoost and regression were most precise.•An ensemble model using all three methods was more precise than individual models.•New Zealand productivity index maps for radiata pine were developed.