BIOLP is an Integer Programming model based on the Balanced Compromise Programming multi-criteria decision method. The aim of BIOLP is to determine how a set of land use types should be distributed over space and time in order to optimize the multi-dimensional land performance of a region. Trajectories were defined as the succession of specific land use types over 30 years, assuming that land use changes can only occur at fixed intervals of 10 years. A database that represents the Tabacay catchment (Ecuador) as a set of land units with associated performance values was used as the input for BIOLP, which was then executed to determine the trajectories distribution that optimizes regional performance. The sensitivity of BIOLP to uncertainty in the input data, simulated through random variations on the performance values, was also tested. BIOLP showed a relative stability on its results under these conditions of stochastic, restricted changes. Additionally, the behaviour of BIOLP under different settings of its balancing and relative importance parameters was studied. Stronger variations on the outcomes were observed in this case, which indicate the influential role that such parameters play. Finally, the inclusion of performance thresholds in BIOLP was tested through the addition of sample constraints that required some of the criteria at stake to exceed predefined values. The outcome of the optimization exercises makes clear that the phenomenon of trade off between the provisioning service of the land (income) and the regulation and maintenance services (runoff, sediment, SOC) is crucial. BIOLP succeeds in accounting for this complex multi-dimensional phenomenon when determining the optimal spatio-temporal distributions of land use types. Despite this complexity, it is confirmed that the weights attributed to the provisioning or to the regulation and maintenance services are the main determinants for having the land use distributions dominated by either agriculture or forest.