Analysing the impact of environmental changes due to climate change, e.g. geographical spread of climate-sensitive infections (CSIs) and agriculture crop modelling, may require land surface modelling (LSM) to predict future land surface conditions. There are multiple LSMs to choose from. The paper proposes a multivariate spatio-temporal data science method to understand the inherent uncertainties in four LSMs and the variations between them in Nordic areas for the net primary production.